Showing 1–50 of 80 results for author: Heggie, D C

Dynamics of star clusters with tangentially anisotropic velocity distribution

Authors: Václav Pavlík, Douglas C. Heggie, Anna Lisa Varri, Enrico Vesperini

Abstract: Recent high-precision observations with HST and Gaia enabled new investigations of the internal kinematics of star clusters (SCs) and the dependence of kinematic properties on the stellar mass. These studies raised new questions about the dynamical evolution of self-gravitating stellar systems. We aim to develop a more complete theoretical understanding of how various kinematical properties of sta… ▽ More Recent high-precision observations with HST and Gaia enabled new investigations of the internal kinematics of star clusters (SCs) and the dependence of kinematic properties on the stellar mass. These studies raised new questions about the dynamical evolution of self-gravitating stellar systems. We aim to develop a more complete theoretical understanding of how various kinematical properties of stars affect the global dynamical development of their host SCs. We perform N-body simulations of SCs with isotropic, radially anisotropic and tangentially anisotropic initial velocity distributions. We also study the effect of an external Galactic tidal field. We find three main results. First, compared to the conventional, isotropic case, the relaxation processes are accelerated in the tangentially anisotropic models and, in agreement with our previous investigations, slower in the radially anisotropic ones. This leads to, e.g., more rapid mass segregation in the central regions of the tangential models or their earlier core collapse. Second, although all SCs become isotropic in the inner regions after several relaxation times, we observe differences in the anisotropy profile evolution in the outer cluster regions - all tidally filling models gain tangential anisotropy there while the underfilling models become radially anisotropic. Third, we observe different rates of evolution towards energy equipartition (EEP). While all SCs evolve towards EEP in their inner regions (regardless of the filling factor), the outer regions of the tangentially anisotropic and isotropic models are evolving to an "inverted" EEP (i.e., the high-mass stars having higher velocity dispersion than the low-mass ones). The extent (both spatial and temporal) of this inversion can be attributed to the initial velocity anisotropy - it grows with increasing tangential anisotropy and decreases as the radial anisotropy rises. △ Less

Submitted 29 May, 2024; originally announced May 2024.

Comments: 13 pages, 13 figures, 1 table, 1 appendix. Accepted for publication in A&A

Punctuated Chaos and Indeterminism in Self-gravitating Many-body Systems

Authors: Tjarda C. N. Boekholt, Simon F. Portegies Zwart, Douglas C. Heggie

Abstract: Dynamical chaos is a fundamental manifestation of gravity in astrophysical, many-body systems. The spectrum of Lyapunov exponents quantifies the associated exponential response to small perturbations. Analytical derivations of these exponents are critical for understanding the stability and predictability of observed systems. This essay presents a new model for chaos in systems with eccentric and… ▽ More Dynamical chaos is a fundamental manifestation of gravity in astrophysical, many-body systems. The spectrum of Lyapunov exponents quantifies the associated exponential response to small perturbations. Analytical derivations of these exponents are critical for understanding the stability and predictability of observed systems. This essay presents a new model for chaos in systems with eccentric and crossing orbits. Here, exponential divergence is not a continuous process but rather the cumulative effect of an ever-increasing linear response driven by discrete events at regular intervals, i.e., punctuated chaos. We show that long-lived systems with punctuated chaos can magnify Planck length perturbations to astronomical scales within their lifetime, rendering them fundamentally indeterministic. △ Less

Submitted 28 August, 2023; originally announced August 2023.

Comments: Essay written for the Gravity Research Foundation 2023 Awards for Essays on Gravitation (Honorable Mention). Accepted by International Journal of Modern Physics D

The Effect of Stars on the Dark Matter Spike Around a Black Hole: A Tale of Two Treatments

Authors: Stuart L. Shapiro, Douglas C. Heggie

Abstract: We revisit the role that gravitational scattering off stars plays in establishing the steady-state distribution of collisionless dark matter (DM) around a massive black hole (BH). This is a physically interesting problem that has potentially observable signatures, such as $γ-$rays from DM annihilation in a density spike. The system serves as a laboratory for comparing two different dynamical appro… ▽ More We revisit the role that gravitational scattering off stars plays in establishing the steady-state distribution of collisionless dark matter (DM) around a massive black hole (BH). This is a physically interesting problem that has potentially observable signatures, such as $γ-$rays from DM annihilation in a density spike. The system serves as a laboratory for comparing two different dynamical approaches, both of which have been widely used: a Fokker-Planck treatment and a two-component conduction fluid treatment. In our Fokker-Planck analysis we extend a previous analytic model to account for a nonzero flux of DM particles into the BH, as well as a cut-off in the distribution function near the BH due to relativistic effects or, further out, possible DM annihilation. In our two-fluid analysis, following an approximate analytic treatment, we recast the equations as a "heated Bondi accretion" problem and solve the equations numerically without approximation. While both the Fokker-Planck and two-fluid methods yield basically the same DM density and velocity dispersion profiles away from the boundaries in the spike interior, there are other differences, especially the determination of the DM accretion rate. We discuss limitations of the two treatments, including the assumption of an isotropic velocity dispersion. △ Less

Submitted 16 September, 2022; originally announced September 2022.

Comments: 12 pages, 6 figures

Journal ref: Phys. Rev. D 106, 043018 (2022)

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

Authors: Philip G. Breen, Simon Rozier, Douglas C. Heggie, Anna Lisa Varri

Abstract: 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 aniso… ▽ More 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. △ Less

Submitted 8 February, 2021; originally announced February 2021.

Comments: 17 pages, 18 figures, accepted for publication in MNRAS

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

Authors: Douglas C. Heggie, Philip G. Breen, Anna Lisa Varri

Abstract: 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… ▽ More 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. △ Less

Submitted 2 March, 2021; v1 submitted 3 February, 2020; originally announced February 2020.

Comments: 9 pages, 6 figures, footnote added in Sec.3.2, no other change to the version published in MNRAS

Spherical models of star clusters with potential escapers

Authors: Ian Claydon, Mark Gieles, Anna Lisa Varri, Douglas C. Heggie, Alice Zocchi

