Understanding the formation of stellar clusters requires following the interplay between gas and newly formed stars accurately. We therefore couple the magnetohydrodynamics code FLASH to the N-body code ph4 and the stellar evolution code SeBa using the Astrophysical Multipurpose Software Environment (AMUSE) to model stellar dynamics, evolution, and collisional N-body dynamics and the formation of binary and higher-order multiple systems, while implementing stellar feedback in the form of radiation, stellar winds and supernovae in FLASH. We use this novel numerical method to simulate the formation and early evolution of open clusters of $\sim 1000$ stars formed from clouds with a mass range of $10^3$-$10^5$ M$_\odot$. Analyzing the effects of stellar feedback on the gas and stars of the natal cluster, we find that our clusters are resilient to disruption, even in the presence of intense feedback. This can even slightly increase the amount of dense, Jeans unstable gas by sweeping up shells; thus, a stellar wind strong enough to trap its own H II region shows modest triggering of star formation. Our clusters are born moderately mass segregated, an effect enhanced by feedback, and retained after the ejection of their natal gas, in agreement with observations.

To fully understand the diverse population of exoplanets, we must study their early lives within open clusters, the birthplace of most stars with masses $>0.5M_\odot$ (including those currently in the field). Indeed, when we observe planets within clustered environments, we notice highly eccentric and odd systems that suggest the importance of dynamical pathways created by interactions with additional bodies (as in the case of HD 285507b). However, it has proven difficult to investigate these effects, as many current numerical solvers for the multi-scale $N$-body problem are simplified and limited in scope. To remedy this, we aim to create a physically complete computational solution to explore the role of stellar close encounters and interplanetary interactions in producing the observed exoplanet populations for both open cluster stars and field stars. We present a new code, TYCHO, which employs a variety of different computational techniques, including multiple $N$-body integration methods, close encounter handling, modified Monte Carlo scattering experiments, and a variety of empirically informed initial conditions. We discuss the methodology in detail, and its implementation within the AMUSE software framework. Approximately 1$\%$ of our systems are promptly disrupted by star-star encounters contributing to the rogue planets occurrence rate. Additionally, we find that close encounters which perturb long-period planets lead to 38.3$\%$ of solar-system-like planetary systems becoming long-term unstable.

Hierarchical code coupling strategies make it possible to combine the results of individual numerical solvers into a self-consistent symplectic solution. We explore the possibility of allowing such a coupling strategy to be non-intrusive. In that case, the underlying numerical implementation is not affected by the coupling itself, but its functionality is carried over in the interface. This method is efficient for solving the equations of motion for a self-gravitating system over a wide range of scales. We adopt a dedicated integrator for solving each particular part of the problem and combine the results to a self-consistent solution. In particular, we explore the possibilities of combining the evolution of one or more microscopic systems that are embedded in a macroscopic system. The here presented generalizations of Bridge include higher-order coupling strategies (from the classic 2nd order up to 10th-order), but we also demonstrate how multiple bridges can be nested and how additional processes can be introduced at the bridge time-step to enrich the physics, for example by incorporating dissipative processes. Such augmentation allows for including additional processes in a classic Newtonian N-body integrator without alterations to the underlying code. These additional processes include for example the Yarkovsky effect, dynamical friction or relativistic dynamics. Some of these processes operate on all particles whereas others apply only to a subset. The presented method is non-intrusive in the sense that the underlying methods remain operational without changes to the code (apart from adding the get- and set-functions to enable the bridge operator). As a result, the fundamental integrators continue to operate with their internal time step and preserve their local optimizations and parallelism. ... abridged ...

We perform a simulation using the Astrophysical Multipurpose Software Environment of the Orion Trapezium star cluster in which the evolution of the stars and the dynamics of planetary systems are taken into account. The initial conditions from earlier simulations were selected in which the size and mass distributions of the observed circumstellar disks in this cluster are satisfactorily reproduced. Four, five, or size planets per star were introduced in orbit around the 500 solar-like stars with a maximum orbital separation of 400au. Our study focuses on the production of free-floating planets. A total of 357 become unbound from a total of 2522 planets in the initial conditions of the simulation. Of these, 281 leave the cluster within the crossing timescale of the star cluster; the others remain bound to the cluster as free-floating intra-cluster planets. Five of these free-floating intra-cluster planets are captured at a later time by another star. The two main mechanisms by which planets are lost from their host star, ejection upon a strong encounter with another star or internal planetary scattering, drive the evaporation independent of planet mass of orbital separation at birth. The effect of small perturbations due to slow changes in the cluster potential are important for the evolution of planetary systems. In addition, the probability of a star to lose a planet is independent of the planet mass and independent of its initial orbital separation. As a consequence, the mass distribution of free-floating planets indistinguishable from the mass distribution of planets bound to their host star.

We describe a star cluster formation model that includes individual star formation from self-gravitating, magnetized gas, coupled to collisional stellar dynamics. The model uses the Astrophysical Multi-purpose Software Environment (AMUSE) to integrate an adaptive-mesh magnetohydrodynamics code (FLASH) with a fourth order Hermite N-body code (ph4), a stellar evolution code (SeBa), and a method for resolving binary evolution (multiples). This combination yields unique star formation simulations that allow us to study binaries formed dynamically from interactions with both other stars and dense, magnetized gas subject to stellar feedback during the birth and early evolution of stellar clusters. We find that for massive stars, our simulations are consistent with the observed dynamical binary fractions and mass ratios. However, our binary fraction drops well below observed values for lower mass stars, presumably due to unincluded binary formation during initial star formation. Further, we observe a build up of binaries near the hard-soft boundary that may be an important mechanism driving early cluster contraction.

The GraphBLAS standard (GraphBlas.org) is being developed to bring the potential of matrix based graph algorithms to the broadest possible audience. Mathematically the Graph- BLAS defines a core set of matrix-based graph operations that can be used to implement a wide class of graph algorithms in a wide range of programming environments. This paper provides an introduction to the mathematics of the GraphBLAS. Graphs represent connections between vertices with edges. Matrices can represent a wide range of graphs using adjacency matrices or incidence matrices. Adjacency matrices are often easier to analyze while incidence matrices are often better for representing data. Fortunately, the two are easily connected by matrix mul- tiplication. A key feature of matrix mathematics is that a very small number of matrix operations can be used to manipulate a very wide range of graphs. This composability of small number of operations is the foundation of the GraphBLAS. A standard such as the GraphBLAS can only be effective if it has low performance overhead. Performance measurements of prototype GraphBLAS implementations indicate that the overhead is low.

We present the results of a survey of N-body simulations aimed at exploring the evolution of compact binaries in multiple-population globular clusters.We show that as a consequence of the initial differences in the structural properties of the first-generation (FG) and the second-generation (SG) populations and the effects of dynamical processes on binary stars, the SG binary fraction decreases more rapidly than that of the FG population. The difference between the FG and SG binary fraction is qualitatively similar to but quantitatively smaller than that found for wider binaries in our previous investigations.The evolution of the radial variation of the binary fraction is driven by the interplay between binary segregation, ionization and ejection. Ionization and ejection counteract in part the effects of mass segregation but for compact binaries the effects of segregation dominate and the inner binary fraction increases during the cluster evolution. We explore the variation of the difference between the FG and the SG binary fraction with the distance from the cluster centre and its dependence on the binary binding energy and cluster structural parameters. The difference between the binary fraction in the FG and the SG populations found in our simulations is consistent with the results of observational studies finding a smaller binary fraction in the SG population.

The discovery of multiple stellar populations in globular clusters has implications for all the aspects of the study of these stellar systems. In this paper, by means of N-body simulations, we study the evolution of binary stars in multiple-population clusters and explore the implications of the initial differences in the spatial distribution of different stellar populations for the evolution and survival of their binary stars. Our simulations show that initial differences between the spatial distribution of first-generation (FG) and second-generation (SG) stars can leave a fingerprint in the current properties of the binary population. SG binaries are disrupted more efficiently than those of the FG population resulting in a global SG binary fraction smaller than that of the FG. As for surviving binaries, dynamical evolution produces a difference between the SG and the FG binary binding energy distribution with the SG population characterized by a larger fraction of high binding energy (more bound) binaries. We have also studied the dependence of the binary properties on the distance from the cluster centre. Although the global binary fraction decreases more rapidly for the SG population, the local binary fraction measured in the cluster inner regions may still be dominated by SG binaries. The extent of the differences between the surviving FG and SG binary binding energy distribution also varies radially within the cluster and is larger in the cluster inner regions.

Globular clusters are among the most congested stellar systems in the Universe. Internal dynamical evolution drives them toward states of high central density, while simultaneously concentrating the most massive stars and binary systems in their cores. As a result, these clusters are expected to be sites of frequent close encounters and physical collisions between stars and binaries, making them efficient factories for the production of interesting and observable astrophysical exotica. I describe some elements of the competition among stellar dynamics, stellar evolution, and other processes that control globular cluster dynamics, with particular emphasis on pathways that may lead to the formation of blue stragglers.

