-
Understanding the Dynamical State of Globular Clusters: Core-Collapsed vs Non Core-Collapsed
Authors:
Sourav Chatterjee,
Stefan Umbreit,
John M. Fregeau,
Frederic A. Rasio
Abstract:
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…
▽ More
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.
△ Less
Submitted 20 November, 2012; v1 submitted 12 July, 2012;
originally announced July 2012.
-
On the Possibility of Tidal Formation of Binary Planets Around Ordinary Stars
Authors:
Philipp Podsiadlowski,
Saul Rappaport,
John M. Fregeau,
Rosemary A. Mardling
Abstract:
The planet formation process and subsequent planet migration may lead to configurations resulting in strong dynamical interactions among the various planets. Well-studied possible outcomes include collisions between planets, scattering events that eject one or more of the planets, and a collision of one or more of the planets with the parent star. In this work we consider one other possibility tha…
▽ More
The planet formation process and subsequent planet migration may lead to configurations resulting in strong dynamical interactions among the various planets. Well-studied possible outcomes include collisions between planets, scattering events that eject one or more of the planets, and a collision of one or more of the planets with the parent star. In this work we consider one other possibility that has seemingly been overlooked in the various scattering calculations presented in the literature: the tidal capture of two planets which leads to the formation of a binary planet (or binary brown dwarf) in orbit about the parent star. We carry out extensive numerical simulations of such dynamical and tidal interactions to explore the parameter space for the formation of such binary planets. We show that tidal formation of binary planets is possible for typical planet masses and distances from the host star. The detection (or lack thereof) of planet-planet binaries can thus be used to constrain the properties of planetary systems, including their mutual spacing during formation, and the fraction of close planets in very eccentric orbits which are believed to form by a closely related process.
△ Less
Submitted 8 July, 2010;
originally announced July 2010.
-
Monte Carlo Simulations of Globular Cluster Evolution. V. Binary Stellar Evolution
Authors:
Sourav Chatterjee,
John M. Fregeau,
Stefan Umbreit,
Frederic A. Rasio
Abstract:
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 co…
▽ More
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.
△ Less
Submitted 23 December, 2009;
originally announced December 2009.
-
Monte Carlo Simulations of Globular Cluster Evolution. VI. The Influence of an Intermediate Mass Black Hole
Authors:
Stefan Umbreit,
John M. Fregeau,
Sourav Chatterjee,
Frederic A. Rasio
Abstract:
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 allowin…
▽ More
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.
△ Less
Submitted 28 February, 2012; v1 submitted 28 October, 2009;
originally announced October 2009.
-
Evolution of the Binary Fraction in Dense Stellar Systems
Authors:
J. M. Fregeau,
N. Ivanova,
F. A. Rasio
Abstract:
Using our recently improved Monte Carlo evolution code, we study the evolution of the binary fraction in globular clusters. In agreement with previous N-body simulations, we find generally that the hard binary fraction in the core tends to increase with time over a range of initial cluster central densities for initial binary fractions <~ 90%. The dominant processes driving the evolution of the…
▽ More
Using our recently improved Monte Carlo evolution code, we study the evolution of the binary fraction in globular clusters. In agreement with previous N-body simulations, we find generally that the hard binary fraction in the core tends to increase with time over a range of initial cluster central densities for initial binary fractions <~ 90%. The dominant processes driving the evolution of the core binary fraction are mass segregation of binaries into the cluster core and preferential destruction of binaries there. On a global scale, these effects and the preferential tidal stripping of single stars tend to roughly balance, leading to overall cluster binary fractions that are roughly constant with time. Our findings suggest that the current hard binary fraction near the half-mass radius is a good indicator of the hard primordial binary fraction. However, the relationship between the true binary fraction and the fraction of main-sequence stars in binaries (which is typically what observers measure) is non-linear and rather complicated. We also consider the importance of soft binaries, which not only modify the evolution of the binary fraction, but can drastically change the evolution of the cluster as a whole. Finally, we describe in some detail the recent addition of single and binary stellar evolution to our cluster evolution code.
△ Less
Submitted 24 July, 2009; v1 submitted 23 July, 2009;
originally announced July 2009.
-
The Dynamical Effects of White Dwarf Birth Kicks in Globular Star Clusters
Authors:
John M. Fregeau,
Harvey B. Richer,
Frederic A. Rasio,
Jarrod R. Hurley
Abstract:
Recent observations of the white dwarf (WD) populations in the Galactic globular cluster NGC 6397 suggest that WDs receive a kick of a few km/s shortly before they are born. Using our Monte Carlo cluster evolution code, which includes accurate treatments of all relevant physical processes operating in globular clusters, we study the effects of the kicks on their host cluster and on the WD popula…
▽ More
Recent observations of the white dwarf (WD) populations in the Galactic globular cluster NGC 6397 suggest that WDs receive a kick of a few km/s shortly before they are born. Using our Monte Carlo cluster evolution code, which includes accurate treatments of all relevant physical processes operating in globular clusters, we study the effects of the kicks on their host cluster and on the WD population itself. We find that in clusters whose velocity dispersion is comparable to the kick speed, WD kicks are a significant energy source for the cluster, prolonging the initial cluster core contraction phase significantly so that at late times the cluster core to half-mass radius ratio is a factor of up to ~ 10 larger than in the no-kick case. WD kicks thus represent a possible resolution of the large discrepancy between observed and theoretically predicted values of this key structural parameter. Our modeling also reproduces the observed trend for younger WDs to be more extended in their radial distribution in the cluster than older WDs.
△ Less
Submitted 2 March, 2009; v1 submitted 6 February, 2009;
originally announced February 2009.
-
The Imprints of IMBHs on the Structure of Globular Clusters: Monte-Carlo Simulations
Authors:
Stefan Umbreit,
John M. Fregeau,
Frederic A. Rasio
Abstract:
We present the first results of a series of Monte-Carlo simulations investigating the imprint of a central black hole on the core structure of a globular cluster. We investigate the three-dimensional and the projected density profile of the inner regions of idealized as well as more realistic globular cluster models, taking into account a stellar mass spectrum, stellar evolution and allowing for…
▽ More
We present the first results of a series of Monte-Carlo simulations investigating the imprint of a central black hole on the core structure of a globular cluster. We investigate the three-dimensional and the projected density profile of the inner regions of idealized as well as more realistic globular cluster models, taking into account a stellar mass spectrum, 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 simulations published previously in the literature.
△ Less
Submitted 5 March, 2008;
originally announced March 2008.
-
X-Ray Binaries and the Dynamical States of Globular Clusters
Authors:
John M. Fregeau
Abstract:
We summarize and discuss recent work (Fregeau 2007) that presents the confluence of three results suggesting that most Galactic globular clusters are still in the process of core contraction, and have not yet reached the thermal equilibrium phase driven by binary scattering interactions: that 1) the three clusters that appear to be overabundant in X-ray binaries per unit encounter frequency are…
▽ More
We summarize and discuss recent work (Fregeau 2007) that presents the confluence of three results suggesting that most Galactic globular clusters are still in the process of core contraction, and have not yet reached the thermal equilibrium phase driven by binary scattering interactions: that 1) the three clusters that appear to be overabundant in X-ray binaries per unit encounter frequency are observationally classified as "core-collapsed," 2) recent numerical simulations of cluster evolution with primordial binaries show that structural parameters of clusters in the binary-burning phase agree only with "core-collapsed" clusters, and 3) a cluster in the binary-burning phase for the last few Gyr should have about 5 times more dynamically formed X-ray sources than if it were in the core contraction phase for the same time.
△ Less
Submitted 7 January, 2008;
originally announced January 2008.
-
X-Ray Binaries and the Current Dynamical States of Galactic Globular Clusters
Authors:
John M. Fregeau
Abstract:
It has been known for over 30 years that Galactic globular clusters (GCs) are overabundant by orders of magnitude in bright X-ray sources per unit mass relative to the disk population. Recently a quantitative understanding of this phenomenon has developed, with a clear correlation between the number of X-ray sources in a cluster, $N_X$, and the cluster's encounter frequency, $Γ$, becoming appare…
▽ More
It has been known for over 30 years that Galactic globular clusters (GCs) are overabundant by orders of magnitude in bright X-ray sources per unit mass relative to the disk population. Recently a quantitative understanding of this phenomenon has developed, with a clear correlation between the number of X-ray sources in a cluster, $N_X$, and the cluster's encounter frequency, $Γ$, becoming apparent. We derive a refined version of $Γ$ that incorporates the finite lifetime of X-ray sources and the dynamical evolution of clusters. With it we find we are able to explain the few clusters that lie off the $N_X$--$Γ$ correlation, and resolve the discrepancy between observed GC core radii and the values predicted by theory. Our results suggest that most GCs are still in the process of core contraction and have not yet reached the thermal equilibrium phase driven by binary scattering interactions.
△ Less
Submitted 3 December, 2007; v1 submitted 24 October, 2007;
originally announced October 2007.
-
Binary Stars and Globular Cluster Dynamics
Authors:
John M. Fregeau
Abstract:
In this brief proceedings article I summarize the review talk I gave at the IAU 246 meeting in Capri, Italy, glossing over the well-known results from the literature, but paying particular attention to new, previously unpublished material. This new material includes a careful comparison of the apparently contradictory results of two independent methods used to simulate the evolution of binary po…
▽ More
In this brief proceedings article I summarize the review talk I gave at the IAU 246 meeting in Capri, Italy, glossing over the well-known results from the literature, but paying particular attention to new, previously unpublished material. This new material includes a careful comparison of the apparently contradictory results of two independent methods used to simulate the evolution of binary populations in dense stellar systems (the direct N-body method of Hurley, et al. 2007 and the approximate Monte Carlo method of Ivanova, et al. 2005), that shows that the two methods may not actually yield contradictory results, and suggests future work to more directly compare the two methods.
△ Less
Submitted 16 October, 2007;
originally announced October 2007.
-
Effects of stellar collisions on star cluster evolution and core collapse
Authors:
Sourav Chatterjee,
John M. Fregeau,
Frederic A. Rasio
Abstract:
We systematically study the effects of collisions on the overall dynamical evolution of dense star clusters using Monte Carlo simulations over many relaxation times. We derive many observable properties of these clusters, including their core radii and the radial distribution of collision products. We also study different aspects of collisions in a cluster taking into account the shorter lifetim…
▽ More
We systematically study the effects of collisions on the overall dynamical evolution of dense star clusters using Monte Carlo simulations over many relaxation times. We derive many observable properties of these clusters, including their core radii and the radial distribution of collision products. We also study different aspects of collisions in a cluster taking into account the shorter lifetimes of more massive stars, which has not been studied in detail before. Depending on the lifetimes of the significantly more massive collision products, observable properties of the cluster can be modified qualitatively; for example, even without binaries, core collapse can sometimes be avoided simply because of stellar collisions.
△ Less
Submitted 16 October, 2007;
originally announced October 2007.
-
Eccentric double white dwarfs as LISA sources in globular clusters
Authors:
B. Willems,
V. Kalogera,
A. Vecchio,
N. Ivanova,
F. A. Rasio,
J. M. Fregeau,
K. Belczynski
Abstract:
We consider the formation of double white dwarfs (DWDs) through dynamical interactions in globular clusters. Such interactions can give rise to eccentric DWDs, in contrast to the exclusively circular population expected to form in the Galactic disk. We show that for a 5-year Laser Interferometer Space Antenna (LISA) mission and distances as far as the Large Magellanic Cloud, multiple harmonics f…
▽ More
We consider the formation of double white dwarfs (DWDs) through dynamical interactions in globular clusters. Such interactions can give rise to eccentric DWDs, in contrast to the exclusively circular population expected to form in the Galactic disk. We show that for a 5-year Laser Interferometer Space Antenna (LISA) mission and distances as far as the Large Magellanic Cloud, multiple harmonics from eccentric DWDs can be detected at a signal-to-noise ratio higher than 8 for at least a handful of eccentric DWDs, given their formation rate and typical lifetimes estimated from current cluster simulations. Consequently the association of eccentricity with stellar-mass LISA sources does not uniquely involve neutron stars, as is usually assumed. Due to the difficulty of detecting (eccentric) DWDs with present and planned electromagnetic observatories, LISA could provide unique dynamical identifications of these systems in globular clusters.
△ Less
Submitted 17 January, 2008; v1 submitted 29 May, 2007;
originally announced May 2007.
-
Neutron Stars and Black Holes in Star Clusters
Authors:
F. A. Rasio,
H. Baumgardt,
A. Corongiu,
F. D'Antona,
G. Fabbiano,
J. M. Fregeau,
K. Gebhardt,
C. O. Heinke,
P. Hut,
N. Ivanova,
T. J. Maccarone,
S. M. Ransom,
N. A. Webb
Abstract:
This article was co-authored by all invited speakers at the Joint Discussion on `Neutron Stars and Black Holes in Star Clusters,' which took place during the IAU General Assembly in Prague, Czech Republic, on August 17 and 18, 2006. Each section presents a short summary of recent developments in a key area of research, incorporating the main ideas expressed during the corresponding panel discuss…
▽ More
This article was co-authored by all invited speakers at the Joint Discussion on `Neutron Stars and Black Holes in Star Clusters,' which took place during the IAU General Assembly in Prague, Czech Republic, on August 17 and 18, 2006. Each section presents a short summary of recent developments in a key area of research, incorporating the main ideas expressed during the corresponding panel discussion at the meeting. Our meeting, which had close to 300 registered participants, was broadly aimed at the large community of astronomers around the world working on the formation and evolution of compact objects and interacting binary systems in dense star clusters, such as globular clusters and galactic nuclei. Great advances have occurred in this field during the past few years, including the introduction of fundamentally new theoretical paradigms for the formation and evolution of compact objects in binaries as well as countless new discoveries by astronomers that have challenged many accepted models. Some of the highlights include: a nearly complete census of all the millisecond pulsars in 47 Tuc; first detections of many new radio pulsars in other clusters, particularly Terzan 5; detailed studies of X-ray binary populations and their luminosity functions in many galaxies and extragalactic globular clusters; increasing evidence for intermediate-mass black holes in clusters and greatly improved theoretical understanding of their possible formation processes.
△ Less
Submitted 20 November, 2006;
originally announced November 2006.
-
Monte Carlo Simulations of Globular Cluster Evolution. IV. Direct Integration of Strong Interactions
Authors:
John M. Fregeau,
Frederic A. Rasio
Abstract:
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…
▽ More
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.
△ Less
Submitted 6 December, 2006; v1 submitted 11 August, 2006;
originally announced August 2006.
-
Observing IMBH-IMBH Binary Coalescences via Gravitational Radiation
Authors:
John M. Fregeau,
Shane L. Larson,
M. Coleman Miller,
Richard O'Shaughnessy,
Frederic A. Rasio
Abstract:
Recent numerical simulations have suggested the possibility of forming double intermediate mass black holes (IMBHs) via the collisional runaway scenario in young dense star clusters. The two IMBHs formed would exchange into a common binary shortly after their birth, and quickly inspiral and merge. Since space-borne gravitational wave (GW) observatories such as LISA will be able to see the late p…
▽ More
Recent numerical simulations have suggested the possibility of forming double intermediate mass black holes (IMBHs) via the collisional runaway scenario in young dense star clusters. The two IMBHs formed would exchange into a common binary shortly after their birth, and quickly inspiral and merge. Since space-borne gravitational wave (GW) observatories such as LISA will be able to see the late phases of their inspiral out to several Gpc, and LIGO will be able to see the merger and ringdown out to similar distances, they represent potentially significant GW sources. In this Letter we estimate the rate at which LISA and LIGO will see their inspiral and merger in young star clusters, and discuss the information that can be extracted from the observations. We find that LISA will likely see tens of IMBH--IMBH inspirals per year, while advanced LIGO could see ~10 merger and ringdown events per year, with both rates strongly dependent on the distribution of cluster masses and densities.
△ Less
Submitted 21 June, 2006; v1 submitted 30 May, 2006;
originally announced May 2006.
-
Massive Black Hole Binaries from Collisional Runaways
Authors:
M. Atakan Gürkan,
John M. Fregeau,
Frederic A. Rasio
Abstract:
Recent theoretical work has solidified the viability of the collisional runaway scenario in young dense star clusters for the formation of very massive stars (VMSs), which may be precursors to intermediate-mass black holes (IMBHs). We present first results from a numerical study of the collisional runaway process in dense star clusters containing primordial binaries. Stellar collisions during bi…
▽ More
Recent theoretical work has solidified the viability of the collisional runaway scenario in young dense star clusters for the formation of very massive stars (VMSs), which may be precursors to intermediate-mass black holes (IMBHs). We present first results from a numerical study of the collisional runaway process in dense star clusters containing primordial binaries. Stellar collisions during binary scattering encounters offer an alternate channel for runaway growth, somewhat independent of direct collisions between single stars. We find that clusters with binary fractions >~10% yield two VMSs via collisional runaways, presenting the exotic possibility of forming IMBH--IMBH binaries in star clusters. We discuss the implications for gravitational wave observations, and the impact on cluster structure.
△ Less
Submitted 11 February, 2006; v1 submitted 29 December, 2005;
originally announced December 2005.
-
Star Cluster Evolution with Primordial Binaries
Authors:
John M. Fregeau,
M. Atakan Gurkan,
Frederic A. Rasio
Abstract:
Observations and theoretical work suggest that globular clusters may be born with initially very large binary fractions. We present first results from our newly modified Monte-Carlo cluster evolution code, which treats binary interactions exactly via direct N-body integration. It is shown that binary scattering interactions generate significantly less energy than predicted by the recipes that ha…
▽ More
Observations and theoretical work suggest that globular clusters may be born with initially very large binary fractions. We present first results from our newly modified Monte-Carlo cluster evolution code, which treats binary interactions exactly via direct N-body integration. It is shown that binary scattering interactions generate significantly less energy than predicted by the recipes that have been used in the past to model them in approximate cluster evolution methods. The new result that the cores of globular clusters in the long-lived binary-burning phase are smaller than previously predicted weakens the agreement with observations, thus implying that more than simply stellar dynamics is at work in shaping the globular clusters we observe today.
△ Less
Submitted 1 December, 2005;
originally announced December 2005.
-
Dynamical Interactions of Planetary Systems in Dense Stellar Environments
Authors:
John M. Fregeau,
Sourav Chatterjee,
Frederic A. Rasio
Abstract:
We study dynamical interactions of star--planet binaries with other single stars. We derive analytical cross sections for all possible outcomes, and confirm them with numerical scattering experiments. We find that a wide mass ratio in the binary introduces a region in parameter space that is inaccessible to comparable-mass systems, in which the nature of the dynamical interaction is fundamentall…
▽ More
We study dynamical interactions of star--planet binaries with other single stars. We derive analytical cross sections for all possible outcomes, and confirm them with numerical scattering experiments. We find that a wide mass ratio in the binary introduces a region in parameter space that is inaccessible to comparable-mass systems, in which the nature of the dynamical interaction is fundamentally different from what has traditionally been considered in the literature on binary scattering. We study the properties of the planetary systems that result from the scattering interactions for all regions of parameter space, paying particular attention to the location of the "hard--soft" boundary. The structure of the parameter space turns out to be significantly richer than a simple statement of the location of the "hard--soft" boundary would imply. We consider the implications of our findings, calculating characteristic lifetimes for planetary systems in dense stellar environments, and applying the results to previous analytical studies, as well as past and future observations. Recognizing that the system PSR B1620-26 in the globular cluster M4 lies in the "new" region of parameter space, we perform a detailed analysis quantifying the likelihood of different scenarios in forming the system we see today.
△ Less
Submitted 23 November, 2005; v1 submitted 26 October, 2005;
originally announced October 2005.
-
Binary Mergers and Growth of Black Holes in Dense Star Clusters
Authors:
Ryan M. O'Leary,
Frederic A. Rasio,
John M. Fregeau,
Natalia Ivanova,
Richard O'Shaughnessy
Abstract:
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 ejectio…
▽ More
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.
△ Less
Submitted 1 October, 2005; v1 submitted 10 August, 2005;
originally announced August 2005.
-
The Evolution of Binary Fractions in Globular Clusters
Authors:
N. Ivanova,
K. Belczynski,
J. M. Fregeau,
F. A. Rasio
Abstract:
We study the evolution of binary stars in globular clusters using a new Monte Carlo approach combining a population synthesis code (StarTrack), and a simple treatment of dynamical interactions in the dense cluster core using a new tool for computing 3-body and 4-body interactions (Fewbody). We find that the combination of stellar evolution and dynamical interactions (binary-single and binary-bin…
▽ More
We study the evolution of binary stars in globular clusters using a new Monte Carlo approach combining a population synthesis code (StarTrack), and a simple treatment of dynamical interactions in the dense cluster core using a new tool for computing 3-body and 4-body interactions (Fewbody). We find that the combination of stellar evolution and dynamical interactions (binary-single and binary-binary) leads to a rapid depletion of the binary population in the cluster core. The maximum binary fraction today in the core of a typical dense cluster like 47 Tuc, assuming an initial binary fraction of 100%, is only about 5-10%. We show that this is in good agreement with recent HST observations of close binaries in the core of 47 Tuc, provided that a realistic distribution of binary periods is used to interpret the results. Our findings also have important consequences for the dynamical modeling of globular clusters, suggesting that ``realistic models'' should incorporate much larger initial binary fractions than has usually been done in the past.
△ Less
Submitted 7 January, 2005;
originally announced January 2005.
-
Binary Evolution and Neutron Stars in Globular Clusters
Authors:
Natalia Ivanova,
John M. Fregeau,
Frederic A. Rasio
Abstract:
We investigate the dynamical formation and evolution of binaries containing neutron stars in dense globular clusters. Our numerical simulations combine a simple Monte Carlo prescription for stellar dynamics, a sophisticated binary population synthesis code, and a small-N-body integrator for computing 3-body and 4-body interactions. Our results suggest that there is no ``retention problem,'' i.e.…
▽ More
We investigate the dynamical formation and evolution of binaries containing neutron stars in dense globular clusters. Our numerical simulations combine a simple Monte Carlo prescription for stellar dynamics, a sophisticated binary population synthesis code, and a small-N-body integrator for computing 3-body and 4-body interactions. Our results suggest that there is no ``retention problem,'' i.e., that, under standard assumptions, globular clusters can retain enough neutron stars to produce the observed numbers of millisecond pulsars. We also identify the dominant evolutionary and dynamical channels through which globular clusters produce their two main types of binary millisecond pulsars
△ Less
Submitted 19 May, 2004;
originally announced May 2004.
-
Stellar Collisions During Binary-Binary and Binary-Single Star Interactions
Authors:
J. M. Fregeau,
P. Cheung,
S. F. Portegies Zwart,
F. A. Rasio
Abstract:
Physical collisions between stars occur frequently in dense star clusters, either via close encounters between two single stars, or during strong dynamical interactions involving binary stars. Here we study stellar collisions that occur during binary-single and binary-binary interactions, by performing numerical scattering experiments. Our results include cross sections, branching ratios, and sa…
▽ More
Physical collisions between stars occur frequently in dense star clusters, either via close encounters between two single stars, or during strong dynamical interactions involving binary stars. Here we study stellar collisions that occur during binary-single and binary-binary interactions, by performing numerical scattering experiments. Our results include cross sections, branching ratios, and sample distributions of parameters for various outcomes. For interactions of hard binaries containing main-sequence stars, we find that the normalized cross section for at least one collision to occur (between any two of the four stars involved) is essentially unity, and that the probability of collisions involving more than two stars is significant. Hydrodynamic calculations have shown that the effective radius of a collision product can be 2-30 times larger than the normal main-sequence radius for a star of the same total mass. We study the effect of this expansion, and find that it increases the probability of further collisions considerably. We discuss these results in the context of recent observations of blue stragglers in globular clusters with masses exceeding twice the main-sequence turnoff mass. We also present Fewbody, a new, freely available numerical toolkit for simulating small-N gravitational dynamics that is particularly suited to performing scattering experiments.
△ Less
Submitted 20 May, 2004; v1 submitted 5 January, 2004;
originally announced January 2004.
-
On the Maximum Binary Fraction in Globular Cluster Cores
Authors:
N. Ivanova,
K. Belczynski,
J. M. Fregeau,
F. A. Rasio
Abstract:
We study the evolution of binary stars in globular clusters using a novel approach combining a state-of-the-art population synthesis code with a simple treatment of dynamical interactions in the dense cluster core. We find that the combination of stellar evolution and dynamical interactions (binary-single and binary-binary) leads to a rapid depletion of the binary population in the cluster core.…
▽ More
We study the evolution of binary stars in globular clusters using a novel approach combining a state-of-the-art population synthesis code with a simple treatment of dynamical interactions in the dense cluster core. We find that the combination of stellar evolution and dynamical interactions (binary-single and binary-binary) leads to a rapid depletion of the binary population in the cluster core. The maximum binary fraction today in the core of a typical dense cluster like 47 Tuc, assuming an initial binary fraction of 100%, is only about 5%. We show that this is in good agreement with recent HST observations of close binaries in the core of 47 Tuc, provided that a realistic distribution of binary periods is used to interpret the results. Our findings also have important consequences for the dynamical modeling of globular clusters, suggesting that ``realistic models'' should incorporate much larger initial binary fractions than has usually been done in the past.
△ Less
Submitted 18 December, 2003;
originally announced December 2003.
-
Accretion Onto Fast X-Ray Pulsars
Authors:
S. A. Rappaport,
J. M. Fregeau,
H. Spruit
Abstract:
The recent emergence of a new class of accretion-powered, transient, millisecond X-ray pulsars presents some difficulties for the conventional picture of accretion onto rapidly rotating magnetized neutron stars and their spin behavior during outbursts. In particular, it is unclear from the standard paradigm how these systems manage to accrete over such a wide range in dM/dt (i.e., >~ a factor of…
▽ More
The recent emergence of a new class of accretion-powered, transient, millisecond X-ray pulsars presents some difficulties for the conventional picture of accretion onto rapidly rotating magnetized neutron stars and their spin behavior during outbursts. In particular, it is unclear from the standard paradigm how these systems manage to accrete over such a wide range in dM/dt (i.e., >~ a factor of 50), and why the neutron stars exhibit a high rate of spindown in at least a number of cases. Following up on prior suggestions, we propose that `fast' X-ray pulsars can continue to accrete, and that their accretion disks terminate at approximately the corotation radius. We demonstrate the existence of such disk solutions by modifying the Shakura-Sunyaev equations with a simple magnetic torque prescription. The solutions are completely analytic, and have the same dependence on dM/dt and alpha (the viscosity parameter) as the original Shakura-Sunyaev solutions; but, the radial profiles can be considerably modified, depending on the degree of fastness. We apply these results to compute the torques expected during the outbursts of the transient millisecond pulsars, and find that we can explain the large spindown rates that are observed for quite plausible surface magnetic fields of ~10^9 G.
△ Less
Submitted 8 October, 2003;
originally announced October 2003.
-
Monte Carlo Simulations of Globular Cluster Evolution. III. Primordial Binary Interactions
Authors:
J. M. Fregeau,
M. A. Gurkan,
K. J. Joshi,
F. A. Rasio
Abstract:
We study the dynamical evolution of globular clusters using our 2D Monte Carlo code with the inclusion of primordial binary interactions for equal-mass stars. We use approximate analytical cross sections for energy generation from binary-binary and binary-single interactions. After a brief period of slight contraction or expansion of the core over the first few relaxation times, all clusters ent…
▽ More
We study the dynamical evolution of globular clusters using our 2D Monte Carlo code with the inclusion of primordial binary interactions for equal-mass stars. We use approximate analytical cross sections for energy generation from binary-binary and binary-single interactions. After a brief period of slight contraction or expansion of the core over the first few relaxation times, all clusters enter a much longer phase of stable "binary burning" lasting many tens of relaxation times. The structural parameters of our models during this phase match well those of most observed globular clusters. At the end of this phase, clusters that have survived tidal disruption undergo deep core collapse, followed by gravothermal oscillations. Our results clearly show that the presence of even a small fraction of binaries in a cluster is sufficient to support the core against collapse significantly beyond the normal core collapse time predicted without the presence of binaries. For tidally truncated systems, collapse is easily delayed sufficiently that the cluster will undergo complete tidal disruption before core collapse. As a first step toward the eventual goal of computing all interactions exactly using dynamical three- and four-body integration, we have incorporated an exact treatment of binary-single interactions in our code. We show that results using analytical cross sections are in good agreement with those using exact three-body integration, even for small binary fractions where binary-single interactions are energetically most important.
△ Less
Submitted 9 May, 2003; v1 submitted 25 January, 2003;
originally announced January 2003.
-
Mass Segregation in Globular Clusters
Authors:
J. M. Fregeau,
K. J. Joshi,
S. F. Portegies Zwart,
F. A. Rasio
Abstract:
We present the results of a new study of mass segregation in two-component star clusters, based on a large number of numerical N-body simulations using our recently developed dynamical Monte Carlo code. Specifically, we follow the dynamical evolution of clusters containing stars with individual masses m_1 as well as a tracer population of objects with individual masses m_2=μm_1, using N=10^5 tot…
▽ More
We present the results of a new study of mass segregation in two-component star clusters, based on a large number of numerical N-body simulations using our recently developed dynamical Monte Carlo code. Specifically, we follow the dynamical evolution of clusters containing stars with individual masses m_1 as well as a tracer population of objects with individual masses m_2=μm_1, using N=10^5 total stars. For heavy tracers, which could represent stellar remnants such as neutron stars or black holes in a globular cluster, we characterize in a variety of ways the tendency for these objects to concentrate in or near the cluster core. In agreement with simple theoretical arguments, we find that the characteristic time for this mass segregation process varies as 1/μ. For models with very light tracers (μ<~ 10^-2), which could represent free-floating planets or brown dwarfs, we find the expected depletion of light objects in the cluster core, but also sometimes a significant enhancement in the halo. Using these results we estimate the optical depth to gravitational microlensing by planetary mass objects or brown dwarfs in typical globular clusters. For some initial conditions, the optical depth in the halo due to very low-mass objects could be much greater than that of luminous stars. If we apply our results to M22, using the recent null detection of Sahu, Anderson, & King (2001), we find an upper limit of ~25% at the 63% confidence level for the current mass fraction of M22 in the form of very low-mass objects.
△ Less
Submitted 28 January, 2002; v1 submitted 2 November, 2001;
originally announced November 2001.
-
Binaries and Globular Cluster Dynamics
Authors:
Frederic A. Rasio,
John M. Fregeau,
Kriten J. Joshi
Abstract:
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 sour…
▽ More
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.
△ Less
Submitted 28 February, 2001;
originally announced March 2001.
-
Binary-induced collapse of a compact, collisionless cluster
Authors:
Matthew D. Duez,
Eric T. Engelhard,
John M. Fregeau,
Kevin M. Huffenberger,
Stuart L. Shapiro
Abstract:
We improve and extend Shapiro's model of a relativistic, compact object which is stable in isolation but is driven dynamically unstable by the tidal field of a binary companion. Our compact object consists of a dense swarm of test particles moving in randomly-oriented, initially circular, relativistic orbits about a nonrotating black hole. The binary companion is a distant, slowly inspiraling po…
▽ More
We improve and extend Shapiro's model of a relativistic, compact object which is stable in isolation but is driven dynamically unstable by the tidal field of a binary companion. Our compact object consists of a dense swarm of test particles moving in randomly-oriented, initially circular, relativistic orbits about a nonrotating black hole. The binary companion is a distant, slowly inspiraling point mass. The tidal field of the companion is treated as a small perturbation on the background Schwarzschild geometry near the hole; the resulting metric is determined by solving the perturbation equations of Regge and Wheeler and Zerilli in the quasi-static limit. The perturbed spacetime supports Bekenstein's conjecture that the horizon area of a near-equilibrium black hole is an adiabatic invariant. We follow the evolution of the system and confirm that gravitational collapse can be induced in a compact collisionless cluster by the tidal field of a binary companion.
△ Less
Submitted 28 June, 1999; v1 submitted 5 April, 1999;
originally announced April 1999.