Abstract: An increasing number of observations of the outer regions of globular clusters (GCs) have shown a flattening of the velocity dispersion profile and an extended surface density profile. Formation scenarios of GCs can lead to different explanations of these peculiarities, therefore the dynamics of stars in the outskirts of GCs are an important tool in tracing back the evolutionary history and format… ▽ More An increasing number of observations of the outer regions of globular clusters (GCs) have shown a flattening of the velocity dispersion profile and an extended surface density profile. Formation scenarios of GCs can lead to different explanations of these peculiarities, therefore the dynamics of stars in the outskirts of GCs are an important tool in tracing back the evolutionary history and formation of star clusters. One possible explanation for these features is that GCs are embedded in dark matter halos. Alternatively, these features are the result of a population of energetically unbound stars that can be spatially trapped within the cluster, known as potential escapers (PEs). We present a prescription for the contribution of these energetically unbound members to a family of self-consistent, distribution function-based models, which, for brevity, we call the Spherical Potential Escapers Stitched (SPES) models. We show that, when fitting to mock data of bound and unbound stars from an N-body model of a tidally-limited star cluster, the SPES models correctly reproduce the density and velocity dispersion profiles up to the Jacobi radius, and they are able to recover the value of the Jacobi radius itself to within 20%. We also provide a comparison to the number density and velocity dispersion profiles of the Galactic cluster 47 Tucanae. Such a case offers a proof of concept that an appropriate modelling of PEs is essential to accurately interpret Gaia data in the outskirts of GCs, and, in turn, to formulate meaningful present-day constraints for GC formation scenarios in the early universe. △ Less

Submitted 15 March, 2019; v1 submitted 14 March, 2019; originally announced March 2019.

Comments: 16 pages, 11 figures, submitted to MNRAS

Mapping the stability of stellar rotating spheres via linear response theory

Authors: Simon Rozier, Jean-Baptiste Fouvry, Philip G. Breen, Anna Lisa Varri, Christophe Pichon, Douglas C. Heggie

Abstract: Rotation is ubiquitous in the Universe, and recent kinematic surveys have shown that early type galaxies and globular clusters are no exception. Yet the linear response of spheroidal rotating stellar systems has seldom been studied. This paper takes a step in this direction by considering the behaviour of spherically symmetric systems with differential rotation. Specifically, the stability of seve… ▽ More Rotation is ubiquitous in the Universe, and recent kinematic surveys have shown that early type galaxies and globular clusters are no exception. Yet the linear response of spheroidal rotating stellar systems has seldom been studied. This paper takes a step in this direction by considering the behaviour of spherically symmetric systems with differential rotation. Specifically, the stability of several sequences of Plummer spheres is investigated, in which the total angular momentum, as well as the degree and flavour of anisotropy in the velocity space are varied. To that end, the response matrix method is customised to spherical rotating equilibria. The shapes, pattern speeds and growth rates of the systems' unstable modes are computed. Detailed comparisons to appropriate N-body measurements are also presented. The marginal stability boundary is charted in the parameter space of velocity anisotropy and rotation rate. When rotation is introduced, two sequences of growing modes are identified corresponding to radially and tangentially-biased anisotropic spheres respectively. For radially anisotropic spheres, growing modes occur on two intersecting surfaces (in the parameter space of anisotropy and rotation), which correspond to fast and slow modes, depending on the net rotation rate. Generalised, approximate stability criteria are finally presented. △ Less

Submitted 25 February, 2019; originally announced February 2019.

Comments: 18 pages, 20 figures. Submitted to MNRAS

The Kinematic Richness of Star Clusters I. Isolated Spherical Models with Primordial Anisotropy

Authors: Philip G. Breen, Anna Lisa Varri, Douglas C. Heggie

Abstract: We investigate the dynamical evolution of isolated equal-mass star cluster models by means of direct N-body simulations, primarily focusing on the effects of the presence of primordial anisotropy in the velocity space. We found evidence of the existence of a monotonic relationship between the moment of core collapse and the amount and flavour of anisotropy in the stellar system. Specifically, equi… ▽ More We investigate the dynamical evolution of isolated equal-mass star cluster models by means of direct N-body simulations, primarily focusing on the effects of the presence of primordial anisotropy in the velocity space. We found evidence of the existence of a monotonic relationship between the moment of core collapse and the amount and flavour of anisotropy in the stellar system. Specifically, equilibria characterised by the same initial structural properties (Plummer density profile) and with different degrees of tangentially-biased (radially-biased) anisotropy, reach core collapse earlier (later) than isotropic models. We interpret this result in light of an accelerated (delayed) phase of the early evolution of collisional stellar systems "anisotropic-response"), which we have characterised both in terms of the evolution of the velocity moments and of a fluid model of two-body relaxation. For the case of the most tangentially anisotropic model the initial phase of evolution involves a catastrophic collapse of the inner part of the system which continues until an isotropic velocity distribution is reached. This study represents a first step towards a comprehensive investigation of the role played by kinematic richness in the long-term dynamical evolution of collisional systems. △ Less

Submitted 10 July, 2017; originally announced July 2017.

Comments: 13 pages, 12 figures, accepted for publication in MNRAS

An Approximate Analytic Model of a Star Cluster with Potential Escapers

Authors: Kathryne J. Daniel, Douglas C. Heggie, Anna Lisa Varri

Abstract: In the context of a star cluster moving on a circular galactic orbit, a "potential escaper" is a cluster star that has orbital energy greater than the escape energy, and yet is confined within the Jacobi radius of the stellar system. On the other hand analytic models of stellar clusters typically have a truncation energy equal to the cluster escape energy, and therefore explicitly exclude these en… ▽ More In the context of a star cluster moving on a circular galactic orbit, a "potential escaper" is a cluster star that has orbital energy greater than the escape energy, and yet is confined within the Jacobi radius of the stellar system. On the other hand analytic models of stellar clusters typically have a truncation energy equal to the cluster escape energy, and therefore explicitly exclude these energetically unbound stars. Starting from the landmark analysis performed by Henon of periodic orbits of the circular Hill equations, we present a numerical exploration of the population of "non-escapers", defined here as those stars which remain within two Jacobi radii for several galactic periods, with energy above the escape energy. We show that they can be characterised by the Jacobi integral and two further approximate integrals, which are based on perturbation theory and ideas drawn from Lidov-Kozai theory. Finally we use these results to construct an approximate analytic model that includes a phase space description of a population resembling that of potential escapers, in addition to the usual bound population. △ Less