We study the effects of the external tidal field on the violent relaxation phase of star clusters dynamical evolution, with particular attention to the kinematical properties of the equilibrium configurations emerging at the end of this phase.We show that star clusters undergoing the process of violent relaxation in the tidal field of their host galaxy can acquire significant internal differential rotation and are characterized by a distinctive radial variation of the velocity anisotropy. These kinematical properties are the result of the symmetry breaking introduced by the external tidal field in the collapse phase and of the action of the Coriolis force on the orbit of the stars. The resulting equilibrium configurations are characterized by differential rotation, with a peak located between one and two half-mass radii. As for the anisotropy, similar to clusters evolving in isolation, the systems explored in this Letter are characterized by an inner isotropic core, followed by a region of increasing radial anisotropy. However for systems evolving in an external tidal field the degree of radial anisotropy reaches a maximum in the cluster intermediate regions and then progressively decreases, with the cluster outermost regions being characterized by isotropy or a mild tangential anisotropy. Young or old but less-relaxed dynamically young star clusters may keep memory of these kinematical fingerprints of their early dynamical evolution.

We perform a series of simulations of evolving star clusters using AMUSE (the Astrophysical Multipurpose Software Environment), a new community-based multi-physics simulation package, and compare our results to existing work. These simulations model a star cluster beginning with a King model distribution and a selection of power-law initial mass functions, and contain a tidal cut-off. They are evolved using collisional stellar dynamics and include mass loss due to stellar evolution. After determining that the differences between AMUSE results and prior publications are understood, we explored the variation in cluster lifetimes due to the random realization noise introduced by transforming a King model to specific initial conditions. This random realization noise can affect the lifetime of a simulated star cluster by up to 30%. Two modes of star cluster dissolution were identified: a mass evolution curve that contains a run-away cluster dissolution with a sudden loss of mass, and a dissolution mode that does not contain this feature. We refer to these dissolution modes as "dynamical" and "relaxation" dominated respectively. For Salpeter-like initial mass functions, we determined the boundary between these two modes in terms of the dynamical and relaxation time scales.

We present a method to couple N-body star cluster simulations to a cosmological tidal field, using the Astrophysical Multipurpose Software Environment. We apply this method to star clusters embedded in the CosmoGrid dark matter-only LambdaCDM simulation. Our star clusters are born at z = 10 (corresponding to an age of the Universe of about 500 Myr) by selecting a dark matter particle and initializing a star cluster with 32,000 stars on its location. We then follow the dynamical evolution of the star cluster within the cosmological environment. We compare the evolution of star clusters in two Milky-Way size haloes with a different accretion history. The mass loss of the star clusters is continuous irrespective of the tidal history of the host halo, but major merger events tend to increase the rate of mass loss. From the selected two dark matter haloes, the halo that experienced the larger number of mergers tends to drive a smaller mass loss rate from the embedded star clusters, even though the final masses of both haloes are similar. We identify two families of star clusters: native clusters, which become part of the main halo before its final major merger event, and the immigrant clusters, which are accreted upon or after this event; native clusters tend to evaporate more quickly than immigrant clusters. Accounting for the evolution of the dark matter halo causes immigrant star clusters to retain more mass than when the z=0 tidal field is taken as a static potential. The reason for this is the weaker tidal field experienced by immigrant star clusters before merging with the larger dark matter halo.

We present the open source Astrophysical Multi-purpose Software Environment (AMUSE, www.amusecode.org), a component library for performing astrophysical simulations involving different physical domains and scales. It couples existing codes within a Python framework based on a communication layer using MPI. The interfaces are standardized for each domain and their implementation based on MPI guarantees that the whole framework is well-suited for distributed computation. It includes facilities for unit handling and data storage. Currently it includes codes for gravitational dynamics, stellar evolution, hydrodynamics and radiative transfer. Within each domain the interfaces to the codes are as similar as possible. We describe the design and implementation of AMUSE, as well as the main components and community codes currently supported and we discuss the code interactions facilitated by the framework. Additionally, we demonstrate how AMUSE can be used to resolve complex astrophysical problems by presenting example applications.

We present the results of the "Cosmogrid" cosmological N-body simulation suites based on the concordance LCDM model. The Cosmogrid simulation was performed in a 30Mpc box with 2048^3 particles. The mass of each particle is 1.28x10^5 Msun, which is sufficient to resolve ultra-faint dwarfs. We found that the halo mass function shows good agreement with the Sheth & Tormen fitting function down to ~10^7 Msun. We have analyzed the spherically averaged density profiles of the three most massive halos which are of galaxy group size and contain at least 170 million particles. The slopes of these density profiles become shallower than -1 at the inner most radius. We also find a clear correlation of halo concentration with mass. The mass dependence of the concentration parameter cannot be expressed by a single power law, however a simple model based on the Press-Schechter theory proposed by Navarro et al. gives reasonable agreement with this dependence. The spin parameter does not show a correlation with the halo mass. The probability distribution functions for both concentration and spin are well fitted by the log-normal distribution for halos with the masses larger than ~10^8 Msun. The subhalo abundance depends on the halo mass. Galaxy-sized halos have 50% more subhalos than ~10^{11} Msun halos have.

We describe the work of the Meta-Institute for Computational Astrophysics (MICA), the first professional scientific organization based in virtual worlds. MICA was an experiment in the use of this technology for science and scholarship, lasting from the early 2008 to June 2012, mainly using the Second Life and OpenSimulator as platforms. We describe its goals and activities, and our future plans. We conducted scientific collaboration meetings, professional seminars, a workshop, classroom instruction, public lectures, informal discussions and gatherings, and experiments in immersive, interactive visualization of high-dimensional scientific data. Perhaps the most successful of these was our program of popular science lectures, illustrating yet again the great potential of immersive VR as an educational and outreach platform. While the members of our research groups and some collaborators found the use of immersive VR as a professional telepresence tool to be very effective, we did not convince a broader astrophysics community to adopt it at this time, despite some efforts; we discuss some possible reasons for this non-uptake. On the whole, we conclude that immersive VR has a great potential as a scientific and educational platform, as the technology matures and becomes more broadly available and accepted.

In this paper we study the long-term dynamical evolution of multiple-population clusters, focusing on the evolution of the spatial distributions of the first- (FG) and second-generation (SG) stars.In previous studies we have suggested that SG stars formed from the ejecta of FG AGB stars are expected initially to be concentrated in the cluster inner regions. Here, by means of N-body simulations, we explore the time scales and the dynamics of the spatial mixing of the FG and the SG populations and their dependence on the SG initial concentration.Our simulations show that, as the evolution proceeds, the radial profile of the SG/FG number ratio, NSG/NFG, is characterized by three regions: 1) a flat inner part; 2) a declining part in which FG stars are increasingly dominant; and 3) an outer region where the NSG/NFG profile flattens again (the NSG/NFG profile may rise slightly again in the outermost cluster regions). The radial variation of NSG/NFG implies that the fraction of SG stars determined by observations covering a limited range of radial distances is not, in general, equal to the SG global fraction, (NSG/NFG)glob. The distance at which NSG/NFG equals (NSG/NFG)glob is approximately between 1 and 2 cluster half-mass radii. The results of our simulations suggest that in many Galactic globular clusters the SG should still be more spatially concentrated than the FG.[abridged]

Several dynamical scenarios have been proposed that can lead to prompt mass segregation on the crossing time scale of a young cluster. They generally rely on cool and/or clumpy initial conditions, and are most relevant to small systems. As a counterpoint, we present a novel dynamical mechanism that can operate in relatively large, homogeneous, cool or cold systems. This mechanism may be important in understanding the assembly of large mass-segregated clusters from smaller clumps.

We introduce a general-purpose framework for interconnecting scientific simulation programs using a homogeneous, unified interface. Our framework is intrinsically parallel, and conveniently separates all component numerical modules in memory. This strict separation allows automatic unit conversion, distributed execution of modules on different cores within a cluster or grid, and orderly recovery from errors. The framework can be efficiently implemented and incurs an acceptable overhead. In practice, we measure the time spent in the framework to be less than 1% of the wall-clock time. Due to the unified structure of the interface, incorporating multiple modules addressing the same physics in different ways is relatively straightforward. Different modules may be advanced serially or in parallel. Despite initial concerns, we have encountered relatively few problems with this strict separation between modules, and the results of our simulations are consistent with earlier results using more traditional monolithic approaches. This framework provides a platform to combine existing simulation codes or develop new physical solver codes within a rich "ecosystem" of interchangeable modules.

