Search

We present the results of Monte Carlo simulations aimed at exploring the evolution towards energy equipartition of first- (1G) and second-generation (2G) stars in multiple-population globular clusters and how this evolution is affected by the initial differences between the spatial distributions of the two populations. Our results show that these initial differences have fundamental implications for the evolution towards energy equipartition of the two populations. We find that 2G stars, which are assumed to be initially more centrally concentrated than 1G stars, are generally characterized by a more rapid evolution towards energy equipartition. The evolution towards energy equipartition depends on the velocity dispersion component and is more rapid for the tangential velocity dispersion. The extent of the present-day differences between the degree of energy equipartition of 2G and 1G stars depends on the cluster's dynamical age and may be more significant in the tangential velocity dispersion and at intermediate distances from the cluster's center around the half-mass radius., Comment: 13 pages, 13 figures, Published in MNRAS

We report a new hot post-asymptotic giant branch (PAGB) star in the Galactic globular cluster (GC) E3, which is one of the first of the identified PAGB stars in a GC to show a binary signature. The star stands out as the brightest source in E3 in the \mbox{{\em Astrosat}}/UVIT images. We confirmed its membership with the cluster E3 using Gaia DR3 kinematics and parallax measurements. We supplemented the photometric observations with radial velocities (RVs) from high-resolution spectroscopic observations at two epochs and with ground- and space-based photometric observations from 0.13 $\mu$m to 22 $\mu$m. We find that the RVs vary over $\sim$6 \kms\ between the two epochs. This is an indication of the star being in a binary orbit. A simulation of possible binary systems with the observed RVs suggests a binary period of either 39.12 days or 17.83 days with mass ratio q$\geq$1.0. The [Fe/H] derived using the high-resolution spectra is $\sim -$0.7 dex, which closely matches the cluster metallicity. The spectroscopic and photometric measurements suggest \Teff\ and $\log g$ of the star as 17\,500$\pm$1\,000~K and 2.37$\pm$0.20~dex, respectively. Various PAGB evolutionary tracks on the Hertzsprung--Russell (H-R) diagram suggest a current mass of the star in the range 0.51$-$0.55 \Msun. The star is enriched with C and O abundances, showing similar CNO abundances compared to the other PAGB stars in GCs with the evidence of the third dredge-up on the AGB phase., Comment: Published in A&A Letters, Published Version

In this work, we investigate the black hole (BH) population of globular clusters (GCs) in Milky Way- (MW) and Andromeda- (M31) like galaxies. We combine the population synthesis code MASinGa and the MOCCA-Survey Database I to infer the properties of GCs harbouring a BH subsystem (BHS), an IMBH, or neither of those. We find that the typical number of GCs with a BHS, an IMBH, or none become comparable in the galactic outskirts, whilst the inner galactic regions are dominated by GCs without a significant dark component. Our models suggest that GCs harbouring a BHS are slightly heavier and with larger half-mass radii compared to the overall population. We retrieve the properties of binary BHs (BBHs) that have either merged in the last 3 Gyr or survived in their parent cluster until present-day. We find that around 80\% of the merging BBHs form due to dynamical interactions while the remaining originate from evolution of primordial binaries. We infer a merger rate for BBHs in the local Universe of $1.0-23\,\,\rm{yr^{-1}\,Gpc^{-3}}$, depending on the adopted assumptions. We find around 100-240 BBHs survive until present-day and are mostly concentrated in the inner few kpc of the galaxy. We estimate also the number of BHs transported into the galactic nucleus by infalling star clusters, finding around 1,000-3,000 BHs and 100-200 BBHs are transported over a time span of 12 Gyr. This enables us to constrains the total amount of BHs and BBHs binaries lurking in nuclear star cluster, i.e. $N_{BHs}=(1.4-2.2)\times10^4$ and $N_{BBHs}=700-1,100$., Comment: 18 pages,15 figures

A comprehensive study of the co-evolution of globular cluster systems (GCS) in galaxies requires the ability to model both the large scale dynamics (0.01 - 10 kpc) regulating their orbital evolution, and the small scale dynamics (sub-pc - AU) regulating the internal dynamics of each globular cluster (GC). In this work we present a novel method that combine semi-analytic models of GCS with fully self-consistent Monte Carlo models to simultaneously evolve large GCSs. We use the population synthesis code MASinGa and the MOCCA-Survey Database I to create synthetic GC populations aimed at representing the observed features of GCs in the Milky Way (MW) and Andromeda (M31). Our procedure enables us to recover the spatial and mass distribution of GCs in such galaxies, and to constrain the amount of mass that GCs left either in the halo as dispersed debris, or in the galactic centre, where they can contribute to the formation of a nuclear star cluster (NSC) and can bring stellar and possibly intermediate mass black holes there. The final masses reported by our simulations are of a few order of magnitudes smaller than the observed values. These differences show that mass build-up of a NSC and central BHs in galaxies like MW and M31 cannot be solely explained by the infalling GC scenario. This build-up is likely to depend on the interplay between interactions and mergers of infalling GCs and gas. The latter can contribute to both in-situ star formation in the NSC and growth of the central BH., Comment: 17 pages, 10 figures