Submitted 13 March, 2017; originally announced March 2017.

Comments: 23 pages, 19 figures, 1 table

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

Authors: Tjarda C. N. Boekholt, F. Inti Pelupessy, Douglas C. Heggie, Simon F. Portegies Zwart

Abstract: 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 Syste… ▽ More 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. △ Less

Submitted 22 June, 2016; originally announced June 2016.

Comments: Accepted for MNRAS, 10 pages, 6 figures

Evolution of star clusters on eccentric orbits

Authors: Maxwell Xu Cai, Mark Gieles, Douglas C. Heggie, Anna Lisa Varri

Abstract: 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… ▽ More 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. △ Less

Submitted 5 October, 2015; originally announced October 2015.

Comments: 7 pages, 8 figures, accepted for publication in the Monthly Notices of the Royal Astronomical Society

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

Authors: Filippo Contenta, Anna Lisa Varri, Douglas C. Heggie

Abstract: 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… ▽ More 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. △ Less

Submitted 9 February, 2015; originally announced February 2015.

Comments: 5 pages, 5 figures, 1 table. Accepted for publication in MNRAS Letters

Michel Henon's contributions to collisional stellar systems

Authors: Douglas C. Heggie

Abstract: 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… ▽ More 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". △ Less

Submitted 21 November, 2014; v1 submitted 17 November, 2014; originally announced November 2014.

Comments: Invited talk at the "hommage a Michel Henon" conference held at the "Institut Henri Poincare" Paris, GRAVASCO trimestre, autumn 2013. Proceeding by Hermann Press, editors: Jean-Michel Alimi, Roya Mohayaee, Jerome Perez. Videos of all talks are available at: https://www.youtube.com/playlist?list=PL9kd4mpdvWcBLrN8u04IZW6-akNgLG7-_ Improved abstract and two figures in this version

Towards an N-body Model for the Globular Cluster M4

Authors: Douglas C. Heggie

Abstract: 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… ▽ More 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. △ Less

Submitted 19 September, 2014; originally announced September 2014.

Comments: 10 pages, 12 figures; MNRAS, in press

M4 Core Project with HST - III. Search for variable stars in the primary field

Authors: V. Nascimbeni, L. R. Bedin, D. C. Heggie, M. van den Berg, M. Giersz, G. Piotto, K. Brogaard, A. Bellini, A. P. Milone, R. M. Rich, D. Pooley, J. Anderson, L. Ubeda, S. Ortolani, L. Malavolta, A. Cunial, A. Pietrinferni

Abstract: We present the results of a photometric search for variable stars in the core of the Galactic globular cluster M4. The input data are a large and unprecedented set of deep Hubble Space Telescope WFC3 images (large program GO-12911; 120 orbits allocated), primarily aimed at probing binaries with massive companions by detecting their astrometric wobbles. Though these data were not optimised to carry… ▽ More We present the results of a photometric search for variable stars in the core of the Galactic globular cluster M4. The input data are a large and unprecedented set of deep Hubble Space Telescope WFC3 images (large program GO-12911; 120 orbits allocated), primarily aimed at probing binaries with massive companions by detecting their astrometric wobbles. Though these data were not optimised to carry out a time-resolved photometric survey, their exquisite precision, spatial resolution and dynamic range enabled us to firmly detect 38 variable stars, of which 20 were previously unpublished. They include 19 cluster-member eclipsing binaries (confirming the large binary fraction of M4), RR Lyrae, and objects with known X-ray counterparts. We improved and revised the parameters of some among published variables. △ Less

Submitted 7 May, 2014; originally announced May 2014.

Comments: 11 pages, 5 figures, 2 tables. Accepted for publication in MNRAS

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

Authors: Douglas C. Heggie, Mirek Giersz

Abstract: 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 distr… ▽ More 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. △ Less

Submitted 15 January, 2014; originally announced January 2014.

Comments: 9 pages, 10 figures, 2 tables; accepted for publication in Monthly Notices of the Royal Astronomical Society

The M4 Core Project with HST -- II. Multiple Stellar Populations at the Bottom of the Main Sequence

Authors: A. P. Milone, A. F. Marino, L. R. Bedin, G. Piotto, S. Cassisi, A. Dieball, J. Anderson, H. Jerjen, M. Asplund, A. Bellini, K. Brogaard, A. Dotter, M. Giersz, D. C. Heggie, C. Knigge, R. M. Rich, M. van den Berg, R. Buonanno

Abstract: The M4 Core Project with HST is designed to exploit the Hubble Space Telescope to investigate the central regions of M4, the Globular Cluster closest to the Sun. In this paper we combine optical and near-infrared photometry to study multiple stellar populations in M4. We detected two sequences of M-dwarfs containing ~38% (MS_I) and ~62% (MS_II) of MS stars below the main-sequence (MS) knee. We com… ▽ More The M4 Core Project with HST is designed to exploit the Hubble Space Telescope to investigate the central regions of M4, the Globular Cluster closest to the Sun. In this paper we combine optical and near-infrared photometry to study multiple stellar populations in M4. We detected two sequences of M-dwarfs containing ~38% (MS_I) and ~62% (MS_II) of MS stars below the main-sequence (MS) knee. We compare our observations with those of NGC2808, which is the only other GCs where multiple MSs of very low-mass stars have been studied to date. We calculate synthetic spectra for M-dwarfs, assuming the chemical composition mixture inferred from spectroscopic studies of stellar populations along the red giant branch, and different Helium abundances, and we compare predicted and observed colors. Observations are consistent with two populations, one with primordial abundance and another with enhanced nitrogen and depleted oxygen. △ Less

Submitted 6 January, 2014; originally announced January 2014.

Comments: 8 pages, 6 figure, accepted for publication in MNRAS

The M4 Core Project with HST --- I. Overview and First-Epoch

Authors: L. R. Bedin, J. Anderson, D. C. Heggie, G. Piotto, A. P. Milone, M. Giersz, V. Nascimbeni, A. Bellini, R. M. Rich, M. van den Berg, D. Pooley, K. Brogaard, S. Ortolani, L. Malavolta, L. Ubeda, A. F. Marino, .

