BZPEERPapers, people, evidence and research paths
Menu

Discover

SearchFind papers, people, topics, institutions and opportunities from one place.GraphSee how papers, authors, topics and evidence connect.

Research tools

HomeContinue from the papers, people and topics that matter now.CollectionsTurn saved papers into reading lists, evidence maps and shareable context.ResearchReconstruct the path that led you from search to paper to question.CollaborationMove from a paper, topic or expert request into a focused thread.PublishCreate full-text research articles with metadata, formulas and versions.

Network

ExpertsFind specialists whose work and reviews explain why they can help.CommunitiesJoin research groups organized around papers, topics and shared evidence.

Opportunities

ListingsBrowse roles, calls and collaborations connected to real research context.

Account

Log inReturn to your saved papers, graph views and active threads.Create accountStart saving papers, building a research trail and joining teams.
Log inCreate account

Source author record

Douglas C. Heggie

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

ResearcherUnclaimed source record
astro-ph.GAastro-phastro-ph.SRastro-ph.EPastro-ph.IMastro-ph.HEphysics.comp-ph

Catalog footprint

What is connected

30works
7topics
4close collaborators

Actions

Connect this record

Log in to claim
Open graphBrowse works

Research graph

See the researcher in context

Open full explorer

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

Building this map preview

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

Published work

30 published item(s)

preprint2021arXiv

$l=1$: Weinberg's weakly damped mode in an $N$-body model of a spherical stellar system

Spherical stellar systems such as King models, in which the distribution function is a decreasing function of energy and depends on no other invariant, are stable in the sense of collisionless dynamics. But Weinberg showed, by a clever application of the matrix method of linear stability, that they may be nearly unstable, in the sense of possessing {\sl weakly} damped modes of oscillation. He also demonstrated the presence of such a mode in an $N$-body model by endowing it with initial conditions generated from his perturbative solution. In the present paper we provide evidence for the presence of this same mode in $N$-body simulations of the King $W_0 = 5$ model, in which the initial conditions are generated by the usual Monte Carlo sampling of the King distribution function. It is shown that the oscillation of the density centre correlates with variations in the structure of the system out to a radius of about 1 virial radius, but anticorrelates with variations beyond that radius. Though the oscillations appear to be continually reexcited (presumably by the motions of the particles) we show by calculation of power spectra that Weinberg's estimate of the period (strictly, $2π$ divided by the real part of the eigenfrequency) lies within the range where the power is largest. In addition, however, the power spectrum displays another very prominent feature at shorter periods, around 5 crossing times.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

The Kinematic Richness of Star Clusters - II. Stability of Spherical Anisotropic Models with Rotation

We study the bar instability in collisionless, rotating, anisotropic, stellar systems, using N-body simulations and also the matrix technique for calculation of modes with the perturbed collisionless Boltzmann equation. These methods are applied to spherical systems with an initial Plummer density distribution, but modified kinematically in two ways: the velocity distribution is tangentially anisotropic, using results of Dejonghe, and the system is set in rotation by reversing the velocities of a fraction of stars in various regions of phase space, a la Lynden-Bell. The aim of the N-body simulations is first to survey the parameter space, and, using those results, to identify regions of phase space (by radius and orbital inclination) which have the most important influence on the bar instability. The matrix method is then used to identify the resonant interactions in the system which have the greatest effect on the growth rate of a bar. Complementary series of N-body simulations examine these processes in relation to the evolving frequency distribution and the pattern speed. Finally, the results are synthesised with an existing theoretical framework, and used to consider the old question of constructing a stability criterion.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

The Origin of Chaos in the Orbit of Comet 1P/Halley

According to Munoz-Gutierrez et al. (2015) the orbit of comet 1P/Halley is chaotic with a surprisingly small Lyapunov time scale of order its orbital period. In this work we analyse the origin of chaos in Halley's orbit and the growth of perturbations, in order to get a better understanding of this unusually short time scale. We perform N-body simulations to model Halley's orbit in the Solar System and measure the separation between neighbouring trajectories. To be able to interpret the numerical results, we use a semi-analytical map to demonstrate different growth modes, i.e. linear, oscillatory or exponential, and transitions between these modes. We find the Lyapunov time scale of Halley's orbit to be of order 300 years, which is significantly longer than previous estimates in the literature. This discrepancy could be due to the different methods used to measure the Lyapunov time scale. A surprising result is that next to Jupiter, also encounters with Venus contribute to the exponential growth in the next 3000 years. Finally, we note an interesting application of the sub-linear, oscillatory growth mode to an ensemble of bodies moving through the Solar System. Whereas in the absence of encounters with a third body the ensemble spreads out linearly in time, the accumulation of weak encounters can increase the lifetime of such systems due to the oscillatory behaviour.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Evolution of star clusters on eccentric orbits

