Interacting galaxies often have complexes of hundreds of young stellar clusters of individual masses $\sim 10^{4-6}~M_\odot$ in regions that are a few hundred parsecs across. These cluster complexes interact dynamically, and their coalescence is a candidate for the origin of some ultracompact dwarf galaxies (UCDs). Individual clusters with short relaxation times are candidates for the production of intermediate-mass black holes of a few hundred solar masses, via runaway stellar collisions prior to the first supernovae in a cluster. It is therefore possible that a cluster complex hosts multiple intermediate-mass black holes that may be ejected from their individual clusters due to mergers or binary processes, but bound to the complex as a whole. Here we explore the dynamical interaction between initially free-flying massive black holes and clusters in an evolving cluster complex. We find that, after hitting some clusters, it is plausible that the massive black hole will be captured in an ultracompact dwarf forming near the center of the complex. In the process, the hole typically triggers electromagnetic flares via stellar disruptions, and is also likely to be a prominent source of gravitational radiation for the advanced ground-based detectors LIGO and VIRGO. We also discuss other implications of this scenario, notably that the central black hole could be considerably larger than expected in other formation scenarios for ultracompact dwarfs.

The last century has seen enormous progress in our understanding of the Universe. We know the life cycles of stars, the structure of galaxies, the remnants of the big bang, and have a general understanding of how the Universe evolved. We have come remarkably far using electromagnetic radiation as our tool for observing the Universe. However, gravity is the engine behind many of the processes in the Universe, and much of its action is dark. Opening a gravitational window on the Universe will let us go further than any alternative. Gravity has its own messenger: Gravitational waves, ripples in the fabric of spacetime. They travel essentially undisturbed and let us peer deep into the formation of the first seed black holes, exploring redshifts as large as z ~ 20, prior to the epoch of cosmic re-ionisation. Exquisite and unprecedented measurements of black hole masses and spins will make it possible to trace the history of black holes across all stages of galaxy evolution, and at the same time constrain any deviation from the Kerr metric of General Relativity. eLISA will be the first ever mission to study the entire Universe with gravitational waves. eLISA is an all-sky monitor and will offer a wide view of a dynamic cosmos using gravitational waves as new and unique messengers to unveil The Gravitational Universe. It provides the closest ever view of the early processes at TeV energies, has guaranteed sources in the form of verification binaries in the Milky Way, and can probe the entire Universe, from its smallest scales around singularities and black holes, all the way to cosmological dimensions.

A common result to all EMRI investigations on rates is that the possibility that a compact object merges with the MBH after only one intense burst of GWs is much more likely than a slow adiabatic inspiral, an EMRI. The later is referred to as a "plunge" because the compact object dives into the MBH. The event rates for plunges are orders of magnitude larger than slow inspirals. On the other hand, nature MBH's are most likely Kerr and the magnitude of the spin has been sized up to be high. We calculate the number of periapsis passages that a compact object set on to an extremely radial orbit goes through before being actually swallowed by the Kerr MBH and we then translate it into an event rate for a LISA-like observatory, such as the proposed ESA mission eLISA/NGO. We prove that a "plunging" compact object is conceptually indistinguishable from an adiabatic, slow inspiral; plunges spend on average up to hundred of thousands of cycles in the bandwidth of the detector for a two years mission. This has an important impact on the event rate, enhancing in some cases significantly, depending on the spin of the MBH and the inclination. Moreover, it has been recently proved that the production of low-eccentricity EMRIs is severely blocked by the presence of a blockade in the rate at which orbital angular momenta change takes place. This is the result of relativistic precession on to the stellar potential torques and hence affects EMRIs originating via resonant relaxation at distances of about $\sim 10^{-2}$ pc from the MBH. Since high-eccentricity EMRIs are a result of two-body relaxation, they are not affected by this phenomenon. Therefore we predict that eLISA EMRI event rates will be dominated by high-eccentricity binaries, as we present here.