Abstract: We present an overview of the ongoing Hubble Space Telescope large program GO-12911. The program is focused on the core of M4, the nearest Galactic globular cluster, and the observations are designed to constrain the number of binaries with massive companions (black holes, neutron stars, or white dwarfs) by measuring the ``wobble'' of the luminous (main-sequence) companion around the center of mas… ▽ More We present an overview of the ongoing Hubble Space Telescope large program GO-12911. The program is focused on the core of M4, the nearest Galactic globular cluster, and the observations are designed to constrain the number of binaries with massive companions (black holes, neutron stars, or white dwarfs) by measuring the ``wobble'' of the luminous (main-sequence) companion around the center of mass of the pair, with an astrometric precision of ~50 micro-arcseconds. The high spatial resolution and stable medium-band PSFs of WFC3/UVIS will make these measurements possible. In this work we describe: (i) the motivation behind this study, (ii) our observing strategy, (iii) the many other investigations enabled by this unique data set, and which of those our team is conducting, and (iv) a preliminary reduction of the first-epoch data-set collected on October 10, 2012. △ Less

Submitted 4 December, 2013; originally announced December 2013.

Comments: 25 pages, 14 figures (9 at low resolution), 3 tables. Published in: Astronomische Nachrichten, Volume 334, Issue 10, pages 1062-1085, December 2013. http://onlinelibrary.wiley.com/doi/10.1002/asna.201311911/abstract

On black hole sub-systems in idealized nuclear star clusters

Authors: Philip G. Breen, Douglas C. Heggie

Abstract: 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 s… ▽ More 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. △ Less

Submitted 21 August, 2013; originally announced August 2013.

Comments: 6 pages, 4 figures, Accepted for publication in MNRAS

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

Authors: Philip G. Breen, Douglas C. Heggie

Abstract: 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… ▽ More 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). △ Less

Submitted 11 April, 2013; originally announced April 2013.

Comments: 20 pages, 17 figures. Accepted for publication in MNRAS

Few-Body Modes of Binary Formation in Core Collapse

Authors: Ataru Tanikawa, Douglas C. Heggie, Piet Hut, Junichiro Makino

Abstract: 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 diss… ▽ More 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. △ Less

Submitted 20 November, 2013; v1 submitted 20 August, 2012; originally announced August 2012.

Comments: 44 pages, 18 figures, 1 table

Gravothermal oscillations in multi-component models of star clusters

Authors: Philip G. Breen, Douglas C. Heggie

Abstract: 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 gravot… ▽ More 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. △ Less

Submitted 11 July, 2012; originally announced July 2012.

Comments: 9 pages, 3 figures, 8 tables. Accepted for publication in MNRAS

The Tidal Tails of 47 Tucanae

Authors: Richard R. Lane, Andreas H. W. Küpper, Douglas C. Heggie

Abstract: 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 uncover… ▽ More 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. △ Less

Submitted 3 July, 2012; v1 submitted 11 April, 2012; originally announced April 2012.

Comments: Original paper has 10 pages, 10 figures and 2 tables. Please note that this now includes an erratum. Erratum has 6 pages, 8 figures and 2 tables. Ignore the exclamation marks in Section 2 of the erratum, these are an artifact of the LaTeX class file used to produce the manuscript

Journal ref: Monthly Notices of the Royal Astronomical Society (2012), 423, 2845

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

Authors: Mirek Giersz, Douglas C. Heggie, Jarrod Hurley, Arkadiusz Hypki

Abstract: 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 t… ▽ More 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. △ Less

Submitted 17 December, 2013; v1 submitted 29 December, 2011; originally announced December 2011.

Comments: 15 pages, 24 figures. Matched to published version in MNRAS

More on the structure of tidal tails

Authors: Andreas H. W. Kuepper, Richard R. Lane, Douglas C. Heggie

Abstract: 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… ▽ More 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. △ Less

Submitted 15 December, 2011; v1 submitted 21 November, 2011; originally announced November 2011.

Comments: 16 pages, 10 figures, accepted for publication in MNRAS

Gravothermal oscillations in two-component models of star clusters

Authors: Philip G. Breen, Douglas C. Heggie

Abstract: 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… ▽ More 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. △ Less

Submitted 24 October, 2011; originally announced October 2011.

Comments: 13 pages, 18 figures and 8 tables. Accepted for publication in MNRAS

Problems of Collisional Stellar Dynamics

Authors: Douglas C. Heggie

Abstract: 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 anal… ▽ More 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. △ Less

Submitted 24 March, 2011; v1 submitted 22 February, 2011; originally announced February 2011.

Comments: To appear in Bulletin of the Astronomical Society of India, Chandra Centennial Issue, eds D.J. Saikia and Virginia Trimble; 20 pages, 3 figures; minor editorial corrections

The life cycle of star cluster in a tidal field

Authors: Mark Gieles, Douglas C. Heggie, HongSheng Zhao

Abstract: 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-simi… ▽ More 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. △ Less

Submitted 16 June, 2011; v1 submitted 10 January, 2011; originally announced January 2011.

Comments: Updated to match final journal style

Journal ref: Mon. Not. R. Astron. Soc. 413, 2509-2524 (2011)

The effects of fly-bys on planetary systems

Authors: Daniel Malmberg, Melvyn B. Davies, Douglas C. Heggie

Abstract: 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 simi… ▽ More 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. △ Less

Submitted 21 September, 2010; originally announced September 2010.

Comments: 22 pages, 15 figures, 5 tables. Accepted for publication in MNRAS

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

Authors: Mirek Giersz, Douglas C. Heggie

Abstract: 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 cl… ▽ More 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. △ Less

Submitted 18 August, 2010; originally announced August 2010.

Comments: 16 pages, 16 figures, revised version submitted to MNRAS

Peculiarities in Velocity Dispersion and Surface Density Profiles of Star Clusters

Authors: A. H. W. Kuepper, P. Kroupa, H. Baumgardt, D. C. Heggie