We study the evolution of star clusters on circular and eccentric orbits using direct $N$-body simulations. We model clusters with initially $N=8{\rm k}$ and $N=16{\rm k}$ single stars of the same mass, orbiting around a point-mass galaxy. For each orbital eccentricity that we consider, we find the apogalactic radius at which the cluster has the same lifetime as the cluster with the same $N$ on a circular orbit. We show that then, the evolution of bound particle number and half-mass radius is approximately independent of eccentricity. Secondly, when we scale our results to orbits with the same semi-major axis, we find that the lifetimes are, to first order, independent of eccentricity. When the results of Baumgardt and Makino for a singular isothermal halo are scaled in the same way, the lifetime is again independent of eccentricity to first order, suggesting that this result is independent of the Galactic mass profile. From both sets of simulations we empirically derive the higher order dependence of the lifetime on eccentricity. Our results serve as benchmark for theoretical studies of the escape rate from clusters on eccentric orbits. Finally, our results can be useful for generative models for cold streams and cluster evolution models that are confined to spherical symmetry and/or time-independent tides, such as Fokker-Planck models, Monte Carlo models, and (fast) semi-analytic models.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Neutron star natal kicks and the long-term survival of star clusters

We investigate the dynamical evolution of a star cluster in an external tidal field by using N-body simulations, with focus on the effects of the presence or absence of neutron star natal velocity kicks.We show that, even if neutron stars typically represent less than 2% of the total bound mass of a star cluster, their primordial kinematic properties may affect the lifetime of the system by up to almost a factor of four. We interpret this result in the light of two known modes of star cluster dissolution, dominated by either early stellar evolution mass loss or two-body relaxation. The competition between these effects shapes the mass loss profile of star clusters, which may either dissolve abruptly ("jumping"), in the pre-core-collapse phase, or gradually ("skiing"), after having reached core collapse.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Michel Henon's contributions to collisional stellar systems

The theory of star cluster dynamics was a major topic in Hénon's early research career. Here we summarise his contributions under three headings: (i) the Monte Carlo method, (ii) homological evolution of star clusters, and (iii) escape from star clusters. In each case we also trace some aspects of how Hénon's contributions have been developed or applied in subsequent decades up to the present. We also propose that Hénon's work be commemorated by adopting the names "Hénon units" and "Hénon's Principle".

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

MOCCA code for star cluster simulations - III. Stellar-mass black holes in the globular cluster M22

Using a Monte Carlo code, we construct a dynamic evolutionary model of the Galactic globular cluster M22 (NGC6656). The initial conditions are chosen so that, after about 12Gyr of stellar and dynamical evolution, the model is an approximate fit to the surface brightness and velocity dispersion profiles of the cluster, to its mass function, and to the current binary fraction. Depending on the distribution of black hole natal kicks, we predict that the present-day population of stellar-mass black holes ranges from about 40 (no kicks) down to essentially zero (kicks distributed like those of neutron stars). Provided that natal kicks do not eject all new black holes, it is suggested that clusters with a present-day half-mass relaxation time above about 1Gyr are the ones that may still retain an appreciable population of black holes.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Towards an N-body Model for the Globular Cluster M4

This paper describes an N-body model for the dynamical evolution of the nearby globular cluster M4. The initial conditions, with N = 484710 particles, were generated from a published study of this cluster with a Monte Carlo code. With the Monte Carlo code, these initial conditions led, after 12 Gyr of dynamical and stellar evolution, to a model which resembles M4 in terms of its surface brightness and velocity dispersion profiles, and its local luminosity function. Though the N-body model reported here is marred by some errors, its evolution can be compared with that of the published Monte Carlo model, with a result from the synthetic evolution code EMACSS, and with M4 itself.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Dynamical evolution of black hole sub-systems in idealised star clusters