We introduce a new prescription for the evolution of globular clusters (GCs) during the initial embedded gas phase into a Monte Carlo method. With a simplified version of the Monte Carlo MOCCA code embedded in the AMUSE framework, we study the survival of GCs after the removal of primordial gas. We first test our code and show that our results for the evolution of mass and Lagrangian radii are in good agreement with those obtained with N-body simulations. The Monte Carlo code enables a more rapid exploration of the evolution of systems with a larger number of stars than N-body simulations. We have carried out a new survey of simulations to explore the evolution of globular clusters with up to $N = 500000$ stars for a range of different star formation efficiencies and half-mass radii. Our study shows the range of initial conditions leading to the clusters' dissolution and those for which the clusters can survive this early evolutionary phase., Comment: 13 pages, 6 figures

We present results from a suite of eight direct N-body simulations, performed with \textsc{Nbody6++GPU}, representing realistic models of rotating star clusters with up to $1.1\times 10^5$ stars. Our models feature primordial (hard) binaries, a continuous mass spectrum, differential rotation, and tidal mass loss induced by the overall gravitational field of the host galaxy. We explore the impact of rotation and stellar evolution on the star cluster dynamics. In all runs for rotating star clusters we detect a previously predicted mechanism: an initial phase of violent relaxation followed by the so-called gravogyro catastrophe. We find that the gravogyro catastrophe reaches a finite amplitude, which depends in strength on the level of the bulk rotation, and then levels off. After this phase the angular momentum is transferred from high-mass to low-mass particles in the cluster (both stars and compact objects). Simultaneously, the system becomes gravothermally unstable and collapses, thus undergoing the so-called gravothermal-gravogyro catastrophe. Comparing models with and without stellar evolution, we find an interesting difference. When stellar evolution is not taken into account, the whole process proceeds at a faster pace. The population of heavy objects tend to form a triaxial structure that rotates in the cluster centre. When stellar evolution is taken into account, we find that such a {\it rotating bar} is populated by stellar black holes and their progenitors. The triaxial structure becomes axisymmetric over time, but we also find that the models without stellar evolution suffer repeated gravogyro catastrophes as sufficient angular momentum and mass are removed by the tidal field., Comment: 19 pages, 9 figures, 2 tables, Accepted by MNRAS

We present the implementation of updated stellar evolution recipes in the codes \texttt{Nbody6++GPU, MOCCA} and \texttt{McLuster}. We test them through numerical simulations of star clusters containing $1.1\times 10^5$ stars (with $2.0\times 10^4$ in primordial hard binaries) performing high-resolution direct $N$-body (\texttt{Nbody6++GPU}) and Monte-Carlo (\texttt{MOCCA}) simulations to an age of 10~Gyr. We compare models implementing either delayed or core-collapse supernovae mechanisms, a different mass ratio distribution for binaries, and white dwarf natal kicks enabled/disabled. Compared to \texttt{Nbody6++GPU}, the \texttt{MOCCA} models appear to be denser, with a larger scatter in the remnant masses, and a lower binary fraction on average. The \texttt{MOCCA} models produce more black holes (BHs) and helium white dwarfs (WDs), whilst \texttt{Nbody6++GPU} models are characterised by a much larger amount of WD-WD binaries. The remnant kick velocity and escape speed distributions are similar for the BHs and neutron stars (NSs), and some NSs formed via electron-capture supernovae, accretion-induced collapse or merger-induced collapse escape the cluster in all simulations. The escape speed distributions for the WDs, on the other hand, are very dissimilar. We categorise the stellar evolution recipes available in \texttt{Nbody6++GPU}, \texttt{MOCCA} and \texttt{Mcluster} into four levels: the one implemented in previous \texttt{Nbody6++GPU} and \texttt{MOCCA} versions (\texttt{level A}), state-of-the-art prescriptions (\texttt{level B}), some in a testing phase (\texttt{level C}), and those that will be added in future versions of our codes., Comment: 31 pages, 15 figures, 5 tables, Accepted by MNRAS