We present a data format for the output of general N-body simulations, allowing the presence of individual time steps. By specifying a standard, different N-body integrators and different visualization and analysis programs can all share the simulation data, independent of the type of programs used to produce the data. Our Particle Stream Data Format, PSDF, is specified in YAML, based on the same approach as XML but with a simpler syntax. Together with a specification of PSDF, we provide background and motivation, as well as specific examples in a variety of computer languages. We also offer a web site from which these examples can be retrieved, in order to make it easy to augment existing codes in order to give them the option to produce PSDF output.

We describe AMUSE, the Astrophysical Multipurpose Software Environment, a programming framework designed to manage multi-scale, multi-physics simulations in a hierarchical, extensible, and internally consistent way. Constructed as a collection of individual modules, AMUSE allows computational tools for different physical domains to be easily combined into a single task. It facilitates the coupling of modules written in different languages by providing inter-language tools and a standard programming interface that represents a balance between generality and computational efficiency. The framework currently incorporates the domains of stellar dynamics, stellar evolution, gas dynamics, and radiative transfer. We present some applications of the framework and outline plans for future development of the package.

We introduce a general-purpose framework for interconnecting scientific simulation programs using a homogeneous, unified software interface. Our framework is intrinsically parallel, and conveniently separates all components in memory. It performs unit conversion between different modules automatically and defines common data structures to communicate across different codes. We use the framework to simulate embedded star clusters. For this purpose we couple solvers for gravitational dynamics, stellar evolution and hydrodynamics to self consistently resolve the dynamical evolution simultaneousy with the internal nuclear evolution of the stars and the hydrodynamic response of the gas. We find, in contrast to earlier studies, that the survival of a young star cluster depends only weakly on the efficiency of star formation. The main reason for this weak dependency is the asymmetric expulsion of the embedding gas from the cluster.

The discovery of multiple stellar populations in globular clusters raises fundamental questions concerning the formation and dynamical history of these systems. In a previous study aimed at exploring the formation of second-generation (SG) stars from the ejecta of first-generation (FG) AGB stars, and the subsequent dynamical evolution of the cluster, we showed that SG stars are expected to form in a dense subsystem concentrated in the inner regions of the FG cluster. In this paper we explore the implications of the structural properties of multiple-population clusters, and in particular the presence of the inner SG subsystem, for the disruption of binary stars. We quantify the enhancement of the binary disruption rate due to the presence of the central SG subsystem for a number of different initial conditions. Our calculations show that SG binaries, which are assumed to be more concentrated in the cluster inner regions, are disrupted at a substantially larger rate than FG binaries. Assuming a similar initial fraction of FG and SG binaries, our dynamical study indicates that the SG population is now expected to contain a significantly smaller binary fraction than the FG population.

Many observational studies have revealed the presence of multiple stellar generations in Galactic globular clusters. These studies suggest that second-generation stars make up a significant fraction of the current mass of globular clusters, with the second-generation mass fraction ranging from ~50 to 80 per cent in individual clusters. In this Letter we carry out hydrodynamical simulations to explore the dependence of the mass of second-generation stars on the initial mass and structural parameters and stellar initial mass function (IMF) of the parent cluster. We then use the results of these simulations to estimate the fraction, f_{SG,H}, of the mass of the Galactic stellar halo composed of second-generation stars that originated in globular clusters. We study the dependence of f_{SG,H} on the parameters of the initial mass function of the Galactic globular cluster system. For a broad range of initial conditions, we find that the fraction of mass of the Galactic stellar halo in second-generation stars is always small, f_{SG,H}<4-6 per cent for a Kroupa-1993 IMF and f_{SG,H}<7-9 per cent for a Kroupa-2001 IMF.

A large number of spectroscopic studies have provided evidence of the presence of multiple populations in globular clusters by revealing patterns in the stellar chemical abundances. This paper is aimed at studying the origin of these abundance patterns. We explore a model in which second generation (SG) stars form out of a mix of pristine gas and ejecta of the first generation of asymptotic giant branch stars. We first study the constraints imposed by the spectroscopic data of SG stars in globular clusters on the chemical properties of the asymptotic and super asymptotic giant branch ejecta. With a simple one-zone chemical model, we then explore the formation of the SG population abundance patterns focussing our attention on the Na-O, Al-Mg anticorrelations and on the helium distribution function. We carry out a survey of models and explore the dependence of the final SG chemical properties on the key parameters affecting the gas dynamics and the SG formation process. Finally, we use our chemical evolution framework to build specific models for NGC 2808 and M4, two Galactic globular clusters which show different patterns in the Na-O and Mg-Al anticorrelation and have different helium distributions. We find that the amount of pristine gas involved in the formation of SG stars is a key parameter to fit the observed O-Na and Mg-Al patterns. The helium distribution function for these models is in general good agreement with the observed one. Our models, by shedding light on the role of different parameters and their interplay in determining the final SG chemical properties, illustrate the basic ingredients, constraints and problems encountered in this self-enrichment scenario which must be addressed by more sophisticated chemical and hydrodynamic simulations.

Spectroscopic and photometric observations show that many globular clusters host multiple stellar populations, challenging the common paradigm that globular clusters are "simple stellar populations" composed of stars of uniform age and chemical composition. The chemical abundances of second-generation (SG) stars constrain the sources of gas out of which these stars must have formed, indicating that the gas must contain matter processed through the high-temperature CNO cycle. First-generation massive Asymptotic Giant Branch (AGB) stars have been proposed as the source of this gas. In a previous study, by means of hydrodynamical and N-body simulations, we have shown that the AGB ejecta collect in a cooling flow in the cluster core, where the gas reaches high densities, ultimately forming a centrally concentrated subsystem of SG stars. In this Letter we show that the high gas density can also lead to significant accretion onto a pre-existing seed black hole. We show that gas accretion can increase the black hole mass by up to a factor of 100. The details of the gas dynamics are important in determining the actual black hole growth. Assuming a near-universal seed black hole mass and small cluster-to-cluster variations in the duration of the SG formation phase, the outcome of our scenario is one in which the present intermediate-mass black hole (IMBH) mass may have only a weak dependence on the current cluster properties. The scenario presented provides a natural mechanism for the formation of an IMBH at the cluster center during the SG star-formation phase.

Young massive clusters are dense aggregates of young stars that form the fundamental building blocks of galaxies. Several examples exist in the Milky Way Galaxy and the Local Group, but they are particularly abundant in starburst and interacting galaxies. The few young massive clusters that are close enough to resolve are of prime interest for studying the stellar mass function and the ecological interplay between stellar evolution and stellar dynamics. The distant unresolved clusters may be effectively used to study the star-cluster mass function, and they provide excellent constraints on the formation mechanisms of young cluster populations. Young massive clusters are expected to be the nurseries for many unusual objects, including a wide range of exotic stars and binaries. So far only a few such objects have been found in young massive clusters, although their older cousins, the globular clusters, are unusually rich in stellar exotica. In this review we focus on star clusters younger than $\sim100$ Myr, more than a few current crossing times old, and more massive than $\sim10^4$ \Msun, irrespective of cluster size or environment. We describe the global properties of the currently known young massive star clusters in the Local Group and beyond, and discuss the state of the art in observations and dynamical modeling of these systems. In order to make this review readable by observers, theorists, and computational astrophysicists, we also review the cross-disciplinary terminology.

Understanding the universe is hampered by the elusiveness of its most common constituent, cold dark matter. Almost impossible to observe, dark matter can be studied effectively by means of simulation and there is probably no other research field where simulation has led to so much progress in the last decade. Cosmological N-body simulations are an essential tool for evolving density perturbations in the nonlinear regime. Simulating the formation of large-scale structures in the universe, however, is still a challenge due to the enormous dynamic range in spatial and temporal coordinates, and due to the enormous computer resources required. The dynamic range is generally dealt with by the hybridization of numerical techniques. We deal with the computational requirements by connecting two supercomputers via an optical network and make them operate as a single machine. This is challenging, if only for the fact that the supercomputers of our choice are separated by half the planet, as one is located in Amsterdam and the other is in Tokyo. The co-scheduling of the two computers and the 'gridification' of the code enables us to achieve a 90% efficiency for this distributed intercontinental supercomputer.

We describe the Meta-Institute for Computational Astrophysics (MICA), the first professional scientific organization based exclusively in virtual worlds (VWs). The goals of MICA are to explore the utility of the emerging VR and VWs technologies for scientific and scholarly work in general, and to facilitate and accelerate their adoption by the scientific research community. MICA itself is an experiment in academic and scientific practices enabled by the immersive VR technologies. We describe the current and planned activities and research directions of MICA, and offer some thoughts as to what the future developments in this arena may be.