In this paper, globular star clusters which contain a sub-system of stellar-mass black holes (BH) are investigated. This is done by considering two-component models, as these are the simplest approximation of more realistic multi-mass systems, where one component represents the BH population and the other represents all the other stars. These systems are found to undergo a long phase of evolution where the centre of the system is dominated by a dense BH sub-system. After mass segregation has driven most of the BH into a compact sub-system, the evolution of the BH sub-system is found to be influenced by the cluster in which it is contained. The BH sub-system evolves in such a way as to satisfy the energy demands of the whole cluster, just as the core of a one component system must satisfy the energy demands of the whole cluster. The BH sub-system is found to exist for a significant amount of time. It takes approximately 10t_{rh,i}, where t_{rh,i} is the initial half-mass relaxation time, from the formation of the compact BH sub-system up until the time when 90% of the sub-system total mass is lost (which is of order 10^{3} times the half-mass relaxation time of the BH sub-system at its time of formation). Based on theoretical arguments the rate of mass loss from the BH sub-system (\dot{M}_2) is predicted to be -(beta*zeta*M)/(alpha*t_{rh}), where M is the total mass, t_{rh} is the half-mass relaxation time, and alpha, beta, zeta are three dimensionless parameters (see Section 2 for details). An interesting consequence of this is that the rate of mass loss from the BH sub-system is approximately independent of the stellar mass ratio (m_2/m_1) and the total mass ratio (M_2/M_1) (in the range m_2/m_1 >~ 10 and M_2/M_1 ~ 10^{-2}, where m_1, m_2 are the masses of individual low-mass and high-mass particles respectively, and M_1, M_2 are the corresponding total masses).

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Few-Body Modes of Binary Formation in Core Collapse

At the moment of deepest core collapse, a star cluster core contains less than ten stars. This small number makes the traditional treatment of hard binary formation, assuming a homogeneous background density, suspect. In a previous paper, we have found that indeed the conventional wisdom of binary formation, based on three-body encounters, is incorrect. Here we refine that insight, by further dissecting the subsequent steps leading to hard binary formation. For this purpose, we add some analysis tools in order to make the study less subjective. We find that the conventional treatment does remain valid for direct three-body scattering, but fails for resonant three-body scattering. Especially democratic resonance scattering, which forms an important part of the analytical theory of three-body binary formation, takes too much space and time to be approximated as being isolated, in the context of a cluster core around core collapse. We conclude that, while three-body encounters can be analytically approximated as isolated, subsequent strong perturbations typically occur whenever those encounters give rise to democratic resonances. We present analytical estimates postdicting our numerical results. If we only had been a bit more clever, we could have predicted this qualitative behaviour.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

MOCCA Code for Star Cluster Simulations - II. Comparison with N-body Simulations

We describe a major upgrade of a Monte Carlo code which has previously been used for many studies of dense star clusters. We outline the steps needed in order to calibrate the results of the new Monte Carlo code against $N$-body simulations for large $N$ systems, up to $N=200000$. The new version of the Monte Carlo code (called MOCCA), in addition to the features of the old version, incorporates the direct Fewbody integrator (Fregeau et al. 2004) for three- and four-body interactions, and a new treatment of the escape process based on Fukushige & Heggie (2000). Now stars which fulfil the escape criterion are not removed immediately, but can stay in the system for a certain time which depends on the excess of the energy of a star above the escape energy. They are called potential escapers. With the addition of the Fewbody integrator the code can follow all interaction channels which are important for the rate of creation of various types of objects observed in star clusters, and ensures that the energy generation by binaries is treated in a manner similar to the $N$-body model. There are at most three new parameters which have to be adjusted against $N$-body simulations for large $N$: two (or one, depending on the chosen approach) connected with the escape process, and one responsible for the determination of the interaction probabilities. The values adopted for the free parameters have at most a weak dependence on $N$. They allow MOCCA to reproduce $N$-body results with reasonable precision, not only for the rate of cluster evolution and the cluster mass distribution, but also for the detailed distributions of mass and binding energy of binaries. Additionally, the code can follow the rate of formation of blue stragglers and black hole - black hole binaries.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

On black hole sub-systems in idealized nuclear star clusters