We have carried out a set of Monte Carlo simulations to study a number of fundamental aspects of the dynamical evolution of multiple stellar populations in globular clusters with different initial masses, fractions of second generation (2G) stars, and structural properties. Our simulations explore and elucidate: 1) the role of early and long-term dynamical processes and stellar escape in the evolution of the fraction of 2G stars and the link between the evolution of the fraction of 2G stars and various dynamical parameters; 2) the link between the fraction of 2G stars inside the cluster and in the population of escaping stars during a cluster's dynamical evolution; 3) the dynamics of the spatial mixing of the first-generation (1G) and 2G stars and the details of the structural properties of the two populations as they evolve toward mixing; 4) the implications of the initial differences between the spatial distribution of 1G and 2G stars for the evolution of the anisotropy in the velocity distribution and the expected radial profile of the 1G and 2G anisotropy for clusters at different stages of their dynamical history; 5) the variation of the degree of energy equipartition of the 1G and the 2G populations as a function of the distance from the cluster's centre and the cluster's evolutionary phase., Comment: 16 pages, 19 figures. Accepted for publication in MNRAS

Anisotropy of gravitational wave (GW) emission results in a net momentum gained by the black hole (BH) merger product, leading to a recoil velocity up to $\sim10^3\text{ km s}^{-1}$, which may kick it out of a globular cluster (GC). We estimate GW kick retention fractions of merger products assuming different models for BH spin magnitude and orientation (MS0 - random, MS1 - spin as a function of mass and metalicity, MS2 - constant value of $0.5$). We check how they depend on BH-BH merger time and properties of the cluster. We analyze the implications of GW kick retention fractions on intermediate massive BH (IMBH) formation by repeated mergers in a GC. We also calculate final spin of the merger product, and investigate how it correlates with effective spin of the binary. We used data from MOCCA (MOnte Carlo Cluster simulAtor) GC simulations to get a realistic sample of BH-BH mergers, assigned each BH spin value according to a studied model, and calculated recoil velocity and final spin based on most recent theoretical formulas. We discovered that for physically motivated models, GW kick retention fractions are about $30\%$ and display small dependence on assumptions about spin, but are much more prone to cluster properties. In particular, we discovered a strong dependence of GW kick retention fractions on cluster density. We also show that GW kick retention fractions are high in final life stages of the cluster, but low at the beginning. Finally, we derive formulas connecting final spin with effective spin for primordial binaries, and with maximal effective spin for dynamical binaries., Comment: 13 pages, 9 figures, accepted for publication in MNRAS

The first neutron star-neutron star (NS-NS) merger was discovered on August 17, 2017 through gravitational waves (GW170817) and followed with electromagnetic observations. This merger was detected in an old elliptical galaxy with no recent star formation. We perform a suite of numerical calculations to understand the formation mechanism of this merger. We probe three leading formation mechanisms of double compact objects: classical isolated binary star evolution, dynamical evolution in globular clusters and nuclear cluster formation to test whether they are likely to produce NS-NS mergers in old host galaxies. Our simulations with optimistic assumptions show current NS-NS merger rates at the level of 10^-2 yr^-1 from binary stars, 5 x 10^-5 yr^-1 from globular clusters and 10^-5 yr^-1 from nuclear clusters for all local elliptical galaxies (within 100 Mpc^3). These models are thus in tension with the detection of GW170817 with an observed rate 1.5 yr^-1 (per 100 Mpc^3; LIGO/Virgo estimate). Our results imply that either (i) the detection of GW170817 by LIGO/Virgo at their current sensitivity in an elliptical galaxy is a statistical coincidence; or that (ii) physics in at least one of our three models is incomplete in the context of the evolution of stars that can form NS-NS mergers; or that (iii) another very efficient (unknown) formation channel with a long delay time between star formation and merger is at play., Comment: A&A: accepted