Abstract: Based on our recent work on tidal tails of star clusters (Kuepper et al. 2009) we investigate star clusters of a few 10^4 Msun by means of velocity dispersion profiles and surface density profiles. We use a comprehensive set of $N$-body computations of star clusters on various orbits within a realistic tidal field to study the evolution of these profiles with time, and ongoing cluster dissolution… ▽ More Based on our recent work on tidal tails of star clusters (Kuepper et al. 2009) we investigate star clusters of a few 10^4 Msun by means of velocity dispersion profiles and surface density profiles. We use a comprehensive set of $N$-body computations of star clusters on various orbits within a realistic tidal field to study the evolution of these profiles with time, and ongoing cluster dissolution From the velocity dispersion profiles we find that the population of potential escapers, i.e. energetically unbound stars inside the Jacobi radius, dominates clusters at radii above about 50% of the Jacobi radius. Beyond 70% of the Jacobi radius nearly all stars are energetically unbound. The velocity dispersion therefore significantly deviates from the predictions of simple equilibrium models in this regime. We furthermore argue that for this reason this part of a cluster cannot be used to detect a dark matter halo or deviations from Newtonian gravity. By fitting templates to the about 10^4 computed surface density profiles we estimate the accuracy which can be achieved in reconstructing the Jacobi radius of a cluster in this way. We find that the template of King (1962) works well for extended clusters on nearly circular orbits, but shows significant flaws in the case of eccentric cluster orbits. This we fix by extending this template with 3 more free parameters. Our template can reconstruct the tidal radius over all fitted ranges with an accuracy of about 10%, and is especially useful in the case of cluster data with a wide radial coverage and for clusters showing significant extra-tidal stellar populations. No other template that we have tried can yield comparable results over this range of cluster conditions. All templates fail to reconstruct tidal parameters of concentrated clusters, however. (abridged) △ Less

Submitted 25 May, 2010; originally announced May 2010.

Comments: 23 pages, 13 figures, accepted for publication in MNRAS

Tidal Tails of Star Clusters

Authors: A. H. W. Kuepper, P. Kroupa, H. Baumgardt, D. C. Heggie

Abstract: Based on recent findings of a formation mechanism of substructure in tidal tails by Kuepper, Macleod & Heggie (2008) we investigate a more comprehensive set of N-body models of star clusters on orbits about a Milky-Way-like potential. We find that the predicted epicyclic overdensities arise in any tidal tail no matter which orbit the cluster follows as long as the cluster lives long enough for t… ▽ More Based on recent findings of a formation mechanism of substructure in tidal tails by Kuepper, Macleod & Heggie (2008) we investigate a more comprehensive set of N-body models of star clusters on orbits about a Milky-Way-like potential. We find that the predicted epicyclic overdensities arise in any tidal tail no matter which orbit the cluster follows as long as the cluster lives long enough for the overdensities to build up. The distance of the overdensities along the tidal tail from the cluster centre depends for circular orbits only on the mass of the cluster and the strength of the tidal field, and therefore decreases monotonically with time, while for eccentric orbits the orbital motion influences the distance, causing a periodic compression and stretching of the tails and making the distance oscillate with time. We provide an approximation for estimating the distance of the overdensities in this case. We describe an additional type of overdensity which arises in extended tidal tails of clusters on eccentric orbits, when the acceleration of the tidal field on the stellar stream is no longer homogeneous. Moreover, we conclude that a pericentre passage or a disk shock is not the direct origin of an overdensity within a tidal tail. Escape due to such tidal perturbations does not take place immediately after the perturbation but is rather delayed and spread over the orbit of the cluster. All observable overdensities are therefore of the mentioned two types. In particular, we note that substructured tidal tails do not imply the existence of dark-matter sub-structures in the haloes of galaxies. △ Less

Submitted 14 October, 2009; v1 submitted 14 September, 2009; originally announced September 2009.

Comments: 20 pages, 17 figures, MNRAS accepted

1Gyr in the Life of the Globular Cluster NGC 6397

Authors: Douglas C. Heggie, Mirek Giersz

Abstract: 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… ▽ More 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. △ Less

Submitted 30 April, 2009; originally announced April 2009.

Comments: 5 pages, 7 figures, to appear in MNRAS Letters

Monte Carlo Simulations of Star Clusters - VI. The globular cluster NGC 6397

Authors: M. Giersz, D. C. Heggie

Abstract: We describe Monte Carlo models for the dynamical evolution of the nearby globular cluster NGC 6397. 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 w… ▽ More We describe Monte Carlo models for the dynamical evolution of the nearby globular cluster NGC 6397. 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 the surface brightness profile, the velocity dispersion profile, and the luminosity function in two fields. We describe in particular those aspects of the evolution which distinguish this cluster from M4, which has a roughly similar mass and Galactocentric distance, but a qualitatively different surface brightness profile. Within the limitations of our modelling, we conclude that the most plausible explanation for the difference is fluctuations: both clusters are post-collapse objects, but sometimes have resolvable cores and sometimes not. △ Less

Submitted 13 February, 2009; v1 submitted 8 January, 2009; originally announced January 2009.

Comments: 11 pages, 12 figures, MNRAS, in press; revised in response to referee report

On the relationship between instability and Lyapunov times for the 3-body problem

Authors: D. J. Urminsky, D. C. Heggie

Abstract: In this study we consider the relationship between the survival time and the Lyapunov time for 3-body systems. It is shown that the Sitnikov problem exhibits a two-part power law relationship as demonstrated previously for the general 3-body problem. Using an approximate Poincare map on an appropriate surface of section, we delineate escape regions in a domain of initial conditions and use these… ▽ More In this study we consider the relationship between the survival time and the Lyapunov time for 3-body systems. It is shown that the Sitnikov problem exhibits a two-part power law relationship as demonstrated previously for the general 3-body problem. Using an approximate Poincare map on an appropriate surface of section, we delineate escape regions in a domain of initial conditions and use these regions to analytically obtain a new functional relationship between the Lyapunov time and the survival time for the 3-body problem. The marginal probability distributions of the Lyapunov and survival times are discussed and we show that the probability density function of Lyapunov times for the Sitnikov problem is similar to that for the general 3-body problem. △ Less

Submitted 30 October, 2008; originally announced October 2008.

Comments: 9 pages, 19 figures, accepted for publication in MNRAS

A search for spectroscopic binaries in the galactic globular cluster M4. Based on 5973 individual spectra collected at VLT

Authors: V. Sommariva, G. Piotto, M. Rejkuba, L. R. Bedin, D. C. Heggie, R. D. Mathieu, S. Villanova