Two coalescing black holes (BHs) represent a conspicuous source of gravitational waves (GWs). The merger involves 17 parameters in the general case of Kerr BHs, so that a successful identification and parameter extraction of the information encoded in the waves will provide us with a detailed description of the physics of BHs. A search based on matched-filtering for characterization and parameter extraction requires the development of some $10^{15}$ waveforms. If a third additional BH perturbed the system, the waveforms would not be applicable, and we would need to increase the number of templates required for a valid detection. In this letter, we calculate the probability that more than two BHs interact in the regime of strong relativity in a dense stellar cluster. We determine the physical properties necessary in a stellar system for three black holes to have a close encounter in this regime and also for an existing binary of two BHs to have a strong interaction with a third hole. In both cases the event rate is negligible. While dense stellar systems such as galactic nuclei, globular clusters and nuclear stellar clusters are the breeding grounds for the sources of gravitational waves that ground-based detectors like Advanced LIGO and Advanced VIRGO will be exploring, the analysis of the waveforms in full general relativity needs only to evaluate the two-body problem. This reduces the number of templates of waveforms to create by orders of magnitude.

We explore the behaviour of accreting protoclusters with a Monte Carlo dynamical code in order to evaluate the relative roles of accretion, two body relaxation and stellar collisions in the cluster evolution. We corroborate the suggestion of Clarke & Bonnell that the number of stellar collisions should scale as $N^{5/3} \dot M^{2/3}$ (independent of other cluster parameters, where N is the number of stars in the cluster and $\dot M$ the rate of mass accretion) and thus strengthen the argument that stellar collisions are more likely in populous (large N) clusters. We however find that the estimates of Clarke & Bonnell were pessimistic in the sense that we find that more than 99 % of the stellar collisions occur within the post-adiabatic regime as the cluster evolves towards core collapse, driven by a combination of accretion and two-body relaxation. We discuss how the inclusion of binaries may reduce the number of collisions through the reversal of core collapse but also note that it opens up another collisional channel involving the merger of stars within hard binaries; future Nbody simulations are however required in order to explore this issue.

We consider whether stellar collisions can explain the observed depletion of red giants in the Galactic center. We model the stellar population with two different IMFs: 1) the Miller-Scalo and 2) a much flatter IMF. In the former case, low-mass main-sequence stars dominate the population, and collisions are unable to remove red giants out to 0.4 pc although brighter red giants much closer in may be depleted via collisions with stellar-mass black holes. For a much flatter IMF, the stellar population is dominated by compact remnants (i.e. black holes, white dwarfs and neutron stars). The most common collisions are then those between main-sequence stars and compact remnants. Such encounters are likely to destroy the main-sequence stars and thus prevent their evolution into red giants. In this way, the red-giant population could be depleted out to 0.4 pc matching observations. If this is the case, it implies the Galactic center contains a much larger population of stellar-mass black holes than would be expected from a regular IMF. This may in turn have implications for the formation and growth of the central supermassive black hole.

Stars on eccentric orbits around a massive black hole (MBH) emit bursts of gravitational waves (GWs) at periapse. Such events may be directly resolvable in the Galactic centre. However, if the star does not spiral in, the emitted GWs are not resolvable for extra-galactic MBHs, but constitute a source of background noise. We estimate the power spectrum of this extreme mass ratio burst background (EMBB) and compare it to the anticipated instrumental noise of the Laser Interferometer Space Antenna (LISA). To this end, we model the regions close to a MBH, accounting for mass-segregation, and for processes that limit the presence of stars close to the MBH, such as GW inspiral and hydrodynamical collisions between stars. We find that the EMBB is dominated by GW bursts from stellar mass black holes, and the magnitude of the noise spectrum (f S_GW)^{1/2} is at least a factor ~10 smaller than the instrumental noise. As an additional result of our analysis, we show that LISA is unlikely to detect relativistic bursts in the Galactic centre.

