Recently, more than 100 Wolf-Rayet and OB stars were identified in the Galactic center. About a third of these sources are not spatially associated with any of the known star clusters in this region. We probe the distribution of drifted sources in numerical models of the massive clusters in the Galactic center and compare it to the observed distribution of isolated massive sources in this region. We find that stars as massive as 100 Msun drift away from the center of each cluster by up to ~60 pc using the cluster models. Our best model reproduces ~60% of the known isolated massive stars out to 80 pc from the center of the Arches cluster. This number increases to 70%-80% when we only consider the region of ~20 pc from the Arches cluster.

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.

Most stars form in a cluster environment. These stars are initially surrounded by discs from which potentially planetary systems form. Of all cluster environments starburst clusters are probably the most hostile for planetary systems in our Galaxy. The intense stellar radiation and extreme density favour rapid destruction of circumstellar discs via photoevaporation and stellar encounters. Evolving a virialized model of the Arches cluster in the Galactic tidal field we investigate the effect of stellar encounters on circumstellar discs in a prototypical starburst cluster. Despite its proximity to the deep gravitational potential of the Galactic centre only a moderate fraction of members escapes to form an extended pair of tidal tails. Our simulations show that encounters destroy one third of the circumstellar discs in the cluster core within the first 2.5 Myr of evolution, preferentially affecting the least and most massive stars. A small fraction of these events causes rapid ejection and the formation of a weaker second pair of tidal tails that is overpopulated by disc-poor stars. Two predictions arise from our study: (i) If not destroyed by photoevaporation protoplanetary discs of massive late B- and early O-type stars represent the most likely hosts of planet formation in starburst clusters. (ii) Multi-epoch K- and L-band photometry of the Arches cluster would provide the kinematically selected membership sample required to detect the additional pair of disc-poor tidal tails.

We have performed a series of N-body simulations to model the Arches cluster. Our aim is to find the best fitting model for the Arches cluster by comparing our simulations with observational data and to constrain the parameters for the initial conditions of the cluster. By neglecting the Galactic potential and stellar evolution, we are able to efficiently search through a large parameter space to determine e.g. the IMF, size, and mass of the cluster. We find, that the cluster's observed present-day mass function can be well explained with an initial Salpeter IMF. The lower mass-limit of the IMF cannot be well constrained from our models. In our best models, the total mass and the virial radius of the cluster are initially (5.1 +/- 0.8) 10^4 Msun and 0.76 +/- 0.12 pc, respectively. The concentration parameter of the initial King model is w0 = 3-5.

The young star clusters we observe today are the building blocks of a new generation of stars and planets in our Galaxy and beyond. Despite their fundamental role we still lack knowledge about the conditions under which star clusters form and the impact of these often harsh environments on the evolution of their stellar and substellar members. We demonstrate the vital role numerical simulations play to uncover both key issues. Using dynamical models of different star cluster environments we show the variety of effects stellar interactions potentially have. Moreover, our significantly improved measure of mass segregation reveals that it can occur rapidly even for star clusters without substructure. This finding is a critical step to resolve the controversial debate on mass segregation in young star clusters and provides strong constraints on their initial conditions.

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.

Galaxy interactions are a common phenomenon in clusters of galaxies. Especially major mergers are of particular importance, because they can change the morphological type of galaxies. They have an impact on the mass function of galaxies and they trigger star formation - the main driver of the Galactic Matter Cycle. Therefore, we conducted a study of major mergers by means of a multi-phase ISM code. This code is based on a TREE-SPH-code combined with a sticky particle method allowing for star formation controlled by the properties of a multi-phase ISM. This is in contrast to the usually implemented Schmidt law depending mainly on the gas density. Previously, this code was used on isolated galaxies. Since our star formation recipe is not restricted to a special type of galaxy, it is interesting to apply it to interacting galaxies, too. Our study on major mergers includes a research of global properties of the interacting system, namely the star formation rate and the star formation efficiency, the evaporation and condensation rates, as well as the mass exchange of distinct components, namely stars, diffuse ISM, and clouds. Investigating these properties provides insight to interrelations between various physical processes. The results indicate that the star formation efficiency as well as the evaporation and condensation rates are influenced by the interaction.