Abstract: We present a large multi-epoch high resolution spectroscopic investigation for the search of binary candidates in the Galactic Globular Cluster (GGC) M4. The aim of our work is the identification of the binary candidates, and the determination of the binary fraction and of the binary radial distribution. We present a large multi-epoch high resolution spectroscopic investigation for the search of… ▽ More We present a large multi-epoch high resolution spectroscopic investigation for the search of binary candidates in the Galactic Globular Cluster (GGC) M4. The aim of our work is the identification of the binary candidates, and the determination of the binary fraction and of the binary radial distribution. We present a large multi-epoch high resolution spectroscopic investigation for the search of binary candidates in the Galactic Globular Cluster (GGC) M4. The aim of our work is the identification of the binary candidates, and the determination of the binary fraction and of the binary radial distribution. The average radial velocity of the observed cluster members is $70.29+/-0.07(+/-0.3)(+/-0.1)km/s. The search for variations in radial velocities among the stars with multi-epoch observations yielded 57 binary star candidates. Our radial velocity measurement accuracy allowed us to identify at a 3sigma level binaries with radial velocity variations larger than ~0.3km/s for the target stars with V<15, and larger than ~0.5km/s for the targets with V>15. We identified 4 binary star candidates out of 97 observed targets inside the core radius, and 53 candidates out of 2372 observed stars outside the core radius. Accounting for the incompleteness affecting our survey, the lower limit for the total binary fraction is f=3.0+/-0.3%. The lower limit for the binary fraction in the cluster core is f=5.1+/-2.3%, while outside the core it decreases to f=3.0+/-0.4%. Similarly, we found f=4.5+/-0.4% and f=1.8+/-0.6% for the binary fraction inside and outside the half mass radius. △ Less

Submitted 10 October, 2008; originally announced October 2008.

Comments: 13 pages, 12 figures, Accepted for publication in A&A

On the structure of tidal tails

Authors: Andreas H. W. Küpper, Andrew Macleod, Douglas C. Heggie

Abstract: We examine the longitudinal distribution of the stars escaping from a cluster along tidal tails. Using both theory and simulations, we show that, even in the case of a star cluster in a circular galactic orbit, when the tide is steady, the distribution exhibits maxima at a distance of many tidal radii from the cluster. We examine the longitudinal distribution of the stars escaping from a cluster along tidal tails. Using both theory and simulations, we show that, even in the case of a star cluster in a circular galactic orbit, when the tide is steady, the distribution exhibits maxima at a distance of many tidal radii from the cluster. △ Less

Submitted 15 April, 2008; originally announced April 2008.

Comments: 6 pages, 4 figures, accepted for publication in MNRAS

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.

Comments: 15 pages, 24 figures

Modelling individual globular clusters

Authors: D. C. Heggie, M. Giersz

Abstract: Astronomers have constructed models of globular clusters for over 100 years. These models mainly fall into two categories: (i) static models, such as King's model and its variants, and (ii) evolutionary models. Most attention has been given to static models, which are used to estimate mass-to-light ratios and mass segregation, and to combine data from proper motions and radial velocities. Evolut… ▽ More Astronomers have constructed models of globular clusters for over 100 years. These models mainly fall into two categories: (i) static models, such as King's model and its variants, and (ii) evolutionary models. Most attention has been given to static models, which are used to estimate mass-to-light ratios and mass segregation, and to combine data from proper motions and radial velocities. Evolutionary models have been developed for a few objects using the gaseous model, the Fokker-Planck model, Monte Carlo models and N-body models. These models have had a significant role in the search for massive black holes in globular clusters, for example. In this presentation the problems associated with these various techniques will be summarised, and then we shall describe new work with Giersz's Monte Carlo code, which has been enhanced recently to include the stellar evolution of single and binary stars. We describe in particular recent attempts to model the nearby globular cluster M4, including predictions on the spatial distribution of binary stars and their semi-major axis distribution, to illustrate the effects of about 12 Gyr of dynamical evolution. We also discuss work on an approximate way of predicting the "initial" conditions for such modelling. △ Less

Submitted 16 November, 2007; originally announced November 2007.

Comments: 10 pages,9 figures,invited review, to appear in the proceedings of IAU Symposium 246, 'Dynamical Evolution of Dense Stellar Systems', eds. E. Vesperini, M. Giersz, A. Sills, Capri, Sept. 2007. Sent to the editors 9 October 2007

Monte Carlo simulations of star clusters with primordial binaries. Comparison with N-body simulations and observations

Authors: Mirek Giersz, Douglas C. Heggie

Abstract: We outline the steps needed in to calibrate the Monte Carlo code in order to perform large scale simulations of real globular clusters. We calibrate the results against $N$-body simulations for $N = 2500$, 10000 and for the old open cluster M67. The calibration is done by choosing appropriate free code parameters. We outline the steps needed in to calibrate the Monte Carlo code in order to perform large scale simulations of real globular clusters. We calibrate the results against $N$-body simulations for $N = 2500$, 10000 and for the old open cluster M67. The calibration is done by choosing appropriate free code parameters. △ Less

Submitted 4 November, 2007; originally announced November 2007.

Comments: 5 pages, 8 figures, to appear in the proceedings of the 246th IAU symposium on 'Dynamical Evolution of Dense Stellar Systems', eds. E. Vesperini, M. Giersz, A. Sills, Capri, Sept. 2007

Modeling Dense Stellar Systems

Authors: Piet Hut, Shin Mineshige, Douglas C. Heggie, Junichiro Makino

Abstract: Black holes and neutron stars present extreme forms of matter that cannot be created as such in a laboratory on Earth. Instead, we have to observe and analyze the experiments that are ongoing in the Universe. The most telling observations of black holes and neutron stars come from dense stellar systems, where stars are crowded close enough to each other to undergo frequent interactions. It is th… ▽ More Black holes and neutron stars present extreme forms of matter that cannot be created as such in a laboratory on Earth. Instead, we have to observe and analyze the experiments that are ongoing in the Universe. The most telling observations of black holes and neutron stars come from dense stellar systems, where stars are crowded close enough to each other to undergo frequent interactions. It is the interplay between black holes, neutron stars and other objects in a dense environment that allows us to use observations to draw firm conclusions about the properties of these extreme forms of matter, through comparisons with simulations. The art of modeling dense stellar systems through computer simulations forms the main topic of this review. △ Less

Submitted 29 July, 2007; originally announced July 2007.

Comments: 23 pages, to appear in Prog. Theor. Phys

Journal ref: Prog.Theor.Phys.118:187-209,2007