In this paper we use N-body simulations to study the effects of primordial mass segregation on the early and long-term evolution of star clusters. Our simulations show that in segregated clusters early mass loss due to stellar evolution triggers a stronger expansion than for unsegregated clusters. Tidally limited, strongly segregated clusters may dissolve rapidly as a consequence of this early expansion, while segregated clusters initially underfilling their Roche lobe can survive the early expansion and have a lifetime similar to that of unsegregated clusters. Long-lived initially segregated clusters tend to have looser structure and reach core collapse later in their evolution than initially unsegregated clusters. We have also compared the effects of dynamical evolution on the global stellar mass function (MF) of low-mass main sequence stars. In all cases the MF flattens as the cluster loses stars. The amount of MF flattening induced by a given amount of mass loss in a rapidly dissolving initially segregated cluster is less than for an unsegregated cluster. The evolution of the MF of a long-lived segregated cluster, on the other hand, is very similar to that of an initially unsegregated cluster.

We present MUSE, a software framework for combining existing computational tools for different astrophysical domains into a single multiphysics, multiscale application. MUSE facilitates the coupling of existing codes written in different languages by providing inter-language tools and by specifying an interface between each module and the framework that represents a balance between generality and computational efficiency. This approach allows scientists to use combinations of codes to solve highly-coupled problems without the need to write new codes for other domains or significantly alter their existing codes. MUSE currently incorporates the domains of stellar dynamics, stellar evolution and stellar hydrodynamics for studying generalized stellar systems. We have now reached a "Noah's Ark" milestone, with (at least) two available numerical solvers for each domain. MUSE can treat multi-scale and multi-physics systems in which the time- and size-scales are well separated, like simulating the evolution of planetary systems, small stellar associations, dense stellar clusters, galaxies and galactic nuclei. In this paper we describe three examples calculated using MUSE: the merger of two galaxies, the merger of two evolving stars, and a hybrid N-body simulation. In addition, we demonstrate an implementation of MUSE on a distributed computer which may also include special-purpose hardware, such as GRAPEs or GPUs, to accelerate computations. The current MUSE code base is publicly available as open source at http://muse.li

Some globular clusters are observed to host a population of second generation (SG) stars which show chemical anomalies and must have formed from gas containing matter processed in the envelopes of first generation (FG) cluster stars. We study the SG formation process by means of 1D hydrodynamical simulations, assuming that the SG is formed by the gas ejected by AGB stars. This gas collects in a cooling flow into the cluster core, where it forms a SG star subsystem strongly concentrated in the cluster innermost regions with structural properties largely independent of the FG initial properties. We also present the results of a model in which pristine gas contributes to the SG formation. In this model a very helium-rich SG population and one with a moderate helium enrichment form; the resulting SG bimodal helium distribution resembles that observed for SG stars in NGC 2808. By means of N-body simulations, we study the two-population cluster dynamical evolution. In our simulations, a large fraction of FG stars are lost early in the cluster evolution due to the expansion and stripping of the cluster outer layers resulting from early mass loss associated with FG SN ejecta. The SG population is largely unscathed by this early mass loss, and this early evolution leads to values of the SG to FG number ratio consistent with observations. We also demonstrate possible evolutionary routes leading to the loss of most of the FG population, leaving an SG-dominated cluster. Until mixing of the two populations is complete, the radial profile of the SG to FG number ratio is characterized by a flat inner part and a declining portion in the outer cluster regions. (abridged)

The advent of affordable parallel computers such as Beowulf PC clusters and, more recently, of multi-core PCs has been highly beneficial for a large number of scientists and smaller institutions that might not otherwise have access to substantial computing facilities. However, there has not been an analogous progress in the development and dissemination of parallel software: scientists need the expertise to develop parallel codes and have to invest a significant amount of time in the development of tools even for the most common data analysis tasks. We describe the Beowulf Analysis Symbolic INterface (BASIN) a multi-user parallel data analysis and visualization framework. BASIN is aimed at providing scientists with a suite of parallel libraries for astrophysical data analysis along with general tools for data distribution and parallel operations on distributed data to allow them to easily develop new parallel libraries for their specific tasks.

Internal dynamical evolution can drive stellar systems into states of high central density. For many star clusters and galactic nuclei, the time scale on which this occurs is significantly less than the age of the universe. As a result, such systems are expected to be sites of frequent interactions among stars, binary systems, and stellar remnants, making them efficient factories for the production of compact binaries, intermediate-mass black holes, and other interesting and eminently observable astrophysical exotica. We describe some elements of the competition among stellar dynamics, stellar evolution, and other mechanisms to control the dynamics of stellar systems, and discuss briefly the techniques by which these systems are modeled and studied. Particular emphasis is placed on pathways leading to massive black holes in present-day globular clusters and other potentially detectable sources of gravitational radiation.

We present performance measurements of direct gravitational N -body simulation on the grid, with and without specialized (GRAPE-6) hardware. Our inter-continental virtual organization consists of three sites, one in Tokyo, one in Philadelphia and one in Amsterdam. We run simulations with up to 196608 particles for a variety of topologies. In many cases, high performance simulations over the entire planet are dominated by network bandwidth rather than latency. With this global grid of GRAPEs our calculation time remains dominated by communication over the entire range of N, which was limited due to the use of three sites. Increasing the number of particles will result in a more efficient execution. Based on these timings we construct and calibrate a model to predict the performance of our simulation on any grid infrastructure with or without GRAPE. We apply this model to predict the simulation performance on the Netherlands DAS-3 wide area computer. Equipping the DAS-3 with GRAPE-6Af hardware would achieve break-even between calculation and communication at a few million particles, resulting in a compute time of just over ten hours for 1 N -body time unit. Key words: high-performance computing, grid, N-body simulation, performance modelling

We describe source code level parallelization for the {\tt kira} direct gravitational $N$-body integrator, the workhorse of the {\tt starlab} production environment for simulating dense stellar systems. The parallelization strategy, called ``j-parallelization'', involves the partition of the computational domain by distributing all particles in the system among the available processors. Partial forces on the particles to be advanced are calculated in parallel by their parent processors, and are then summed in a final global operation. Once total forces are obtained, the computing elements proceed to the computation of their particle trajectories. We report the results of timing measurements on four different parallel computers, and compare them with theoretical predictions. The computers employ either a high-speed interconnect, a NUMA architecture to minimize the communication overhead or are distributed in a grid. The code scales well in the domain tested, which ranges from 1024 - 65536 stars on 1 - 128 processors, providing satisfactory speedup. Running the production environment on a grid becomes inefficient for more than 60 processors distributed across three sites.

We simulate open clusters containing up to 182 stars initially in the form of singles, binaries and triples. Due to the high interaction rate a large number of stable quadruples, quintuples, sextuples, and higher-order hierarchies form during the course of the simulations. For our choice of initial conditions, the formation rate of quadruple systems after about 2 Myr is roughly constant with time at $\sim 0.008$ per cluster per Myr. The formation rate of quintuple and sextuple systems are about half and one quarter, respectively, of the quadruple formation rate, and both rates are also approximately constant with time. We present reaction channels and relative probabilities for the formation of persistent systems containing up to six stars. The reaction networks for the formation and destruction of quintuple and sextuple systems can become quite complicated, although the branching ratios remain largely unchanged during the course of the cluster evolution. The total number of quadruples is about a factor of three smaller than observed in the solar neighbourhood.

We study the first 100Myr of the evolution of isolated star clusters initially containing 144179 stars, including 13107 (10%) primordial hard binaries. Our calculations include the effects of both stellar and binary evolution. Gravitational interactions among the stars are computed by direct N-body integration using high precision GRAPE-6 hardware. The evolution of the core radii and central concentrations of our simulated clusters are compared with the observed sample of young (about 100Myr) star clusters in the large Magellanic cloud. Even though our simulations start with a rich population of primordial binaries, core collapse during the early phase of the cluster evolution is not prevented. Throughout the simulations, the fraction of binaries remains roughly constant (about 10%). Due to the effects of mass segregation the mass function of intermediate-mass main-sequence stars becomes as flat as $\alpha=-1.8$ in the central part of the cluster (where the initial Salpeter mass function had $\alpha=-2.35$). About 6--12% of the neutron stars were retained in our simulations; the fraction of retained black holes is 40--70%. In each simulation about three neutron stars become members of close binaries with a main-sequence companion. Such a binary will eventually become an x-ray binary, when the main-sequence star starts to fill its Roche lobe. Black holes are found more frequently in binaries; in each simulated cluster we find about 11 potential x-ray binaries containing a black hole. Abstract abbreviated....