Recent observational evidence, numerical simulations and theoretical arguments seem to indicate that stellar mass black holes (BH) persist in stellar systems such as globular star clusters for much longer than previously thought. Previously, theory implied that the BH would segregate into a compact system with short dynamical time scales, and that the BH would escape long before the present. But stellar systems can exist in a state of balanced evolution, where the energy generated in the core is regulated by the process of two-body relaxation in the bulk of the system. If the system has a centrally concentrated BH subsystem and there is no massive central BH, then the energy is generated by three-body encounters in the core of this subsystem. Therefore, in balanced evolution, the evolution of the BH subsystem is regulated by the much longer time scales of the host system. In the present paper the implications of these results for idealized nuclear star clusters (NSC) are discussed. Though previous theory implied that BH would be almost absent from many NSC -- those with relatively short dynamical time scales -- it is argued here that, based on the results of idealized models, many such NSC could still be host to substantial BH subsystems.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Gravothermal oscillations in multi-component models of star clusters

In this paper, gravothermal oscillations are investigated in multi-component star clusters which have power law initial mass functions (IMF). For the power law IMFs, the minimum masses ($m_{min}$) were fixed and three different maximum stellar masses ($m_{max}$) were used along with different power-law exponents ($α$) ranging from 0 to -2.35 (Salpeter). The critical number of stars at which gravothermal oscillations first appear with increasing $N$ was found using the multi-component gas code SPEDI. The total mass ($M_{tot}$) is seen to give an approximate stability condition for power law IMFs with fixed values of $m_{max}$ and $m_{min}$ independent of $α$. The value $M_{tot}/m_{max} \simeq 12000$ is shown to give an approximate stability condition which is also independent of $m_{max}$, though the critical value is somewhat higher for the steepest IMF that was studied. For appropriately chosen cases, direct N-body runs were carried out in order to check the results obtained from SPEDI. Finally, evidence of the gravothermal nature of the oscillations found in the N-body runs is presented.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

The Tidal Tails of 47 Tucanae

The Galactic globular cluster 47 Tucanae (47 Tuc) shows a rare increase in its velocity dispersion profile at large radii, indicative of energetic, yet bound, stars at large radii dominating the velocity dispersion and, potentially, of ongoing evaporation. Escaping stars will form tidal tails, as seen with several Galactic globular clusters, however, the tidal tails of 47 Tuc are yet to be uncovered. We model these tails of 47 Tuc using the most accurate input data available, with the specific aim of determining their locations, as well as the densities of the epicyclic overdensities within the tails. The overdensities from our models show an increase of 3-4% above the Galactic background and, therefore, should be easily detectable using matched filtering techniques. We find that the most influential parameter with regard to both the locations and densities of the epicyclic overdensities is the Heliocentric distance to the cluster. Hence, uncovering these tidal features observationally will contribute greatly to the ongoing problem of determining the distance to 47 Tuc, tightly constraining the distance of the cluster independent of other methods. Using our streakline method for determining the locations of the tidal tails and their overdensities, we show how, in principle, the shape and extent of the tidal tails of any Galactic globular cluster can be determined without resorting to computationally expensive N-body simulations.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Gravothermal oscillations in two-component models of star clusters

In this paper, gravothermal oscillations are investigated in two-component clusters with a range of different stellar mass ratios and total component mass ratios. The critical number of stars at which gravothermal oscillations first appeared is found using a gas code. The nature of the oscillations is investigated and it is shown that the oscillations can be understood by focusing on the behaviour of the heavier component, because of mass segregation. It is argued that, during each oscillation, the re-collapse of the cluster begins at larger radii while the core is still expanding. This re-collapse can halt and reverse a gravothermally driven expansion. This material outside the core contracts because it is losing energy both to the cool expanding core and to the material at larger radii. The core collapse times for each model are also found and discussed. For an appropriately chosen case, direct N -body runs were carried out, in order to check the results obtained from the gas model, including evidence of the gravothermal nature of the oscillations and the temperature inversion that drives the expansion.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

More on the structure of tidal tails