Predictions for Triple Stars with and without a Pulsar in Star Clusters

Authors: M. Trenti, S. Ransom, P. Hut, D. C. Heggie

Abstract: Though about 80 pulsar binaries have been detected in globular clusters so far, no pulsar has been found in a triple system in which all three objects are of comparable mass. Here we present predictions for the abundance of such triple systems, and for the most likely characteristics of these systems. Our predictions are based on an extensive set of more than 500 direct simulations of star clust… ▽ More Though about 80 pulsar binaries have been detected in globular clusters so far, no pulsar has been found in a triple system in which all three objects are of comparable mass. Here we present predictions for the abundance of such triple systems, and for the most likely characteristics of these systems. Our predictions are based on an extensive set of more than 500 direct simulations of star clusters with primordial binaries, and a number of additional runs containing primordial triples. Our simulations employ a number N_{tot} of equal mass stars from N_{tot}=512 to N_{tot}=19661 and a primordial binary fraction from 0-50%. In addition, we validate our results against simulations with N=19661 that include a mass spectrum with a turn-off mass at 0.8 M_{sun}, appropriate to describe the old stellar populations of galactic globular clusters. Based on our simulations, we expect that typical triple abundances in the core of a dense cluster are two orders of magnitude lower than the binary abundances, which in itself already suggests that we don't have to wait too long for the first comparable-mass triple with a pulsar to be detected. △ Less

Submitted 3 April, 2008; v1 submitted 29 May, 2007; originally announced May 2007.

Comments: 11 pages, minor changes to match MNRAS accepted version

Exact equilibria of a stellar system in a linearised tidal field

Authors: David G. M. Mitchell, Douglas C. Heggie

Abstract: We study the motion of stars in a star cluster which revolves in a circular orbit about its parent galaxy. The star cluster is modelled as an ellipsoid of uniform spatial density. We exhibit two 2-parameter families of self-consistent equilibrium models in which the velocity at each point is confined to a line in velocity space. We exhibit the link between this problem and that of a uniform rota… ▽ More We study the motion of stars in a star cluster which revolves in a circular orbit about its parent galaxy. The star cluster is modelled as an ellipsoid of uniform spatial density. We exhibit two 2-parameter families of self-consistent equilibrium models in which the velocity at each point is confined to a line in velocity space. We exhibit the link between this problem and that of a uniform rotating ellipsoidal galaxy. With minimal adaptation, Freeman's bar models yield a third family. △ Less

Submitted 2 January, 2007; originally announced January 2007.

Comments: 4 pages, 5 figures, accepted by MNRAS

Journal ref: Mon.Not.Roy.Astron.Soc.376:705-708,2007

Dynamics of Planetary Systems in Star Clusters

Authors: R. Spurzem, M. Giersz, D. C. Heggie, D. N. C. Lin

Abstract: At least 10-15% of nearby sun-like stars have known Jupiter-mass planets. In contrast, very few planets are found in mature open and globular clusters such as the Hyades and 47 Tuc. We explore here the possibility that this dichotomy is due to the post-formation disruption of planetary systems associated with the stellar encounters in long-lived clusters. One supporting piece of evidence for thi… ▽ More At least 10-15% of nearby sun-like stars have known Jupiter-mass planets. In contrast, very few planets are found in mature open and globular clusters such as the Hyades and 47 Tuc. We explore here the possibility that this dichotomy is due to the post-formation disruption of planetary systems associated with the stellar encounters in long-lived clusters. One supporting piece of evidence for this scenario is the discovery of freely floating low-mass objects in star forming regions. We use two independent numerical approaches, a hybrid Monte Carlo and a direct $N$-body method, to simulate the impact of the encounters. We show that the results of numerical simulations are in reasonable agreement with analytical determinations in the adiabatic and impulsive limits. They indicate that distant stellar encounters generally do not significantly modify the compact and nearly circular orbits. However, moderately close stellar encounters, which are likely to occur in dense clusters, can excite planets' orbital eccentricity and induce dynamical instability in systems which are closely packed with multiple planets. We compute effective cross sections for the dissolution of planetary systems and show that, for all initial eccentricities, dissolution occurs on time scales which are longer than the dispersion of small stellar associations, but shorter than the age of typical open and globular clusters. Although it is much more difficult to disrupt short-period planets, close encounters can excite modest eccentricity among them, such that subsequent tidal dissipation leads to orbital decay, tidal inflation, and even disruption of the close-in planets. △ Less

Submitted 22 December, 2008; v1 submitted 28 December, 2006; originally announced December 2006.

Comments: 57 pages, 14 figures, 4 tables, major revision by authors, accepted for publication at the Astrophysical Journal (ApJ)

Report number: KITP Preprint NSF-KITP-09-14

Journal ref: Astrophys.J.697:458-482,2009

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

Authors: Douglas C. Heggie, Piet Hut, Shin Mineshige, Jun Makino, Holger Baumgardt