Some young star clusters show a degree of mass segregation that is inconsistent with the effects of standard two-body relaxation from an initially unsegregated system without substructure, in virial equilibrium, and it is unclear whether current cluster formation models can account for this degree of initial segregation in clusters of significant mass. In this Letter we demonstrate that mergers of small clumps that are either initially mass segregated, or in which mass segregation can be produced by two-body relaxation before they merge, generically lead to larger systems which inherit the progenitor clumps' segregation. We conclude that clusters formed in this way are naturally mass segregated, accounting for the anomalous observations and suggesting that this process of prompt mass segregation due to initial clumping should be taken fully into account in constructing cluster dynamical models.

We report the discovery of an X-ray pulsar in the young, massive Galactic star cluster Westerlund 1. We detected a coherent signal from the brightest X-ray source in the cluster, CXO J164710.2-455216, during two Chandra observations on 2005 May 22 and June 18. The period of the pulsar is 10.6107(1) s. We place an upper limit to the period derivative of Pdot<2e-10 s/s, which implies that the spin-down luminosity is Edot<3e33 erg/s. The X-ray luminosity of the pulsar is L_X = 3(+10,-2)e33 (D/5 kpc)^2 erg/s, and the spectrum can be described by a kT = 0.61+/-0.02 keV blackbody with a radius of R_bb = 0.27+/-0.03 (D/5 kpc}) km. Deep infrared observations reveal no counterpart with K<18.5, which rules out binary companions with M>1 Msun. Taken together, the properties of the pulsar indicate that it is a magnetar. The rarity of slow X-ray pulsars and the position of CXO J164710.2-455216 only 1.6' from the core of Westerlund 1 indicates that it is a member of the cluster with >99.97% confidence. Westerlund 1 contains 07V stars with initial masses M_i=35 Msun and >50 post-main-sequence stars that indicate the cluster is 4+/-1 Myr old. Therefore, the progenitor to this pulsar had an initial mass M_i>40 Msun. This is the most secure result among a handful of observational limits to the masses of the progenitors to neutron stars.

The first locations of short gamma-ray bursts (GRBs) in elliptical galaxies suggest they are produced by the mergers of double neutron star (DNS) binaries in old stellar populations. Globular clusters, where the extreme densities of very old stars in cluster cores create and exchange compact binaries efficiently, are a natural environment to produce merging NSs. They also allow some short GRBs to be offset from their host galaxies, as opposed to DNS systems formed from massive binary stars which appear to remain in galactic disks. Starting with a simple scaling from the first DNS observed in a galactic globular, which will produce a short GRB in ~300My, we present numerical simulations which show that ~10-30% of short GRBs may be produced in globular clusters vs. the much more numerous DNS mergers and short GRBs predicted for galactic disks. Reconciling the rates suggests the disk short GRBs are more beamed, perhaps by both the increased merger angular momentum from the DNS spin-orbit alignment (random for the DNS systems in globulars) and a larger magnetic field on the secondary NS.

We simulate the inner 100pc of the Milky-Way Galaxy to study the formation and evolution of the population of star clusters and intermediate mass black holes. For this study we perform extensive direct N-body simulations of the star clusters which reside in the bulge, and of the inner few tenth of parsecs of the super massive black hole in the Galactic center. In our N-body simulations the dynamical friction of the star cluster in the tidal field of the bulge are taken into account via (semi)analytic soluations. The N-body calculations are used to calibrate a (semi)analytic model of the formation and evolution of the bulge. We find that about 10% of the clusters born within 100pc of the Galactic center undergo core collapse during their inward migration and form intermediate-mass black holes (IMBHs) via runaway stellar merging. After the clusters dissolve, these IMBHs continue their inward drift, carrying a few of the most massive stars with them. We predict that region within about 10 parsec of the SMBH is populated by about 50IMBHs of some 1000Msun. Several of these are expected to be accompanied still by some of the most massive stars from the star cluster. We also find that within a few milliparsec of the SMBH there is a steady population of several IMBHs. This population drives the merger rate between IMBHs and the SMBH at a rate of about one per 10Myr, sufficient to build the accumulate majority of mass of the SMBH. Mergers of IMBHs with SMBHs throughout the universe are detectable by LISA, at a rate of about two per week.

We consider the formation of the recently discovered ``hot Jupiter'' planet orbiting the primary component of the triple star system HD188753. Although the current outer orbit of the triple is too tight for a Jupiter-like planet to have formed and migrated to its current location, the binary may have been much wider in the past. We assume here that the planetary system formed in an open star cluster, the dynamical evolution of which subsequently led to changes in the system's orbital parameters and binary configuration. We calculate cross sections for various scenarios that could have led to the multiple system currently observed, and conclude that component A of HD188753 with its planet were most likely formed in isolation to be swapped in a triple star system by a dynamical encounter in an open star cluster. We estimate that within 500pc of the Sun there are about 1200 planetary systems which, like Hd188753, have orbital parameters unfavorable for forming planets but still having a planet, making it quite possible that the HD188753 system was indeed formed by a dynamical encounter in an open star cluster.

Close encounters and physical collisions between stars in young dense clusters may lead to the formation of very massive stars and black holes via runaway merging. We examine critically some details of this process, using N-body simulations and simple analytical estimates to place limits on the cluster parameters for which it expected to occur. For small clusters, the mass of the runaway is effectively limited by the total number of high-mass stars in the system. For sufficiently dense larger clusters, the runaway mass is determined by the fraction of stars that can mass segregate to the cluster core while still on the main sequence. The result is in the range commonly cited for intermediate-mass black holes, such as that recently reported in the Galactic center.

We simulate a cluster of 144179 stars including 13107 primordial hard binaries (10% of the total number of single stars and binary centers of mass), using direct integration of the equations of motion of all stars and binaries and incorporating the effects of stellar and binary evolution. The initial conditions are representative of young dense star clusters in the Local Group and other nearby galaxies like the Antennae and M82. We find that the early phase of core collapse, driven by mass segregation, is not appreciably delayed by the presence of a large number of hard binaries. By the end of the simulation, at an age of 115Myr, the cluster radius has expanded by about a factor of two. This may be explained as adiabatic expansion driven by the loss (via stellar evolution) of about 40% of the initial total mass. Binary dynamics apparently has little effect on the early cluster expansion. During the evolution, the total binary fraction drops at a roughly constant rate of about 0.01% per Myr. The fraction of very hard binaries, however increases at about 0.025% per Myr. By the end of the simulation the cluster contains 37 binaries containing at least one black hole; roughly half (17) of these contain two black holes.

Close encounters and physical collisions between stars in young dense clusters can result in new channels for stellar evolution, and may lead to the formation of very massive stars and black holes via runaway merging. We present some details of this process, using the results of N-body simulations and simple analytical estimates to place limits on the cluster parameters for which it expected to occur. For small clusters, the mass of the runaway is effectively limited by the total number of high-mass stars in the system. For larger clusters, the runaway mass is determined by the fraction of stars that can mass-segregate to the cluster core while still on the main sequence. In typical cases, the result is in the range commonly cited for intermediate-mass black holes. This mechanism may therefore have important implications for the formation of massive black holes and black-hole binaries in dense cluster cores.

We have followed the evolution of multi-mass star clusters containing massive central black holes by N-body simulations on the GRAPE6 computers of Tokyo University. We find a strong cluster expansion and significant structural changes of the clusters. Star clusters with IMBHs have power-law density profiles $rho ~ r^{-alpha}$ with slopes $alpha=1.55$ inside the influence sphere of the central black hole. This leads to a constant density profile of bright stars in projection, which rules out the presence of intermediate mass black holes in core collapse clusters. If the star clusters are surrounded by a tidal field, a central IMBH speeds up the destruction of the cluster until a remnant of a few hundred stars remains, which stays bound to the IMBH for a long time. We also discuss the efficiency of different detection mechanisms for finding IMBHs in star clusters.

We discuss the structural change and degree of mass segregation of young dense star clusters within about 100pc of the Galactic center. In our calculations, which are performed with GRAPE-6, the equations of motion of all stars and binaries are calculated accurately but the external potential of the Galaxy is solved (semi)analytically. The simulations are preformed to model the Arches star cluster. We find that star clusters with are less strongly perturbed by the tidal field and dynamical friction are much stronger affected by mass segregation; resulting in a significant pile-up of massive stars in the cluster center. At an age of about 3.5Myr more than 90 per cent of the stars more massive than ~10Msun are concentrated within the half-mass radius of the surviving cluster. Star clusters which are strongly perturbed by the tidal field of the parent Galaxy are much less affected by mass segregation.

A luminous X-ray source is associated with a cluster (MGG-11) of young stars \~200pc from the center of the starburst galaxy M82. The properties of the X-ray source are best explained by a black hole with a mass of at least 350Msun, which is intermediate between stellar-mass and supermassive black holes. A nearby but somewhat more massive star cluster (MGG-9) shows no evidence of such an intermediate mass black hole, raising the issue of just what physical characteristics of the clusters can account for this difference. Here we report numerical simulations of the evolution and the motions of stars within the clusters, where stars are allowed to mergers with each other. We find that for MGG-11 dynamical friction leads to the massive stars sinking rapidly to the center of the cluster to participate in a runaway collision, thereby producing a star of 800-3000Msun, which ultimately collapses to an black hole of intermediate mass. No such runaway occurs in the cluster MGG-9 because the larger cluster radius leads to a mass-segregation timescale a factor of five longer than for MGG-11.

