Search

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., Comment: 44 pages, 18 figures, 1 table

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 ($\alpha$) 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 $\alpha$. 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., Comment: 9 pages, 3 figures, 8 tables. Accepted for publication in MNRAS

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., Comment: 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

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., Comment: 15 pages, 24 figures. Matched to published version in MNRAS

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/(\Omega^2-\partial^2\Phi/\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., Comment: 16 pages, 10 figures, accepted for publication in MNRAS

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., Comment: 13 pages, 18 figures and 8 tables. Accepted for publication in MNRAS

We demonstrate that free graphene sheet edges can curl back on themselves,reconstructing as nanotubes. This results in lower formation energies than any other non-functionalised edge structure reported to date in the literature. We determine the critical tube size and formation barrier and compare with density functional simulations of other edge terminations including a new reconstructed Klein edge. Simulated high resolution electron microscopy images show why such rolled edges may be difficult to detect. Rolled zigzag edges serve as metallic conduction channels, separated from the neighbouring bulk graphene by a chain of insulating sp$^3$-carbon atoms, and introduce Van Hove singularities into the graphene density of states., Comment: To appear in Phys. Rev. Lett

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., Comment: 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 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., Comment: Updated to match final journal style

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., Comment: 22 pages, 15 figures, 5 tables. Accepted for publication in MNRAS

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., Comment: 16 pages, 16 figures, revised version submitted to MNRAS

Most globular clusters have half-mass radii of a few pc with no apparent correlation with their masses. This is different from elliptical galaxies, for which the Faber-Jackson relation suggests a strong positive correlation between mass and radius. Objects that are somewhat in between globular clusters and low-mass galaxies, such as ultra-compact dwarf galaxies, have a mass-radius relation consistent with the extension of the relation for bright ellipticals. Here we show that at an age of 10 Gyr a break in the mass-radius relation at ~10^6 Msun is established because objects below this mass, i.e. globular clusters, have undergone expansion driven by stellar evolution and hard binaries. From numerical simulations we find that the combined energy production of these two effects in the core comes into balance with the flux of energy that is conducted across the half-mass radius by relaxation. An important property of this `balanced' evolution is that the cluster half-mass radius is independent of its initial value and is a function of the number of bound stars and the age only. It is therefore not possible to infer the initial mass-radius relation of globular clusters and we can only conclude that the present day properties are consistent with the hypothesis that all hot stellar systems formed with the same mass-radius relation and that globular clusters have moved away from this relation because of a Hubble time of stellar and dynamical evolution., Comment: 5 pages, 3 figures, MNRAS Letters (accepted)

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), Comment: 23 pages, 13 figures, accepted for publication in MNRAS

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., Comment: 20 pages, 17 figures, MNRAS accepted

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., Comment: 5 pages, 7 figures, to appear in MNRAS Letters

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., Comment: 11 pages, 12 figures, MNRAS, in press; revised in response to referee report

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., Comment: 9 pages, 19 figures, accepted for publication in MNRAS

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., Comment: 13 pages, 12 figures, Accepted for publication in A&A

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., Comment: 6 pages, 4 figures, accepted for publication in MNRAS

We describe Monte Carlo models for the dynamical evolution of the nearby globular cluster M4. The code includes treatments of two-body relaxation, 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 satisfactory match to the surface brightness profile, the velocity dispersion profile, and the luminosity function in two fields. We describe in particular the evolution of the core, and find that M4 (which has a classic King profile) is actually a post-collapse cluster, its core radius being sustained by binary burning. We also consider the distribution of its binaries, including those which would be observed as photometric binaries and as radial-velocity binaries. We also consider the populations of white dwarfs, neutron stars, black holes and blue stragglers, though not all channels for blue straggler formation are represented yet in our simulations., Comment: 10 pages, 13 figures. Some numbers and figures may change slightly before journal submission

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., Comment: 15 pages, 24 figures

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., Comment: 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

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., Comment: 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

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., Comment: 23 pages, to appear in Prog. Theor. Phys

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., Comment: 11 pages, minor changes to match MNRAS accepted version