Abstract: (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,… ▽ More (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. △ Less

Submitted 12 April, 2007; v1 submitted 30 November, 2006; originally announced November 2006.

Comments: 4 pages, 2 figures, to appear in PASJ Letters, PASJ Vol. 59, No. 3. (to be issued on June 25, 2007). Minor revisions

McScatter: a Simple Three-Body Scattering Package with Stellar Evolution

Authors: Douglas C. Heggie, Simon Portegies Zwart, Jarrod Hurley

Abstract: We describe a simple computer package which illustrates a method of combining stellar dynamics with stellar evolution. Though the method is intended for elaborate applications (especially the dynamical evolution of rich star clusters) it is illustrated here in the context of three-body scattering, i.e. interactions between a binary star and a field of single stars. We describe the interface betw… ▽ More We describe a simple computer package which illustrates a method of combining stellar dynamics with stellar evolution. Though the method is intended for elaborate applications (especially the dynamical evolution of rich star clusters) it is illustrated here in the context of three-body scattering, i.e. interactions between a binary star and a field of single stars. We describe the interface between the dynamics and the two independent packages which describe the internal evolution of single stars and binaries. We also give an example application, and introduce a stand alone utility for the visual presentation of simulation results. △ Less

Submitted 13 April, 2006; originally announced April 2006.

Comments: 16 pages, Accepted for publication in New Astronomy. Source codes available at: http://manybody.org/manybody/McScatter.html and http://www.manybody.org/manybody/roche.html

Journal ref: New Astron. 12 (2006) 20-28

Star Clusters with Primordial Binaries: II. Dynamical Evolution of Models in a Tidal Field

Authors: M. Trenti, D. C. Heggie, P. Hut

Abstract: [abridged] We extend our analysis of the dynamical evolution of simple star cluster models, in order to provide comparison standards that will aid in interpreting the results of more complex realistic simulations. We augment our previous primordial-binary simulations by introducing a tidal field, and starting with King models of different central concentrations. We present the results of N-body… ▽ More [abridged] We extend our analysis of the dynamical evolution of simple star cluster models, in order to provide comparison standards that will aid in interpreting the results of more complex realistic simulations. We augment our previous primordial-binary simulations by introducing a tidal field, and starting with King models of different central concentrations. We present the results of N-body calculations of the evolution of equal-mass models, starting with primordial binary fractions of 0 - 100 %, and N values from 512 to 16384. We also attempt to extrapolate some of our results to the larger number of particles that are necessary to model globular clusters. We characterize the steady-state `deuterium main sequence' phase in which primordial binaries are depleted in the core in the process of `gravitationally burning'. In this phase we find that the ratio of the core to half-mass radius, r_c/r_h, is similar to that measured for isolated systems. In addition to the generation of energy due to hardening and depletion of the primordial binary population, the overall evolution of the star clusters is driven by a competing process: the tidal disruption of the system. We find that the depletion of primordial binaries before tidal dissolution of the system is possible only if the initial number is below 0.05 N, in the case of a King model with W_0=7 and N=4096 (which is one of our longest living models). We compare our findings, obtained by means of direct N-body simulations but scaled, where possible, to larger N, with similar studies carried out by means of Monte Carlo methods. △ Less

Submitted 6 October, 2006; v1 submitted 17 February, 2006; originally announced February 2006.

Comments: 15 pages, 18 figures, matches MNRAS accepted version, some sections reorganized but no major changes

Journal ref: Mon.Not.Roy.Astron.Soc.374:344-356,2007

Star Clusters with Primordial Binaries: I. Dynamical Evolution of Isolated Models

Authors: D. C. Heggie, M. Trenti, P. Hut

Abstract: In order to interpret the results of complex realistic star cluster simulations, which rely on many simplifying approximations and assumptions, it is essential to study the behavior of even more idealized models, which can highlight the essential physical effects and are amenable to more exact methods. With this aim, we present the results of N-body calculations of the evolution of equal-mass mo… ▽ More In order to interpret the results of complex realistic star cluster simulations, which rely on many simplifying approximations and assumptions, it is essential to study the behavior of even more idealized models, which can highlight the essential physical effects and are amenable to more exact methods. With this aim, we present the results of N-body calculations of the evolution of equal-mass models, starting with primordial binary fractions of 0 - 100 %, with values of N ranging from 256 to 16384. This allows us to extrapolate the main features of the evolution to systems comparable in particle number with globular clusters. In this range, we find that the steady-state `deuterium main sequence' is characterized by a ratio of the core radius to half-mass radius that follows qualitatively the analytical estimate by Vesperini & Chernoff (1994), although the N dependence is steeper than expected. Interestingly, for an initial binary fraction f greater than 10%, the binary heating in the core during the post collapse phase almost saturates (becoming nearly independent of f), and so little variation in the structural properties is observed. Thus, although we observe a significantly lower binary abundance in the core with respect to the Fokker-Planck simulations by Gao et al. (1991), this is of little dynamical consequence. At variance with the study of Gao et al. (1991), we see no sign of gravothermal oscillations before 150 halfmass relaxation times. At later times, however, oscillations become prominent. We demonstrate the gravothermal nature of these oscillations. △ Less

Submitted 17 February, 2006; originally announced February 2006.

Comments: 14 pages, 22 figures, MNRAS accepted

Journal ref: Mon.Not.Roy.Astron.Soc.368:677-689,2006

Gravitational Scattering

Authors: Douglas C. Heggie

Abstract: We review some modern applications of the theory of few-body encounters between binaries and single stars. Then we focus on a new treatment of adiabatic encounters, in a regime which is of importance in encounters between a star and a planetary system in a star cluster. We review some modern applications of the theory of few-body encounters between binaries and single stars. Then we focus on a new treatment of adiabatic encounters, in a regime which is of importance in encounters between a star and a planetary system in a star cluster. △ Less

Submitted 5 January, 2007; v1 submitted 20 December, 2005; originally announced December 2005.

Comments: 8 pages, 2 figures. To appear in Few-Body Problem, Annales Universitatis Turkuensis, C. Flynn, ed. Small correction to Fig.1

The Classical Gravitational N-Body Problem

Authors: Douglas C. Heggie

Abstract: Let a number, N, of particles interact classically through Newton's Laws of Motion and Newton's inverse square Law of Gravitation. The resulting equations of motion provide an approximate mathematical model with numerous applications in astrophysics, including the motion of the moon and other bodies in the Solar System (planets, asteroids, comets and meteor particles); stars in stellar systems r… ▽ More Let a number, N, of particles interact classically through Newton's Laws of Motion and Newton's inverse square Law of Gravitation. The resulting equations of motion provide an approximate mathematical model with numerous applications in astrophysics, including the motion of the moon and other bodies in the Solar System (planets, asteroids, comets and meteor particles); stars in stellar systems ranging from binary and other multiple stars to star clusters and galaxies; and the motion of dark matter particles in cosmology. For N=1 and N=2 the equations can be solved analytically. The case N=3 provides one of the richest of all unsolved dynamical problems -- the general three-body problem. For problems dominated by one massive body, as in many planetary problems, approximate methods based on perturbation expansions have been developed. In stellar dynamics, astrophysicists have developed numerous numerical and theoretical approaches to the problem for larger values of N, including treatments based on the Boltzmann equation and the Fokker-Planck equation; such N-body systems can also be modelled as self-gravitating gases, and thermodynamic insights underpin much of our qualitative understanding. △ Less

Submitted 11 August, 2005; v1 submitted 28 March, 2005; originally announced March 2005.

Comments: 29 pages, 4 figures, to appear in Encyclopaedia of Mathematical Physics, Elsevier, 2006; proof corrections and one new reference inserted