We investigate the epicyclic motion of stars escaping from star clusters. Using streaklines, we visualise the path of escaping stars and show how epicyclic motion leads to over- and underdensities in tidal tails of star clusters moving on circular and eccentric orbits about a galaxy. Additionally, we investigate the effect of the cluster mass on the tidal tails, by showing that their structure is better matched when the perturbing effect of the cluster mass is included. By adjusting streaklines to results of N-body computations we can accurately and quickly reproduce all observed substructure, especially the streaky features often found in simulations which may be interpreted in observations as multiple tidal tails. Hence, we can rule out tidal shocks as the origin of such substructures. Finally, from the adjusted streakline parameters we can verify that for the star clusters we studied escape mainly happens from the tidal radius of the cluster, given by x_L = (GM/(Ω^2-\partial^2Φ/\partial R^2))^{1/3}. We find, however, that there is another limiting radius, the "edge" radius, which gives the smallest radius from which a star can escape during one cluster orbit about the galaxy. For eccentric cluster orbits the edge radius shrinks with increasing orbital eccentricity (for fixed apocentric distance) but is always significantly larger than the respective perigalactic tidal radius. In fact, the edge radii of the clusters we investigated, which are extended and tidally filling, agree well with their (fitted) King radii, which may indicate a fundamental connection between these two quantities.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Problems of Collisional Stellar Dynamics

The discovery of dynamical friction was Chandrasekhar's best known contribution to the theory of stellar dynamics, but his work ranged from the few-body problem to the limit of large N (in effect, galaxies). Much of this work was summarised in the text "Principles of Stellar Dynamics" (Chandrasekhar 1942, 1960), which ranges from a precise calculation of the time of relaxation, through a long analysis of galaxy models, to the behaviour of star clusters in tidal fields. The later edition also includes the work on dynamical friction and related issues. In this review we focus on progress in the collisional aspects of these problems, i.e. those where few-body interactions play a dominant role, and so we omit further discussion of galaxy dynamics. But we try to link Chandrasekhar's fundamental discoveries in collisional problems with the progress that has been made in the 50 years since the publication of the enlarged edition.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

The life cycle of star cluster in a tidal field

The evolution of globular clusters due to 2-body relaxation results in an outward flow of energy and at some stage all clusters need a central energy source to sustain their evolution. Henon provided the insight that we do not need to know the details of the energy production in order to understand the relaxation-driven evolution of the cluster, at least outside the core. He provided two self-similar solutions for the evolution of clusters based on the view that the cluster as a whole determines the amount of energy that is produced in the core: steady expansion for isolated clusters, and homologous contraction for clusters evaporating in a tidal field. We combine these models: the half-mass radius increases during the first half of the evolution, and decreases in the second half; while the escape rate approaches a constant value set by the tidal field. We refer to these phases as `expansion dominated' and `evaporation dominated'. These simple analytical solutions immediately allow us to construct evolutionary tracks and isochrones in terms of cluster half-mass density, cluster mass and galacto-centric radius. From a comparison to the Milky Way globular clusters we find that roughly 1/3 of them are in the second, evaporation-dominated phase and for these clusters the density inside the half-mass radius varies with the galactocentric distance R as rho_h ~ 1/R^2. The remaining 2/3 are still in the first, expansion-dominated phase and their isochrones follow the environment-independent scaling rho_h ~ M^2; that is, a constant relaxation time-scale. We find substantial agreement between Milky Way globular cluster parameters and the isochrones, which suggests that there is, as Henon suggested, a balance between the flow of energy and the central energy production for almost all globular clusters.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Monte Carlo Simulations of Star Clusters - VII. The globular cluster 47 Tuc

We describe Monte Carlo models for the dynamical evolution of the massive globular cluster 47 Tuc (NGC 104). The code includes treatments of two-body relaxation, most kinds of three- and four-body interactions involving primordial binaries and those formed dynamically, the Galactic tide, and the internal evolution of both single and binary stars. We arrive at a set of initial parameters for the cluster which, after 12Gyr of evolution, gives a model with a fairly satisfactory match to surface brightness and density profiles, the velocity dispersion profile, the luminosity function in two fields, and the acceleration of pulsars. Our models appear to require a relatively steep initial mass function for stars above about turnoff, with an index of about 2.8 (where the Salpeter mass function has an index of 2.35), and a relatively flat initial mass function (index about 0.4) for the lower main sequence. According to the model, the current mass is estimated at 0.9 million solar masses, of which about 34% consists of remnants. We find that primordial binaries are gradually taking over from mass loss by stellar evolution as the main dynamical driver of the core. Despite the high concentration of the cluster, core collapse will take at least another 20Gyr.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

The effects of fly-bys on planetary systems