If binary intermediate-mass black holes (IMBHs; with masses between 100 and $10^4 \Msun$) form in dense stellar clusters, their inspiral will be detectable with the planned Laser Interferometer Space Antenna (LISA) out to several Gpc. Here we present a study of the dynamical evolution of such binaries using a combination of direct $N$-body techniques (when the binaries are well separated) and three-body relativistic scattering experiments (when the binaries are tight enough that interactions with stars occur one at a time). We find that for reasonable IMBH masses there is only a mild effect on the structure of the surrounding cluster even though the binary binding energy can exceed the binding energy of the cluster. We demonstrate that, contrary to standard assumptions, the eccentricity in the LISA band can be in {\em some} cases as large as $\sim 0.2 - 0.3$ and that it induces a measurable phase difference from circular binaries in the last year before merger. We also show that, even though energy input from the binary decreases the density of the core and slows down interactions, the total time to coalescence is short enough (typically less than a hundred million years) that such mergers will be unique snapshots of clustered star formation.

We show that collisions with stellar--mass black holes can partially explain the absence of bright giant stars in the Galactic Centre, first noted by Genzel et al, 1996. We show that the missing objects are low--mass giants and AGB stars in the range 1-3 M$_{\odot}$. Using detailed stellar evolution calculations, we find that to prevent these objects from evolving to become visible in the depleted K bands, we require that they suffer collisions on the red giant branch, and we calculate the fractional envelope mass losses required. Using a combination of Smoothed Particle Hydrodynamic calculations, restricted three--body analysis and Monte Carlo simulations, we compute the expected collision rates between giants and black holes, and between giants and main--sequence stars in the Galactic Centre. We show that collisions can plausibly explain the missing giants in the $10.5<K<12$ band. However, depleting the brighter ($K<10.5$) objects out to the required radius would require a large population of black hole impactors which would in turn deplete the $10.5<K<12$ giants in a region much larger than is observed. We conclude that collisions with stellar--mass black holes cannot account for the depletion of the very brightest giants, and we use our results to place limits on the population of stellar--mass black holes in the Galactic Centre.

In young star clusters, the density can be high enough and the velocity dispersion low enough for stars to collide and merge with a significant probability. This has been suggested as a possible way to build up the high-mass portion of the stellar mass function and as a mechanism leading to the formation of one or two very massive stars (M > 150 Msun) through a collisional runaway. I quickly review the standard theory of stellar collisions, covering both the stellar dynamics of dense clusters and the hydrodynamics of encounters between stars. The conditions for collisions to take place at a significant rate are relatively well understood for idealised spherical cluster models without initial mass segregation, devoid of gas and composed of main-sequence (MS) stars. In this simplified situation, 2-body relaxation drives core collapse through mass segregation and a collisional phase ensues if the core collapse time is shorter than the MS lifetime of the most massive stars initially present. The outcome of this phase is still highly uncertain. A more realistic situation is that of a cluster still containing large amounts of interstellar gas from which stars are accreting. As stellar masses increase, the central regions of the cluster contracts. This little-explored mechanism can potentially lead to very high stellar densities but it is likely that, except for very rich systems, the contraction is halted by few-body interactions before collisions set in. A complete picture, combining both scenarios, will need to address many uncertainties, including the role of cluster sub-structure, the dynamical effect of interstellar gas, non-MS stars and the structure and evolution of merged stars.

Massive black holes (MBHs) with a mass below ~ 1e7 Msun are likely to reside at the centre of dense stellar nuclei shaped by 2-body relaxation, close interactions with the MBH and direct collisions. In this contribution, we stress the role of mass segregation of stellar-mass black holes into the innermost tenths of a parsec and point to the importance of hydrodynamical collisions between stars. At the Galactic centre, collisions must affect giant stars and some of the S-stars.

Black hole binaries with extreme ($\gtrsim 10^4:1$) or intermediate ($\sim 10^2-10^4:1$) mass ratios are among the most interesting gravitational wave sources that are expected to be detected by the proposed Laser Interferometer Space Antenna. These sources have the potential to tell us much about astrophysics, but are also of unique importance for testing aspects of the general theory of relativity in the strong field regime. Here we discuss these sources from the perspectives of astrophysics, data analysis, and applications to testing general relativity, providing both a description of the current state of knowledge and an outline of some of the outstanding questions that still need to be addressed. This review grew out of discussions at a workshop in September 2006 hosted by the Albert Einstein Institute in Golm, Germany.

