-
Modelling Neutron Star-Black Hole Binaries: Future Pulsar Surveys and Gravitational Wave Detectors
Authors:
Debatri Chattopadhyay,
Simon Stevenson,
Jarrod R. Hurley,
Matthew Bailes,
Floor Broekgaarden
Abstract:
Binaries comprised of a neutron star (NS) and a black hole (BH) have so far eluded observations as pulsars and with gravitational waves (GWs). We model the formation and evolution of these NS+BH binaries - including pulsar evolution - using the binary population synthesis code COMPAS. We predict the presence of a total of 50-2000 binaries containing a pulsar and a BH (PSR+BHs) in the Galactic fiel…
▽ More
Binaries comprised of a neutron star (NS) and a black hole (BH) have so far eluded observations as pulsars and with gravitational waves (GWs). We model the formation and evolution of these NS+BH binaries - including pulsar evolution - using the binary population synthesis code COMPAS. We predict the presence of a total of 50-2000 binaries containing a pulsar and a BH (PSR+BHs) in the Galactic field. We find the population observable by the next generation of radio telescopes, represented by the SKA and MeerKAT, current (LIGO/Virgo) and future (LISA) GW detectors. We conclude that the SKA will observe 1-80 PSR+BHs, with 0-4 binaries containing millisecond pulsars. MeerKAT is expected to observe 0-40 PSR+BH systems. Future radio detections of NS+BHs will constrain uncertain binary evolution processes such as BH natal kicks. We show that systems in which the NS formed first (NSBH) can be distinguished from those where the BH formed first (BHNS) by their pulsar and binary properties. We find 40% of the LIGO/Virgo observed NS+BHs from a Milky-Way like field population will have a chirp mass $\geq 3.0$ M$_\odot$. We estimate the spin distributions of NS+BHs with two models for the spins of BHs. The remnants of BHNS mergers will have a spin of $\sim$0.4, whilst NSBH merger remnants can have a spin of $\sim$0.6 or $\sim$0.9 depending on the model for BH spins. We estimate that approximately 25-1400 PSR+BHs will be radio alive whilst emitting GWs in the LISA frequency band, raising the possibility of joint observation by the SKA and LISA.
△ Less
Submitted 19 April, 2021; v1 submitted 26 November, 2020;
originally announced November 2020.
-
Modelling Double Neutron Stars: Radio and Gravitational Waves
Authors:
Debatri Chattopadhyay,
Simon Stevenson,
Jarrod R. Hurley,
Luca J. Rossi,
Chris Flynn
Abstract:
We have implemented prescriptions for modelling pulsars in the rapid binary population synthesis code COMPAS. We perform a detailed analysis of the double neutron star (DNS) population, accounting for radio survey selection effects. The surface magnetic field decay timescale (${\sim}1000$\,Myr) and mass scale (${\sim}0.02$\,M$_\odot$) are the dominant uncertainties in our model. Mass accretion dur…
▽ More
We have implemented prescriptions for modelling pulsars in the rapid binary population synthesis code COMPAS. We perform a detailed analysis of the double neutron star (DNS) population, accounting for radio survey selection effects. The surface magnetic field decay timescale (${\sim}1000$\,Myr) and mass scale (${\sim}0.02$\,M$_\odot$) are the dominant uncertainties in our model. Mass accretion during common envelope evolution plays a non-trivial role in recycling pulsars. We find a best-fit model that is in broad agreement with the observed Galactic DNS population. Though the pulsar parameters (period and period derivative) are strongly biased by radio selection effects, the observed orbital parameters (orbital period and eccentricity) closely represent the intrinsic distributions. The number of radio observable DNSs in the Milky Way at present is about 2500 in our model, corresponding to approximately 10\% of the predicted total number of DNSs in the galaxy. Using our model calibrated to the Galactic DNS population, we make predictions for DNS mergers observed in gravitational waves. The DNS chirp mass distribution varies from 1.1M$_\odot$ to 2.1M$_\odot$ and the median is found to be 1.14\,M$_\mathrm{\odot}$. The expected effective spin $χ_\mathrm{eff}$ for isolated DNSs is $\lesssim$0.03 from our model. We predict that 34\% of the current Galactic isolated DNSs will merge within a Hubble time, and have a median total mass of 2.7\,M$_\mathrm{\odot}$. Finally, we discuss implications for fast radio bursts and post-merger remnant gravitational-waves.
△ Less
Submitted 25 March, 2020; v1 submitted 5 December, 2019;
originally announced December 2019.
-
BSE versus StarTrack: implementations of new wind, remnant-formation, and natal-kick schemes in NBODY7 and their astrophysical consequences
Authors:
Sambaran Banerjee,
Krzysztof Belczynski,
Christopher L. Fryer,
Peter Berczik,
Jarrod R. Hurley,
Rainer Spurzem,
Long Wang
Abstract:
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…
▽ More
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.
△ Less
Submitted 9 April, 2020; v1 submitted 20 February, 2019;
originally announced February 2019.
-
No cataclysmic variables missing: higher merger rate brings into agreement observed and predicted space densities
Authors:
Diogo Belloni,
Matthias R. Schreiber,
Mónica Zorotovic,
Krystian Iłkiewicz,
Jarrod R. Hurley,
Mirek Giersz,
Felipe Lagos
Abstract:
The predicted and observed space density of cataclysmic variables (CVs) have been for a long time discrepant by at least an order of magnitude. The standard model of CV evolution predicts that the vast majority of CVs should be period bouncers, whose space density has been recently measured to be $ρ\lesssim 2 \times 10^{-5}$ pc$^{-3}$. We performed population synthesis of CVs using an updated vers…
▽ More
The predicted and observed space density of cataclysmic variables (CVs) have been for a long time discrepant by at least an order of magnitude. The standard model of CV evolution predicts that the vast majority of CVs should be period bouncers, whose space density has been recently measured to be $ρ\lesssim 2 \times 10^{-5}$ pc$^{-3}$. We performed population synthesis of CVs using an updated version of the Binary Stellar Evolution (BSE) code for single and binary star evolution. We find that the recently suggested empirical prescription of consequential angular momentum loss (CAML) brings into agreement predicted and observed space densities of CVs and period bouncers. To progress with our understanding of CV evolution it is crucial to understand the physical mechanism behind empirical CAML. Our changes to the BSE code are also provided in details, which will allow the community to accurately model mass transfer in interacting binaries in which degenerate objects accrete from low-mass main-sequence donor stars.
△ Less
Submitted 14 June, 2018;
originally announced June 2018.
-
Mass Evaporation Rate of Globular Clusters in a Strong Tidal Field
Authors:
Juan P. Madrid,
Nathan W. C. Leigh,
Jarrod R. Hurley,
Mirek Giersz
Abstract:
The mass evaporation rate of globular clusters evolving in a strong Galactic tidal field is derived through the analysis of large, multi-mass $N$-body simulations. For comparison, we also study the same evaporation rates using MOCCA Monte Carlo models for globular cluster evolution. Our results show that the mass evaporation rate is a dynamical value; that is, far from a constant single number fou…
▽ More
The mass evaporation rate of globular clusters evolving in a strong Galactic tidal field is derived through the analysis of large, multi-mass $N$-body simulations. For comparison, we also study the same evaporation rates using MOCCA Monte Carlo models for globular cluster evolution. Our results show that the mass evaporation rate is a dynamical value; that is, far from a constant single number found in earlier analytical work and commonly used in the literature. Moreover, the evaporation rate derived with these simulations is higher than values previously published. These models also show that the value of the mass evaporation rate depends on the strength of the tidal field. We give an analytical estimate of the mass evaporation rate as a function of time and galactocentric distance xi(R_{GC},t). Upon extrapolating this formula to smaller R_{GC} values, our results provide tentative evidence for a very high xi value at small R_{GC}. Our results suggest that the corresponding mass loss in the inner Galactic potential could be high and it should be accounted for when star clusters pass within it. This has direct relevance to nuclear cluster formation/growth via the infall of globular clusters through dynamical friction. As an illustrative example, we estimate how the evaporation rate increases for a ~10^5 solar masses globular cluster that decays through dynamical friction into the galactic centre. We discuss the findings of this work in relation to the formation of nuclear star clusters by inspiraling globular clusters.
△ Less
Submitted 20 June, 2017;
originally announced June 2017.
-
Evolution of star cluster systems in isolated galaxies: first results from direct $N$-body simulations
Authors:
L. J. Rossi,
K. Bekki,
J. R. Hurley
Abstract:
The evolution of star clusters is largely affected by the tidal field generated by the host galaxy. It is thus in principle expected that under the assumption of an "universal" initial cluster mass function the properties of the evolved present-day mass function of star cluster systems should show a dependency on the properties of the galactic environment in which they evolve. To explore this expe…
▽ More
The evolution of star clusters is largely affected by the tidal field generated by the host galaxy. It is thus in principle expected that under the assumption of an "universal" initial cluster mass function the properties of the evolved present-day mass function of star cluster systems should show a dependency on the properties of the galactic environment in which they evolve. To explore this expectation a sophisticated model of the tidal field is required in order to study the evolution of star cluster systems in realistic galaxies. Along these lines, in the present work we first describe a method developed for coupling $N$-body simulations of galaxies and star clusters. We then generate a database of galaxy models along the Hubble sequence and calibrate evolutionary equations to the results of direct $N$-body simulations of star clusters in order to predict the clusters' mass evolution as function of the galactic environment. We finally apply our methods to explore the properties of evolved "universal" initial cluster mass functions and any dependence on the host galaxy morphology and mass distribution. The preliminary results show that an initial power-law distribution of the masses "universally" evolves into a log-normal distribution, with the properties correlated with the stellar mass and stellar mass density density of the host galaxy.
△ Less
Submitted 8 August, 2016;
originally announced August 2016.
-
A dynamical gravitational wave source in a dense cluster
Authors:
Jarrod R. Hurley,
Anna C. Sippel,
Christopher A. Tout,
Sverre J. Aarseth
Abstract:
Making use of a new N-body model to describe the evolution of a moderate-size globular cluster we investigate the characteristics of the population of black holes within such a cluster. This model reaches core-collapse and achieves a peak central density typical of the dense globular clusters of the Milky Way. Within this high-density environment we see direct confirmation of the merging of two st…
▽ More
Making use of a new N-body model to describe the evolution of a moderate-size globular cluster we investigate the characteristics of the population of black holes within such a cluster. This model reaches core-collapse and achieves a peak central density typical of the dense globular clusters of the Milky Way. Within this high-density environment we see direct confirmation of the merging of two stellar remnant black-holes in a dynamically-formed binary, a gravitational wave source. We describe how the formation, evolution and ultimate ejection/destruction of binary systems containing black holes impacts the evolution of the cluster core. Also, through comparison with previous models of lower density, we show that the period distribution of black hole binaries formed through dynamical interactions in this high-density model favours the production of gravitational wave sources. We confirm that the number of black holes remaining in a star cluster at late times and the characteristics of the binary black hole population depend on the nature of the star cluster, critically on the number density of stars and by extension the relaxation timescale.
△ Less
Submitted 3 July, 2016;
originally announced July 2016.
-
The Globular Cluster System of NGC 6822
Authors:
J. Veljanoski,
A. M. N. Ferguson,
A. D. Mackey,
A. P. Huxor,
J. R. Hurley,
E. J. Bernard,
P. Cote,
M. J. Irwin,
N. F. Martin,
W. S. Burgett,
K. C. Chambers,
H. Flewelling,
R. Kudritzki,
C. Waters
Abstract:
We present a comprehensive analysis of the globular cluster (GC) system of the Local Group dwarf irregular galaxy NGC 6822. Our study is based on homogeneous optical and near-IR photometry, as well as long-slit spectroscopic observations which are used to determine new radial velocities for 6 GCs, two of which had no previous spectroscopic information. We construct optical-near IR colour-colour di…
▽ More
We present a comprehensive analysis of the globular cluster (GC) system of the Local Group dwarf irregular galaxy NGC 6822. Our study is based on homogeneous optical and near-IR photometry, as well as long-slit spectroscopic observations which are used to determine new radial velocities for 6 GCs, two of which had no previous spectroscopic information. We construct optical-near IR colour-colour diagrams and through comparison to simple stellar population models infer that the GCs have old ages consistent with being 9 Gyr or older, while their metallicities are in the range between -1.6 < [Fe/H] < -0.4. We conduct a kinematic analysis of the GC population and find tentative evidence for weak net rotation of the GC system, in the same sense as that exhibited by the underlying spheroid. The most likely amplitude of rotation is ~10 km/s, approximately half the magnitude of the observed velocity dispersion. Finally, we use the GCs to estimate the dynamical mass of NGC 6822 within 11 kpc and we formally find it to be in the range between (3-4)10^9 Msun. This implies an overall mass-to-light ratio in the range of ~ 30-40 and indicates that NGC 6822 is highly dark matter dominated. The mass and the corresponding mass-to-light ratio estimates are affected by various additional systematic effects due to limitations of the data and the model that are not necessary reflected in the formal uncertainties.
△ Less
Submitted 2 June, 2015;
originally announced June 2015.
-
Different Dynamical Ages for the Two Young and Coeval LMC Star Clusters NGC 1805 and NGC 1818 Imprinted on Their Binary Populations
Authors:
Aaron M. Geller,
Richard de Grijs,
Chengyuan Li,
Jarrod R. Hurley
Abstract:
The two Large Magellanic Cloud star clusters NGC 1805 and NGC 1818 are approximately the same chronological age (~30 Myr), but show different radial trends in binary frequency. The F-type stars (1.3 - 2.2 MSun) in NGC 1818 have a binary frequency that decreases towards the core, while the binary frequency for stars of similar mass in NGC 1805 is flat with radius, or perhaps bimodal (with a peak in…
▽ More
The two Large Magellanic Cloud star clusters NGC 1805 and NGC 1818 are approximately the same chronological age (~30 Myr), but show different radial trends in binary frequency. The F-type stars (1.3 - 2.2 MSun) in NGC 1818 have a binary frequency that decreases towards the core, while the binary frequency for stars of similar mass in NGC 1805 is flat with radius, or perhaps bimodal (with a peak in the core). We show here, through detailed N-body modeling, that both clusters could have formed with the same primordial binary frequency and with binary orbital elements and masses drawn from the same distributions (defined from observations of open clusters and the field of our Galaxy). The observed radial trends in binary frequency for both clusters are best matched with models that have initial substructure. Furthermore, both clusters may be evolving along a very similar dynamical sequence, with the key difference that NGC 1805 is dynamically older than NGC 1818. The F-type binaries in NGC 1818 still show evidence of an initial period of rapid dynamical disruptions (which occur preferentially in the core), while NGC 1805 has already begun to recover a higher core binary frequency, owing to mass segregation (which will eventually produce a distribution in binary frequency that rises only towards the core, as is observed in old Milky Way star clusters). This recovery rate increases for higher-mass binaries, and therefore even at one age in one cluster, we predict a similar dynamical sequence in the radial distribution of the binary frequency as a function of binary primary mass.
△ Less
Submitted 17 March, 2015;
originally announced March 2015.
-
Merging Binary Stars and the magnetic white dwarfs
Authors:
Gordon P. Briggs,
Lilia Ferrario,
Christopher A. Tout,
Dayal T. Wickramasinghe,
Jarrod R. Hurley
Abstract:
A magnetic dynamo driven by differential rotation generated when stars merge can explain strong fields in certain classes of magnetic stars, including the high field magnetic white dwarfs (HFMWDs). In their case the site of the differential rotation has been variously proposed to be within a common envelope, the massive hot outer regions of a merged degenerate core or an accretion disc formed by a…
▽ More
A magnetic dynamo driven by differential rotation generated when stars merge can explain strong fields in certain classes of magnetic stars, including the high field magnetic white dwarfs (HFMWDs). In their case the site of the differential rotation has been variously proposed to be within a common envelope, the massive hot outer regions of a merged degenerate core or an accretion disc formed by a tidally disrupted companion that is subsequently incorporated into a degenerate core. We synthesize a population of binary systems to investigate the stellar merging hypothesis for observed single HFMWDs. Our calculations provide mass distribution and the fractions of white dwarfs that merge during a common envelope phase or as double degenerate systems in a post common envelope phase. We vary the common envelope efficiency parameter alpha and compare with observations. We find that this hypothesis can explain both the observed incidence of magnetism and the mass distribution of HFMWDs for a wide range of alpha. In this model, the majority of the HFMWDs are of the Carbon Oxygen type and merge within a common envelope. Less than about a quarter of a per cent of HFMWDs originate from double degenerate stars that merge after common envelope evolution and these populate the high-mass tail of the HFMWD mass distribution.
△ Less
Submitted 17 December, 2014;
originally announced December 2014.
-
Tracking Cluster Debris (TraCD)- I. Dissolution of clusters and searching the solar cradle
Authors:
Guido R. I. Moyano Loyola,
Chris Flynn,
Jarrod R. Hurley,
Brad K. Gibson
Abstract:
The capability to reconstruct dissolved stellar systems in dynamical and chemical space is a key factor in improving our understanding of the evolution of the Milky Way. Here we concentrate on the dynamical aspect and given that a significant portion of the stars in the Milky Way have been born in stellar associations or clusters that have lived a few Myr up to several Gyr, we further restrict our…
▽ More
The capability to reconstruct dissolved stellar systems in dynamical and chemical space is a key factor in improving our understanding of the evolution of the Milky Way. Here we concentrate on the dynamical aspect and given that a significant portion of the stars in the Milky Way have been born in stellar associations or clusters that have lived a few Myr up to several Gyr, we further restrict our attention to the evolution of star clusters. We have carried out our simulations in two steps: (1) we create a simulation of dissolution and mixing processes which yields a close fit to the present-day Milky Way dynamics and (2) we have evolved three sets of stellar clusters with masses of 400, 1000 and 15000 M$_{\odot}$ to dissolution. The birth location of these sets was 4, 6, 8 and 10 kpc for the 400 and 1000 M$_{\odot}$ clusters and 4, 6, 8, 10 and 12 kpc for the 15000 M$_{\odot}$. We have focused our efforts on studying the state of the escapers from these clusters after 4.5 Gyr of evolution with particular attention to stars that reach the Solar annulus, i.e. 7.5 $\le$ R$_{\rm{gc}}$ $\le$ 8.5 kpc. We give results for Solar twins and siblings over a wide range of radii and cluster masses for two dissolution mechanisms. From kinematics alone, we conclude that the Sun was $\sim$50 per cent more likely to have been born near its current Galactocentric radius, rather than have migrated (radially) $\sim$2 kpc since birth. We conclude our analysis by calculating magnitudes and colours of our single stars for comparison with the samples that the $Gaia$, $Gaia$-ESO and GALAH-AAO surveys will obtain. In terms of reconstructing dissolved star clusters we find that on short time-scales we cannot rely on kinematic evolution alone and thus it will be necessary to extend our study to include information on chemical space.
△ Less
Submitted 30 March, 2015; v1 submitted 26 November, 2014;
originally announced November 2014.
-
Dynamical Interactions Make Hot Jupiters in Open Star Clusters
Authors:
Michael M. Shara,
Jarrod R. Hurley,
Rosemary A. Mardling
Abstract:
Explaining the origin and evolution of exoplanetary "hot Jupiters" remains a significant challenge. One possible mechanism for their production is planet-planet interactions, which produces hot Jupiters from planets born far from their host stars but near their dynamical stability limits. In the much more likely case of planets born far from their dynamical stability limits, can hot Jupiters can b…
▽ More
Explaining the origin and evolution of exoplanetary "hot Jupiters" remains a significant challenge. One possible mechanism for their production is planet-planet interactions, which produces hot Jupiters from planets born far from their host stars but near their dynamical stability limits. In the much more likely case of planets born far from their dynamical stability limits, can hot Jupiters can be formed in star clusters? Our N-body simulations of planetary systems inside star clusters answer this question in the affirmative, and show that hot Jupiter formation is not a rare event. We detail three case studies of the dynamics-induced births of hot Jupiters on highly eccentric orbits that can only occur inside star clusters. The hot Jupiters' orbits bear remarkable similarities to those of some of the most extreme exoplanets known: HAT-P-32 b, HAT-P-2 b, HD 80606 b and GJ 876 d. If stellar perturbations formed these hot Jupiters then our simulations predict that these very hot, inner planets are often accompanied by much more distant gas giants in highly eccentric orbits.
△ Less
Submitted 3 November, 2015; v1 submitted 25 November, 2014;
originally announced November 2014.
-
Reconstructing the initial mass function of disc-bulge Galactic globular clusters from N-body simulations
Authors:
L. J. Rossi,
J. R. Hurley
Abstract:
We propose an evolutionary model to describe the dynamical evolution of star cluster systems in tidal fields, in which we calibrated the parametric equations defining the model by running direct N-body simulations of star clusters with a wide range of initial masses and set of orbital parameters, living within the external tidal field generated by a disc-like galaxy. We derived a new method to sol…
▽ More
We propose an evolutionary model to describe the dynamical evolution of star cluster systems in tidal fields, in which we calibrated the parametric equations defining the model by running direct N-body simulations of star clusters with a wide range of initial masses and set of orbital parameters, living within the external tidal field generated by a disc-like galaxy. We derived a new method to solve numerically the evolutionary equations, allowing us to infer constraints on the mass of a star cluster from its age, present-day mass, orbital parameters and external gravitational potential. The result has been applied to the metal-rich subsample of Galactic globular clusters, being a good representation of a disc-bulge population. We reconstructed the initial mass function of these objects from the present-day mass function, finding that a lognormal distribution is well preserved during the evolution of the globular cluster system. The evolution of a power-law initial mass function has been evaluated, confirming that it transforms into a lognormal distribution of the cluster masses within an Hubble time. Our results are consistent with a formation scenario in which metal-rich Galactic globular clusters formed from giant molecular clouds in high-pressure regions during the early phases of the evolution of the Galactic disc and bulge.
△ Less
Submitted 7 October, 2014;
originally announced October 2014.
-
The Effects of Orbital Inclination on the Scale Size and Evolution of Tidally Filling Star Clusters
Authors:
Jeremy J. Webb,
Alison Sills,
William E. Harris,
Jarrod R. Hurley
Abstract:
We have performed N-body simulations of tidally filling star clusters with a range of orbits in a Milky Way-like potential to study the effects of orbital inclination and eccentricity on their structure and evolution. At small galactocentric distances Rgc, a non-zero inclination results in increased mass loss rates. Tidal heating and disk shocking, the latter sometimes consisting of two shocking e…
▽ More
We have performed N-body simulations of tidally filling star clusters with a range of orbits in a Milky Way-like potential to study the effects of orbital inclination and eccentricity on their structure and evolution. At small galactocentric distances Rgc, a non-zero inclination results in increased mass loss rates. Tidal heating and disk shocking, the latter sometimes consisting of two shocking events as the cluster moves towards and away from the disk, help remove stars from the cluster. Clusters with inclined orbits at large Rgc have decreased mass loss rates than the non-inclined case, since the strength the disk potential decreases with Rgc. Clusters with inclined and eccentric orbits experience increased tidal heating due to a constantly changing potential, weaker disk shocks since passages occur at higher Rgc, and an additional tidal shock at perigalacticon. The effects of orbital inclination decrease with orbital eccentricity, as a highly eccentric cluster spends the majority of its lifetime at a large Rgc. The limiting radii of clusters with inclined orbits are best represented by the rt of the cluster when at its maximum height above the disk, where the cluster spends the majority of its lifetime and the rate of change in rt is a minimum. Conversely, the effective radius is independent of inclination in all cases.
△ Less
Submitted 2 September, 2014;
originally announced September 2014.
-
The Effect of Orbital Eccentricity on the Dynamical Evolution of Star Clusters
Authors:
Jeremy J. Webb,
Nathan Leigh,
Alison Sills,
William E. Harris,
Jarrod R. Hurley
Abstract:
We use N-body simulations to explore the influence of orbital eccentricity on the dynamical evolution of star clusters. Specifically we compare the mass loss rate, velocity dispersion, relaxation time, and the mass function of star clusters on circular and eccentric orbits. For a given perigalactic distance, increasing orbital eccentricity slows the dynamical evolution of a cluster due to a weaker…
▽ More
We use N-body simulations to explore the influence of orbital eccentricity on the dynamical evolution of star clusters. Specifically we compare the mass loss rate, velocity dispersion, relaxation time, and the mass function of star clusters on circular and eccentric orbits. For a given perigalactic distance, increasing orbital eccentricity slows the dynamical evolution of a cluster due to a weaker mean tidal field. However, we find that perigalactic passes and tidal heating due to an eccentric orbit can partially compensate for the decreased mean tidal field by energizing stars to higher velocities and stripping additional stars from the cluster, accelerating the relaxation process. We find that the corresponding circular orbit which best describes the evolution of a cluster on an eccentric orbit is much less than its semi-major axis or time averaged galactocentric distance. Since clusters spend the majority of their lifetimes near apogalacticon, the properties of clusters which appear very dynamically evolved for a given galactocentric distance can be explained by an eccentric orbit. Additionally we find that the evolution of the slope of the mass function within the core radius is roughly orbit-independent, so it could place additional constraints on the initial mass and initial size of globular clusters with solved orbits. We use our results to demonstrate how the orbit of Milky Way globular clusters can be constrained given standard observable parameters like galactocentric distance and the slope of the mass function. We then place constraints on the unsolved orbits of NGC 1261,NGC 6352, NGC 6496, and NGC 6304 based on their positions and mass functions.
△ Less
Submitted 4 July, 2014; v1 submitted 13 May, 2014;
originally announced May 2014.
-
The impact of galaxy geometry and mass evolution on the survival of star clusters
Authors:
Juan P. Madrid,
Jarrod R. Hurley,
Marie Martig
Abstract:
Direct N-body simulations of globular clusters in a realistic Milky Way-like potential are carried out using the code NBODY6 to determine the impact of the host galaxy disk mass and geometry on the survival of star clusters. A relationship between disk mass and star cluster dissolution timescale is derived. These N-body models show that doubling the mass of the disk from 5x10^10 solar masses to 10…
▽ More
Direct N-body simulations of globular clusters in a realistic Milky Way-like potential are carried out using the code NBODY6 to determine the impact of the host galaxy disk mass and geometry on the survival of star clusters. A relationship between disk mass and star cluster dissolution timescale is derived. These N-body models show that doubling the mass of the disk from 5x10^10 solar masses to 10x10^10 solar masses halves the dissolution time of a satellite star cluster orbiting the host galaxy at 6 kpc from the galactic center. Different geometries in a disk of identical mass can determine either the survival or dissolution of a star cluster orbiting within the inner 6 kpc of the galactic center. Furthermore, disk geometry has measurable effects on the mass loss of star clusters up to 15 kpc from the galactic center. N-body simulations performed with a fine output time step show that at each disk crossing the outer layers of star clusters experience an increase in velocity dispersion of ~5% of the average velocity dispersion in the outer section of star clusters. This leads to an enhancement of mass-loss -- a clearly discernable effect of disk shocking. By running models with different inclinations we determine that star clusters with an orbit perpendicular to the Galactic plane have larger mass loss rates than both clusters evolving in the Galactic plane or in an inclined orbit.
△ Less
Submitted 10 February, 2014;
originally announced February 2014.
-
Consequences of Dynamical Disruption and Mass Segregation for the Binary Frequencies of Star Clusters
Authors:
Aaron M. Geller,
Richard de Grijs,
Chengyuan Li,
Jarrod R. Hurley
Abstract:
The massive (13,000-26,000 Msolar), young (15-30 Myr) Large Magellanic Cloud star cluster NGC 1818 reveals an unexpected increasing binary frequency with radius for F-type stars (1.3-2.2 Msolar). This is in contrast to many older star clusters that show a decreasing binary frequency with radius. We study this phenomenon with sophisticated N-body modeling, exploring a range of initial conditions, f…
▽ More
The massive (13,000-26,000 Msolar), young (15-30 Myr) Large Magellanic Cloud star cluster NGC 1818 reveals an unexpected increasing binary frequency with radius for F-type stars (1.3-2.2 Msolar). This is in contrast to many older star clusters that show a decreasing binary frequency with radius. We study this phenomenon with sophisticated N-body modeling, exploring a range of initial conditions, from smooth virialized density distributions to highly substructured and collapsing configurations. We find that many of these models can reproduce the cluster's observed properties, although with a modest preference for substructured initial conditions. Our models produce the observed radial trend in binary frequency through disruption of soft binaries (with semi-major axes, a > 3000 AU), on approximately a crossing time (~ 5.4 Myr), preferentially in the cluster core. Mass segregation subsequently causes the binaries to sink towards the core. After roughly one initial half-mass relaxation time (t_rh(0) ~ 340 Myr) the radial binary frequency distribution becomes bimodal, the innermost binaries having already segregated towards the core, leaving a minimum in the radial binary frequency distribution that marches outwards with time. After 4-6 t_rh(0), the rising distribution in the halo disappears, leaving a radial distribution that rises only towards the core. Thus, both a radial binary frequency distribution that falls towards the core (as observed for NGC 1818) and one that rises towards the core (as for older star clusters) can arise naturally from the same evolutionary sequence owing to binary disruption and mass segregation in rich star clusters.
△ Less
Submitted 3 October, 2013;
originally announced October 2013.
-
An Age Difference of 2 Gyr between a Metal-Rich and a Metal-Poor Globular Cluster
Authors:
B. M. S. Hansen,
J. S. Kalirai,
J. Anderson,
A. Dotter,
H. B. Richer,
R. M. Rich,
M. M. Shara,
G. G. Fahlman,
J. R. Hurley,
I. R. King,
D. Reitzel,
P. B. Stetson
Abstract:
Globular clusters trace the formation history of the spheroidal components of both our Galaxy and others, which represent the bulk of star formation over the history of the universe. They also exhibit a range of metallicities, with metal-poor clusters dominating the stellar halo of the Galaxy, and higher metallicity clusters found within the inner Galaxy, associated with the stellar bulge, or the…
▽ More
Globular clusters trace the formation history of the spheroidal components of both our Galaxy and others, which represent the bulk of star formation over the history of the universe. They also exhibit a range of metallicities, with metal-poor clusters dominating the stellar halo of the Galaxy, and higher metallicity clusters found within the inner Galaxy, associated with the stellar bulge, or the thick disk. Age differences between these clusters can indicate the sequence in which the components of the Galaxy formed, and in particular which clusters were formed outside the Galaxy and later swallowed along with their original host galaxies, and which were formed in situ. Here we present an age determination of the metal-rich globular cluster 47 Tucanae by fitting the properties of the cluster white dwarf population, which implies an absolute age of 9.9 (0.7) Gyr at 95% confidence. This is about 2.0 Gyr younger than inferred for the metal-poor cluster NGC 6397 from the same models, and provides quantitative evidence that metal-rich clusters like 47 Tucanae formed later than the metal-poor halo clusters like NGC 6397.
△ Less
Submitted 31 July, 2013;
originally announced August 2013.
-
Stars on the run: escaping from stellar clusters
Authors:
G. R. I. Moyano Loyola,
J. R. Hurley
Abstract:
A significant proportion of Milky Way stars are born in stellar clusters, which dissolve over time so that the members become part of the disc and halo populations of the Galaxy. In the present work we will assume that these young stellar clusters live mainly within the disc of the Galaxy and that they can have primordial binary percentages ranging from 0% to as high as 70%. We have evolved models…
▽ More
A significant proportion of Milky Way stars are born in stellar clusters, which dissolve over time so that the members become part of the disc and halo populations of the Galaxy. In the present work we will assume that these young stellar clusters live mainly within the disc of the Galaxy and that they can have primordial binary percentages ranging from 0% to as high as 70%. We have evolved models of such clusters to an age of 4 Gyr through N-body simulations, paying attention to the stars and binaries that escape in the process. We have quantified the contribution of these escaping stars to the Galaxy population by analysing their escape velocity and evolutionary stage at the moment of escape. In this way we could analyse the mechanisms that produced these escapers, whether evaporation through weak two- body encounters, energetic close encounters or stellar evolution events, e.g. supernovae. In our models we found that the percentage of primordial binaries in a star cluster does not produce significant variations in the velocities of the stars that escape in the velocity range of 0-20 km/s. However, in the high-velocity 20-100 km/s range the number of escapers increased markedly as the primordial binary percentage increased. We could also infer that dissolving stellar clusters such as those that we have modelled can populate the Galactic halo with giant stars for which the progenitors were stars of up to 2.4 Msun. Furthermore, choices made for the velocity kicks of remnants do influence the production of hyper-velocity stars - and to a lesser extent stars in the high-velocity range - but once again the difference for the 99% of stars in the 0-20 km/s range is not significant.
△ Less
Submitted 4 July, 2013;
originally announced July 2013.
-
The Influence of Orbital Eccentricity on Tidal Radii of Star Clusters
Authors:
Jeremy J. Webb,
William E. Harris,
Alison Sills,
Jarrod R. Hurley
Abstract:
We have performed N-body simulations of star clusters orbiting in a spherically symmetric smooth galactic potential. The model clusters cover a range of initial half-mass radii and orbital eccentricities in order to test the historical assumption that the tidal radius of a cluster is imposed at perigalacticon. The traditional assumption for globular clusters is that since the internal relaxation t…
▽ More
We have performed N-body simulations of star clusters orbiting in a spherically symmetric smooth galactic potential. The model clusters cover a range of initial half-mass radii and orbital eccentricities in order to test the historical assumption that the tidal radius of a cluster is imposed at perigalacticon. The traditional assumption for globular clusters is that since the internal relaxation time is larger than its orbital period, the cluster is tidally stripped at perigalacticon. Instead, our simulations show that a cluster with an eccentric orbit does not need to fully relax in order to expand. After a perigalactic pass, a cluster re-captures previously unbound stars, and the tidal shock at perigalacticon has the effect of energizing inner region stars to larger orbits. Therefore, instead of the limiting radius being imposed at perigalacticon, it more nearly traces the instantaneous tidal radius of the cluster at any point in the orbit. We present a numerical correction factor to theoretical tidal radii calculated at perigalacticon which takes into consideration both the orbital eccentricity and current orbital phase of the cluster.
△ Less
Submitted 3 January, 2013;
originally announced January 2013.
-
Multiple stellar-mass black holes in globular clusters: theoretical confirmation
Authors:
Anna C. Sippel,
Jarrod R. Hurley
Abstract:
While tens or hundreds of stellar-remnant black holes are expected to form in globular star clusters, it is still unclear how many of those will be retained upon formation, and how many will be ejected through subsequent dynamical interactions. No such black holes have been found in any Milky Way globular cluster until the recent discovery of stellar-mass black holes in the globular cluster M22 (N…
▽ More
While tens or hundreds of stellar-remnant black holes are expected to form in globular star clusters, it is still unclear how many of those will be retained upon formation, and how many will be ejected through subsequent dynamical interactions. No such black holes have been found in any Milky Way globular cluster until the recent discovery of stellar-mass black holes in the globular cluster M22 (NGC 6656) with now an estimated population of 5-100 black holes. We present a direct N-body model of a star cluster of the same absolute and dynamical age as M22. Imposing an initial retention fraction of approx. 10% for black holes, 16 stellar-remnant black holes are retained at a cluster age of 12 Gyr, in agreement with the estimate for M22. Of those 16 BHs, two are in a binary system with a main sequence star each while also one pure black hole binary is present. We argue that multiple black holes can be present in any Milky Way cluster with an extended core radius, such as M22 or the model presented here.
△ Less
Submitted 28 November, 2012;
originally announced November 2012.
-
Direct N-Body Modeling of the Old Open Cluster NGC 188: A Detailed Comparison of Theoretical and Observed Binary Star and Blue Straggler Populations
Authors:
Aaron M. Geller,
Jarrod R. Hurley,
Robert D. Mathieu
Abstract:
(Abridged) Following on from a recently completed radial-velocity survey of the old (7 Gyr) open cluster NGC 188 in which we study in detail the solar-type hard binaries and blue stragglers of the cluster, here we investigate the dynamical evolution of NGC 188 through a sophisticated N-body model. We employ the observed binary properties of the young (150 Myr) open cluster M35, where possible, to…
▽ More
(Abridged) Following on from a recently completed radial-velocity survey of the old (7 Gyr) open cluster NGC 188 in which we study in detail the solar-type hard binaries and blue stragglers of the cluster, here we investigate the dynamical evolution of NGC 188 through a sophisticated N-body model. We employ the observed binary properties of the young (150 Myr) open cluster M35, where possible, to guide our choices for parameters of the initial binary population. At 7 Gyr the main-sequence solar-type hard-binary population in the model matches that of NGC 188 in both binary frequency and distributions of orbital parameters. This agreement between the model and observations is in a large part due to the similarities between the NGC 188 and M35 solar-type binaries. Indeed, among the 7 Gyr main-sequence binaries in the model, only those with P>1000 days show potentially observable evidence for modifications by dynamical encounters. This emphasizes the importance of defining accurate initial conditions for star cluster models, which we propose is best accomplished through comparisons with observations of young open clusters like M35. Furthermore, this suggests that observations of the present-day binaries in even old open clusters can provide valuable information on their primordial binary populations. However, despite the model's successes at matching the true cluster, the model underproduces blue stragglers and produces an overabundance of long-period circular main-sequence--white-dwarf binaries as compared to NGC 188. We conclude that improvements in the physics of mass transfer and common envelope may in fact solve both discrepancies with the observations. This project highlights the unique accessibility of open clusters to both comprehensive observational surveys and full-scale N-body simulations, and underscores the importance of open clusters to the study of star cluster dynamics.
△ Less
Submitted 2 November, 2012; v1 submitted 4 October, 2012;
originally announced October 2012.
-
A direct N-body model of core-collapse and core oscillations
Authors:
Jarrod R. Hurley,
Michael M. Shara
Abstract:
We report on the results of a direct N-body simulation of a star cluster that started with N = 200 000, comprising 195 000 single stars and 5 000 primordial binaries. The code used for the simulation includes stellar evolution, binary evolution, an external tidal field and the effects of two-body relaxation. The model cluster is evolved to 12 Gyr, losing more than 80% of its stars in the process.…
▽ More
We report on the results of a direct N-body simulation of a star cluster that started with N = 200 000, comprising 195 000 single stars and 5 000 primordial binaries. The code used for the simulation includes stellar evolution, binary evolution, an external tidal field and the effects of two-body relaxation. The model cluster is evolved to 12 Gyr, losing more than 80% of its stars in the process. It reaches the end of the main core-collapse phase at 10.5 Gyr and experiences core oscillations from that point onwards -- direct numerical confirmation of this phenomenon. However, we find that after a further 1 Gyr the core oscillations are halted by the ejection of a massive binary comprised of two black holes from the core, producing a core that shows no signature of the prior core-collapse. We also show that the results of previous studies with N ranging from 500 to 100 000 scale well to this new model with larger N. In particular, the timescale to core-collapse (in units of the relaxation timescale), mass segregation, velocity dispersion, and the energies of the binary population all show similar behaviour at different N.
△ Less
Submitted 23 August, 2012;
originally announced August 2012.
-
N-body models of globular clusters: metallicity, half-light radii and mass-to-light ratios
Authors:
Anna C. Sippel,
Jarrod R. Hurley,
Juan P. Madrid,
William E. Harris
Abstract:
Size differences of approx. 20% between red (metal-rich) and blue (metal-poor) sub-populations of globular clusters have been observed, generating an ongoing debate as to weather these originate from projection effects or the difference in metallicity. We present direct N-body simulations of metal-rich and metal-poor stellar populations evolved to study the effects of metallicity on cluster evolut…
▽ More
Size differences of approx. 20% between red (metal-rich) and blue (metal-poor) sub-populations of globular clusters have been observed, generating an ongoing debate as to weather these originate from projection effects or the difference in metallicity. We present direct N-body simulations of metal-rich and metal-poor stellar populations evolved to study the effects of metallicity on cluster evolution. The models start with N = 100000 stars and include primordial binaries. We also take metallicity dependent stellar evolution and an external tidal field into account. We find no significant difference for the half-mass radii of those models, indicating that the clusters are structurally similar. However, utilizing observational tools to fit half-light (or effective) radii confirms that metallicity effects related to stellar evolution combined with dynamical effects such as mass segregation produce an apparent size difference of 17% on average. The metallicity effect on the overall cluster luminosity also leads to higher mass-to-light ratios for metal-rich clusters.
△ Less
Submitted 23 August, 2012;
originally announced August 2012.
-
The Size Scale of Star Clusters
Authors:
Juan P. Madrid,
Jarrod R. Hurley,
Anna C. Sippel
Abstract:
Direct N-body simulations of star clusters in a realistic Milky Way-like potential are carried out using the code NBODY6. Based on these simulations a new relationship between scale size and galactocentric distance is derived: the scale size of star clusters is proportional to the hyperbolic tangent of the galactocentric distance. The half-mass radius of star clusters increases systematically with…
▽ More
Direct N-body simulations of star clusters in a realistic Milky Way-like potential are carried out using the code NBODY6. Based on these simulations a new relationship between scale size and galactocentric distance is derived: the scale size of star clusters is proportional to the hyperbolic tangent of the galactocentric distance. The half-mass radius of star clusters increases systematically with galactocentric distance but levels off when star clusters orbit the galaxy beyond ~40 kpc. These simulations show that the half-mass radius of individual star clusters varies significantly as they evolve over a Hubble time, more so for clusters with shorter relaxation times, and remains constant through several relaxation times only in certain situations when expansion driven by the internal dynamics of the star cluster and the influence of the host galaxy tidal field balance each other. Indeed, the radius of a star cluster evolving within the inner 20 kpc of a realistic galactic gravitational potential is severely truncated by tidal interactions and does not remain constant over a Hubble time. Furthermore, the half-mass radius of star clusters measured with present day observations bears no memory of the original cluster size. Stellar evolution and tidal stripping are the two competing physical mechanisms that determine the present day size of globular clusters. These simulations also show that extended star clusters can form at large galactocentric distances while remaining fully bound to the host galaxy. There is thus no need to invoke accretion from an external galaxy to explain the presence of extended clusters at large galactocentric distances in a Milky Way-type galaxy.
△ Less
Submitted 1 August, 2012;
originally announced August 2012.
-
N-body Models of Extended Clusters
Authors:
Jarrod R. Hurley,
A. Dougal Mackey
Abstract:
We use direct N-body simulations to investigate the evolution of star clusters with large size-scales with the particular goal of understanding the so-called extended clusters observed in various Local Group galaxies, including M31 and NGC6822. The N-body models incorporate a stellar mass function, stellar evolution and the tidal field of a host galaxy. We find that extended clusters can arise nat…
▽ More
We use direct N-body simulations to investigate the evolution of star clusters with large size-scales with the particular goal of understanding the so-called extended clusters observed in various Local Group galaxies, including M31 and NGC6822. The N-body models incorporate a stellar mass function, stellar evolution and the tidal field of a host galaxy. We find that extended clusters can arise naturally within a weak tidal field provided that the tidal radius is filled at the start of the evolution. Differences in the initial tidal filling-factor can produce marked differences in the subsequent evolution of clusters and the size-scales that would be observed. These differences are more marked than any produced by internal evolution processes linked to the properties of cluster binary stars or the action of an intermediate-mass black hole, based on models performed in this work and previous work to date. Models evolved in a stronger tidal field show that extended clusters cannot form and evolve within the inner regions of a galaxy such as M31. Instead our results support the suggestion many extended clusters found in large galaxies were accreted as members of dwarf galaxies that were subsequently disrupted. Our results also enhance the recent suggestion that star clusters evolve to a common sequence in terms of their size and mass.
△ Less
Submitted 29 August, 2010;
originally announced August 2010.
-
The Progeny of Stellar Dynamics and Stellar Evolution within an N-body model of NGC 188
Authors:
A. M. Geller,
J. R. Hurley,
R. D. Mathieu
Abstract:
We present a direct N-body simulation modeling the evolution of the old (7 Gyr) open cluster NGC 188. This is the first N-body open cluster simulation whose initial binary population is directly defined by observations of a specific open cluster: M35 (150 Myr). We compare the simulated color-magnitude diagram at 7 Gyr to that of NGC 188, and discuss the blue stragglers produced in the simulation…
▽ More
We present a direct N-body simulation modeling the evolution of the old (7 Gyr) open cluster NGC 188. This is the first N-body open cluster simulation whose initial binary population is directly defined by observations of a specific open cluster: M35 (150 Myr). We compare the simulated color-magnitude diagram at 7 Gyr to that of NGC 188, and discuss the blue stragglers produced in the simulation. We compare the solar-type main sequence binary period and eccentricity distributions of the simulation to detailed observations of similar binaries in NGC 188. These results demonstrate the importance of detailed observations in guiding N-body open cluster simulations. Finally, we discuss the implications of a few discrepancies between the NGC 188 model and observations and suggest a few methods for bringing N-body open cluster simulations into better agreement with observations.
△ Less
Submitted 23 November, 2009;
originally announced November 2009.
-
Star Cluster Evolution in Dark Matter Dominated Galaxies
Authors:
A. K. Praagman,
J. R. Hurley,
C. Power
Abstract:
We investigate the influence of the external tidal field of a dark matter halo on the dynamical evolution of star clusters using direct N-body simulations, where we assume that the halo is described by a Navarro, Frenk & White mass profile which has an inner density cusp. We assess how varying the mass and concentration of the halo affects the rate at which the star cluster loses mass and we fin…
▽ More
We investigate the influence of the external tidal field of a dark matter halo on the dynamical evolution of star clusters using direct N-body simulations, where we assume that the halo is described by a Navarro, Frenk & White mass profile which has an inner density cusp. We assess how varying the mass and concentration of the halo affects the rate at which the star cluster loses mass and we find that increasing halo mass and concentration drives enhanced mass loss rates and in principle shorter cluster disruption timescales. In addition, we examine disruption timescales in a three-component model of a galaxy (bulge, disk and dark matter halo) and find good agreement with results based on an empirical model of the Galactic potential if we assume a halo mass of ~1e12 solar masses. In general, dark matter halos are expected to contribute significantly to the masses of galaxies and should not be ignored when modelling the evolution of star clusters. We extend our results to discuss how this can have a potentially profound effect on the disruption timescales of globular clusters, suggesting that we may underestimate the rate at which primordial globular clusters are disrupted.
△ Less
Submitted 11 May, 2009;
originally announced May 2009.
-
On The Maximum Mass of Stellar Black Holes
Authors:
Krzysztof Belczynski,
Tomasz Bulik,
Chris L. Fryer,
Ashley Ruiter,
Francesca Valsecchi,
Jorick S. Vink,
Jarrod R. Hurley
Abstract:
We present the spectrum of compact object masses: neutron stars and black holes that originate from single stars in different environments. In particular, we calculate the dependence of maximum black hole mass on metallicity and on some specific wind mass loss rates (e.g., Hurley et al. and Vink et al.). Our calculations show that the highest mass black holes observed in the Galaxy M_bh = 15 Msun…
▽ More
We present the spectrum of compact object masses: neutron stars and black holes that originate from single stars in different environments. In particular, we calculate the dependence of maximum black hole mass on metallicity and on some specific wind mass loss rates (e.g., Hurley et al. and Vink et al.). Our calculations show that the highest mass black holes observed in the Galaxy M_bh = 15 Msun in the high metallicity environment (Z=Zsun=0.02) can be explained with stellar models and the wind mass loss rates adopted here. To reach this result we had to set Luminous Blue Variable mass loss rates at the level of about 0.0001 Msun/yr and to employ metallicity dependent Wolf-Rayet winds. With such winds, calibrated on Galactic black hole mass measurements, the maximum black hole mass obtained for moderate metallicity (Z=0.3 Zsun=0.006) is M_bh,max = 30 Msun. This is a rather striking finding as the mass of the most massive known stellar black hole is M_bh = 23-34 Msun and, in fact, it is located in a small star forming galaxy with moderate metallicity. We find that in the very low (globular cluster-like) metallicity environment the maximum black hole mass can be as high as M_bh,max = 80 Msun (Z=0.01 Zsun=0.0002). It is interesting to note that X-ray luminosity from Eddington limited accretion onto an 80 Msun black hole is of the order of about 10^40 erg/s and is comparable to luminosities of some known ULXs. We emphasize that our results were obtained for single stars only and that binary interactions may alter these maximum black hole masses (e.g., accretion from a close companion). This is strictly a proof-of-principle study which demonstrates that stellar models can naturally explain even the most massive known stellar black holes.
△ Less
Submitted 9 March, 2010; v1 submitted 17 April, 2009;
originally announced April 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.
-
N-body Simulations with Live Stellar Evolution
Authors:
Ross P. Church,
Christopher A. Tout,
Jarrod R. Hurley
Abstract:
An N-body code containing live stellar evolution through combination of the software packages NBODY6 and STARS is presented. Operational details of the two codes are outlined and the changes that have been made to combine them discussed. We have computed the evolution of clusters of 10 000 stars using the combined code and we compare the results with those obtained using NBODY6 and the synthetic…
▽ More
An N-body code containing live stellar evolution through combination of the software packages NBODY6 and STARS is presented. Operational details of the two codes are outlined and the changes that have been made to combine them discussed. We have computed the evolution of clusters of 10 000 stars using the combined code and we compare the results with those obtained using NBODY6 and the synthetic stellar evolution code SSE. We find that, providing the physics package within STARS is set up correctly to match the parameters of the models used to construct SSE, the results are very similar. This provides a good indication that the new code is working well. We also demonstrate how this physics can be changed simply in the new code with convective overshooting as an example. Similar changes in SSE would require considerable reworking of the model fits. We conclude by outlining proposed future development of the code to include more complete models of single stars and binary star systems.
△ Less
Submitted 4 February, 2009;
originally announced February 2009.
-
Self-consistent simulations of star cluster formation from gas clouds under the influence of galaxy-scale tidal fields
Authors:
Jarrod R. Hurley,
Kenji Bekki
Abstract:
We present the first results of a project aimed at following the formation and long-term dynamical evolution of star clusters within the potential of a host galaxy. Here we focus on a model evolved within a simplified potential representing the Large Magellanic Cloud. This demonstrates for the first time the self-consistent formation of a bound star cluster from a giant molecular cloud. The mode…
▽ More
We present the first results of a project aimed at following the formation and long-term dynamical evolution of star clusters within the potential of a host galaxy. Here we focus on a model evolved within a simplified potential representing the Large Magellanic Cloud. This demonstrates for the first time the self-consistent formation of a bound star cluster from a giant molecular cloud. The model cluster reproduces the density profiles and structural characteristics of observed star clusters.
△ Less
Submitted 8 July, 2008;
originally announced July 2008.
-
Evolution of stellar collision products in open clusters. I. Blue stragglers in N-body models of M67
Authors:
E. Glebbeek,
O. R. Pols,
J. R. Hurley
Abstract:
Stellar collisions are an important formation channel for blue straggler stars in globular and old open clusters. Hydrodynamical simulations have shown that the remnants of such collisions are out of thermal equilibrium, are not strongly mixed and can rotate very rapidly.
Detailed evolution models of collision products are needed to interpret observed blue straggler populations and to use them…
▽ More
Stellar collisions are an important formation channel for blue straggler stars in globular and old open clusters. Hydrodynamical simulations have shown that the remnants of such collisions are out of thermal equilibrium, are not strongly mixed and can rotate very rapidly.
Detailed evolution models of collision products are needed to interpret observed blue straggler populations and to use them to probe the dynamical history of a star cluster. We expand on previous studies by presenting an efficient procedure to import the results of detailed collision simulations into a fully implicit stellar evolution code. Our code is able to evolve stellar collision products in a fairly robust manner and allows for a systematic study of their evolution.
Using our code we have constructed detailed models of the collisional blue stragglers produced in the $N$-body simulation of M67 performed by Hurley \emph{et al.} in 2005. We assume the collisions are head-on and thus ignore the effects of rotation in this paper.
Our detailed models are more luminous than normal stars of the same mass and in the same stage of evolution, but cooler than homogeneously mixed versions of the collision products. The increased luminosity and inefficient mixing decrease the remaining main-sequence lifetimes of the collision products, which are much shorter than predicted by the simple prescription commonly used in $N$-body simulations.
△ Less
Submitted 4 June, 2008;
originally announced June 2008.
-
Deep ACS Imaging in the Globular Cluster NGC6397: Dynamical Models
Authors:
Jarrod R. Hurley,
Michael M. Shara,
Harvey B. Richer,
Ivan R. King,
D. Saul Davis,
Jason S. Kalirai,
Brad M. S. Hansen,
Aaron Dotter,
Jay Anderson,
Gregory G. Fahlman,
R. Michael Rich
Abstract:
We present N-body models to complement deep imaging of the metal-poor core-collapsed cluster NGC6397 obtained with the Hubble Space Telescope. All simulations include stellar and binary evolution in-step with the stellar dynamics and account for the tidal field of the Galaxy. We focus on the results of a simulation that began with 100000 objects (stars and binaries), 5% primordial binaries and P…
▽ More
We present N-body models to complement deep imaging of the metal-poor core-collapsed cluster NGC6397 obtained with the Hubble Space Telescope. All simulations include stellar and binary evolution in-step with the stellar dynamics and account for the tidal field of the Galaxy. We focus on the results of a simulation that began with 100000 objects (stars and binaries), 5% primordial binaries and Population II metallicity. After 16 Gyr of evolution the model cluster has about 20% of the stars remaining and has reached core-collapse. We compare the color-magnitude diagrams of the model at this age for the central region and an outer region corresponding to the observed field of NGC6397 (about 2-3 half-light radii from the cluster centre). This demonstrates that the white dwarf population in the outer region has suffered little modification from dynamical processes - contamination of the luminosity function by binaries and white dwarfs with non-standard evolution histories is minimal and should not significantly affect measurement of the cluster age. We also show that the binary fraction of main-sequence stars observed in the NGC6397 field can be taken as representative of the primordial binary fraction of the cluster. For the mass function of the main-sequence stars we find that although this has been altered significantly by dynamics over the cluster lifetime, especially in the central and outer regions, that the position of the observed field is close to optimal for recovering the initial mass function of the cluster stars (below the current turn-off mass). More generally we look at how the mass function changes with radius in a dynamically evolved stellar cluster and suggest where the best radial position to observe the initial mass function is for clusters of any age.
△ Less
Submitted 28 February, 2008;
originally announced February 2008.
-
Monte Carlo Simulations of Star Clusters - IV. Calibration of the Monte Carlo Code and Comparison with Observations for the Open Cluster M67
Authors:
Mirek Giersz,
Douglas C. Heggie,
Jarrod R. Hurley
Abstract:
We outline the steps needed in order to incorporate the evolution of single and binary stars into a particular Monte Carlo code for the dynamical evolution of a star cluster. We calibrate the results against N-body simulations, and present models for the evolution of the old open cluster M67 (which has been studied thoroughly in the literature with N-body techniques). The calibration is done by…
▽ More
We outline the steps needed in order to incorporate the evolution of single and binary stars into a particular Monte Carlo code for the dynamical evolution of a star cluster. We calibrate the results against N-body simulations, and present models for the evolution of the old open cluster M67 (which has been studied thoroughly in the literature with N-body techniques). The calibration is done by choosing appropriate free code parameters. We describe in particular the evolution of the binary, white dwarf and blue straggler populations, though not all channels for blue straggler formation are represented yet in our simulations. Calibrated Monte Carlo runs show good agreement with results of N-body simulations not only for global cluster parameters, but also for e.g. binary fraction, luminosity function and surface brightness. Comparison of Monte Carlo simulations with observational data for M67 shows that is possible to get reasonably good agreement between them. Unfortunately, because of the large statistical fluctuations of the numerical data and uncertainties in the observational data the inferred conclusions about the cluster initial conditions are not firm.
△ Less
Submitted 24 January, 2008;
originally announced January 2008.
-
Ratios of star cluster core and half-mass radii: a cautionary note on intermediate-mass black holes in star clusters
Authors:
Jarrod R. Hurley
Abstract:
There is currently much interest in the possible presence of intermediate-mass black holes in the cores of globular clusters. Based on theoretical arguments and simulation results it has previously been suggested that a large core radius -- or particularly a large ratio of the core radius to half-mass radius -- is a promising indicator for finding such a black hole in a star cluster. In this stu…
▽ More
There is currently much interest in the possible presence of intermediate-mass black holes in the cores of globular clusters. Based on theoretical arguments and simulation results it has previously been suggested that a large core radius -- or particularly a large ratio of the core radius to half-mass radius -- is a promising indicator for finding such a black hole in a star cluster. In this study N-body models of 100000 stars with and without primordial binaries are used to investigate the long-term structural evolution of star clusters. Importantly, the simulation data is analysed using the same processes by which structural parameters are extracted from observed star clusters. This gives a ratio of the core and half-mass (or half-light) radii that is directly comparable to the Galactic globular cluster sample. As a result, it is shown that the ratios observed for the bulk of this sample can be explained without the need for an intermediate-mass black hole. Furthermore, it is possible that clusters with large core to half-light radius ratios harbour a black-hole binary (comprised of stellar mass black holes) rather than a single massive black hole. This work does not rule out the existence of intermediate-mass black holes in the cores of at least some star clusters.
△ Less
Submitted 5 May, 2007;
originally announced May 2007.
-
The core binary fractions of star clusters from realistic simulations
Authors:
Jarrod R. Hurley,
Sverre J. Aarseth,
Michael M. Shara
Abstract:
We investigate the evolution of binary fractions in star clusters using N-body models of up to 100000 stars. Primordial binary frequencies in these models range from 5% to 50%. Simulations are performed with the NBODY4 code and include a full mass spectrum of stars, stellar evolution, binary evolution and the tidal field of the Galaxy. We find that the overall binary fraction of a cluster almost…
▽ More
We investigate the evolution of binary fractions in star clusters using N-body models of up to 100000 stars. Primordial binary frequencies in these models range from 5% to 50%. Simulations are performed with the NBODY4 code and include a full mass spectrum of stars, stellar evolution, binary evolution and the tidal field of the Galaxy. We find that the overall binary fraction of a cluster almost always remains close to the primordial value, except at late times when a cluster is near dissolution. A critical exception occurs in the central regions where we observe a marked increase in binary fraction with time -- a simulation starting with 100000 stars and 5% binaries reached a core binary frequency as high as 40% at the end of the core-collapse phase (occurring at 16 Gyr with ~20000 stars remaining). Binaries are destroyed in the core by a variety of processes as a cluster evolves, but the combination of mass-segregation and creation of new binaries in exchange interactions produces the observed increase in relative number. We also find that binaries are cycled into and out of cluster cores in a manner that is analogous to convection in stars. For models of 100000 stars we show that the evolution of the core-radius up to the end of the initial phase of core-collapse is not affected by the exact value of the primordial binary frequency (for frequencies of 10% or less). We discuss the ramifications of our results for the likely primordial binary content of globular clusters.
△ Less
Submitted 2 April, 2007;
originally announced April 2007.
-
Populating the Galaxy with close Be binaries
Authors:
P. D. Kiel,
J. R. Hurley,
J. R. Murray,
K. Hayasaki
Abstract:
Be/X-ray binaries comprise roughly two-thirds of the high-mass X-ray binaries (HMXBs), which is a class of X-ray binaries that results from the high mass of the companion or donor star (> 10 solar masses). Currently the formation and evolution of X-ray producing Be binaries is a matter of great debate. Modelling of these systems requires knowledge of Be star evolution and also consideration of h…
▽ More
Be/X-ray binaries comprise roughly two-thirds of the high-mass X-ray binaries (HMXBs), which is a class of X-ray binaries that results from the high mass of the companion or donor star (> 10 solar masses). Currently the formation and evolution of X-ray producing Be binaries is a matter of great debate. Modelling of these systems requires knowledge of Be star evolution and also consideration of how the evolution changes when the star is in close proximity to a companion. Within this work we complete a full population synthesis study of Be binaries for the Galaxy. The results for the first time match aspects of the observational data, most notably the observed upper limit to the period distribution. We conclude that greater detailed studies on the evolution of Be stars within X-ray binaries needs to be completed, so that rapid binary evolution population synthesis packages may best evolve these systems.
△ Less
Submitted 7 March, 2007;
originally announced March 2007.
-
The Space Motion of the Globular Cluster NGC 6397
Authors:
Jasonjot S. Kalirai,
Jay Anderson,
Harvey B. Richer,
Ivan R. King,
James P. Brewer,
Giovanni Carraro,
Saul D. Davis,
Gregory G. Fahlman,
Brad M. S. Hansen,
Jarrod R. Hurley,
Sebastien Lepine,
David B. Reitzel,
R. Michael Rich,
Michael M. Shara,
Peter B. Stetson
Abstract:
As a by-product of high-precision, ultra-deep stellar photometry in the Galactic globular cluster NGC 6397 with the Hubble Space Telescope, we are able to measure a large population of background galaxies whose images are nearly point-like. These provide an extragalactic reference frame of unprecedented accuracy, relative to which we measure the most accurate absolute proper motion ever determin…
▽ More
As a by-product of high-precision, ultra-deep stellar photometry in the Galactic globular cluster NGC 6397 with the Hubble Space Telescope, we are able to measure a large population of background galaxies whose images are nearly point-like. These provide an extragalactic reference frame of unprecedented accuracy, relative to which we measure the most accurate absolute proper motion ever determined for a globular cluster. We find mu_alpha = 3.56 +/- 0.04 mas/yr and mu_delta = -17.34 +/- 0.04 mas/yr. We note that the formal statistical errors quoted for the proper motion of NGC 6397 do not include possible unavoidable sources of systematic errors, such as cluster rotation. These are very unlikely to exceed a few percent. We use this new proper motion to calculate NGC 6397's UVW space velocity and its orbit around the Milky Way, and find that the cluster has made frequent passages through the Galactic disk.
△ Less
Submitted 6 February, 2007; v1 submitted 26 January, 2007;
originally announced January 2007.
-
Populating the Galaxy with low-mass X-ray binaries
Authors:
Paul D. Kiel,
Jarrod R. Hurley
Abstract:
We perform binary population synthesis calculations to investigate the incidence of low-mass X-ray binaries and their birth rate in the Galaxy. We use a binary evolution algorithm that models all the relevant processes including tidal circularization and synchronization. Parameters in the evolution algorithm that are uncertain and may affect X-ray binary formation are allowed to vary during the…
▽ More
We perform binary population synthesis calculations to investigate the incidence of low-mass X-ray binaries and their birth rate in the Galaxy. We use a binary evolution algorithm that models all the relevant processes including tidal circularization and synchronization. Parameters in the evolution algorithm that are uncertain and may affect X-ray binary formation are allowed to vary during the investigation. We agree with previous studies that under standard assumptions of binary evolution the formation rate and number of black-hole low-mass X-ray binaries predicted by the model are more than an order of magnitude less than what is indicated by observations. We find that the common-envelope process cannot be manipulated to produce significant numbers of black-hole low-mass X-ray binaries. However, by simply reducing the mass-loss rate from helium stars adopted in the standard model, to a rate that agrees with the latest data, we produce a good match to the observations. Including low-mass X-ray binaries that evolve from intermediate-mass systems also leads to favourable results. We stress that constraints on the X-ray binary population provided by observations are used here merely as a guide as surveys suffer from incompleteness and much uncertainty is involved in the interpretation of results.
△ Less
Submitted 2 May, 2006;
originally announced May 2006.
-
A Complete N-body Model of the Old Open Cluster M67
Authors:
Jarrod R. Hurley,
Onno R. Pols,
Sverre J. Aarseth,
Christopher A. Tout
Abstract:
The old open cluster M67 is an ideal testbed for current cluster evolution models because of its dynamically evolved structure and rich stellar populations that show clear signs of interaction between stellar, binary and cluster evolution. Here we present the first truly direct N-body model for M67, evolved from zero age to 4 Gyr taking full account of cluster dynamics as well as stellar and bin…
▽ More
The old open cluster M67 is an ideal testbed for current cluster evolution models because of its dynamically evolved structure and rich stellar populations that show clear signs of interaction between stellar, binary and cluster evolution. Here we present the first truly direct N-body model for M67, evolved from zero age to 4 Gyr taking full account of cluster dynamics as well as stellar and binary evolution. Our preferred model starts with 12000 single stars and 12000 binaries placed in a Galactic tidal field at 8.0 kpc from the Galactic Centre. Our choices for the initial conditions and for the primordial binary population are explained in detail. At 4 Gyr, the age of M67, the total mass has reduced by 90% as a result of mass loss and stellar escapes. The mass and half-mass radius of luminous stars in the cluster are a good match to observations although the model is more centrally concentrated than observations indicate. The stellar mass and luminosity functions are significantly flattened by preferential escape of low-mass stars. We find that M67 is dynamically old enough that information about the initial mass function is lost, both from the current luminosity function and from the current mass fraction in white dwarfs. The model contains 20 blue stragglers at 4 Gyr which is slightly less than the 28 observed in M67. Nine are in binaries. The blue stragglers were formed by a variety of means and we find formation paths for the whole variety observed in M67. Both the primordial binary population and the dynamical cluster environment play an essential role in shaping the population. A substantial population of short-period primordial binaries (with periods less than a few days) is needed to explain the observed number of blue stragglers in M67.
△ Less
Submitted 10 July, 2005;
originally announced July 2005.
-
Inclusion of Binaries in Evolutionary Population Synthesis
Authors:
Fenghui Zhang,
Zhanwen Han,
Lifang Li,
Jarrod R. Hurley
Abstract:
Using evolutionary population synthesis we present integrated colours, integrated spectral energy distributions and absorption-line indices defined by the Lick Observatory image dissector scanner (referred to as Lick/IDS) system, for an extensive set of instantaneous burst binary stellar populations with and without binary interactions. The ages of the populations are in the range $1-15 $Gyr and…
▽ More
Using evolutionary population synthesis we present integrated colours, integrated spectral energy distributions and absorption-line indices defined by the Lick Observatory image dissector scanner (referred to as Lick/IDS) system, for an extensive set of instantaneous burst binary stellar populations with and without binary interactions. The ages of the populations are in the range $1-15 $Gyr and the metallicities are in the range $0.0001-0.03$. By comparing the results for populations with and without binary interactions we show that the inclusion of binary interactions makes the integrated $\rm U-B$, $\rm B-V$, $\rm V-R$ and $\rm R-I$ colours and all Lick/IDS spectral absorption indices (except for $\rm H_β$) substantially smaller. In other words binary evolution makes a population appear bluer. This effect raises the derived age and metallicity of the population.
We calculate several sets of additional solar-metallicity binary stellar populations to explore the influence of input parameters to the binary evolution algorithm (the common-envelope ejection efficiency and the stellar wind mass-loss rate) on the resulting integrated colours. We also look at the dependence on the choice of distribution functions used to generate the initial binary population. The results show that variations in the choice of input model parameters and distributions can significantly affect the results. However, comparing the discrepancies that exist between the colours of various models, we find that the differences are less than those produced between the models with and those without binary interactions. Therefore it is very necessary to consider binary interactions in order to draw accurate conclusions from evolutionary population synthesis work.
△ Less
Submitted 30 November, 2004;
originally announced December 2004.
-
Integrated Spectral Energy Distributions and Absorption Feature Indices of Single Stellar Populations
Authors:
Fenghui Zhang,
Zhanwen Han,
Lifang Li,
Jarrod R. Hurley
Abstract:
Using evolutionary population synthesis, we present integrated spectral energy distributions and absorption-line indices defined by the Lick Observatory image dissector scanner (referred to as Lick/IDS) system, for an extensive set of instantaneous burst single stellar populations (SSPs). The ages of the SSPs are in the range 1-19 Gyr and the metallicities [Fe/H] are in the range -2.3 - 0.2. Our…
▽ More
Using evolutionary population synthesis, we present integrated spectral energy distributions and absorption-line indices defined by the Lick Observatory image dissector scanner (referred to as Lick/IDS) system, for an extensive set of instantaneous burst single stellar populations (SSPs). The ages of the SSPs are in the range 1-19 Gyr and the metallicities [Fe/H] are in the range -2.3 - 0.2. Our models use the rapid single stellar evolution algorithm of Hurley, Pols and Tout for the stellar evolutionary tracks, the empirical and semi-empirical calibrated BaSeL-2.0 model of Lejeune, Cuisinier and Buser for the library of stellar spectra and the empirical fitting functions of Worthey, Faber, Gonzalez and Burstein for the Lick/IDS spectral absorption-feature indices.
Applying our synthetic Lick/IDS absorption-line indices to the merit function, we obtain the age and the metallicity of the central region of M32, which can be well explained by an instantaneous SSP with an age of 6.5 Gyr and a metallicity similar to solar. Applying the derived age and the metallicity from the merit function to a number of index-index diagrams, we find that the plots of Hbeta-Fe5015 and Hbeta-Fe5782 are the best index-index diagrams from which we can directly obtain reasonable age and metallicity.
△ Less
Submitted 11 March, 2004;
originally announced March 2004.
-
Evolutionary Population Synthesis for Binary Stellar Populations
Authors:
Fenghui Zhang,
Zhanwen Han,
Lifang Li,
Jarrod R. Hurley
Abstract:
We present integrated colours, integrated spectral energy distributions, and absorption-line indices, for instantaneous burst solar-metallicity binary stellar populations with ages in the range 1 - 15Gyr. By comparing the results for populations with and without binary interactions we show that the inclusion of binary interactions makes the appearance of the population substantially bluer -- thi…
▽ More
We present integrated colours, integrated spectral energy distributions, and absorption-line indices, for instantaneous burst solar-metallicity binary stellar populations with ages in the range 1 - 15Gyr. By comparing the results for populations with and without binary interactions we show that the inclusion of binary interactions makes the appearance of the population substantially bluer -- this is the case for each of the quantities we have considered. This effect raises the derived age and metallicity of the population. Therefore it is necessary to consider binary interactions in order to draw accurate conclusions from evolutionary population synthesis work.
△ Less
Submitted 11 March, 2004;
originally announced March 2004.
-
Core radius evolution of star clusters
Authors:
Mark I. Wilkinson,
Jarrod R. Hurley,
A. Dougal Mackey,
Gerard F. Gilmore,
Christopher A. Tout
Abstract:
We use N-body simulations of star clusters to investigate the possible dynamical origins of the observed spread in core radius among intermediate-age and old star clusters in the Large Magellanic Cloud (LMC). Two effects are considered, a time-varying external tidal field and variations in primordial hard binary fraction. Simulations of clusters orbiting a point-mass galaxy show similar core rad…
▽ More
We use N-body simulations of star clusters to investigate the possible dynamical origins of the observed spread in core radius among intermediate-age and old star clusters in the Large Magellanic Cloud (LMC). Two effects are considered, a time-varying external tidal field and variations in primordial hard binary fraction. Simulations of clusters orbiting a point-mass galaxy show similar core radius evolution for clusters on both circular and elliptical orbits and we therefore conclude that the tidal field of the LMC has not yet significantly influenced the evolution of the intermediate-age clusters. The presence of large numbers of hard primordial binaries in a cluster leads to core radius expansion; however, the magnitude of the effect is insufficient to explain the observations. Further, the range of binary fractions required to produce significant core radius growth is inconsistent with the observational evidence that all the LMC clusters have similar stellar luminosity functions.
△ Less
Submitted 29 April, 2003;
originally announced April 2003.
-
White Dwarf Sequences in Dense Star Clusters
Authors:
Jarrod R. Hurley,
Michael M. Shara
Abstract:
We use the results of realistic N-body simulations to investigate the appearance of the white dwarf population in dense star clusters. We show that the presence of a substantial binary population in a star cluster, and the interaction of this population with the cluster environment, has serious consequences for the morphology of the observed white dwarf sequence and the derived white dwarf cooli…
▽ More
We use the results of realistic N-body simulations to investigate the appearance of the white dwarf population in dense star clusters. We show that the presence of a substantial binary population in a star cluster, and the interaction of this population with the cluster environment, has serious consequences for the morphology of the observed white dwarf sequence and the derived white dwarf cooling age of the cluster. We find that over time the dynamical evolution of the cluster -- mass-segregation, stellar interactions and tidal stripping -- hampers the use of white dwarfs as tracers of the initial mass function, and also leads to a significant enhancement of the white dwarf mass fraction. Future observations of star clusters should be conducted slightly interior to the half-mass radius of the cluster in order to best obtain information about the cluster age and initial mass function from the white dwarf luminosity function. The evolution of binary stars and the cluster environment must necessarily be accounted for when studying the white dwarf populations of dynamically evolved star clusters.
△ Less
Submitted 6 February, 2003;
originally announced February 2003.
-
Star Clusters as Type Ia Supernova Factories
Authors:
Michael M. Shara,
Jarrod R. Hurley
Abstract:
We find a remarkably enhanced production rate in star clusters (relative to the field) of very short period, massive double-white-dwarf stars and of giant-white dwarf binaries. These results are based on N-body simulations performed with the new GRAPE-6 special purpose hardware and are important in identifying and characterizing the progenitors of type Ia supernovae. The high incidence of very c…
▽ More
We find a remarkably enhanced production rate in star clusters (relative to the field) of very short period, massive double-white-dwarf stars and of giant-white dwarf binaries. These results are based on N-body simulations performed with the new GRAPE-6 special purpose hardware and are important in identifying and characterizing the progenitors of type Ia supernovae. The high incidence of very close double-white-dwarf systems is the result of dynamical encounters between (mostly) primordial binaries and other cluster stars. Orbital hardening rapidly drives these degenerate binaries to periods under ~10 hours. Gravitational radiation emission and mergers producing supra-Chandrasekhar objects follow in less than a Hubble time. If most stars are born in clusters then estimates of the double white dwarf merger rates in galaxies (due to cluster dynamical interaction) must be increased more than tenfold. A majority of the Roche lobe overflow giant-white dwarf binaries are not primordial; they are produced in exchange reactions. Most cases resulted in a common-envelope and formation of a double-white-dwarf binary rather than Supersoft X-ray sources leading possibly to a type Ia supernova.
△ Less
Submitted 8 February, 2002;
originally announced February 2002.
-
Evolution of binary stars and the effect of tides on binary populations
Authors:
Jarrod R. Hurley,
Christopher A. Tout,
Onno R. Pols
Abstract:
We present a rapid binary evolution algorithm that enables modelling of even the most complex binary systems. In addition to all aspects of single star evolution, features such as mass transfer, mass accretion, common-envelope evolution, collisions, supernova kicks and angular momentum loss mechanisms are included. In particular, circularization and synchronization of the orbit by tidal interact…
▽ More
We present a rapid binary evolution algorithm that enables modelling of even the most complex binary systems. In addition to all aspects of single star evolution, features such as mass transfer, mass accretion, common-envelope evolution, collisions, supernova kicks and angular momentum loss mechanisms are included. In particular, circularization and synchronization of the orbit by tidal interactions are calculated for convective, radiative and degenerate damping mechanisms. We use this algorithm to study the formation and evolution of various binary systems. We also investigate the effect that tidal friction has on the outcome of binary evolution. Using the rapid binary code, we generate a series of large binary populations and evaluate the formation rate of interesting individual species and events. By comparing the results for populations with and without tidal friction we quantify the hitherto ignored systematic effect of tides and show that modelling of tidal evolution in binary systems is necessary in order to draw accurate conclusions from population synthesis work. Tidal synchronism is important but because orbits generally circularize before Roche-lobe overflow the outcome of the interactions of systems with the same semi-latus rectum is almost independent of eccentricity. It is not necessary to include a distribution of eccentricities in population synthesis of interacting binaries, however, the initial separations should be distributed according to the observed distribution of semi-latera recta rather than periods or semi-major axes.
△ Less
Submitted 14 January, 2002;
originally announced January 2002.
-
The Promiscuous Nature of Stars in Clusters
Authors:
Jarrod R. Hurley,
Michael M. Shara
Abstract:
The recent availability of special purpose computers designed for calculating gravitational interactions of N-bodies at extremely high speed has provided the means to model globular clusters on a star-by-star basis for the first time. By endeavouring to make the N-body codes that operate on these machines as realistic as possible, the addition of stellar evolution being one example, much is bein…
▽ More
The recent availability of special purpose computers designed for calculating gravitational interactions of N-bodies at extremely high speed has provided the means to model globular clusters on a star-by-star basis for the first time. By endeavouring to make the N-body codes that operate on these machines as realistic as possible, the addition of stellar evolution being one example, much is being learnt about the interaction between the star cluster itself and the stars it contains. A fascinating aspect of this research is the ability to follow the orbits of individual stars in detail and to document the formation of observed exotic systems. This has revealed that many stars within a star cluster lead wildly promiscuous lives, interacting, often intimately and in rapid succession, with a variety of neighbours.
△ Less
Submitted 14 January, 2002;
originally announced January 2002.
-
Could Edge-Lit Type Ia Supernovae be Standard Candles
Authors:
Eniko Regos,
Christopher A. Tout,
Dayal Wickramasinghe,
Jarrod R. Hurley,
Onno R. Pols
Abstract:
The progenitors of Type Ia supernovae (SNe Ia) have not been identified. Though they are no longer fashionable we investigate the consequences if a significant number of SNe Ia were edge-lit detonations (ELDs) of carbon/oxygen white dwarfs that have accreted a critical mass of helium. Our best understanding of the Phillips relation between light curve speed and peak luminosity assigns both these…
▽ More
The progenitors of Type Ia supernovae (SNe Ia) have not been identified. Though they are no longer fashionable we investigate the consequences if a significant number of SNe Ia were edge-lit detonations (ELDs) of carbon/oxygen white dwarfs that have accreted a critical mass of helium. Our best understanding of the Phillips relation between light curve speed and peak luminosity assigns both these phenomena to the amount of $^{56}$Ni produced. In ELDs there are two sites of $^{56}$Ni synthesis. If the peak luminosity is determined primarily by the C/O ratio in the core it is primarily a function of its progenitor's initial mass. If the light curve decay speed is determined by the total mass of iron group elements ejected this is a function of the total mass of the ELD at the time of explosion. In general these two masses are correlated and an empirical relation between peak luminosity and light curve shape can be expected. However when we perform population synthesis for progenitors of different metallicities we find a systematic shift in this relation that would make distant SNe Ia fainter than those nearby. The abundances of alpha-rich isotopes, such as $^{44}$Ca, in the solar system indicate that only about 40 per cent of SNe Ia are edge-lit so any systematic effect that could be present would be correspondingly diluted. If we examine only the small subset of ELDs that accrete from a naked helium star, rather than a He white dwarf, the systematic effect disappears.
△ Less
Submitted 23 October, 2002; v1 submitted 14 December, 2001;
originally announced December 2001.