Most of the observed extrasolar planets are found on tight and often eccentric orbits. The high eccentricities are not easily explained by planet-formation models, which predict that planets should be on rather circular orbits. Here we explore whether fly-bys involving planetary systems with properties similar to those of the gas giants in the solar system, can produce planets with properties similar to the observed planets. Using numerical simulations, we show that fly-bys can cause the immediate ejection of planets, and sometimes also lead to the capture of one or more planets by the intruder. More common, however, is that fly-bys only perturb the orbits of planets, sometimes leaving the system in an unstable state. Over time-scales of a few million to several hundred million years after the fly-by, this perturbation can trigger planet-planet scatterings, leading to the ejection of one or more planets. For example, in the case of the four gas giants of the solar system, the fraction of systems from which at least one planet is ejected more than doubles in 10^8 years after the fly-by. The remaining planets are often left on more eccentric orbits, similar to the eccentricities of the observed extrasolar planets. We combine our results of how fly-bys effect solar-system-like planetary systems, with the rate at which encounters in young stellar clusters occur. For example, we measure the effects of fly-bys on the four gas giants in the solar system. We find, that for such systems, between 5 and 15 per cent suffer ejections of planets in 10^8 years after fly-bys in typical open clusters. Thus, encounters in young stellar clusters can significantly alter the properties of any planets orbiting stars in clusters. As a large fraction of stars which populate the solar neighbourhood form in stellar clusters, encounters can significantly affect the properties of the observed extrasolar planets.

0 citations0 reviewsNo code linkedOpen access
preprint2009arXiv

1Gyr in the Life of the Globular Cluster NGC 6397

M4 and NGC 6397 are two very similar galactic globular clusters, which differ mainly in their surface brightness profile. M4 has a classic King-like profile, whereas NGC 6397 has a more concentrated profile, which is often interpreted as that of a post-core collapse cluster. Heggie & Giersz (2008), however, found that M4 is also a post-core collapse cluster, and Giersz & Heggie (2009) concluded that the main reason for the difference between the two surface brightness profiles is fluctuations. This conclusion was reached on the basis of Monte Carlo models, however, and in the present Letter we verify that similar fluctuations occur in N-body models. The models were initialised by generating initial conditions from the Monte Carlo model of NGC 6397 at the simulated age of 12Gyr, and one was followed for 1Gyr. The new models help to clarify the nature of the fluctuations, which have the nature of semi-regular oscillations with a time scale of order 10^8yr. They are influenced by the dynamical role which is played by primordial binaries in the evolution of the core.

0 citations0 reviewsNo code linkedOpen access
preprint2007arXiv

The core radius of a star cluster containing a massive black hole

(Simplified) We present a theoretical framework which establishes how the core radius of a star cluster varies with the mass of an assumed central black hole. Our result is that the ratio of core to half-mass radius varies as the 3/4 power of the ratio of the black hole to cluster mass, when the system is well relaxed. The theory compares favourably with a number of simulations of this problem, which extend to black hole masses of order 10% of the cluster mass. Though strictly limited as yet to clusters with stars of equal mass, our conclusion strengthens the view that clusters with large core radii are the most promising candidates in which to find a massive black hole.

0 citations0 reviewsNo code linkedOpen access
preprint1996arXiv

The Effect of Encounters on the Eccentricity of Binaries in Clusters

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

0 citations0 reviewsNo code linkedOpen access
preprint1995arXiv

Statistics of N-Body Simulations. III. Unequal Masses

We describe results from large numbers of $N$-body simulations containing from $250$ to $1000$ stars each. The distribution of stellar masses is a power law, and the systems are isolated. While the collapse of the core exhibits the expected segregation of different masses, we find that the post-collapse evolution is, at a first approximation, homologous. This is quite surprising because there is no reason for supposing that mass segregation should not continue to have a substantial effect on the evolution of the cluster. In fact the spatial distribution of the mean stellar mass is nearly static throughout the post-collapse regime, except for the overall expansion of the systems, and this helps to explain why the post-collapse evolution is nearly self-similar. Self-similarity is also exhibited by the distribution of anisotropy and the profile of departures from equipartition, which show little change during the post-collapse phase. The departures from energy equipartition and isotropy are small in the core and increase with radius. During post-collapse evolution massive stars (mainly) are removed from the system by binary activity. This effect dominates the preferential escape of low-mass stars due to standard two-body relaxation processes.