The Galactic massive black hole (MBH), with a mass of Mbh=3.6\times10^6 Solar masses, is the closest known MBH, at a distance of only 8 kpc. The proximity of this MBH makes it possible to observe gravitational waves from stars with periapse in the observational frequency window of the Laser Interferometer Space Antenna (LISA). This is possible even if the orbit of the star is very eccentric, so that the orbital frequency is many orders of magnitude below the LISA frequency window, as suggested by Rubbo et al. (2006). Here we give an analytical estimate of the detection rate of such gravitational wave bursts. The burst rate is critically sensitive to the inner cut-off of the stellar density profile. Our model accounts for mass-segregation and for the physics determining the inner radius of the cusp, such as stellar collisions, energy dissipation by gravitational wave emission, and consequences of the finite number of stars. We find that stellar black holes have a burst rate of the order of 1 per year, while the rate is of order 0.1 per year for main sequence stars and white dwarfs. These analytical estimates are supported by a series of Monte Carlo samplings of the expected distribution of stars around the Galactic MBH, which yield the full probability distribution for the rates. We estimate that no burst will be observable from the Virgo cluster.

Observations suggest that star clusters often form in binaries or larger bound groups. Therefore, mergers between two clusters are likely to occur. If these clusters both harbor an intermediate-mass black hole (IMBH; 10^{2-4} Msun) in their center, they can become a strong source of gravitational waves when the black holes merge with each other. In order to understand the dynamical processes that operate in such a scenario, one has to study the evolution of the merger of two such young massive star clusters, and more specifically, their respective IMBHs. We employ the direct-summation Nbody4 numerical tool on special-purpose GRAPE6 hardware to simulate a merger of two stellar clusters each containing 63,000 particles and a central IMBH. This allows us to study accurately the orbital evolution of the colliding clusters and the embedded massive black holes. Within ~7 Myr the clusters have merged and the IMBHs constitute a hard binary. The final coalescence happens in ~10^8 yrs. The implication of our analysis is that intermediate-mass black holes merging as the result of coalescence of young dense clusters could provide a source for the Laser Interferometer Space Antenna (LISA) space-based gravitational wave detector mission. We find that interactions with stars increase the eccentricity of the IMBH binary to about 0.8. Although the binary later circularizes by emission of gravitational waves, the residual eccentricity can be detectable through its influence on the phase of the waves if the last few years of inspiral are observed.

We present a review of some of the current major challenges in stellar cluster research, including young clusters, globular clusters, and galactic nuclei. Topics considered include: primordial mass segregation and runaway mergers, expulsion of gas from clusters, the production of stellar exotica seen in some clusters (eg blue stragglers and extreme horizontal--branch stars), binary populations within clusters, the black--hole population within stellar clusters, the final parsec problem, stellar dynamics around a massive black hole, and stellar collisions. The Modest Questions posed here are the outcome of discussions which took place at the Modest-6A workshop held in Lund, Sweden, in December, 2005. Modest-6A was organised as part of the activities of the Modest Collaboration (see www.manybody.org for further details)

We study the process of mass segregation through 2-body relaxation in galactic nuclei with a central massive black hole (MBH). This study has bearing on a variety of astrophysical questions, from the distribution of X-ray binaries at the Galactic centre, to tidal disruptions of main-sequence and giant stars, to inspirals of compact objects into the MBH, an important category of events for the future space borne gravitational wave interferometer LISA. In relatively small galactic nuclei, typical hosts of MBHs with masses in the range 1e4-1e7 Msun, the relaxation induces the formation of a steep density cusp around the MBH and strong mass segregation. Using a spherical stellar dynamical Monte-Carlo code, we simulate the long-term relaxational evolution of galactic nucleus models with a spectrum of stellar masses. Our focus is the concentration of stellar black holes to the immediate vicinity of the MBH. Special attention is given to models developed to match the conditions in the Milky Way nucleus.