We know from observations that globular clusters are very efficient catalysts in forming unusual short-period binary systems or their offspring, such as low-mass X-ray binaries (LMXBs; neutron stars accreting matter from low-mass stellar companions), cataclysmic variables (CVs; white dwarfs accreting matter from stellar companions), and millisecond pulsars (MSPs; rotating neutron stars with spin periods of a few ms). Although there has been little direct evidence, the overabundance of these objects in globular clusters has been attributed by numerous authors to the high densities in the cores, which leads to an increase in the formation rate of exotic binary systems through close stellar encounters. Many such close binary systems emit X-radiation at low luminosities (L_x < 10^{34} erg/s) and are being found in large numbers through observations with the Chandra X-ray Observatory. Here we present conclusive observational evidence for a link between the number of close binaries observed in X-rays in a globular cluster and the stellar encounter rate of the cluster. We also make an estimate of the total number of LMXBs in globular clusters in our Galaxy.

We present a comparison between the observational data on the kinematical structure of G1 in M31, obtained with the Hubble WFPC2 and STIS instruments, and the results of dynamical simulations carried out using the special-purpose computer GRAPE-6. We have obtained good fits for models starting with single cluster King-model initial conditions and even better fits when starting our simulations with a dynamically constructed merger remnant of two star clusters. In the latter case, the results from our simulations are in excellent agreement with the observed profiles of luminosity, velocity dispersion, rotation, and ellipticity. We obtain a mass-to-light ratio of $M/L = 4.0 \pm 0.4$ and a total cluster mass of $M=(8\pm 1)\times 10^6 M_\odot$. Given that our dynamical model can fit all available observational data very well, there seems to be no need to invoke the presence of an intermediate-mass black hole in the center of G1.

A star cluster in a galactic nucleus sinks toward the galactic center due to dynamical friction. As it spirals inward, the cluster loses mass due to stellar evolution, relaxation driven evaporation, and tidal stripping, eventually dissolving in the galactic tidal field. We model the inspiral of dense young star clusters near the center of our Galaxy to study the extent of the region of parameter space in which the cluster can reach the inner parsec of the Galaxy within a few million years. Since we neglect changes in cluster structure due to internal evolution, the present study is most applicable to star clusters less than about one initial relaxation time old. We find that only star clusters with initial masses $\apgt 10^5$ {\msun} can reach the Galactic center from an initial distance of $\apgt 60$ pc within one initial relaxation time or a few million years, whichever is smaller.

We use direct N-body simulations to study the inspiral and internal evolution of dense star clusters near the Galactic center. These clusters sink toward the center due to dynamical friction with the stellar background, and may go into core collapse before being disrupted by the Galactic tidal field. If a cluster reaches core collapse before disruption, its dense core, which has become rich in massive stars, survives to reach close to the Galactic center. When it eventually dissolves, the cluster deposits a disproportionate number of massive stars in the innermost parsec of the Galactic nucleus. Comparing the spatial distribution and kinematics of the massive stars with observations of IRS 16, a group of young He I stars near the Galactic center, we argue that this association may have formed in this way.

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.

We have modeled in detail the evolution of rich open star clusters such as the Pleiades, Praesepe and Hyades, using simulations that include stellar dynamics as well as the effects of stellar evolution. The dynamics is modeled via direct N-body integration, while the evolution of single stars and binaries is followed through the use of fitting formulae and recipes. The feedback of stellar and binary evolution on the dynamical evolution of the stellar system is taken into account self-consistently. Our model clusters dissolve in the tidal field of the Galaxy in a time span on the order of a billion years. The rate of mass loss is rather constant, about 1Msun per million years. The binary fraction at first is nearly constant in time, then increases slowly near the end of a cluster's lifetime. For clusters which are more than about 100 million years old the fractions of stars in the form of binaries, giants and collision products in the inner few core radii are considerably higher than in the outer regions, beyond the cluster's half mass radius. When stars with masses >2Msun escape from the cluster, they tend to do so with velocities higher than average.... <Abstract abbreviated, read manuscript for rest>

We present a detailed comparison between the latest observational data on the kinematical structure of the core of M15, obtained with the Hubble STIS and WFPC2 instruments, and the results of dynamical simulations carried out using the special-purpose GRAPE-6 computer. The observations imply the presence of a significant amount of dark matter in the cluster core. In our dynamical simulations, neutron stars and/or massive white dwarfs concentrate to the center through mass segregation, resulting in a sharp increase in $M/L$ toward the center. While consistent with the presence of a central black hole, the Hubble data can also be explained by this central concentration of stellar-mass compact objects. The latter interpretation is more conservative, since such remnants result naturally from stellar evolution, although runaway merging leading to the formation of a black hole may also occur for some range of initial conditions. We conclude that no central massive object is required to explain the observational data, although we cannot conclusively exclude such an object at the level of $\sim500-1000$ solar masses. Our findings are unchanged when we reduce the assumed neutron-star retention fraction in our simulations from 100% to 0%.

We summarize the main results from MODEST-1, the first workshop on MOdeling DEnse STellar systems. Our goal is to go beyond traditional population synthesis models, by introducing dynamical interactions between single stars, binaries, and multiple systems. The challenge is to define and develop a software framework to enable us to combine in one simulation existing computer codes in stellar evolution, stellar dynamics, and stellar hydrodynamics. With this objective, the workshop brought together experts in these three fields, as well as other interested astrophysicists and computer scientists. We report here our main conclusions, questions and suggestions for further steps toward integrating stellar evolution and stellar (hydro)dynamics.

We study the growth rate of stars via stellar collisions in dense star clusters, calibrating our analytic calculations with direct N-body simulations of up to 65536 stars, performed on the GRAPE family of special-purpose computers. We find that star clusters with initial half-mass relaxation times of about 25 Myr are dominated by stellar collisions, the first collisions occurring at or near the point of core collapse, which is driven by the segregation of the most massive stars to the cluster center, where they end up in hard binaries. The majority of collisions occur with the same star, resulting in the runaway growth of a supermassive object. This object can grow up to about 0.1% of the mass of the entire star cluster and could manifest itself as an intermediate-mass black hole (IMBH). The phase of runaway growth lasts until mass loss by stellar evolution arrests core collapse. Star clusters older than about 5 Myr and with present-day half-mass relaxation times less than 100 Myr are expected to contain an IMBH.

The past few years have seen dramatic improvements in the scope and realism of star cluster simulations. Accurate treatments of stellar evolution, coupled with robust descriptions of all phases of binary evolution, have been incorporated self-consistently into several dynamical codes, allowing for the first time detailed study of the interplay between stellar dynamics and stellar physics. The coupling between evolution, dynamics, and the observational appearance of the cluster is particularly strong in young systems and those containing large numbers of primordial binary systems, and important inroads have been made in these areas, particularly in N-body simulations. I discuss some technical aspects of the current generation of N-body integrators, and describe some recent results obtained using these codes.

During the last couple of years, observers have started to make plans for a Virtual Observatory, as a federation of existing data bases, connected through levels of software that enable rapid searches, correlations, and various forms of data mining. We propose to extend the notion of a Virtual Observatory by adding archives of simulations, together with interactive query and visualization capabilities, as well as ways to simulate observations of simulations in order to compare them with observations. For this purpose, we have already organized two small workshops, earlier in 2001, in Tucson and Aspen. We have also provided concrete examples of theory data, designed to be federated with a Virtual Observatory. These data stem from a project to construct an archive for our large-scale simulations using the GRAPE-6 (a 32-Teraflops special purpose computer for stellar dynamics). We are constructing interfaces by which remote observers can observe these simulations. In addition, these data will enable detailed comparisons between different simulations.

We study the evolution and observability of young, compact star clusters near the Galactic center, such as the Arches and Quintuplet systems. The clusters are modeled by integrating the equations of motion of all stars while accounting for the internal evolution of stars and binaries, as well as the effect of the Galactic tidal field. We find that clusters within 150 pc of the Galactic center dissolve within ~55 Myr, but their projected densities drop below the background density in the direction of the Galactic center within only a few Myr, effectively making these clusters undetectable after that time. Detailed observations of the Arches cluster, taken at face value, suggest that its mass function is unusually flat and that the cluster contains an overabundance of stars more massive than 20 Msun. Our dynamical analysis, however, shows that the observed characteristics of the Arches cluster are consistent with a perfectly normal initial mass function. The observed anomalies are then caused by a combination of observational selection effects and the dynamical evolution of the cluster. We calibrate the current parameters of the Arches cluster using a normal initial mass function and conclude that the cluster is more massive than 40000 Msun, has a half mass radius of about 0.35 pc and is located between 50 and 90 pc from the Galactic center.