0 citations0 reviewsNo code linkedOpen access
preprint1995arXiv

THE ORBITAL ECCENTRICITIES OF BINARY MILLISECOND PULSARS IN GLOBULAR CLUSTERS

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

0 citations0 reviewsNo code linkedOpen access
preprint1994arXiv

Pre-Collapse Evolution of Galactic Globular Clusters

This paper is concerned with collisionless aspects of the early evolution of model star clusters. The effects of mass loss through stellar evolution and of a steady tidal field are modelled using $N$-body simulations. Our results (which depend on the assumed initial structure and the mass spectrum) agree qualitatively with those of Chernoff \& Weinberg (1990), who used a Fokker-Planck model with a spherically symmetric tidal cutoff. For those systems which are disrupted, the lifetime to disruption generally exceeds that found by Chernoff \& Weinberg, sometimes by as much as an order of magnitude. Because we do not model collisional effects correctly we cannot establish the fate of the survivors. In terms of theoretical interpretation, we find that tidal disruption must be understood as a loss of {\sl equilibrium}, and not a loss of {\sl stability}, as is sometimes stated.

0 citations0 reviewsNo code linkedOpen access
preprint1994arXiv

Statistics of N-Body Simulations. II. Equal Masses after Core Collapse

This paper presents and analyses statistical results from a large number of $N$-body simulations of isolated systems with equal masses, in which $250\le N\le 2000$. It concentrates on the phase starting around the end of core collapse, when binaries start to play a crucial role. Interactions of hard binaries are consistent with the Spitzer (1987) cross section. The evolution of the half-mass radius after core collapse nearly follows classical theory, The evolution of the inner parts of the system around the time of core bounce is consistent with these simplified models provided that the continuous production of energy, as usually assumed, is replaced by a model of {\sl stochastic} energy production. Similarly, post-collapse evolution of the core requires a modest recalibration of the coefficient of energy generation, especially for small $N$. The distributions of velocity dispersion and anisotropy become remarkably homologous soon after core bounce. The bound mass of the systems very nearly follows a power law with time. A small number of escapers, presumed to be those associated with binary activity, dominate the energy which is carried off. The ``internal" energy of escaping binaries is consistent with theoretical expectations.

0 citations0 reviewsNo code linkedOpen access
preprint1993arXiv

Gravothermal Expansion in an $N$-Body System

This paper describes the numerical evolution of an $N$-body system with a slight ``temperature inversion''; i.e. the maximum velocity dispersion occurs not at the centre but further out. Fluid models predict that the core of such a system expands on a time-scale of thousands of central relaxation times, and here this behaviour is qualitatively confirmed for an $N$-body system of over 3000 bodies. With certain qualifications, this demonstrates the existence in N-body systems of one of the fundamental mechanisms which, in fluid models, drive the gravothermal oscillations discovered by Bettwieser & Sugimoto.

0 citations0 reviewsNo code linkedOpen access
preprint1993arXiv

Star Clusters

This review concentrates almost entirely on globular star clusters. It emphasises the increasing realisation that few of the traditional problems of star cluster astronomy can be studied in isolation: the influence of the Galaxy affects dynamical evolution deep in the core, and the spectrum of stellar masses; in turn the evolution of the core determines the highest stellar densities, and the rate of encounters. In this way external tidal effects indirectly influence the formation and evolution of blue stragglers, binary pulsars, X-ray sources, etc. More controversially, the stellar density appears to influence the relative distribution of normal stars. In the opposite sense, the evolution of individual stars governs much of the early dynamics of a globular cluster, and the existence of large numbers of primordial binary stars has changed important details of our picture of the dynamical evolution. New computational tools which will become available in the next few years will help dynamical theorists to address these questions.

0 citations0 reviewsNo code linkedOpen access
preprint1993arXiv

Statistics of N-body Simulations I. Equal Masess Before Core Collapse

We study the dynamical evolution of idealised stellar systems by averaging results from many $N$-body simulations, each having modest numbers of stars. For isolated systems with stars of uniform mass, we discuss aspects of evolution up to the point of core collapse: relaxation and its $N$-dependence, the evolution of the density profile, the development of the velocity dispersion and anisotropy, and the rate of stellar escape.

0 citations0 reviewsNo code linkedOpen access