All ten LIGO/Virgo binary black hole (BH-BH) coalescences reported from the O1/O2 runs have near zero effective spins. There are only three potential explanations of this fact. If the BH spin magnitudes are large then (i) either both BH spin vectors must be nearly in the orbital plane or (ii) the spin angular momenta of the BHs must be oppositely directed and similar in magnitude. Or, (iii) the BH spin magnitudes are small. We test the third hypothesis within the framework of the classical isolated binary evolution scenario of the BH-BH merger formation. We test three models of angular momentum transport in massive stars: a mildly efficient transport by meridional currents (as employed in the Geneva code), an efficient transport by the Tayler-Spruit magnetic dynamo (as implemented in the MESA code), and a very-efficient transport (as proposed by Fuller et al.) to calculate natal BH spins. We allow for binary evolution to increase the BH spins through accretion and account for the potential spin-up of stars through tidal interactions. Additionally, we update the calculations of the stellar-origin BH masses, include revisions to the history of star formation and to the chemical evolution across cosmic time. We find that we can match simultaneously the observed BH-BH merger rate density, BH masses, and effective spins. Models with efficient angular momentum transport are favored. The updated stellar-mass weighted gas-phase metallicity evolution now used in our models appears to be a key in better reproducing the LIGO/Virgo merger rate estimate. Mass losses during the pair-instability pulsation supernova phase are likely overestimated if the merger GW170729 hosts a BH more massive than 50 Msun. We also estimate rate of BH-NS mergers from recent LIGO/Virgo observations. Our updated models of BH-BH, BH-NS and NS-NS mergers are now publicly available at www.syntheticuniverse.org., Comment: A&A accepted: revised comparison with LIGO/Virgo (41 pages)

We discuss a new scenario for the formation of intermediate mass black holes in dense star clusters. In this scenario, intermediate mass black holes are formed as a result of dynamical interactions of hard binaries containing a stellar mass black hole, with other stars and binaries. We discuss the necessary conditions to initiate the process of intermediate mass black hole formation and the influence of an intermediate mass black hole on the host global globular cluster properties. We discuss two scenarios for intermediate mass black hole formation. The SLOW and FAST scenarios. They occur later or earlier in the cluster evolution and require smaller or extremely large central densities, respectively. In our simulations, the formation of intermediate mass black holes is highly stochastic. In general, higher formation probabilities follow from larger cluster concentrations (i.e. central densities). We further discuss possible observational signatures of the presence of intermediate mass black holes in globular clusters that follow from our simulations. These include the spatial and kinematic structure of the host cluster, possible radio, X-ray and gravitational wave emissions due to dynamical collisions or mass-transfer and the creation of hypervelocity main sequence escapers during strong dynamical interactions between binaries and an intermediate mass black hole. All simulations discussed in this paper were performed with the MOCCA Monte Carlo code. MOCCA accurately follows most of the important physical processes that occur during the dynamical evolution of star clusters but, as with other dynamical codes, it approximates the dissipative processes connected with stellar collisions and binary mergers., Comment: 16 pages, 1 table, 24 figures. printed in MNRAS

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

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

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

We use a self-consistent Monte Carlo treatment of stellar dynamics to investigate black hole binaries that are dynamically ejected from globular clusters to determine if they will be gravitational wave sources. We find that many of the ejected binaries have initially short periods and will merge within a Hubble time due to gravitational wave radiation. Thus they are potential sources for ground-based gravitational wave detectors. We estimate the yearly detection rate for current and advanced ground-based detectors and find a modest enhancement over the rate predicted for binaries produced by pure stellar evolution in galactic fields. We also find that many of the ejected binaries will pass through the longer wavelength Laser Interferometer Space Antenna (LISA) band and may be individually resolvable. We find a low probability that the Galaxy will contain a binary in the LISA band during its three-year mission. Some such binaries may, however, be detectable at Mpc distances implying that there may be resolvable stellar-mass LISA sources beyond our Galaxy. We conclude that globular clusters have a significant effect on the detection rate of ground-based detectors and may produce interesting LISA sources in local group galaxies., Comment: 19 pages, 16 figures, 2 tables, submitted to MNRAS

We study the compact binary population in star clusters, focusing on binaries containing black holes, using a self-consistent Monte Carlo treatment of dynamics and full stellar evolution. We find that the black holes experience strong mass segregation and become centrally concentrated. In the core the black holes interact strongly with each other and black hole-black hole binaries are formed very efficiently. The strong interactions, however, also destroy or eject the black hole-black hole binaries. We find no black hole-black hole mergers within our simulations but produce many hard escapers that will merge in the galactic field within a Hubble time. We also find several highly eccentric black hole-black hole binaries that are potential LISA sources, suggesting that star clusters are interesting targets for space-based detectors. We conclude that star clusters must be taken into account when predicting compact binary population statistics., Comment: 19 pages, 5 Tables, 12 Figures, updated in response to referee report, accepted for publication in MNRAS

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

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

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., Comment: 57 pages, 14 figures, 4 tables, major revision by authors, accepted for publication at the Astrophysical Journal (ApJ)