We present the living application, a method to autonomously manage applications on the grid. During its execution on the grid, the living application makes choices on the resources to use in order to complete its tasks. These choices can be based on the internal state, or on autonomously acquired knowledge from external sensors. By giving limited user capabilities to a living application, the living application is able to port itself from one resource topology to another. The application performs these actions at run-time without depending on users or external workflow tools. We demonstrate this new concept in a special case of a living application: the living simulation. Today, many simulations require a wide range of numerical solvers and run most efficiently if specialized nodes are matched to the solvers. The idea of the living simulation is that it decides itself which grid machines to use based on the numerical solver currently in use. In this paper we apply the living simulation to modelling the collision between two galaxies in a test setup with two specialized computers. This simulation switces at run-time between a GPU-enabled computer in the Netherlands and a GRAPE-enabled machine that resides in the United States, using an oct-tree N-body code whenever it runs in the Netherlands and a direct N-body solver in the United States.

We present Sapporo, a library for performing high-precision gravitational N-body simulations on NVIDIA Graphical Processing Units (GPUs). Our library mimics the GRAPE-6 library, and N-body codes currently running on GRAPE-6 can switch to Sapporo by a simple relinking of the library. The precision of our library is comparable to that of GRAPE-6, even though internally the GPU hardware is limited to single precision arithmetics. This limitation is effectively overcome by emulating double precision for calculating the distance between particles. The performance loss of this operation is small (< 20%) compared to the advantage of being able to run at high precision. We tested the library using several GRAPE-6-enabled N-body codes, in particular with Starlab and phiGRAPE. We measured peak performance of 800 Gflop/s for running with 10^6 particles on a PC with four commercial G92 architecture GPUs (two GeForce 9800GX2). As a production test, we simulated a 32k Plummer model with equal mass stars well beyond core collapse. The simulation took 41 days, during which the mean performance was 113 Gflop/s. The GPU did not show any problems from running in a production environment for such an extended period of time.

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

We describe a novel N-body code designed for simulations of the central regions of galaxies containing massive black holes. The code incorporates Mikkola's 'algorithmic' chain regularization scheme including post-Newtonian terms up to PN2.5 order. Stars moving beyond the chain are advanced using a fourth-order integrator with forces computed on a GRAPE board. Performance tests confirm that the hybrid code achieves better energy conservation, in less elapsed time, than the standard scheme and that it reproduces the orbits of stars tightly bound to the black hole with high precision. The hybrid code is applied to two sample problems: the effect of finite-N gravitational fluctuations on the orbits of the S-stars; and inspiral of an intermediate-mass black hole into the galactic center.

We show that the presence of a rho~1/r^{3/2} dark matter overdensity can be robustly predicted at the center of any galaxy old enough to have grown a power-law density cusp in the stars via the Bahcall-Wolf mechanism. Using both Fokker-Planck and direct N-body integrations, we demonstrate collisional generation of these dark matter "crests" (Collisionally REgenerated STtructures) even in the extreme case that the density of both stars and dark matter were previously lowered by slingshot ejection from a binary supermassive black hole. The time scale for collisional growth of the crest is approximately the two-body relaxation time as defined by the stars, which is < 10 Gyr at the centers of stellar spheroids with luminosities comparable to that of the Milky Way bulge or less. The presence of crests can robustly be predicted in such galaxies, unlike the steeper enhancements, called "spikes," produced by the adiabatic growth of black holes. We discuss special cases where the prospects for detecting dark matter annihilations from the centers of galaxy haloes are significantly affected by the formation of crests.

Direct-summation N-body algorithms compute the gravitational interaction between stars in an exact way and have a computational complexity of O(N^2). Performance can be greatly enhanced via the use of special-purpose accelerator boards like the GRAPE-6A. However the memory of the GRAPE boards is limited. Here, we present a performance analysis of direct N-body codes on two parallel supercomputers that incorporate special-purpose boards, allowing as many as four million particles to be integrated. Both computers employ high-speed, Infiniband interconnects to minimize communication overhead, which can otherwise become significant due to the small number of "active" particles at each time step. We find that the computation time scales well with processor number; for 2*10^6 particles, efficiencies greater than 50% and speeds in excess of 2 TFlops are reached.