Observations of stellar kinematics, gas dynamics and masers around galactic nuclei have now firmly established that many galaxies host central supermassive black holes (SMBHs) with masses in the range $10^6 \sim 10^9$M$_{\odot}$. However, how these SMBHs formed is not well understood. One reason for this situation is the lack of observations of intermediate-mass BHs (IMBHs), which could bridge the gap between stellar-mass BHs and SMBHs. Recently, this missing link (i.e., an IMBH) has been found in observations made by the ASCA and the Chandra of the central region of the starburst galaxy M82 \citep{MT99, PG99, MT01, Ka01}. Subsequent observations by SUBARU have revealed that this IMBH apparently coincides with a young compact star cluster. Based on these findings, we suggest a new formation scenario for SMBHs. In this scenario, IMBHs first form in young compact star clusters through runaway merging of massive stars. While these IMBHs are forming, the host star clusters sink toward the galactic nucleus through dynamical friction, and upon evaporation deposit their IMBHs near the galactic center. The IMBHs then form binaries and eventually merge via gravitational radiation, forming an SMBH.

N-body simulations of star cluster evolution have reached a high degree of realism, by incorporating more and more elements of stellar dynamics, stellar evolution, and hydrodynamics. At the end of this conference, six participants discussed how to present the increasingly realistic data from star cluster simulations in a way that is most useful for a direct comparison with observations.