We present new models of the evolution and dissolution of star clusters evolving under the combined influence of internal relaxation and external tidal fields, using the anisotropic gaseous model based on the Fokker-Planck approximation, and a new escaper loss cone model. This model borrows ideas from loss cones of stellar distributions near massive black holes, and describes physical processes related to escaping stars by a simple model based on two timescales and a diffusion process. We compare our results with those of direct $N$-body models and of direct numerical solutions of the orbit-averaged Fokker-Planck equation. For this comparative study we limit ourselves to idealized single point mass star clusters, in order to present a detailed study of the physical processes determining the rate of mass loss, core collapse and other features of the system's evolution. With the positive results of our study the path is now open in the future to use the computationally efficient gaseous models for future studies with more realism (mass spectrum, stellar evolution)., Comment: 13 pages, 14 figures, added two figures, updated conclusion and discussion, accepted for publication in MNRAS

Spherically symmetric equal mass star clusters containing a large amount of primordial binaries are studied using a hybrid method, consisting of a gas dynamical model for single stars and a Monte Carlo treatment for relaxation of binaries and the setup of close resonant and fly-by encounters of single stars with binaries and binaries with each other (three- and four-body encounters). What differs from our previous work is that each encounter is being integrated using a highly accurate direct few-body integrator which uses regularized variables. Hence we can study the systematic evolution of individual binary orbital parameters (eccentricity, semi-major axis) and differential and total cross sections for hardening, dissolution or merging of binaries (minimum distance) from a sampling of several ten thousands of scattering events as they occur in real cluster evolution including mass segregation of binaries, gravothermal collapse and reexpansion, binary burning phase and ultimately gravothermal oscillations. For the first time we are able to present empirical cross sections for eccentricity variation of binaries in close three- and four-body encounters. It is found that a large fraction of three-body and four-body encounters results in merging. Previous cross sections obtained by Spitzer and Gao for strong encounters can be reproduced, while for weak encounters non-standard processes like formation of hierarchical triples occur., Comment: 16 pages, 19 figures, Latex in the MN style, submitted to MNRAS

This is a summary paper of MODEST-2, a workshop held at the Astronomical Institute ``Anton Pannekoek'' in Amsterdam, 16-17 December 2002. MODEST is a loose collaboration of people interested in MOdelling DEnse STellar systems, particularly those interested in modelling these systems using all the available physics (stellar dynamics, stellar evolution, hydrodynamics and the interplay between the three) by defining interfaces between different codes. In this paper, we summarize 1) the main advances in this endeavour since MODEST-1; 2) the main science goals which can be and should be addressed by these types of simulations; and 3) the most pressing theoretical and modelling advances that we identified., Comment: Accepted by New Astronomy

The new approach outlined in Paper I (Spurzem \& Giersz 1996) to follow the individual formation and evolution of binaries in an evolving, equal point-mass star cluster is extended for the self-consistent treatment of relaxation and close three- and four-body encounters for many binaries (typically a few percent of the initial number of stars in the cluster). The distribution of single stars is treated as a conducting gas sphere with a standard anisotropic gaseous model. A Monte Carlo technique is used to model the motion of binaries, their formation and subsequent hardening by close encounters, and their relaxation (dynamical friction) with single stars and other binaries. The results are a further approach towards a realistic model of globular clusters with primordial binaries without using special hardware. We present, as our main result, the self-consistent evolution of a cluster consisting of 300.000 equal point-mass stars, plus 30.000 equal mass binaries over several hundred half-mass relaxation times, well into the phase where most of the binaries have been dissolved and evacuated from the core. In a self-consistent model it is the first time that such a realistically large number of binaries is evolving in a cluster with an even ten times larger number of single stars. Due to the Monte Carlo treatment of the binaries we can at every moment analyze their external and internal parameters in the cluster as in an N-body simulation., Comment: LaTeX MNRAS Style 21 pages, 34 figures, submitted to MNRAS Nov. 1999, for preprint, see ftp://ftp.ari.uni-heidelberg.de/pub/spurzem/warspaper-98.ps.gz for associated mpeg-files (20 MB and 13 MB, respectively), see ftp://ftp.ari.uni-heidelberg.de/pub/spurzem/movie1.mpg and ftp://ftp.ari.uni-heidelberg.de/pub/spurzem/movie2.mpg

We review recent progress in computational methods for studying the dynamical evolution of large star clusters. After describing some highlights of a new comparison of different models, including N-body models, we mention in detail recent developments in three codes which allow models with anisotropic distributions of velocity to be studied: (i) a Fokker-Planck code using finite differences, (ii) a Monte Carlo code and (iii) a gaseous model code., Comment: 6 pages, postscript, etc., also available at http://www.maths.ed.ac.uk/people/douglas/preprints.html To appear in Andersen J., ed, Highlights of Astronomy Vol. 11. Kluwer, Dordrecht. Pagination differs from that in published version