The study of how stars distribute themselves around a massive black hole (MBH) in the center of a galaxy is an important prerequisite for the understanding of many galactic-center processes. These include the observed overabundance of point X-ray sources at the Galactic center, the prediction of rates and characteristics of tidal disruptions of extended stars by the MBH and of inspirals of compact stars into the MBH, the latter being events of high importance for the future space borne gravitational wave interferometer LISA. In relatively small galactic nuclei, hosting MBHs with masses in the range 10^5-10^7 Msun, the single most important dynamical process is 2-body relaxation. It induces the formation of a steep density cusp around the MBH and strong mass segregation, as more massive stars lose energy to lighter ones and drift to the central regions. Using a spherical stellar dynamical Monte-Carlo code, we simulate the long-term relaxational evolution of galactic nucleus models with a spectrum of stellar masses. Our focus is the concentration of stellar black holes to the immediate vicinity of the MBH. We quantify this mass segregation for a variety of galactic nucleus models and discuss its astrophysical implications. Special attention is given to models developed to match the conditions in the Milky Way nucleus; we examine the presence of compact objects in connection to recent high-resolution X-ray observations.

We present the methods and preparatory work for our study of the collisional runaway scenario to form a very massive star (VMS, M > 400 Msun) at the centre of a young, compact stellar cluster. In the first phase of the process, a very dense central core of massive stars (M =~ 30-120 Msun) forms through mass segregation and gravothermal collapse. This leads to a collisional stage, likely to result in the formation of a VMS (itself a possible progenitor for an intermediate-mass black hole) through a runaway sequence of mergers between the massive stars. In this paper we present the runaway scenario in a general astrophysical context. We then explain the numerical method used to investigate it. Our approach is based on a Monte Carlo code to simulate the stellar dynamics of spherical star clusters using a very large number of particles (a few 1e5 to several 1e6). Finally, we report on test computations carried out to ensure that our implementation of the important physics is sound. In a second paper, we present results from more than 100 cluster simulations realized to determine the conditions leading to the collisional formation of a VMS and the characteristics of the runaway sequences.

We present a new study of the collisional runaway scenario to form an intermediate-mass black hole (IMBH, MBH > 100 Msun) at the centre of a young, compact stellar cluster. The first phase is the formation of a very dense central core of massive stars (Mstar =~ 30-120 Msun) through mass segregation and gravothermal collapse. Previous work established the conditions for this to happen before the massive stars evolve off the main sequence (MS). In this and a companion paper, we investigate the next stage by implementing direct collisions between stars. Using a Monte Carlo stellar dynamics code, we follow the core collapse and subsequent collisional phase in more than 100 models with varying cluster mass, size, and initial concentration. Collisions are treated either as ideal, ``sticky-sphere'' mergers or using realistic prescriptions derived from 3-D hydrodynamics computations. In all cases for which the core collapse happens in less than the MS lifetime of massive stars (~3 Myr), we obtain the growth of a single very massive star (VMS, Mstar =~ 400-4000 Msun) through a runaway sequence of mergers. Mass loss from collisions, even for velocity dispersions as high as sigma1D ~ 1000 km/s, does not prevent the runaway. The region of cluster parameter space leading to runaway is even more extended than predicted in previous work because, in clusters with sigma1D > 300 km/s, collisions accelerate (and, in extreme cases, drive) core collapse. Although the VMS grows rapidly to > 1000 Msun in models exhibiting runaway, we cannot predict accurately its final mass. This is because the termination of the runaway process must eventually be determined by a complex interplay between stellar dynamics, hydrodynamics, and the stellar evolution of the VMS. [abridged]

Current simulations of the rate at which stellar-mass compact objects merge with supermassive black holes (called extreme mass ratio inspirals, or EMRIs) focus on two-body capture by emission of gravitational radiation. The gravitational wave signal of such events will likely involve a significant eccentricity in the sensitivity range of the Laser Interferometer Space Antenna (LISA). We show that tidal separation of stellar-mass compact object binaries by supermassive black holes will instead produce events whose eccentricity is nearly zero in the LISA band. Compared to two-body capture events, tidal separations have a high cross section and result in orbits that have a large pericenter and small apocenter. Therefore, the rate of interactions per binary is high and the resulting systems are very unlikely to be perturbed by other stars into nearly radial plunges. Depending on the fraction of compact objects that are in binaries within a few parsecs of the center, the rate of low-eccentricity LISA events could be comparable to or larger than the rate of high-eccentricity events.