We present a new particle code for modelling the evolution of galaxies. The code is based on a multi-phase description for the interstellar medium (ISM). We included star formation (SF), stellar feedback by massive stars and planetary nebulae, phase transitions and interactions between gas clouds and ambient diffuse gas, namely condensation, evaporation, drag and energy dissipation. The latter is realised by radiative cooling and inelastic cloud-cloud collisions. We present new schemes for SF and stellar feedback. They include a consistent calculation of the star formation efficiency (SFE) based on ISM properties as well as a detailed redistribution of the feedback energy into the different ISM phases. As a first test example we show a model of the evolution of a present day Milky-Way-type galaxy. Though the model exhibits a quasi-stationary behaviour in global properties like mass fractions or surface densities, the evolution of the ISM is locally strongly variable depending on the local SF and stellar feedback. We start only with two distinct phases, but a three-phase ISM is formed soon consisting of cold molecular clouds, a warm gas disk and a hot gaseous halo. Hot gas is also found in bubbles in the disk accompanied by type II supernovae explosions. The star formation rate (SFR) is ~1.6 M_sun/year on average decreasing slowly with time due to gas consumption. In order to maintain a constant SFR gas replenishment, e.g. by infall, of the order 1 M_sun/year is required. Our model is in fair agreement with Kennicutt's (1998) SF law including the cut-off at ~10 M_sun/pc^2. Models with a constant SFE, i.e. no feedback on the SF, fail to reproduce Kennicutt's law.

We present a modified TREESPH code to model galaxies in 3d. The model includes a multi-phase description of the interstellar medium which combines two numerical techniques. A diffuse warm/hot gas phase is modelled by SPH while a sticky particle scheme is used to represent a cloudy medium. Interaction processes, such as star formation and feedback, cooling and mixing by condensation and evaporation, are taken into account. Here we apply our model to the evolution of a Milky Way type galaxy. After an initial stage, a quasi-equilibrium state is reached. It is characterised by a star formation rate of ~1 M_sun/year. Condensation and evaporation rates are in balance at 0.1-1 M_sun/year.

Chemo-dynamical models have been introduced in the late eighties and are a generally accepted tool for understanding galaxy evolution. They have been successfully applied to one-dimensional problems, e.g. the evolution of non-rotating galaxies, and two-dimensional problems, e.g. the evolution of disk galaxies. Recently, also three-dimensional chemo-dynamical models have become available. In these models the dynamics of different components, i.e. dark matter, stars and a multi-phase interstellar medium, are treated in a self-consistent way and several processes allow for an exchange of matter, energy and momentum between the components or different gas phases. Some results of chemo-dynamical models and their comparison with observations of chemical abundances or star formation histories will be reviewed.

We present a 3d code for the dynamical evolution of a multi-phase interstellar medium (ISM) coupled to stars via star formation (SF) and feedback processes. The multi-phase ISM consists of clouds (sticky particles) and diffuse gas (SPH): exchange of matter, energy and momentum is achieved by drag (due to ram pressure) and condensation or evaporation processes. The cycle of matter is completed by SF and feedback by SNe and PNe. A SF scheme based on a variable SF efficiency as proposed by Elmegreen & Efremov (1997) is presented. For a Milky Way type galaxy we get a SF rate of ~1 M_sun/yr with an average SF efficiency of ~5%.

We present a new particle based code with a multi-phase description of the ISM implemented in order to follow the chemo-dynamical evolution of galaxies. The multi-phase ISM consists of clouds (sticky particles) and diffuse gas (SPH): Exchange of matter, energy and momentum is achieved by drag (due to ram pressure) and condensation or evaporation. Based on time scales we show that in Milky-Way-like galaxies the drag force is for molecular clouds only important, if their relative velocities exceed 100 km/s. For the mass exchange we find that clouds evaporate only if the temperature of the ambient gas is higher than one million Kelvin. At lower temperatures condensation takes place at time scales of the order of 1-10 Gyr.

Modeling of interacting galaxies suffers from an extended parameter space prohibiting traditional grid based search strategies. As an alternative approach a combination of a Genetic Algorithm (GA) with fast restricted N-body simulations can be applied. A typical fit takes about 3-6 CPU-hours on a PentiumII processor. Here we present a parallel implementation of our GA which reduces the CPU-requirement of a parameter determination to a few minutes on 100 nodes of a CRAY T3E.

We present optical data which shows that G", the optical counterpart of the gamma-ray pulsar Geminga, pulses in B with a period of 0.237 seconds. The similarity between the optical pulse shape and the gamma-ray light curve indicates that a large fraction of the optical emission is non-thermal in origin - contrary to recent suggestions based upon the total optical flux. The derived magnitude of the pulsed emission is m_B = 26.0 +/- 0.4. Whilst it is not possible to give an accurate figure for the pulsed fraction (due to variations in the sky background) we can give an upper limit of m_B ~ 27 for the unpulsed fraction.