In dense clusters a bewildering variety of interactions between stars can be observed, ranging from simple encounters to collisions and other mass-transfer encounters. With faster and special-purpose computers like GRAPE, the amount of data per simulation is now exceeding 1TB. Visualization of such data has now become a complex 4D data-mining problem, combining space and time, and finding interesting events in these large datasets. We have recently starting using the virtual reality simulator, installed in the Hayden Planetarium in the American Museum for Natural History, to tackle some of these problem. This work (http://www.astro.umd.edu/nemo/amnh/) reports on our first ``observations'', modifications needed for our specific experiments, and perhaps field ideas for other fields in science which can benefit from such immersion. We also discuss how our normal analysis programs can be interfaced with this kind of visualization.

We study the occurrence of physical collisions between stars in young and compact star cluster. The calculations are performed on the GRAPE-4 with the starlab software environment which include the dynamical evolution and the nuclear evolution of all stars and binaries. The selection of the initial conditions is based on existing and well observed star clusters, such as R136 in the 30 Doradus region in the Large Magellanic Cloud and the Arches and Quintuplet star clusters in the vicinity of the Galactic center. Collisions between stars occurred rather frequently in our models. At any time a single star dominates the collision history of the system. The collision rate of this runaway merger scales with the initial relaxation time of the cluster and is independent on other cluster parameters, such as the initial mass function or the initial density profile of the cluster. Subsequent encounters result in a steady grow in mass of the coagulating star, until it escapes or explodes in a supernova. The collision rate in these models is about 0.00022 collisions per star per Myr for a cluster with an initial relaxation time of 1 Myr.

We study the evolution and observability of young compact star clusters within about 200pc of the Galactic center. Calculations are performed using direct N-body integration on the GRAPE-4, including the effects of both stellar and binary evolution and the external influence of the Galaxy. The results of these detailed calculations are used to calibrate a simplified model applicable over a wider range of cluster initial conditions. We find that clusters within 200 pc from the Galactic center dissolve within about 70 Myr. However, their projected densities drop below the background density in the direction of the Galactic center within 20 Myr, effectively making these clusters undetectable after that time. Clusters farther from the Galactic center but at the same projected distance are more strongly affected by this selection effect, and may go undetected for their entire lifetimes. Based on these findings, we conclude that the region within 200 pc of the Galactic center could easily harbor some 50 clusters with properties similar to those of the Arches or the Quintuplet systems.

The evolution of star clusters is studied using N-body simulations in which the evolution of single stars and binaries are taken self-consistently into account. Initial conditions are chosen to represent relatively young Galactic open clusters, such as the Pleiades, Praesepe and the Hyades. The calculations include a realistic mass function, primordial binaries and the external potential of the parent Galaxy. Our model clusters are generally significantly flattened in the Galactic tidal field, and dissolve before deep core collapse occurs. The binary fraction decreases initially due to the destruction of soft binaries, but increases later because lower mass single stars escape more easily than the more massive binaries. At late times, the cluster core is quite rich in giants and white dwarfs. There is no evidence for preferential evaporation of old white dwarfs, on the contrary the formed white dwarfs are likely to remain in the cluster. Stars tend to escape from the cluster through the first and second Lagrange points, in the direction of and away from the Galactic center. Mass segregation manifests itself in our models well within an initial relaxation time. As expected, giants and white dwarfs are much more strongly affected by mass segregation than main-sequence stars. Open clusters are dynamically rather inactive. However, the combined effect of stellar mass loss and evaporation of stars from the cluster potential drives its dissolution on a much shorter timescale than if these effects are neglected. The often-used argument that a star cluster is barely older than its relaxation time and therefore cannot be dynamically evolved is clearly in error for the majority of star clusters.

Mergers of black-hole binaries are expected to release large amounts of energy in the form of gravitational radiation. However, binary evolution models predict merger rates too low to be of observational interest. In this paper we explore the possibility that black holes become members of close binaries via dynamical interactions with other stars in dense stellar systems. In star clusters, black holes become the most massive objects within a few tens of millions of years; dynamical relaxation then causes them to sink to the cluster core, where they form binaries. These black-hole binaries become more tightly bound by superelastic encounters with other cluster members, and are ultimately ejected from the cluster. The majority of escaping black-hole binaries have orbital periods short enough and eccentricities high enough that the emission of gravitational waves causes them to coalesce within a few billion years. We predict a black-hole merger rate of 10^-8 to 10^-7 per year per cubic megaparsec, implying gravity-wave detection rates substantially greater than the corresponding rates from neutron star mergers. For the first generation Laser Interferometer Gravitational-Wave Observatory (LIGO-I), we expect about one detection during the first two years of operation. For its successor LIGO-II, the rate rises to roughly one detection per day. There is about an order of magnitude uncertainty in these numbers.

Mergers of black-hole binaries are expected to release large amounts of energy in the form of gravitational radiation. However, binary evolution models predict merger rates too low to be of observational interest. In this paper we explore the possibility that black holes become members of close binaries via dynamical interactions with other stars in dense stellar systems. In star clusters, black holes become the most massive objects within a few tens of millions of years; dynamical relaxation then causes them to sink to the cluster core, where they form binaries. These black-hole binaries become more tightly bound by superelastic encounters with other cluster members, and are ultimately ejected from the cluster. The majority of escaping black-hole binaries have orbital periods short enough and eccentricities high enough that the emission of gravitational radiation causes them to coalesce within a few billion years. We predict a black-hole merger rate of about $1.6 \times 10^{-7}$ per year per cubic megaparsec, implying gravity wave detection rates substantially greater than the corresponding rates from neutron star mergers. For the first generation Laser Interferometer Gravitational-Wave Observatory (LIGO-I), we expect about one detection during the first two years of operation. For its successor LIGO-II, the rate rises to roughly one detection per day. The uncertainties in these numbers are large. Event rates may drop by about an order of magnitude if the most massive clusters eject their black hole binaries early in their evolution.

The evolution of young compact star clusters is studied using N-body simulations in which both stellar evolution and physical collisions between stars are taken into account. The initial conditions are chosen to represent R136, a compact star cluster in the 30 Doradus region of the Large Magellanic Cloud. The present runs do not include the effects of primordial binaries. We find that physical collisions between stars in these models are frequent, and that the evolution of the most massive stars and the dynamical evolution of the cluster are closely coupled. In all cases, a single star grows steadily in mass through mergers with other stars, forming a very massive (>100 Msun) star in less than 3-4 Myr. The growth rate of this runaway merger is much larger than estimates based on simple cross-section arguments, mainly because the star is typically found in the core and tends to form binaries with other massive stars there. The runaway is ``rejuvenated'' by each new collision, and its lifetime is extended considerably as a consequence. Observationally, such a star will appear in the Hertzsprung-Russell diagram as a blue straggler. When the runaway forms a black hole, the binary in which it is found is usually dissociated.

Using direct N-body simulations which include both the evolution of single stars and the tidal field of the parent galaxy, we study the dynamical evolution of globular clusters and rich open clusters. We compare our results with other N-body simulations and Fokker-Planck calculations. Our simulations, performed on the GRAPE-4, employ up to 32,768 stars. The results are not in agreement with Fokker-Planck models, in the sense that the lifetimes of stellar systems derived using the latter are an order of magnitude smaller than those obtained in our simulations. For our standard run, Fokker-Plank calculations obtained a lifetime of 0.28 Gyr, while our equivalent N-body calculations find about 4 Gyr. The principal reason for the discrepancy is that a basic assumption of the Fokker-Plank approach is not valid for typical cluster parameters. The stellar evolution timescale is comparable to the dynamical timescale, and therefore the assumption of dynamical equilibrium leads to an overestimate of the dynamical effects of mass loss. Our results suggest that the region in parameter space for which Fokker-Planck studies of globular cluster evolution, including the effects of both stellar evolution and the galactic tidal field, are valid is limited. The discrepancy is largest for clusters with short lifetimes.

Binary stars play a critical role in determining both the dynamical evolution and the observational appearance of globular star clusters. We present a summary of three talks presented at IAU JD15, reviewing the current observational and theoretical status of binaries and related ``exotic'' stars in the Galactic globular cluster system.

We present the results from a two-day study in which we discussed various implementations of Smooth Particle Hydrodynamics (SPH), one of the leading methods used across a variety of areas of large-scale astrophysical simulations. In particular, we evaluated the suitability of designing special hardware extensions, to further boost the performance of the high-end general purpose computers currently used for those simulations. We considered a range of hybrid architectures, consisting of a mix of custom LSI and reconfigurable logic, combining the extremely high throughput of Special-Purpose Devices (SPDs) with the flexibility of reconfigurable structures, based on Field Programmable Gate Arrays (FPGAs). The main findings of our workshop consist of a clarification of the decomposition of the computational requirements, together with specific estimates for cost/performance improvements that can be obtained at each stage in this decomposition, by using enabling hardware technology to accelerate the performance of general purpose computers.

Three-body effects greatly complicate stellar evolution. We model the effects of encounters of binaries with single stars, based on parameters chosen from conditions prevalent in the cores of globular clusters. For our three-body encounters, we start with a population of primordial binaries, and choose incoming stars from the evolving single star populations of Portgies Zwart, Hut and Verbunt (1997). In addition, we study the formation of new binaries through tidal capture among the single stars. In subsequent papers in this series, we will combine stellar evolution with the dynamics of a full N-body system.

We present two types of meta-algorithm that can greatly improve the accuracy of existing algorithms for integrating the equations of motion of dynamical systems. The first meta-algorithm takes an integrator that is time-symmetric only for constant time steps, and ensures time-symmetry even in the case of varying time steps. The second meta-algorithm can be applied to any self-starting integration scheme to create time-symmetry, both for constant and for variable time steps, even if the original scheme was not time-symmetric. Our meta-algorithms are most effective for Hamilton systems or systems with periodic solutions. If the system is not Hamiltonian (for example, if some dissipative force exists), our methods are still useful so long as the dissipation is small.

We present the outline of a research project aimed at designing and constructing a hybrid computing system that can be easily scaled up to petaflops speeds. As a first step, we envision building a prototype which will consist of three main components: a general-purpose, programmable front end, a special-purpose, fully hardwired computing engine, and a multi-purpose, reconfigurable system. The driving application will be a suite of particle-based large-scale simulations in various areas of physics. The prototype system will achieve performance in the $\sim 50 - 100$ teraflops range for a broad class of applications in this area. The combination of a hardwired petaflops-class computational engine and a front end with sustained speed on the order of 10 gigaflops can produce extremely high performance, but only for the limited class of problems in which there exists a single bottleneck with computing cost dominating the total. While the calculation for which the Grape-4 (our system's immediate predecessor) was designed is a prime example of such a problem, in many other applications the primary computational bottleneck, while still related to an inverse-square (gravitational, Coulomb, etc.) force, requires less than 99% of the computing power. Although the remainder of the CPU time is typically dominated by just one secondary bottleneck, its nature varies greatly from problem to problem. It is not cost-effective to attempt to design custom chips for each new problem that arises. FPGA-based systems can restore the balance, guaranteeing scalability from the teraflops to the petaflops domain, while still retaining significant flexibility. (abbreviated abstract)

We calculate the gravitational radiation produced by the coalescence of inspiraling binary neutron stars in the Newtonian regime using 3-dimensional numerical simulations. The stars are modeled as polytropes and start out in the point-mass regime at wide separation. The hydrodynamic integration is performed using smooth particle hydrodynamics (SPH) with Newtonian gravity, and the gravitational radiation is calculated using the quadrupole approximation. We have run a number of simulations varying the neutron star radii, equations of state, spins, and mass ratio. The resulting gravitational waveforms and spectra are rich in information about the hydrodynamics of coalescence, and show characteristic dependence on GM/Rc^2, the equation of state, and the mass ratio.

In this paper we describe a new algorithm for the long-term numerical integration of the two-body problem, in which two particles interact under a Newtonian gravitational potential. Although analytical solutions exist in the unperturbed and weakly perturbed cases, numerical integration is necessary in situations where the perturbation is relatively strong. Kustaanheimo--Stiefel (KS) regularization is widely used to remove the singularity in the equations of motion, making it possible to integrate orbits having very high eccentricity. However, even with KS regularization, long-term integration is difficult, simply because the required accuracy is usually very high. We present a new time-integration algorithm which has no secular error in either the binding energy or the eccentricity, while allowing variable stepsize. The basic approach is to take a time-symmetric algorithm, then apply an implicit criterion for the stepsize to ensure strict time reversibility. We describe the algorithm in detail and present the results of numerical tests involving long-term integration of binaries and hierarchical triples. In all cases studied, we found no systematic error in either the energy or the angular momentum. We also found that its calculation cost does not become higher than those of existing algorithms. By contrast, the stabilization technique, which has been widely used in the field of collisional stellar dynamics, conserves energy very well but does not conserve angular momentum.

We present the first comprehensive fitting formula for exchange reactions of arbitrary mass ratios. In a comparison with numerical results, this expression is shown to be accurate in the hard binary limit to within 25\% for most mass ratios. The result will be useful in forming quantitative estimates for the branching ratios of various exchange reactions in astrophysical applications. For example, it can be used to construct quantitative formation scenarios for unusual objects in globular clusters, such as binaries containing a pulsar.

Scattering encounters between binaries and single stars play a central role in determining the dynamical evolution of a star cluster. In addition, three-body scattering can give rise to many interesting exceptional objects: merging can produce blue stragglers; exchange can produce binaries containing millisecond pulsars in environments quite different from those in which the pulsars were spun up; various types of X-ray binaries can be formed, and their activity can be either shut off or triggered as a result of triple interactions. To date, all published results on three-body scattering have relied on human guidance for determining the correct parameter range for the envelope within which to perform Monte--Carlo scattering experiments. In this paper, we describe the first fully automatic determination of cross sections and reaction rates for binary--single-star scattering. Rather than relying on human inspection of pilot calculations, we have constructed a feedback system that ensures near-optimal coverage of parameter space while guaranteeing completeness. We illustrate our approach with a particular example, in which we describe the results of a three-body encounter between three main-sequence stars of different masses. We provide total cross sections, as well as branching ratios for the various different types of two-body mergers, three-body mergers, and exchanges, both non-resonant and resonant. The companion paper in this series, Paper VII, provides a full survey of unequal-mass three-body scattering for hard binaries in the point-mass limit.

We present results on the exchange cross section for the interaction between a hard binary and a field of single stars, for arbitrary masses. The results are based partly on extensive numerical scattering experiments, and partly on analytic estimates of the mass-dependence of the cross section. They can be used to estimate the rate of exchange in an arbitrary mixture of masses, provided that the binary is hard.

We calculate the gravitational radiation produced by the merger and coalescence of inspiraling binary neutron stars using 3-dimensional numerical simulations. The stars are modeled as polytropes and start out in the point-mass limit at wide separation. The hydrodynamic integration is performed using smooth particle hydrodynamics (SPH) with Newtonian gravity, and the gravitational radiation is calculated using the quadrupole approximation. We have run several simulations, varying both the neutron star radius and the equation of state. The resulting gravitational wave energy spectra $dE/df$ are rich in information about the hydrodynamics of merger and coalescence. In particular, our results demonstrate that detailed information on both $GM/Rc^2$ and the equation of state can in principle be extracted from the spectrum.

The hierarchical triple system containing the millisecond pulsar PSR B1620-26 in M4 is the first triple star system ever detected in a globular cluster. Such systems should form in globular clusters as a result of dynamical interactions between binaries. We propose that the triple system containing PSR B1620-26 formed through an exchange interaction between a wide primordial binary and a {\it pre-existing\/} binary millisecond pulsar. This scenario would have the advantage of reconciling the $\sim10^9\,$yr timing age of the pulsar with the much shorter lifetime of the triple system in the core of M4.

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