We report on a very large set of simulations of collisions between two main sequence (MS) stars. These computations were done with the ``Smoothed Particle Hydrodynamics'' method. Realistic stellar structure models for evolved MS stars were used. In order to sample an extended domain of initial parameters space (masses of the stars, relative velocity and impact parameter), more than 15000 simulations were carried out. We considered stellar masses ranging between 0.1 and 75 Msun and relative velocities up to a few thousands km/s. To limit the computational burden, a resolution of 2000-30000 particles per star was used. The primary goal of this study was to build a complete database from which the result of any collision can be interpolated. This allows us to incorporate the effects of stellar collisions with an unprecedented level of realism into dynamical simulations of galactic nuclei and other dense stellar clusters. We make the data describing the initial condition and outcome (mass and energy loss, angle of deflection) of all our simulations freely available on the Internet. We find that the outcome of collisions depends sensitively on the stellar structure and that, in most cases, using polytropic models is inappropriate. Published fitting formulas for the collision outcomes, established from a limited set of collisions, prove of limited use because they do not allow robust extrapolation to other stellar structures or relative velocities.

Using Monte Carlo codes, we follow the collisional evolution of clusters in a variety of scenarios. We consider the conditions under which a cluster of main sequence stars may undergo rapid core collapse due to mass segregation, thus entering a phase of runaway collisions, forming a very massive star (VMS, M >= 1000 Msun) through repeated collisions between single stars. Although collisional mass loss is accounted for realistically, we find that a VMS forms even in proto-galactic nuclei models with a high velocity dispersion (many 100 km/s). Such a VMS may be a progenitor for an intermediate-mass black hole (M >= 100 Msun). In contrast, in galactic nuclei hosting a central massive black hole, collisions are found to be disruptive. The stars which are subject to collisions are progressively ground down by high-velocity collisions and a merger sequence appears impossible.

To be presented is a study of the secular evolution of a spherical stellar system with a central star-accreting black hole (BH) using the anisotropic gaseous model. This method solves numerically moment equations of the full Fokker-Planck equation, with Boltzmann-Vlasov terms on the left- and collisional terms on the right-hand sides. We study the growth of the central BH due to star accretion at its tidal radius and the feedback of this process on to the core collapse as well as the post-collapse evolution of the surrounding stellar cluster in a self-consistent manner. Diffusion in velocity space into the loss-cone is approximated by a simple model. The results show that the self-regulated growth of the BH reaches a certain fraction of the total mass cluster and agree with other methods. Our approach is much faster than competing ones (Monte Carlo, $N$--body) and provides detailed informations about the time and space dependent evolution of all relevant properties of the system. In this work we present the method and study simple models (equal stellar masses, no stellar evolution or collisions). Nonetheless, a generalisation to include such effects is conceptually simple and under way.

We consider spherical stellar clusters with a broad mass function and a relaxation time short enough so that the segregation of massive stars toward the centre occurs before they have time to evolve off the main sequence. The relaxational and collisional dynamics of model clusters is followed with a high-resolution Monte Carlo code. Stellar collisions are treated in a realistic way, through the use of a outcome of a very large set of SPH simulations. We find that, even in proto-galactic nuclei models with high velocity dispersions, run-away growth of a very massive star (VMS, M>100 M_sun) occurs in all cases when the core collapse time is shorter than the MS life time of massive stars, i.e. 3 Myrs. The VMS is a likely progenitor for an intermediate-mass or massive black hole (IMBH/MBH).

We study the early dynamical evolution of young, dense star clusters using Monte Carlo simulations for systems with up to N~10^7 stars. Rapid mass segregation of massive main-sequence stars and the development of the Spitzer instability can drive these systems to core collapse in a small fraction of the initial half-mass relaxation time. If the core collapse time is less than the lifetime of the massive stars, all stars in the collapsing core may then undergo a runaway collision process leading to the formation of a massive black hole. Here we study in detail the first step in this process, up to the occurrence of core collapse. We have performed about 100 simulations for clusters with a wide variety of initial conditions, varying systematically the cluster density profile, stellar IMF, and number of stars. We also considered the effects of initial mass segregation and stellar evolution mass loss. Our results show that, for clusters with a moderate initial central concentration and any realistic IMF, the ratio of core collapse time to initial half-mass relaxation time is typically ~0.1, in agreement with the value previously found by direct N-body simulations for much smaller systems. Models with even higher central concentration initially, or with initial mass segregation (from star formation) have even shorter core-collapse times. Remarkably, we find that, for all realistic initial conditions, the mass of the collapsing core is always close to ~10^-3 of the total cluster mass, very similar to the observed correlation between central black hole mass and total cluster mass in a variety of environments. We discuss the implications of our results for the formation of intermediate-mass black holes in globular clusters and super star clusters, ultraluminous X-ray sources, and seed black holes in proto-galactic nuclei.

First we present a recently developed 3D chemodynamical code for galaxy evolution from the K**2 collaboration. It follows the evolution of all components of a galaxy such as dark matter, stars, molecular clouds and diffuse interstellar matter (ISM). Dark matter and stars are treated as collisionless N-body systems. The ISM is numerically described by a smoothed particle hydrodynamics (SPH) approach for the diffuse (hot) gas and a sticky particle scheme for the (cool) molecular clouds. Physical processs such as star formation, stellar death or condensation and evaporation processes of clouds interacting with the ISM are described locally. An example application of the model to a star forming dwarf galaxy will be shown for comparison with other codes. Secondly we will discuss new kinds of exotic chemodynamical processes, as they occur in dense gas-star systems in galactic nuclei, such as non-standard ``drag''-force interactions, destructive and gas producing stellar collisions. Their implementation in 1D dynamical models of galactic nuclei is presented. Future prospects to generalize these to 3D are work in progress and will be discussed.

Among the astrophysical systems targeted by LISA, stars on relativistic orbits around massive black holes (MBHs) are particularly promising sources. Unfortunately, the prediction for the number and characteristics of such sources suffers from many uncertainties. Stellar dynamical Monte Carlo simulations of the evolution of galactic nucleus models allow more realistic estimates of these quantities. The computations presented here strongly suggest that the closest such extreme mass-ratio binary to be detected by LISA could be a low-mass MS star (MSS) orbiting the MBH at the center of our Milky Way. Only compact stars contribute to the expected detections from other galaxies because MSSs are disrupted by tidal forces too early.

We review possible dynamical formation processes for central massive black holes in dense star clusters. We focus on the early dynamical evolution of young clusters containing a few thousand to a few million stars. One natural formation path for a central seed black hole in these systems involves the development of the Spitzer instability, through which the most massive stars can drive the cluster to core collapse in a very short time. The sudden increase in the core density then leads to a runaway collision process and the formation of a very massive merger remnant, which must then collapse to a black hole. Alternatively, if the most massive stars end their lives before core collapse, a central cluster of stellar-mass black holes is formed. This cluster will likely evaporate before reaching the highly relativistic state necessary to drive a runaway merger process through gravitational radiation, thereby avoiding the formation of a central massive black hole. We summarize the conditions under which these different paths will be followed, and present the results of recent numerical simulations demonstrating the process of rapid core collapse and runaway collisions between massive stars.

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.

One of the main astrophysical processes leading to strong emission of gravitational waves to be detected by the future space-borne interferometer LISA is the capture of a compact star by a black hole with a mass of a few million solar masses in the center of a galaxy. In previous studies, main sequence stars were thought not to contribute because they suffer from early tidal disruption. Here we show that, according to our simulations of the stellar dynamics of the Sgr A* cluster, there must be one to a few low-mass main sequence stars sufficiently bound to the central Galactic black hole to be conspicuous sources in LISA observations. The probability that a white dwarf may be detectable is lower than 0.5 and, in spite of mass segregation, detection of a captured neutron star or stellar black hole in the center of the Milky Way is highly unlikely.

We have recently written a new code to simulate the long term evolution of spherical clusters of stars. It is based on the pioneering Monte Carlo scheme proposed by Henon in the 70's. Our code has been devised in the specific goal to treat dense galactic nuclei. After having described how we treat relaxation in a first paper, we go on and include further physical ingredients that are mostly pertinent to galactic nuclei, namely the presence of a central (growing) black hole (BH) and collisions between MS stars. Stars that venture too close to the BH are destroyed by the tidal field. This process is a channel to feed the BH and a way to produce accretion flares. Collisions between stars have often been proposed as another mechanism to drive stellar matter into the central BH. To get the best handle on the role of this process in galactic nuclei, we include it with unpreceded realism through the use of a set of more than 10000 collision simulations carried out with a SPH (Smoothed Particle Hydrodynamics) code. Stellar evolution has also been introduced in a simple way, similar to what has been done in previous dynamical simulations of galactic nuclei. To ensure that this physics is correctly simulated, we realized a variety of tests whose results are reported here. This unique code, featuring most important physical processes, allows million particle simulations, spanning a Hubble time, in a few CPU days on standard personal computers and provides a wealth of data only rivalized by N-body simulations.

The capture and gradual inspiral of stellar mass objects by a massive black hole at the centre of a galaxy has been proposed as one of the most promising source of gravitational radiation to be detected by LISA. Unfortunately rate estimates for this process suffer from many uncertainties. Here we report on the use of our newly developed Monte Carlo stellar dynamics code to tackle this problem. We present results from simple galactic nuclei models that demonstrate the high potential of our approach and point out the aspects of the problem where an improved treatment seems desirable.

We have developed a new simulation code aimed at studying the stellar dynamics of a galactic central star cluster surrounding a massive black hole. In order to include all the relevant physical ingredients (2-body relaxation, stellar mass spectrum, collisions, tidal disruption,...), we chose to revive a numerical scheme pioneered by Henon in the 70's. It is basically a Monte Carlo resolution of the Fokker-Planck equation. It can cope with any stellar mass spectrum or velocity distribution. Being a particle-based method, it also allows one to take stellar collisions into account in a very realistic way. This first paper covers the basic version of our code which treats the relaxation-driven evolution of stellar cluster without a central BH. A technical description of the code is presented, as well as the results of test computations. Thanks to the use of a binary tree to store potential and rank information and of variable time steps, cluster models with up to 2 million particles can be simulated on a standard personal computer and the CPU time required scales as N*ln(N) with the particle number N. Furthermore, the number of simulated stars needs not be equal to N but can be arbitrarily larger. A companion paper will treat further physical ingredients, mostly relevant to galactic nuclei.

Stellar collisions have long been envisioned to be of great importance in the center of galaxies where densities of 1e6 stars per cubic pc or higher are attained. Not only can they play a unique dynamical role by modifying stellar masses and orbits, but high velocity disruptive encounters occurring in the vicinity of a massive black hole can also be an occasional source of fuel for the starved central engine. In the past few years, we have been building a comprehensive table of SPH (Smoothed Particle Hydrodynamics) collision simulations for main sequence stars. This database is now integrated as a module into our Henon-like Monte Carlo code. The combination of SPH collision simulations with a Monte Carlo cluster evolution code seems ideally suited to study the frequency, characteristics and effects of stellar collisions during the long term evolution of galactic nuclei.

We have developed a new numerical code to simulate the joint evolution of a massive black hole (MBH) and a surrounding stellar cluster at the center of a galaxy. The physics treated in the simulations include: 2-body relaxation, stellar collisions (using a large set of "SPH" hydrodynamical simulations) and tidal disruption of stars by the MBH. In particular, we investigate the rates of star disruptive events that provide the MBH with gas to accrete.