Many stars do not live alone, but instead have one or more stellar companions. Observations show that these binaries, triples and higher-order multiples are common. Whereas the evolution of single stars and binaries have been studied extensively, the same is not true for the evolution of stellar triples. To fill this gap in our general understanding of stellar lives, we aim to systematically explore the long-term evolution of triples and to map out the most common evolutionary pathways that triples go through. We quantitatively study how triples evolve, which processes are most relevant, and how this differs from binary evolution. We simulate the evolution of several large populations of triples with a population synthesis approach. We make use of the triple evolution code TRES to simulate the evolution of each triple in a consistent way; including three-body dynamics (based on the secular approach), stellar evolution and their mutual influences. We simulate the evolution of the system up until mass transfer starts, the system becomes dynamically unstable, or a Hubble time has passed. We find that stellar interactions are common in triples. Compared to a binary population, we find that the fraction of systems that can undergo mass transfer is about 2 to 3 times larger in triples. Moreover, whereas in binaries the orbits typically reach circularisation before Roche-lobe overflow, this is not true anymore in triples. In our simulations, about 40% of systems retain an eccentric orbit. Additionally, we discuss various channels of triple evolution in detail such as those where the secondary or the tertiary is the first star to initiate a mass transfer event.

It has been recently discovered that the transition of X-ray pulsars to the low luminosity state ($L\lesssim 10^{35}\,{\rm erg\,s^{-1}}$) is accompanied by a dramatic spectral changes. Namely, the typical power-law-like spectrum with high energy cutoff transforms into a two-component structure with a possible cyclotron absorption feature on top of it. It was proposed that these spectral characteristics can be explained qualitatively by the emission of cyclotron photons in the atmosphere of the neutron star caused by collisional excitation of electrons to upper Landau levels and further comptonization of the photons by electron gas. The latter is expected to be overheated in a thin top layer of the atmosphere. In this paper, we perform Monte Carlo simulations of the radiative transfer in the atmosphere of an accreting neutron star while accounting for a resonant scattering of polarized X-ray photons by thermally distributed electrons. The spectral shape is shown to be strongly polarization-dependent in soft X-rays ($\lesssim 10\,{\rm keV}$) and near the cyclotron scattering feature. The results of our numerical simulations are tested against the observational data of X-ray pulsar A 0535+262 in the low luminosity state. We show that the spectral shape of the pulsar can be reproduced by the proposed theoretical model. The applications of the discovery to the observational studies of accreting neutron stars are discussed.

In this paper, we consider triple systems composed of main-sequence (MS) stars, and their internal evolution due to stellar and binary evolution. Our focus is on triples that produce white dwarfs (WDs), where Roche lobe overflow of an evolving tertiary triggers accretion onto the inner binary via a circumbinary disk (CBD) driving it toward a mass ratio of unity. We present a combination of analytic- and population synthesis-based calculations performed using the \texttt{SeBa} code to constrain the expected frequency of such systems, given a realistic initial population of MS triples, and provide the predicted distributions of orbital periods. We identify the parameter space for triples that can accommodate a CBD, to inform future numerical simulations of suitable initial conditions. We find that $\lesssim$ 10\% of all MS triples should be able to accommodate a CBD around the inner binary, and compute lower limits for the production rates. This scenario broadly predicts mergers of near equal-mass binaries, producing blue stragglers (BSs), Type Ia supernovae, gamma ray bursts and gravitational wave-induced mergers, along with the presence of an outer WD tertiary companion. We compare our predicted distributions to a sample of field BS binaries, and argue that our proposed mechanism explains the observed range of orbital periods. Finally, the mechanism considered here could produce hypervelocity MS stars, WDs and even millisecond pulsars with masses close to the Chandrasekhar mass limit, and be used to constrain the maximum remnant masses at the time of any supernova explosion.

Circumstellar discs are the precursors of planetary systems and develop shortly after their host star has formed. In their early stages these discs are immersed in an environment rich in gas and neighbouring stars, which can be hostile for their survival. There are several environmental processes that can affect the evolution of circumstellar discs, and external photoevaporation is arguably one of the most important ones. Theoretical and observational evidence point to circumstellar discs losing mass quickly when in the vicinity of massive, bright stars. External photoevaporation can then constrain the time and material available to form planets. The stellar density of the region seems to influence the extent of the effects of external photoevaporation. In this work we perform simulations of star clusters with a range of stellar densities. Stars with masses $\mathrm{M}_* \leq 1.9$ \MSun are surrounded by a disc, and stars with masses $\mathrm{M}_* > 1.9$ \MSun are considered as emitting radiation. Our results show that external photoevaporation is efficient in depleting disc masses and that the degree of its effect is related to stellar density. Dense clusters have $\sim10\%$ surviving discs by 2 Myr, whereas sparse ones have $\sim50\%$. Surviving discs in sparser regions also span a larger range of masses. We compare our results to observations of the Lupus clouds, the Orion Nebula Cluster, the Orion Molecular Cloud-2, Taurus, and NGC 2024, and find that the trends observed between region density and disc masses are similar to those in our simulations.

Gaia Data Release 2 revealed detailed structures of nearby stars in phase space. These include the Hercules stream, whose origin is still debated. Earlier numerical studies conjectured that the observed structures originate from orbits in resonance with the bar, based on static potential models for the Milky Way. We, in contrast, approach the problem via a self-consistent, dynamic, and morphologically well-resolved model, namely a full $N$-body simulation of the Milky Way. Our simulation comprises about 5.1 billion particles in the galactic stellar bulge, bar, disk, and dark-matter halo and is evolved to 10 Gyr. Our model's disk component is composed of 200 million particles, and its simulation snapshots are stored every 10 Myr, enabling us to resolve and classify resonant orbits of representative samples of stars. After identifying the Sun's position in the simulation, we compare the distribution of stars in its neighborhood with Gaia's astrometric data, thereby establishing the role of identified resonantly trapped stars in the formation of Hercules-like structures. From our orbital spectral-analysis we identify multiple resonances and conclude that the Hercules stream is dominated by the 4:1 and 5:1 outer Lindblad and corotation resonances. In total, this yields a trimodal structure of the Hercules stream. From the relation between resonances and ridges in phase space, we conclude that the pattern speed of the Milky-Way bar is 40$-$45 km s$^{-1}$ kpc$^{-1}$.

We make use of APOGEE and $Gaia$ data to identify stars that are consistent with being born in the same association or star cluster as the Sun. We limit our analysis to stars that match solar abundances within their uncertainties, as they could have formed from the same Giant Molecular Cloud (GMC) as the Sun. We constrain the range of orbital actions that solar siblings can have with a suite of simulations of solar birth clusters evolved in static and time-dependent tidal fields. The static components of each galaxy model are the bulge, disk, and halo, while the various time-dependent components include a bar, spiral arms, and GMCs. In galaxy models without GMCs, simulated solar siblings all have $J_R < 122$ km $\rm s^{-1}$ kpc, $990 < L_z < 1986$ km $\rm s^{-1}$ kpc, and $0.15 < J_z < 0.58$ km $\rm s^{-1}$ kpc. Given the actions of stars in APOGEE and $Gaia$, we find 104 stars that fall within this range. One candidate in particular, Solar Sibling 1, has both chemistry and actions similar enough to the solar values that strong interactions with the bar or spiral arms are not required for it to be dynamically associated with the Sun. Adding GMCs to the potential can eject solar siblings out of the plane of the disk and increase their $J_z$, resulting in a final candidate list of 296 stars. The entire suite of simulations indicate that solar siblings should have $J_R < 122$ km $\rm s^{-1}$ kpc, $353 < L_z < 2110$ km $\rm s^{-1}$ kpc, and $J_z < 0.8$ km $\rm s^{-1}$ kpc. Given these criteria, it is most likely that the association or cluster that the Sun was born in has reached dissolution and is not the commonly cited open cluster M67.

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.

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 ...

Escalating observations of exo-minor planets and their destroyed remnants both passing through the solar system and within white dwarf planetary systems motivate an understanding of the orbital history and fate of exo-Kuiper belts and planetesimal discs. Here we explore how the structure of a 40-1000 au annulus of planetesimals orbiting inside of a solar system analogue that is itself initially embedded within a stellar cluster environment varies as the star evolves through all of its stellar phases. We attempt this computationally challenging link in four parts: (1) by performing stellar cluster simulations lasting 100 Myr, (2) by making assumptions about the subsequent quiescent 11 Gyr main-sequence evolution, (3) by performing simulations throughout the giant branch phases of evolution, and (4) by making assumptions about the belt's evolution during the white dwarf phase. Throughout these stages, we estimate the planetesimals' gravitational responses to analogues of the four solar system giant planets, as well as to collisional grinding, Galactic tides, stellar flybys, and stellar radiation. We find that the imprint of stellar cluster dynamics on the architecture of $\gtrsim 100$ km-sized exo-Kuiper belt planetesimals is retained throughout all phases of stellar evolution unless violent gravitational instabilities are triggered either (1) amongst the giant planets, or (2) due to a close ($\ll 10^3$ au) stellar flyby. In the absence of these instabilities, these minor planets simply double their semimajor axis while retaining their primordial post-cluster eccentricity and inclination distributions, with implications for the free-floating planetesimal population and metal-polluted white dwarfs.

We present the smoothed-particle hydrodynamics simulation code, Bonsai-SPH, which is a continuation of our previously developed gravity-only hierarchical $N$-body code (called Bonsai). The code is optimized for Graphics Processing Unit (GPU) accelerators which enables researchers to take advantage of these powerful computational resources. Bonsa-SPH produces simulation results comparable with state-of-the-art, CPU based, codes, but using an order of magnitude less computation time. The code is freely available online and the details are described in this work.

Chaos is present in most stellar dynamical systems and manifests itself through the exponential growth of small perturbations. Exponential divergence drives time irreversibility and increases the entropy in the system. A numerical consequence is that integrations of the N-body problem unavoidably magnify truncation and rounding errors to macroscopic scales. Hitherto, a quantitative relation between chaos in stellar dynamical systems and the level of irreversibility remained undetermined. In this work we study chaotic three-body systems in free fall initially using the accurate and precise N-body code Brutus, which goes beyond standard double-precision arithmetic. We demonstrate that the fraction of irreversible solutions decreases as a power law with numerical accuracy. This can be derived from the distribution of amplification factors of small initial perturbations. Applying this result to systems consisting of three massive black holes with zero total angular momentum, we conclude that up to five percent of such triples would require an accuracy of smaller than the Planck length in order to produce a time-reversible solution, thus rendering them fundamentally unpredictable.

By means of a fully connected artificial neural network, we identified asteroids with the potential to impact Earth. The resulting instrument, named the Hazardous Object Identifier (HOI), was trained on the basis of an artificial set of known impactors which were generated by launching objects from Earth's surface and integrating them backward in time. HOI was able to identify 95.25% of the known impactors simulated that were present in the test set as potential impactors. In addition, HOI was able to identify 90.99% of the potentially hazardous objects identified by NASA, without being trained on them directly.

Planet-forming circumstellar disks are a fundamental part of the star formation process. Since stars form in a hierarchical fashion in groups of up to hundreds or thousands, the UV radiation environment that these disks are exposed to can vary in strength by at least six orders of magnitude. This radiation can limit the masses and sizes of the disks. Diversity in star forming environments can have long lasting effects in disk evolution and in the resulting planetary populations. We perform simulations to explore the evolution of circumstellar disks in young star clusters. We include viscous evolution, as well as the impact of dynamical encounters and external photoevaporation. We find that photoevaporation is an important process in destroying circumstellar disks: in regions of stellar density $\rho \sim 100 \mathrm{\ M}_\odot \mathrm{\ pc}^{-3}\mathrm{\ }$ around 80% of disks are destroyed before 2 Myr of cluster evolution. Our findings are in agreement with observed disk fractions in young star forming regions and support previous estimations that planet formation must start in timescales < 0.1 - 1 Myr.

Since its formulation by Sir Isaac Newton, the problem of solving the equations of motion for three bodies under their own gravitational force has remained practically unsolved. Currently, the solution for a given initialization can only be found by performing laborious iterative calculations that have unpredictable and potentially infinite computational cost, due to the system's chaotic nature. We show that an ensemble of solutions obtained using an arbitrarily precise numerical integrator can be used to train a deep artificial neural network (ANN) that, over a bounded time interval, provides accurate solutions at fixed computational cost and up to 100 million times faster than a state-of-the-art solver. Our results provide evidence that, for computationally challenging regions of phase-space, a trained ANN can replace existing numerical solvers, enabling fast and scalable simulations of many-body systems to shed light on outstanding phenomena such as the formation of black-hole binary systems or the origin of the core collapse in dense star clusters.

As of August 2019, among the more than 4000 confirmed exoplanets, only one has been detected in a globular cluster (GC) M4. The scarce of exoplanet detections motivates us to employ direct $N$-body simulations to investigate the dynamical stability of planets in young massive clusters (YMCs), which are potentially the progenitors of GCs. In an $N=128{\rm k}$ cluster of virial radius 1.7 pc (comparable to Westerlund-1), our simulations show that most wide-orbit planets ($a\geq 20$~au) will be ejected within a timescale of 10 Myr. Interestingly, more than $70\%$ of planets with $a<5$~au survive in the 100 Myr simulations. Ignoring planet-planet scattering and tidal damping, the survivability at $t$ Myr as a function of initial semi-major axis $a_0$ in au in such a YMC can be described as $f_{\rm surv}(a_0, t)=-0.33 \log_{10}(a_0) \left(1 - e^{-0.0482t} \right) + 1$. Upon ejection, about $28.8\%$ of free-floating planets (FFPs) have sufficient speeds to escape from the host cluster at a crossing timescale. The other FFPs will remain bound to the cluster potential, but the subsequent dynamical evolution of the stellar system can result in the delayed ejection of FFPs from the host cluster. Although a full investigation of planet population in GCs requires extending the simulations to multi-Gyr, our results suggest that wide-orbit planets and free-floating planets are unlikely to be found in GCs.

Comets in the Oort cloud evolve under the influence of internal and external perturbations, such as giant planets, stellar passages, and the galactic tidal field. We aim to study the dynamical evolution of the comets in the Oort cloud, accounting for external perturbations (passing stars and the galactic tide). We first construct an analytical model of stellar encounters. We find that individual perturbations do not modify the dynamics of the comets in the cloud unless very close (< 0.5pc) encounters occur. Using proper motions, parallaxes, and radial velocities from Gaia DR2, we construct an astrometric catalogue of 14,659 stars that are within 50pc from the Sun. For all these stars we calculate the time and the closest distance to the Sun. We find that the cumulative effect of relatively distant ($\leq1$ pc) passing stars can perturb the comets in the Oort cloud. Finally, we study the dynamical evolution of the comets in the Oort cloud under the influence of multiple stellar encounters within 2.5pc from the Sun and the galactic tidal field over $\pm10$Myr. We considered two models for the Oort cloud, compact (a $\leq$0.25 pc) and extended (a$ \leq0.5$ pc). We find that the cumulative effect of stellar encounters is the major perturber of the Oort cloud for a compact configuration while for the extended, the galactic tide is the major perturber. In both cases, the effect of passing stars and the galactic tide raises the semi-major axis of $\sim1.1$\% of the comets at the edge of the cloud up to interstellar regions ($a >0.5$pc). This leads to the creation of transitional interstellar comets, which might become interstellar objects due to external perturbations. This raises the question about the existence of a cloud of objects in the interstellar space which might overlap with our Oort cloud if we consider that other planetary systems face similar processes for the ejection of comets.

Many white dwarf stars show signs of having accreted smaller bodies, implying that they may host planetary systems. A small number of these systems contain gaseous debris discs, visible through emission lines. We report a stable 123.4min periodic variation in the strength and shape of the CaII emission line profiles originating from the debris disc around the white dwarf SDSSJ122859.93+104032.9. We interpret this short-period signal as the signature of a solid body held together by its internal strength.

We present stellar_wind.py, a module that provides multiple methods of simulating stellar winds using smoothed particle hydrodynamics codes (SPH) within the astrophysical multipurpose software environment (AMUSE) framework. With the simple wind mode, we create SPH wind particles in a spherically symmetric shell. We inject the wind particles with a velocity equal to their terminal velocity. The accelerating wind mode is similar, but with this method particles can be injected with a lower initial velocity than the terminal velocity and they are accelerated away from the star according to an acceleration function. With the heating wind mode, SPH particles are created with zero initial velocity with respect to the star, but instead wind particles are given an internal energy based on the integrated mechanical luminosity of the star. This mode is designed to be used on longer timescales and larger spatial scales compared to the other two modes and assumes that the star is embedded in a gas cloud. For fast winds, we find that both the simple and accelerating mode can reproduce the desired velocity, density and temperature profiles. For slow winds, the simple wind mode is insufficient due to dominant hydrodynamical effects that change the wind velocities. The accelerating mode, with additional options to account for these hydrodynamical effects, can still reproduce the desired wind profiles. We test the heating mode by simulating both a normal wind and a supernova explosion of a single star in a uniform density medium. The stellar wind simulation results matches the analytical solution for an expanding wind bubble. The supernova simulation gives qualitatively correct results, but the simulated bubble expands faster than the analytical solution predicts. We conclude with an example of a triple star system which includes the colliding winds of all three stars.

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.

Gamma-ray astronomy explores the most energetic photons in nature to address some of the most pressing puzzles in contemporary astrophysics. It encompasses a wide range of objects and phenomena: stars, supernovae, novae, neutron stars, stellar-mass black holes, nucleosynthesis, the interstellar medium, cosmic rays and relativistic-particle acceleration, and the evolution of galaxies. MeV gamma-rays provide a unique probe of nuclear processes in astronomy, directly measuring radioactive decay, nuclear de-excitation, and positron annihilation. The substantial information carried by gamma-ray photons allows us to see deeper into these objects, the bulk of the power is often emitted at gamma-ray energies, and radioactivity provides a natural physical clock that adds unique information. New science will be driven by time-domain population studies at gamma-ray energies. This science is enabled by next-generation gamma-ray instruments with one to two orders of magnitude better sensitivity, larger sky coverage, and faster cadence than all previous gamma-ray instruments. This transformative capability permits: (a) the accurate identification of the gamma-ray emitting objects and correlations with observations taken at other wavelengths and with other messengers; (b) construction of new gamma-ray maps of the Milky Way and other nearby galaxies where extended regions are distinguished from point sources; and (c) considerable serendipitous science of scarce events -- nearby neutron star mergers, for example. Advances in technology push the performance of new gamma-ray instruments to address a wide set of astrophysical questions.

We propose a formation mechanism for twin blue stragglers (BSs) in compact binaries that involves mass transfer from an evolved outer tertiary companion on to the inner binary via a circumbinary disk. We apply this scenario to the observed double BS system Binary 7782 in the old open cluster NGC 188, and show that its observed properties are naturally reproduced within the context of the proposed model. We predict the following properties for twin BSs: (1) For the outer tertiary orbit, the initial orbital period should lie between 220 days $\lesssim$ P$_{\rm out}$ $\lesssim$ 1100 days, assuming initial masses for the inner binary components of $m_{\rm 1} = 1.1$ M$_{\odot}$ and $m_{\rm 2} =$ 0.9 M$_{\odot}$ and an outer tertiary mass of $m_{\rm 3} = 1.4$ M$_{\odot}$. After Roche-lobe overflow, the outer star turns into a white dwarf (WD) of mass 0.43 to 0.54\,\MSun. There is a correlation between the mass of this WD and the outer orbital period: more massive WDs will be on wider orbits. (3) The rotational axes of both BSs will be aligned with each other and the orbital plane of the outer tertiary WD. (4) The BSs will have roughly equal masses, independent of their initial masses (since the lower mass star accretes the most). The dominant accretor should, therefore, be more enriched by the accreted material. Hence, one of the BSs will appear to be more enriched by either He, C and O or by s-process elements, if the donor started Roche lobe overflow on, respectively, the red giant or asymptotic giant branch. (5) Relative to old clusters, twin BSs in close binaries formed from the proposed mechanism should be more frequent in the Galactic field and open clusters with ages $\lesssim$ 4-6 Gyr, since then the donor will have a radiative envelope. (6) The orbit of the binary BS will have a small semi-major axis (typically $\aplt 0.3$\,au) and be close to circular ($e \aplt 0.2$).

Type II supernova progenitors are expected to emit copious amounts of mass in a dense stellar wind prior to the explosion. When the progenitor is a member of a binary, the orbital motion modulates the density of this wind. When the progenitor explodes, the high-velocity ejecta collides with the modulated wind, which in turn produces a modulated radio signal. In this work we derive general analytic relations between the parameters of the radio signal modulations and binary parameter in the limit of large member mass ratio. We use these relations to infer the semi major axis of SN1979c and a lower bound for the mass of the companion. We further constrain the analytic estimates by numerical simulations using the AMUSE framework. In these calculations we simulate the progenitor binary system including the wind and the gravitational effect of a companion star. The simulation output is compared to the observed radio signal in supernova SN1979C. We find that it must have been a binary with an orbital period of about 2000 year. If the exploding star evolved from a $\sim 18 M_{\odot}$ zero-age main-sequence at solar metalicity, we derive a companion mass of $5$ to $12 M_{\odot}$ in an orbit with an eccentricity lower than about 0.8.

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 galactic halo likely grew over time in part by assembling smaller galaxies, the so-called building blocks. We investigate if the properties of these building blocks are reflected in the halo white dwarf (WD) population in the Solar neighborhood. Furthermore, we compute the halo WD luminosity functions (WDLFs) for four major building blocks of five cosmologically motivated stellar haloes. We couple the SeBa binary population synthesis model to the Munich-Groningen semi-analytic galaxy formation model, applied to the high-resolution Aquarius dark matter simulations. Although the semi-analytic model assumes an instantaneous recycling approximation, we model the evolution of zero-age main sequence stars to WDs, taking age and metallicity variations of the population into account. Although the majority of halo stars is old and metal-poor and therefore the WDs in the different building blocks have similar properties (including present-day luminosity), we find in our models that the WDs originating from building blocks that have young and/or metal-rich stars can be distinguished from WDs that were born in other building blocks. In practice however, it will be hard to prove that these WDs really originate from different building blocks, as the variations in the halo WD population due to binary WD mergers result in similar effects. The five joined stellar halo WD populations that we modelled result in WDLFs that are very similar to each other. We find that simple models with a Kroupa or Salpeter initial mass function (IMF) fit the observed luminosity function slightly better, since the Chabrier IMF is more top-heavy, although this result is dependent on our choice of the stellar halo mass density in the Solar neighborhood.

The (yet-to-be confirmed) discovery of a Neptune-sized moon around the ~3.2 Jupiter-mass planet in Kepler 1625 puts interesting constraints on the formation of the system. In particular, the relatively wide orbit of the moon around the planet, at ~40 planetary radii, is hard to reconcile with planet formation theories. We demonstrate that the observed characteristics of the system can be explained from the tidal capture of a secondary planet in the young system. After a quick phase of tidal circularization, the lunar orbit, initially much tighter than 40 planetary radii, subsequently gradually widened due to tidal synchronization of the spin of the planet with the orbit, resulting in a synchronous planet-moon system. Interestingly, in our scenario the captured object was originally a Neptune-like planet, turned into a moon by its capture.

Stars with circumstellar disks may form in environments with high stellar and gas densities which affects the disks through processes like truncation from dynamical encounters, ram pressure stripping, and external photoevaporation. Circumstellar disks also undergo viscous evolution which leads to disk expansion. Previous work indicates that dynamical truncation and viscous evolution play a major role in determining circumstellar disk size and mass distributions. However, it remains unclear under what circumstances each of these two processes dominates. Here we present results of simulations of young stellar clusters taking viscous evolution and dynamical truncations into account. We model the embedded phase of the clusters by adding leftover gas as a background potential which can be present through the whole evolution of the cluster, or expelled after 1 Myr. We compare our simulation results to actual observations of disk sizes, disk masses, and accretion rates in star forming regions. We argue that the relative importance of dynamical truncations and the viscous evolution of the disks changes with time and cluster density. Viscous evolution causes the importance of dynamical encounters to increase in time, but the encounters cease soon after the expulsion of the leftover gas. For the clusters simulated in this work, viscous growth dominates the evolution of the disks.

We analyze the position of the two populations of blue stragglers in the globular cluster M30 in the Hertzsprung-Russell diagram. Both populations of blue stragglers are brighter than the cluster's turn-off, but one population (the blue blue-stragglers) align along the zero-age main-sequence whereas the (red) population is elevated in brightness (or colour) by $\sim 0.75$ mag. Based on stellar evolution and merger simulations we argue that the red population, which composes about 40\% of the blue stragglers in M 30, is formed at a constant rate of $\sim 2.8$ blue stragglers per Gyr over the last $\sim 10$ Gyr. The blue population is formed in a burst that started $\sim 3.2$ Gyr ago at a peak rate of $30$ blue stragglers per Gyr$^{-1}$ with an e-folding time scale of $0.93$ Gyr. We speculate that the burst resulted from the core collapse of the cluster at an age of about 9.8 Gyr, whereas the constantly formed population is the result of mass transfer and mergers through binary evolution. In that case about half the binaries in the cluster effectively result in a blue straggler.

The Solar system was once rich in the short-lived radionuclide (SLR) $^{26}$Al\, but deprived in $^{60}$Fe. Several models have been proposed to explain these anomalous abundances in SLRs, but none has been set within a self-consistent framework of the evolution of the Solar system and its birth environment. The anomalous abundance in $^{26}$Al may have originated from the accreted material in the wind of a massive $\apgt 20$\,$M_\odot$ Wolf-Rayet star, but the star could also have been a member of the parental star-cluster instead of an interloper or an older generation that enriched the proto-solar nebula. The protoplanetary disk at that time was already truncated around the Kuiper-cliff (at $45$ au) by encounters with another cluster members before it was enriched by the wind of the nearby Wolf-Rayet star. The supernova explosion of a nearby star, possibly but not necessarily the exploding Wolf-Rayet star, heated the disk to $\apgt 1500$K, melting small dust grains and causing the encapsulation and preservation of $^{26}$Al into vitreous droplets. This supernova, and possibly several others, caused a further abrasion of the disk and led to its observed tilt of $5.6\pm1.2^\circ$ with respect to the Sun's equatorial plane. The abundance of $^{60}$Fe originates from a supernova shell, but its preservation results from a subsequent supernova. At least two supernovae are needed (one to deliver $^{60}$Fe\, and one to preserve it in the disk) to explain the observed characteristics of the Solar system. The most probable birth cluster then has $N = 2500\pm300$ stars and a radius of $r_{\rm vir} = 0.75\pm0.25$ pc. We conclude that Solar systems equivalent systems form in the Milky Way Galaxy at a rate of about 30 per Myr, in which case approximately 36,000 Solar system analogues roam the Milky Way.

We construct a model for the Milky Way Galaxy composed of a stellar disc and bulge embedded in a dark-matter halo. All components are modelled as $N$-body systems with up to 8 billion equal-mass particles and integrated up to an age of 10\,Gyr. We find that net angular-momentum of the dark-matter halo with a spin parameter of $\lambda=0.06$ is required to form a relatively short bar ($\sim 4$\,kpc) with a high pattern speed (40--50\,km\,s$^{-1}$). By comparing our model with observations of the Milky Way Galaxy, we conclude that a disc mass of $\sim 3.7\times10^{10}M_{\odot}$ and an initial bulge scale length and velocity of $\sim 1$\,kpc and $\sim 300$\,km\,s$^{-1}$, respectively, fit best to the observations. The disc-to-total mass fraction ($f_{\rm d}$) appears to be an important parameter for the evolution of the Galaxy and models with $f_{\rm d}\sim 0.45$ are most similar to the Milky Way Galaxy. In addition, we compare the velocity distribution in the solar neighbourhood in our simulations with observations in the Milky Way Galaxy. In our simulations the observed gap in the velocity distribution, which is expected to be caused by the outer Lindblad resonance (the so-called Hercules stream), appears to be a time-dependent structure. The velocity distribution changes on a time scale of 20--30\,Myr and therefore it is difficult to estimate the pattern speed of the bar from the shape of the local velocity distribution alone.

Scientific discovery is mediated by ideas that, after being formulated in hypotheses, can be tested, validated, and quantified before they eventually lead to accepted concepts. Computer-mediated discovery in astrophysics is no exception, but antiquated code that is only intelligible to scientists who were involved in writing it is holding up scientific discovery in the field. A bold initiative is needed to modernize astrophysics code and make it transparent and useful beyond a small group of scientists. (abridged)

WD 1145+017 is currently the only white dwarf known to exhibit periodic transits of planetary debris as well as absorption lines from circumstellar gas. We present the first simultaneous fast optical spectrophotometry and broad-band photometry of the system, obtained with the Gran Telescopio Canarias (GTC) and the Liverpool Telescope (LT), respectively. The observations spanned $5.5$ h, somewhat longer than the $4.5$-h orbital period of the debris. Dividing the GTC spectrophotometry into five wavelength bands reveals no significant colour differences, confirming grey transits in the optical. We argue that absorption by an optically thick structure is a plausible alternative explanation for the achromatic nature of the transits that can allow the presence of small-sized ($\sim\mu$m) particles. The longest ($87$ min) and deepest ($50$ per cent attenuation) transit recorded in our data exhibits a complex structure around minimum light that can be well modelled by multiple overlapping dust clouds. The strongest circumstellar absorption line, Fe II $\lambda$5169, significantly weakens during this transit, with its equivalent width reducing from a mean out-of-transit value of $2$ \AA\ to $1$ \AA\ in-transit, supporting spatial correlation between the circumstellar gas and dust. Finally, we made use of the Gaia Data Release 2 and archival photometry to determine the white dwarf parameters. Adopting a helium-dominated atmosphere containing traces of hydrogen and metals, and a reddening $E(B-V)=0.01$ we find $T_\mathrm{eff}=15\,020 \pm 520$ K, $\log g=8.07\pm0.07$, corresponding to $M_\mathrm{WD}=0.63\pm0.05\ \mbox{$\mathrm{M}_{\odot}$}$ and a cooling age of $224\pm30$ Myr.

We study the origin of the interstellar object 1I/2017 U1 'Oumuamua by juxtaposing estimates based on the observations with simulations. We speculate that objects like 'Oumuamua are formed in the debris disc as left over from the star and planet formation process, and subsequently liberated. The liberation process is mediated either by interaction with other stars in the parental star-cluster, by resonant interactions within the planetesimal disc or by the relatively sudden mass loss when the host star becomes a compact object. Integrating backward in time in the Galactic potential together with stars from the Gaia-TGAS catalogue we find that about 1.3Myr ago 'Oumuamua passed the nearby star HIP 17288 within a mean distance of $1.3$pc. By comparing nearby observed L-dwarfs with simulations of the Galaxy we conclude that the kinematics of 'Oumuamua is consistent with relatively young objects of $1.1$--$1.7$Gyr. We just met 'Oumuamua by chance, and with a derived mean Galactic density of $\sim 3\times 10^{5}$ similarly sized objects within 100\,au from the Sun or $\sim 10^{14}$ per cubic parsec we expect about 2 to 12 such visitors per year within 1au from the Sun.

In low-mass binary systems, mass transfer is likely to occur via a slow and dense stellar wind when one of the stars is in the AGB phase. Observations show that many binaries that have undergone AGB mass transfer have orbital periods of 1-10 yr, at odds with the predictions of binary population synthesis models. We investigate the mass-accretion efficiency and angular-momentum loss via wind mass transfer in AGB binary systems. We use these quantities to predict the evolution of the orbit. We perform 3D hydrodynamical simulations of the stellar wind lost by an AGB star using the AMUSE framework. We approximate the thermal evolution of the gas by imposing a simple effective cooling balance and we vary the orbital separation and the velocity of the stellar wind. We find that for wind velocities $v_{\infty}$ larger than the relative orbital velocity of the system $v_\mathrm{orb}$ the flow is described by the Bondi-Hoyle-Lyttleton approximation and the angular-momentum loss is modest, leading to an expansion of the orbit. For low wind velocities an accretion disk is formed around the companion and the accretion efficiency as well as the angular-momentum loss are enhanced, implying that the orbit will shrink. We find that the transfer of angular momentum from the orbit to the outflowing gas occurs within a few orbital separations from the center of mass of the binary. Our results suggest that the orbital evolution of AGB binaries can be predicted as a function of the ratio $v_{\infty}/v_\mathrm{orb}$. Our results can provide insight into the puzzling orbital periods of post-AGB binaries and suggest that the number of stars entering into the common-envelope phase will increase. The latter can have significant implications for the expected formation rates of the end products of low-mass binary evolution, such as cataclysmic binaries, type Ia supernova and double white-dwarf mergers. [ABRIDGED]

If the Sun was born in a relatively compact open cluster, it is quite likely that a massive (10MSun) star was nearby when it exploded in a supernova. The repercussions of a supernova can be rather profound, and the current Solar System may still bear the memory of this traumatic event. The truncation of the Kuiper belt and the tilt of the ecliptic plane with respect to the Sun's rotation axis could be such signatures. We simulated the effect of a nearby supernova on the young Solar System using the Astronomical Multipurpose Software Environment. Our calculations are realized in two subsequent steps in which we study the effect of the supernova irradiation on the circumstellar disk and the effect of the impact of the nuclear blast-wave which arrives a few decades later. We find that the blastwave of our adopted supernova exploding at a distance of $0.15$--$0.40$\,pc and at an angle of $35^\circ$--$65^\circ$ with respect to the angular-momentum axis of the circumsolar disk would induce a misalignment between the Sun's equator and its disk to $5^\circ.6\pm1^\circ.2$, consistent with the current value. The blast of a supernova truncates the disk at a radius between $42$ and $55$\,au, which is consistent with the current edge of the Kuiper belt. For the most favored parameters, the irradiation by the supernova as well as the blast wave heat the majority of the disk to $\sim 1200$\,K, which is sufficiently hot to melt chondrules in the circumstellar disk. The majority of planetary system may have been affected by a nearby supernova, some of its repercussions, such as truncation and tilting of the disk, may still be visible in their current planetary system's topology. The amount of material from the supernova blast wave that is accreted by the circumstellar disk is too small by several orders of magnitude to explain the current abundance of the short live radionuclide $^{26}$Al.

Recent developments in computer hardware and software enable researchers to simulate the self-gravitating evolution of galaxies at a resolution comparable to the actual number of stars. Here we present the results of a series of such simulations. We performed $N$-body simulations of disk galaxies with between 100 and 500 million particles over a wide range of initial conditions. Our calculations include a live bulge, disk, and dark matter halo, each of which is represented by self-gravitating particles in the $N$-body code. The simulations are performed using the gravitational $N$-body tree-code Bonsai running on the Piz Daint supercomputer. We find that the time scale over which the bar forms increases exponentially with decreasing disk-mass fraction and that the bar formation epoch exceeds a Hubble time when the disk-mass fraction is $\sim0.35$. These results can be explained with the swing-amplification theory. The condition for the formation of $m=2$ spirals is consistent with that for the formation of the bar, which is also an $m=2$ phenomenon. We further argue that the non-barred grand-design spiral galaxies are transitional, and that they evolve to barred galaxies on a dynamical timescale. We also confirm that the disk-mass fraction and shear rate are important parameters for the morphology of disk galaxies. The former affects the number of spiral arms and the bar formation epoch, and the latter determines the pitch angle of the spiral arms.

Numerical solutions to Newton's equations of motion for chaotic self gravitating systems of more than 2 bodies are often regarded to be irreversible. This is due to the exponential growth of errors introduced by the integration scheme and the numerical round-off in the least significant figure. This secular growth of error is sometimes attributed to the increase in entropy of the system even though Newton's equations of motion are strictly time reversible. We demonstrate that when numerical errors are reduced to below the physical perturbation and its exponential growth during integration the microscopic reversibility is retrieved. Time reversibility itself is not a guarantee for a definitive solution to the chaotic N-body problem. However, time reversible algorithms may be used to find initial conditions for which perturbed trajectories converge rather than diverge. The ability to calculate such a converging pair of solutions is a striking illustration which shows that it is possible to compute a definitive solution to a highly unstable problem. This works as follows: If you (i) use a code which is capable of producing a definitive solution (and which will therefore handle converging pairs of solutions correctly), (ii) use it to study the statistical result of some other problem, and then (iii) find that some other code produces a solution S with statistical properties which are indistinguishable from those of the definitive solution, then solution S may be deemed veracious.

Due to the high stellar densities in young clusters, planetary systems formed in these environments are likely to have experienced perturbations from encounters with other stars. We carry out direct $N$-body simulations of multi-planet systems in star clusters to study the combined effects of stellar encounters and internal planetary dynamics. These planetary systems eventually become part of the Galactic field population the parental cluster dissolves, which is where most presently-known exoplanets are observed. We show that perturbations induced by stellar encounters lead to distinct signatures in the field planetary systems, most prominently, the excited orbital inclinations and eccentricities. Planetary systems that form within the cluster's half-mass radius are more prone to such perturbations. The orbital elements are most strongly excited in the outermost orbit, but the effect propagates to the entire planetary system through secular evolution. Planet ejections may occur long after a stellar encounter. The surviving planets in these reduced systems tend to have, on average, higher inclinations and larger eccentricities compared to systems that were perturbed less strongly. As soon as the parental star cluster dissolves, external perturbations stop affecting the escaped planetary systems, and further evolution proceeds on a relaxation time scale. The outer regions of these ejected planetary systems tend to relax so slowly that their state carries the memory of their last strong encounter in the star cluster. Regardless of the stellar density, we observe a robust anticorrelation between multiplicity and mean inclination/eccentricity. We speculate that the "Kepler dichotomy" observed in field planetary systems is a natural consequence of their early evolution in the parental cluster.

The Oort cloud (OC) probably formed more than 4$\,$Gyr ago and has been moving with the Sun in the Galaxy since, exposed to external influences, most prominently to the Galactic tide and passing field stars. Theories suggest that other stars might posses exocomets distributed similarly to our OC. We study the erosion of the OC and the possibility for capturing exocomets during the encounters with such field stars. We carry out simulations of flybys, where both stars are surrounded by a cloud of comets. We measure how many exocomets are transferred to the OC, how many OC's comets are lost, and how this depends on the other star's mass, velocity and impact parameter. Exocomets are transferred to the OC only during relatively slow ($\lesssim0.5\,$km$\,$s$^{-1}$) and close ($\lesssim10^5\,$AU) flybys and these are expected to be extremely rare. Assuming that all passing stars are surrounded by a cloud of exocomets, we derive that the fraction of exocomets in the OC has been about $10^{-5}$--$10^{-4}$. Finally we simulate the OC for the whole lifetime of the Sun, taking into account the encounters and the tidal effects. The OC has lost 25--65% of its mass, mainly due to stellar encounters, and at most 10% (and usually much less) of its mass can be captured. However, exocomets are often lost shortly after the encounter that delivers them, due to the Galactic tide and consecutive encounters.

Misalignments between the orbital planes of planets and the equatorial planes of their host stars have been observed in our solar system, in transiting exoplanets, and in the orbital planes of debris discs. We present a mechanism that causes such a spin-orbit misalignment for a protoplanetary disc due to its movement through an ambient medium. Our physical explanation of the mechanism is based on the theoretical solutions to the Stark problem. We test this idea by performing self-consistent hydrodynamical simulations and simplified gravitational $N$-body simulations. The $N$-body model reduces the mechanism to the relevant physical processes. The hydrodynamical simulations show the mechanism in its full extent, including gas-dynamical and viscous processes in the disc which are not included in the theoretical framework. We find that a protoplanetary disc embedded in a flow changes its orientation as its angular momentum vector tends to align parallel to the relative velocity vector. Due to the force exerted by the flow, orbits in the disc become eccentric, which produces a net torque and consequentially changes the orbital inclination. The tilting of the disc causes it to contract. Apart from becoming lopsided, the gaseous disc also forms a spiral arm even if the inclination does not change substantially. The process is most effective at high velocities and observational signatures are therefore mostly expected in massive star-forming regions and around winds or supernova ejecta. Our $N$-body model indicates that the interaction with supernova ejecta is a viable explanation for the observed spin-orbit misalignment in our solar system.

Scientific discovery via numerical simulations is important in modern astrophysics. This relatively new branch of astrophysics has become possible due to the development of reliable numerical algorithms and the high performance of modern computing technologies. These enable the analysis of large collections of observational data and the acquisition of new data via simulations at unprecedented accuracy and resolution. Ideally, simulations run until they reach some pre-determined termination condition, but often other factors cause extensive numerical approaches to break down at an earlier stage. In those cases, processes tend to be interrupted due to unexpected events in the software or the hardware. In those cases, the scientist handles the interrupt manually, which is time-consuming and prone to errors. We present the Simulation Monitor (SiMon) to automatize the farming of large and extensive simulation processes. Our method is light-weight, it fully automates the entire workflow management, operates concurrently across multiple platforms and can be installed in user space. Inspired by the process of crop farming, we perceive each simulation as a crop in the field and running simulation becomes analogous to growing crops. With the development of SiMon we relax the technical aspects of simulation management. The initial package was developed for extensive parameter searchers in numerical simulations, but it turns out to work equally well for automating the computational processing and reduction of observational data reduction.

We show that black-hole High-Mass X-ray Binaries (HMXBs) with O- or B-type donor stars and relatively short orbital periods, of order one week to several months may survive spiral in, to then form Wolf-Rayet (WR) X-ray binaries with orbital periods of order a day to a few days; while in systems where the compact star is a neutron star, HMXBs with these orbital periods never survive spiral-in. We therefore predict that WR X-ray binaries can only harbor black holes. The reason why black-hole HMXBs with these orbital periods may survive spiral in is: the combination of a radiative envelope of the donor star, and a high mass of the compact star. In this case, when the donor begins to overflow its Roche lobe, the systems are able to spiral in slowly with stable Roche-lobe overflow, as is shown by the system SS433. In this case the transferred mass is ejected from the vicinity of the compact star (so-called "isotropic re-emission" mass loss mode, or "SS433-like mass loss"), leading to gradual spiral-in. If the mass ratio of donor and black hole is $>3.5$, these systems will go into CE evolution and are less likely to survive. If they survive, they produce WR X-ray binaries with orbital periods of a few hours to one day. Several of the well-known WR+O binaries in our Galaxy and the Magellanic Clouds, with orbital periods in the range between a week and several months, are expected to evolve into close WR-Black-Hole binaries,which may later produce close double black holes. The galactic formation rate of double black holes resulting from such systems is still uncertain, as it depends on several poorly known factors in this evolutionary picture. It might possibly be as high as $\sim 10^{-5}$ per year.

Observations indicate that the dispersal of protoplanetary discs in star clusters occurs on time scales of about 5 Myr. Several processes are thought to be responsible for this disc dispersal. Here we compare two of these processes: dynamical encounters and interaction with the interstellar medium, which includes face-on accretion and ram pressure stripping. We perform simulations of embedded star clusters with parameterisations for both processes to determine the environment in which either of these processes is dominant. We find that face-on accretion, including ram pressure stripping, is the dominant disc truncation process if the fraction of the total cluster mass in stars is $\lesssim 30\,\%$ regardless of the cluster mass and radius. Dynamical encounters require stellar densities $\gtrsim 10^4$ pc$^{-3}$ combined with a mass fraction in stars of $\approx 90\,\%$ to become the dominant process. Our results show that during the embedded phase of the cluster, the truncation of the discs is dominated by face-on accretion and dynamical encounters become dominant when the intra-cluster gas has been expelled. As a result of face-on accretion the protoplanetary discs become compact and their surface density increases. In contrast, dynamical encounters lead to discs that are less massive and remain larger.

Most stars form in star clusters and stellar associated. To understand the roles of star cluster environments in shaping the dynamical evolution of planetary systems, we carry out direct $N$-body simulations of four planetary systems models in three different star cluster environments with respectively N=2k, 8k and 32k stars. In each cluster, an ensemble of initially identical planetary systems are assigned to solar-type stars with $\sim 1 M_{\odot}$ and evolved for 50~Myr. We found that following the depletion of protoplanetary disks, external perturbations and planet-planet interactions are two driving mechanisms responsible for the destabilization of planetary systems. The planet survival rate varies from $\sim 95\%$ in the N=2k cluster to $\sim 60\%$ in the N=32k cluster, which suggests that most planetary systems can indeed survive in low-mass clusters, except in the central regions. We also find that planet ejections through stellar encounters are cumulative processes, as only $\sim 3\%$ of encounters are strong enough to excite the eccentricity by $\Delta e \geq 0.5$. Short-period planets can be perturbed through orbit crossings with long-period planets. When taking into account planet-planet interactions, the planet ejection rate nearly doubles, and therefore multiplicity contributes to the vulnerability of planetary systems. In each ensemble, $\sim 0.2\%$ of planetary orbits become retrograde due to random directions of stellar encounters. Our results predict that young low-mass star clusters are promising sites for next-generation planet surveys, yet low planet detection rates are expected in dense globular clusters such as 47 Tuc. Nevertheless, planets in denser stellar environments are likely to have shorter orbital periods, which enhances their detectability.

Observations indicate that stars generally lose their protoplanetary discs on a timescale of about 5 Myr. Which mechanisms are responsible for the disc dissipation is still debated. Here we investigate the movement through an ambient medium as a possible cause of disc dispersal. The ram pressure exerted by the flow can truncate the disc and the accretion of material with no azimuthal angular momentum leads to further disc contraction. We derive a theoretical model from accretion disc theory that describes the evolution of the disc radius, mass, and surface density profile as a function of the density and velocity of the ambient medium. We test our model by performing hydrodynamical simulations of a protoplanetary disc embedded in a flow with different velocities and densities. We find that our model gives an adequate description of the evolution of the disc radius and accretion rate onto the disc. The total disc mass in the simulations follows the theoretically expected trend, except at the lowest density where our simulated discs lose mass owing to continuous stripping. This stripping may be a numerical rather than a physical effect. Some quantitative differences exist between the model predictions and the simulations. These are at least partly caused by numerical viscous effects in the disc and depend on the resolution of the simulation. Our model can be used as a conservative estimate for the process of face-on accretion onto protoplanetary discs, as long as viscous processes in the disc can be neglected. The model predicts that in dense gaseous environments, discs can shrink substantially in size and can, in theory, sweep up an amount of gas of the order of their initial mass. This process could be relevant for planet formation in dense environments.

Field stars are frequently formed in pairs, and many of these binaries are part of triples or even higher-order systems. Even though, the principles of single stellar evolution and binary evolution, have been accepted for a long time, the long-term evolution of stellar triples is poorly understood. The presence of a third star in an orbit around a binary system can significantly alter the evolution of those stars and the binary system. The rich dynamical behavior in three-body systems can give rise to Lidov-Kozai cycles, in which the eccentricity of the inner orbit and the inclination between the inner and outer orbit vary periodically. In turn, this can lead to an enhancement of tidal effects (tidal friction), gravitational-wave emission and stellar interactions such as mass transfer and collisions. The lack of a self-consistent treatment of triple evolution, including both three-body dynamics as well as stellar evolution, hinders the systematic study and general understanding of the long-term evolution of triple systems. In this paper, we aim to address some of these hiatus, by discussing the dominant physical processes of hierarchical triple evolution, and presenting heuristic recipes for these processes. To improve our understanding on hierarchical stellar triples, these descriptions are implemented in a public source code TrES which combines three-body dynamics (based on the secular approach) with stellar evolution and their mutual influences. Note that modeling through a phase of stable mass transfer in an eccentric orbit is currently not implemented in TrES , but can be implemented with the appropriate methodology at a later stage.

Hot Jupiters (HJs) are Jupiter-like planets that reside very closely to their host star, within $\sim 0.1\,\mathrm{AU}$. Their formation is not well understood. It is generally believed that they cannot have formed in situ, implying that some form of migration must have occurred after their initial formation. We study the production of HJs through secular evolution in multiplanet systems with three to five planets. In this variant of high-$e$ migration, the eccentricity of the orbit of the innermost planet is excited on secular time-scales, triggering orbital migration due to tidal dissipation. We use a secular dynamics code and carry out a population synthesis study. We find that HJs are only produced if the viscous time-scale is short ($\approx 0.014$ yr). In contrast, in up to $\approx 0.3$ of systems, the innermost planet is tidally disrupted. The orbital period distribution is peaked around 5 d, consistent with observations. The median HJ mass is $1\,M_\mathrm{J}$ with a maximum of $\approx 2 \, M_\mathrm{J}$, similar to observed HJs. Approximately 0.1 of the HJs have retrograde orbits with respect to the stellar spin. We do not find any warm Jupiters in our simulations, i.e. planets with semimajor axes between 0.1 and 1 AU.

The frequency of Galactic stellar encounters the Solar system experienced depends on the local density and velocity dispersion along the orbit of the Sun in the Milky Way galaxy. We aim at determining the effect of the radial migration of the solar orbit on the rate of stellar encounters. As a first step we integrate the orbit of the Sun backwards in time in an analytical potential of the Milky Way. We use the present-day phase-space coordinates of the Sun, according to the measured uncertainties. The resulting orbits are inserted in an N-body simulation of the Galaxy, where the stellar velocity dispersion is calculated at each position along the orbit of the Sun. We compute the rate of Galactic stellar encounters by employing three different solar orbits ---migrating from the inner disk, without any substantial migration, and migrating from the outer disk. We find that the rate for encounters within $4\times10^5$ AU from the Sun is approximately 21, 39 and 63 Myr$^{-1}$, respectively. The stronger encounters establish the outer limit of the so-called parking zone, which is the region in the plane of the orbital eccentricities and semi-major axes where the planetesimals of the Solar system have been perturbed only by interactions with stars belonging to the Sun's birth cluster. We estimate the outer edge of the parking zone at semi-major axes of 250--1300 AU (the outward and inward migrating orbits reaching the smallest and largest values, respectively), which is one order of magnitude smaller than the determination made by Portegies Zwart & J\'ilkov\'a (2015). We further discuss the effect of stellar encounters on the stability of the hypothetical Planet 9.

Approximately $0.2 \pm 0.2$ of white dwarfs (WDs) show signs of pollution by metals, which is likely due to the accretion of tidally disrupted planetary material. Models invoking planet-planet interactions after WD formation generally cannot explain pollution at cooling times of several Gyr. We consider a scenario in which a planet is perturbed by Lidov-Kozai oscillations induced by a binary companion and exacerbated by stellar mass loss, explaining pollution at long cooling times. Our computed accretion rates are consistent with observations assuming planetary masses between $\sim 0.01$ and $1\,M_\mathrm{Mars}$, although nongravitational effects may already be important for masses $\lesssim 0.3 \, M_\mathrm{Mars}$. The fraction of polluted WDs in our simulations, $\sim 0.05$, is consistent with observations of WDs with intermediate cooling times between $\sim 0.1$ and 1 Gyr. For cooling times $\lesssim 0.1$ Gyr and $\gtrsim 1$ Gyr, our scenario cannot explain the high observed pollution fractions of up to 0.7. Nevertheless, our results motivate searches for companions around polluted WDs.

Globular clusters (GCs) are known to harbor multiple stellar populations. To explain these observations Bastian et al. suggested a scenario in which a second population is formed by the accretion of enriched material onto the low-mass stars in the initial GC population. The idea is that the low-mass, pre-main sequence stars sweep up gas expelled by the massive stars of the same generation into their protoplanetary disc as they move through the GC core. We perform simulations with 2 different smoothed particle hydrodynamics codes to investigate if a low-mass star surrounded by a protoplanetary disc can accrete the amount of enriched material required in this scenario. We focus on the gas loading rate onto the disc and star as well as on the lifetime of the disc. We find that the gas loading rate is a factor of 2 smaller than the geometric rate, because the effective cross section of the disc is smaller than its surface area. The loading rate is consistent for both codes, irrespective of resolution. The disc gains mass in the high resolution runs, but loses angular momentum on a time scale of 10^4 yrs. Two effects determine the loss of (specific) angular momentum in our simulations: 1) continuous ram pressure stripping and 2) accretion of material with no azimuthal angular momentum. Our study and previous work suggest that the former, dominant process is mainly caused by numerical rather than physical effects, while the latter is not. The latter process causes the disc to become more compact, increasing the surface density profile at smaller radii. The disc size is determined in the first place by the ram pressure when the flow first hits the disc. Further evolution is governed by the decrease in the specific angular momentum of the disc. We conclude that the size and lifetime of the disc are probably not sufficient to accrete the amount of mass required in Bastian et al.'s scenario.

According to Munoz-Gutierrez et al. (2015) the orbit of comet 1P/Halley is chaotic with a surprisingly small Lyapunov time scale of order its orbital period. In this work we analyse the origin of chaos in Halley's orbit and the growth of perturbations, in order to get a better understanding of this unusually short time scale. We perform N-body simulations to model Halley's orbit in the Solar System and measure the separation between neighbouring trajectories. To be able to interpret the numerical results, we use a semi-analytical map to demonstrate different growth modes, i.e. linear, oscillatory or exponential, and transitions between these modes. We find the Lyapunov time scale of Halley's orbit to be of order 300 years, which is significantly longer than previous estimates in the literature. This discrepancy could be due to the different methods used to measure the Lyapunov time scale. A surprising result is that next to Jupiter, also encounters with Venus contribute to the exponential growth in the next 3000 years. Finally, we note an interesting application of the sub-linear, oscillatory growth mode to an ensemble of bodies moving through the Solar System. Whereas in the absence of encounters with a third body the ensemble spreads out linearly in time, the accumulation of weak encounters can increase the lifetime of such systems due to the oscillatory behaviour.

Transiting circumbinary planets are more easily detected around short-period than long-period binaries, but none have yet been observed by {\it Kepler} orbiting binaries with periods shorter than seven days. In triple systems, secular Kozai-Lidov cycles and tidal friction (KLCTF) have been shown to reduce the inner orbital period from $\sim 10^4$ to a few days. Indeed, the majority of short-period binaries are observed to possess a third stellar companion. Using secular evolution analysis and population synthesis, we show that KLCTF makes it unlikely for circumbinary transiting planets to exist around short-period binaries. We find the following outcomes. (1) Sufficiently massive planets in tight and/or coplanar orbits around the inner binary can quench the KL evolution because they induce precession in the inner binary. The KLCTF process does not take place, preventing the formation of a short-period binary. (2) Secular evolution is not quenched and it drives the planetary orbit into a high eccentricity, giving rise to an unstable configuration, in which the planet is most likely ejected from the system. (3) Secular evolution is not quenched but the planet survives the KLCTF evolution. Its orbit is likely to be much wider than the currently observed inner binary orbit, and is likely to be eccentric and inclined with respect to the inner binary. These outcomes lead to two main conclusions: (1) it is unlikely to find a massive planet on a tight and coplanar orbit around a short-period binary, and (2) the properties of circumbinary planets in short-period binaries are constrained by secular evolution.

We present a method for studying the secular gravitational dynamics of hierarchical multiple systems consisting of nested binaries, which is valid for an arbitrary number of bodies and arbitrary hierarchical structure. We derive the Hamiltonian of the system and expand it in terms of the -- assumed to be -- small ratios $x_i$ of binary separations. At the lowest nontrivial expansion order (quadrupole order, second order in $x_i$), the Hamiltonian consists of terms which, individually, depend on binary pairs. At higher orders, in addition to terms depending on binary pairs, we also find terms which, individually, depend on more than two binaries. In general, at order $n$ in $x_i$, individual terms depend on at most $n-1$ binaries. We explicitly derive the Hamiltonian including all terms up and including third order in $x_i$ (octupole order), and including the binary pairwise terms up and including fifth order in $x_i$. These terms are orbit averaged, and we present a new algorithm for efficiently solving the equations of motion. This algorithm is highly suitable for studying the secular evolution of hierarchical systems with complex hierarchies, making long-term integrations of such systems feasible. We show that accurate results are obtained for multiplanet systems with semimajor axis ratios as large as $\approx 0.4$, provided that high-order terms are included. In addition to multiplanet systems with a single star, we apply our results to multistar systems with multiple planets.

We study the secular gravitational dynamics of quadruple systems consisting of a hierarchical triple system orbited by a fourth body. These systems can be decomposed into three binary systems with increasing semimajor axes, binaries A, B and C. The Hamiltonian of the system is expanded in ratios of the three binary separations, and orbit-averaged. Subsequently, we numerically solve the equations of motion. We study highly hierarchical systems that are well described by the lowest-order terms in the Hamiltonian. We find that the qualitative behaviour is determined by the ratio $\mathcal{R}_0$ of the initial Kozai-Lidov (KL) time-scales of the binary pairs AB and BC. If $\mathcal{R}_0\ll 1$, binaries AB remain coplanar if this is initially the case, and KL eccentricity oscillations in binary B are efficiently quenched. If $\mathcal{R}_0\gg 1$, binaries AB become inclined, even if initially coplanar. However, there are no induced KL eccentricity oscillations in binary A. Lastly, if $\mathcal{R}_0\sim 1$, complex KL eccentricity oscillations can occur in binary A that are coupled with the KL eccentricity oscillations in B. Even if binaries A and B are initially coplanar, the induced inclination can result in very high eccentricity oscillations in binary A. These extreme eccentricities could have significant implications for strong interactions such as tidal interactions, gravitational wave dissipation, and collisions and mergers of stars and compact objects. As an example, we apply our results to a planet+moon system orbiting a central star, which in turn is orbited by a distant and inclined stellar companion or planet, and to observed stellar quadruples.

We study mass transfers between debris discs during stellar encounters. We carried out numerical simulations of close flybys of two stars, one of which has a disc of planetesimals represented by test particles. We explored the parameter space of the encounters, varying the mass ratio of the two stars, their pericentre and eccentricity of the encounter, and its geometry. We find that particles are transferred to the other star from a restricted radial range in the disc and the limiting radii of this transfer region depend on the parameters of the encounter. We derive an approximate analytic description of the inner radius of the region. The efficiency of the mass transfer generally decreases with increasing encounter pericentre and increasing mass of the star initially possessing the disc. Depending on the parameters of the encounter, the transfer particles have a specific distributions in the space of orbital elements (semimajor axis, eccentricity, inclination, and argument of pericentre) around their new host star. The population of the transferred particles can be used to constrain the encounter through which it was delivered. We expect that many stars experienced transfer among their debris discs and planetary systems in their birth environment. This mechanism presents a formation channel for objects on wide orbits of arbitrary inclinations, typically having high eccentricity but possibly also close-to-circular (eccentricities of about 0.1). Depending on the geometry, such orbital elements can be distinct from those of the objects formed around the star.

Appreciable star formation, and, therefore, numerous massive stars, are frequently found near supermassive black holes (SMBHs). As a result, core-collapse supernovae in these regions should also be expected. In this paper, we consider the observational consequences of predicting the fate of supernova remnants (SNRs) in the sphere of influence of quiescent SMBHs. We present these results in the context of `autarkic' nuclei, a model that describes quiescent nuclei as steady-state and self-sufficient environments where the SMBH accretes stellar winds with no appreciable inflow of material from beyond the sphere of influence. These regions have properties such as gas density that scale with the mass of the SMBH. Using predictions of the X-ray lifetimes of SNRs originating in the sphere of influence, we make estimates of the number of core collapse SNRs present at a given time. With the knowledge of lifetimes of SNRs and their association with young stars, we predict a number of core-collapse SNRs that grows from ~1 around Milky Way-like (4.3 x 10^6 Msun) SMBHs to ~100 around the highest mass (10^10 Msun) SMBHs. The presence of young SNRs will amplify the X-ray emission near quiescent SMBHs, and we show that the total core-collapse SNR emission has the potential to influence soft X-ray searches for very low-luminosity SMBHs. Our SNR lifetime estimates also allow us to predict star formation rates in these regions. Assuming a steady-state replenishment of massive stars, we estimate a star formation rate density of 2 x 10^-4 Msun/yr/pc^2 around the Milky Way SMBH, and a similar value around other SMBHs due to a weak dependence on SMBH mass. This value is consistent with currently available observations.

We use self-consistent numerical simulations of the evolution and disruption of the Sun's birth cluster in the Milky Way potential to investigate the present-day phase space distribution of the Sun's siblings. The simulations include the gravitational N-body forces within the cluster and the effects of stellar evolution on the cluster population. In addition the gravitational forces due to the Milky Way potential are accounted for in a self-consistent manner. Our aim is to understand how the astrometric and radial velocity data from the Gaia mission can be used to pre-select solar sibling candidates. We vary the initial conditions of the Sun's birth cluster, as well as the parameters of the Galactic potential. We show that the disruption time-scales of the cluster are insensitive to the details of the non-axisymmetric components of the Milky Way model and we make predictions, averaged over the different simulated possibilities, about the number of solar siblings that should appear in surveys such as Gaia or GALAH. We find a large variety of present-day phase space distributions of solar siblings, which depend on the cluster initial conditions and the Milky Way model parameters. We show that nevertheless robust predictions can be made about the location of the solar siblings in the space of parallaxes ($\varpi$), proper motions ($\mu$) and radial velocities ($V_\mathrm{r}$). By calculating the ratio of the number of simulated solar siblings to that of the number of stars in a model Galactic disk, we find that this ratio is above 0.5 in the region given by: $\varpi \geq 5$mas, $4 \leq \mu \leq 6$masyr$^{-1}$, and $-2\leq V_\mathrm{r} \leq 0$kms$^{-1}$. Selecting stars from this region should increase the probability of success in identifying solar siblings through follow up observations [Abridged].

The Galactic center provides a unique laboratory to study the interaction of a supermassive black hole (SMBH) with its gaseous and stellar environment. Simulations to determine the accretion of stellar winds from the surrounding O-stars onto the black hole have been performed earlier, but in those the presence of the S-star system was ignored. The S-stars are a group of young massive B-stars in relatively close orbits around the black hole. Here we simulate those stars in order to study their contribution to the accretion rate, without taking the more distant and massive O-stars into account. We use the Astrophysical Multi-purpose Software Environment (AMUSE) to combine gravitational physics, stellar evolution and hydrodynamics in a single simulation of the S-stars orbiting the supermassive black hole, and use this framework to determine the amount of gas that is accreted onto the black hole. We find that the accretion rate is sensitive to the wind properties of the S-stars (rate of mass-loss and terminal velocity). Our simulations are consistent with the observed accretion rate of the black hole only if the stars exhibit high wind massloss rates that are comparable with those of evolved 7-10 Myr old stars with masses of M=19-25 M_SUN. This is in contrast with observations that have shown that these stars are rather young, main-sequence B-stars. We therefore conclude that the S-stars cannot account for the accretion rate alone.

The Galactic center hosts several hundred early-type stars, about 20% of which lie in the so-called clockwise disk, while the remaining 80% do not belong to any disks. The circumnuclear ring (CNR), a ring of molecular gas that orbits the supermassive black hole (SMBH) with a radius of 1.5 pc, has been claimed to induce precession and Kozai-Lidov oscillations onto the orbits of stars in the innermost parsec. We investigate the perturbations exerted by a gas ring on a nearly-Keplerian stellar disk orbiting a SMBH by means of combined direct N-body and smoothed particle hydrodynamics simulations. We simulate the formation of gas rings through the infall and disruption of a molecular gas cloud, adopting different inclinations between the infalling gas cloud and the stellar disk. We find that a CNR-like ring is not efficient in affecting the stellar disk on a timescale of 3 Myr. In contrast, a gas ring in the innermost 0.5 pc induces precession of the longitude of the ascending node Omega, significantly affecting the stellar disk inclination. Furthermore, the combined effect of two-body relaxation and Omega-precession drives the stellar disk dismembering, displacing the stars from the disk. The impact of precession on the star orbits is stronger when the stellar disk and the inner gas ring are nearly coplanar. We speculate that the warm gas in the inner cavity might have played a major role in the evolution of the clockwise disk.

Recent observations have revealed a variety of young star clusters, including embedded systems, young massive clusters, and associations. We study the formation and dynamical evolution of these clusters using a combination of simulations and theoretical models. Our simulations start with a turbulent molecular cloud that collapses under its own gravity. The stars are assumed to form in the densest regions in the collapsing cloud after an initial free-fall times of the molecular cloud. The dynamical evolution of these stellar distributions are continued by means of direct $N$-body simulations. The molecular clouds typical for the Milky Way Galaxy tend to form embedded clusters which evolve to resemble open clusters. The associations were initially considerably more clumpy, but lost their irregularity in about a dynamical time scale due to the relaxation process. The densest molecular clouds, which are absent in the Milky Way but are typical in starburst galaxies, form massive young star clusters. They indeed are rare in the Milky Way. Our models indicate a distinct evolutionary path from molecular clouds to open clusters and associations or to massive star clusters. The mass-radius relation for both types of evolutionary tracks excellently matches the observations. According to our calculations the time evolution of the half-mass radius for open clusters and associations follows $r_{\rm h}/{\rm pc}=2.7(t_{\rm age}/{\rm pc})^{2/3}$, whereas for massive star clusters $r_{\rm h}/{\rm pc}=0.34(t_{\rm age}/{\rm Myr})^{2/3}$. Both trends are consistent with the observed age-mass-radius relation for clusters in the Milky Way.

We compare the observed size distribution of circum stellar disks in the Orion Trapezium cluster with the results of $N$-body simulations in which we incorporated an heuristic prescription for the evolution of these disks. In our simulations, the sizes of stellar disks are affected by close encounters with other stars (with disks). We find that the observed distribution of disk sizes in the Orion Trapezium cluster is excellently reproduced by truncation due to dynamical encounters alone. The observed distribution appears to be a sensitive measure of the past dynamical history of the cluster, and therewith on the conditions of the cluster at birth. The best comparison between the observed disk size distribution and the simulated distribution is realized with a cluster of $N = 2500\pm500$ stars with a half-mass radius of about 0.5\,pc in virial equilibrium (with a virial ratio of $Q = 0.5$, or somewhat colder $Q \simeq 0.3$), and with a density structure according to a fractal dimension of $F \simeq 1.6$. Simulations with these parameters reproduce the observed distribution of circum stellar disks in about 0.2--0.5\,Myr.

We discuss the events that led to the giant eruption of Eta Carinae, and find that the mid-nineteenth century (in 1838-1843) giant mass-loss outburst has the characteristics of being produced by the merger event of a massive close binary, triggered by the gravitational interaction with a massive third companion star, which is the current binary companion in the Eta Carinae system. We come to this conclusion by a combination of theoretical arguments supported by computer simulations using the Astrophysical Multipurpose Software Environment. According to this model the $\sim 90$\,\MSun\, present primary star of the highly eccentric Eta Carinae binary system is the product of this merger, and its $\sim 30$\,\MSun\, companion originally was the third star in the system. In our model the Homunculus nebula was produced by an extremely enhanced stellar wind, energized by tidal energy dissipation prior to the merger, which enormously boosted the radiation-driven wind mass-loss. The current orbital plane is then aligned with the equatorial plane of the Homunculus, and the symmetric lobes are roughly aligned with the argument of periastron of the current Eta Carina binary. The merger itself then occurred in 1838, which resulted in a massive asymmetric outflow in the equatorial plane of the Homunculus. The 1843 outburst can in our model be attributed to the subsequent encounter when the companion star (once the outer most star in the triple system) plunges through the bloated envelope of the merger product, once when it passed periastron again. We predict that the system has an excess space velocity of order 50\,km/s in the equatorial plane of the Homunculus. Our triple model gives a viable explanation for the high runaway velocities typically observed in LBVs \citep{2015MNRAS.447..598S}.

Astrophysical direct $N$-body methods have been one of the first production algorithms to be implemented using NVIDIA's CUDA architecture. Now, almost seven years later, the GPU is the most used accelerator device in astronomy for simulating stellar systems. In this paper we present the implementation of the Sapporo2 $N$-body library, which allows researchers to use the GPU for $N$-body simulations with little to no effort. The first version, released five years ago, is actively used, but lacks advanced features and versatility in numerical precision and support for higher order integrators. In this updated version we have rebuilt the code from scratch and added support for OpenCL, multi-precision and higher order integrators. We show how to tune these codes for different GPU architectures and present how to continue utilizing the GPU optimal even when only a small number of particles ($N < 100$) is integrated. This careful tuning allows Sapporo2 to be faster than Sapporo1 even with the added options and double precision data loads. The code runs on a range of NVIDIA and AMD GPUs in single and double precision accuracy. With the addition of OpenCL support the library is also able to run on CPUs and other accelerators that support OpenCL.

We perform smoothed-particle hydrodynamical simulations of the explosion of a helium star in a close binary system, and study the effects of the explosion on the companion star as well as the effect of the presence of the companion on the supernova remnant. By simulating the mechanism of the supernova from just after core bounce until the remnant shell passes the stellar companion, we are able to separate the various effects leading to the final system parameters. In the final system, we measure the mass stripping and ablation from, and the velocity kick imparted to, the companion star, as well as the structure of the supernova shell. The presence of the companion star produces a conical cavity in the expanding supernova remnant, and loss of material from the companion causes the supernova remnant to be more metal-rich on one side and more hydrogen-rich (from the companion material) around the cavity. Following the removal of mass from the companion, we study its subsequent evolution and compare it with a single star not subjected to a supernova impact.

The X-ray source HLX-1 near the spiral galaxy ESO 243-49 is currently the best intermediate-mass black hole candidate. It has a peak bolometric luminosity of $10^{42}$ erg s$^{-1}$, which implies a mass inflow rate of $\sim10^{-4}$ MSun yr$^{-1}$, but the origin of this mass is unknown. It has been proposed that there is a star on an eccentric orbit around the black hole which transfers mass at pericentre. To investigate the orbital evolution of this system, we perform stellar evolution simulations using mesa and SPH simulations of a stellar orbit around an intermediate-mass black hole using fi. We run and couple these simulations using the amuse framework. We find that mass is lost through both the first and second Lagrange points and that there is a delay of up to 10 days between the pericentre passage and the peak mass loss event. The orbital evolution timescales we find in our simulations are larger than what is predicted by analytical models, but these models fall within the errors of our results. Despite the fast orbital evolution, we are unable to reproduce the observed change in outburst period. We conclude that the change in the stellar orbit with the system parameters investigated here is unable to account for all observed features of HLX-1.

The discovery of 2012VP113 initiated the debate on the origin of the Sedna family of planetesimals in orbit around the Sun. Sednitos roam the outer regions of the Solar System between the Egeworth--Kuiper belt and the Oort cloud, in extraordinary wide (a>150au) orbits with a large perihelion distance of q>30au compared to the Earth's (a=1au and eccentricity e=(1-q/a) ~ 0.0167 or q=1au). This population is composed of a dozen objects, which we consider a family because they have similar perihelion distance and inclination with respect to the ecliptic i=10--30deg. They also have similar argument of perihelion omega=340+/-55deg. There is no ready explanation for their origin. Here we show that these orbital parameters are typical for a captured population from the planetesimal disk of another star.Assuming the orbital elements of Sednitos have not changed since they acquired their orbits, we reconstruct the encounter that led to their capture. We conclude that they might have been captured in a near miss with a 1.8MSun star that impacted the Sun at ~340au at an inclination with respect to the ecliptic of 17--34deg with a relative velocity at infinity of ~4.3km/s. We predict that the Sednitos-region is populated by 930 planetesimals and the inner Oort cloud acquired ~440 planetesimals through the same encounter.

We have coupled a fast, parametrized star cluster evolution code to a Markov Chain Monte Carlo code to determine the distribution of probable initial conditions of observed star clusters, which may serve as a starting point for future $N$-body calculations. In this paper we validate our method by applying it to a set of star clusters which have been studied in detail numerically with $N$-body simulations and Monte Carlo methods: the Galactic globular clusters M4, 47 Tucanae, NGC 6397, M22, $\omega$ Centauri, Palomar 14 and Palomar 4, the Galactic open cluster M67, and the M31 globular cluster G1. For each cluster we derive a distribution of initial conditions that, after evolution up to the cluster's current age, evolves to the currently observed conditions. We find that there is a connection between the morphology of the distribution of initial conditions and the dynamical age of a cluster and that a degeneracy in the initial half-mass radius towards small radii is present for clusters which have undergone a core collapse during their evolution. We find that the results of our method are in agreement with $N$-body and Monte Carlo studies for the majority of clusters. We conclude that our method is able to find reliable posteriors for the determined initial mass and half-mass radius for observed star clusters, and thus forms an suitable starting point for modeling an observed cluster\rq{}s evolution.

We describe the political and technical complications encountered during the astronomical CosmoGrid project. CosmoGrid is a numerical study on the formation of large scale structure in the universe. The simulations are challenging due to the enormous dynamic range in spatial and temporal coordinates, as well as the enormous computer resources required. In CosmoGrid we dealt with the computational requirements by connecting up to four supercomputers via an optical network and make them operate as a single machine. This was challenging, if only for the fact that the supercomputers of our choice are separated by half the planet, as three of them are located scattered across Europe and fourth one is in Tokyo. The co-scheduling of multiple computers and the 'gridification' of the code enabled us to achieve an efficiency of up to $93\%$ for this distributed intercontinental supercomputer. In this work, we find that high-performance computing on a grid can be done much more effectively if the sites involved are willing to be flexible about their user policies, and that having facilities to provide such flexibility could be key to strengthening the position of the HPC community in an increasingly Cloud-dominated computing landscape. Given that smaller computer clusters owned by research groups or university departments usually have flexible user policies, we argue that it could be easier to instead realize distributed supercomputing by combining tens, hundreds or even thousands of these resources.

We describe the implementation and performance of the ${\rm P^3T}$ (Particle-Particle Particle-Tree) scheme for simulating dense stellar systems. In ${\rm P^3T}$, the force experienced by a particle is split into short-range and long-range contributions. Short-range forces are evaluated by direct summation and integrated with the fourth order Hermite predictor-corrector method with the block timesteps. For long-range forces, we use a combination of the Barnes-Hut tree code and the leapfrog integrator. The tree part of our simulation environment is accelerated using graphical processing units (GPU), whereas the direct summation is carried out on the host CPU. Our code gives excellent performance and accuracy for star cluster simulations with a large number of particles even when the core size of the star cluster is small.

We specify the range to which perturbations penetrate a planetesimal system. Such perturbations can originate from massive planets or from encounters with other stars. The latter can have an origin in the star cluster in which the planetary system was born, or from random encounters once the planetary system has escaped its parental cluster. The probability of a random encounter, either in a star cluster or in the Galactic field depends on the local stellar density, the velocity dispersion and the time spend in that environment. By adopting order of magnitude estimates we argue that the majority of planetary systems born in open clusters will have a {\em Parking zone}, in which planetesimals are affected by encounters in their parental star cluster but remain unperturbed after the star has left the cluster. Objects found in this range of semi-major axis and eccentricity preserve the memory of the encounter that last affected their orbits, and they can therefore be used to reconstruct this encounter. Planetary systems born in a denser environment, such as in a globular cluster are unlikely to have a Parking zone. We further argue that some planetary systems may have a {\em Frozen zone}, in which orbits are not affected either by the more inner massive planets or by external influences. Objects discovered in this zone will have preserved information about their formation in their orbital parameters.

The 13 Myr old star HD106906 is orbited by a debris disk of at least 0.067 M_Moon with an inner and outer radius of 20 AU and 120 AU, respectively, and by a planet at a distance of 650 AU. We use this curious combination of a close low-mass disk and a wide planet to motivate our simulations of this system. We study the parameter space of the initial conditions to quantify the mass loss from the debris disk and its lifetime under the influence of the planet. We find that when the planet orbits closer to the star than about 50 AU and with low inclination relative to the disk (less than about 10 degrees), more disk material is perturbed outside than inside the region constrained by observations on timescales shorter than 1 Myr. Considering the age of the system, such a short lifetime of the disk is incompatible with the timescale for planet--planet scattering which is one of the scenarios suggested to explain the wide separation of the planet. For some configurations when the planet's orbit is inclined with respect to the disk, the latter will start to wobble. We argue that this wobbling is caused by a mechanism similar to the Kozai--Lidov oscillations. We also observe various resonant structures (such as rings and spiral arms) induced in the disk by the planet.

We simulate the formation and evolution of young star clusters using the combination of SPH simulations and direct N-body simulations. We start by performing SPH simulations of the giant molecular cloud with a turbulent velocity field, a mass of $4\times10^4$ to $5\times10^6M_{\odot}$, and a density between $1.7\times10^3$ and $170cm^{-3}$. We continue the SPH simulations for a free-fall time scale, and analyze the resulting structure of the collapsed cloud. We subsequently replace a density-selected subset of SPH particles with stars by adopting a local star-formation efficiency proportional to $\rho^{1/2}$. As a consequence, the local star formation efficiency exceeds 30 %, whereas globally only a few % of the gas is converted to stars. The stellar distribution by the time gas is converted to stars is very clumpy, with typically a dozen bound conglomerates that consist of 100 to $10^4$ stars. We continue to evolve the stars dynamically using the collisional N-body method, which accurately treats all pairwise interactions, stellar collisions and stellar evolution. We analyze the results of the N-body simulations when the stars have an age of 2 Myr and 10 Myr. During the dynamical simulations, massive clusters grow via hierarchical merging of smaller clusters. The shape of the cluster mass function that originates from an individual molecular cloud is consistent with a Schechter function with a power-law slope of -1.73 at 2 Myr and -1.67 at 10 Myr, which fits to observed cluster mass function of the Carina region. The superposition of mass functions have a power-law slope of < -2, which fits the observed mass function of star clusters in the Milky Way, M31 and M83. We further find that the mass of the most massive cluster formed in a single molecular cloud with a mass of $M_g$ scales with $6.1M_g^{0.51}$ which also agrees with recent observation of the GMC and young clusters in M51.

We perform realistic simulations of the Sun's birth cluster in order to predict the current distribution of solar siblings in the Galaxy. We study the possibility of finding the solar siblings in the Gaia catalogue by using only positional and kinematic information. We find that the number of solar siblings predicted to be observed by Gaia will be around 100 in the most optimistic case, and that a phase space only search in the Gaia catalogue will be extremely difficult. It is therefore mandatory to combine the chemical tagging technique with phase space selection criteria in order to have any hope of finding the solar siblings.

The study of intermediate-mass black holes (IMBHs) is a young and promising field of research. Formed by runaway collisions of massive stars in young and dense stellar clusters, intermediate-mass black holes could still be present in the centers of globular clusters, today. Our group investigated the presence of intermediate-mass black holes for a sample of 10 Galactic globular clusters. We measured the inner kinematic profiles with integral-field spectroscopy and determined masses or upper limits of central black holes in each cluster. In combination with literature data we further studied the positions of our results on known black-hole scaling relations (such as M_bh - sigma) and found a similar but flatter correlation for IMBHs. Applying cluster evolution codes, the change in the slope could be explained with the stellar mass loss occurring in clusters in a tidal field over its life time. Furthermore, we present results from several numerical simulations on the topic of IMBHs and integral field units (IFUs). We ran N-body simulations of globular clusters containing IMBHs in a tidal field and studied their effects on mass-loss rates and remnant fractions and showed that an IMBH in the center prevents core collapse and ejects massive objects more rapidly. These simulations were further used to simulate IFU data cubes. For the specific case of NGC 6388 we simulated two different IFU techniques and found that velocity dispersion measurements from individual velocities are strongly biased towards lower values due to blends of neighbouring stars and background light. In addition, we use the Astrophysical Multipurpose Software Environment (AMUSE) to combine gravitational physics, stellar evolution and hydrodynamics to simulate the accretion of stellar winds onto a black hole.

The solar system started to form about 4.56 Gyr ago and despite the long intervening time span, there still exist several clues about its formation. The three major sources for this information are meteorites, the present solar system structure and the planet-forming systems around young stars. In this introduction we give an overview of the current understanding of the solar system formation from all these different research fields. This includes the question of the lifetime of the solar protoplanetary disc, the different stages of planet formation, their duration, and their relative importance. We consider whether meteorite evidence and observations of protoplanetary discs point in the same direction. This will tell us whether our solar system had a typical formation history or an exceptional one. There are also many indications that the solar system formed as part of a star cluster. Here we examine the types of cluster the Sun could have formed in, especially whether its stellar density was at any stage high enough to influence the properties of today's solar system. The likelihood of identifying siblings of the Sun is discussed. Finally, the possible dynamical evolution of the solar system since its formation and its future are considered.

There is mounting observational evidence that most galactic nuclei host both supermassive black holes (SMBHs) and young populations of stars. With an abundance of massive stars, core-collapse supernovae are expected in SMBH spheres of influence. We develop a novel numerical method, based on the Kompaneets approximation, to trace supernova remnant (SNR) evolution in these hostile environments, where radial gas gradients and SMBH tides are present. We trace the adiabatic evolution of the SNR shock until 50% of the remnant is either in the radiative phase or is slowed down below the SMBH Keplerian velocity and is sheared apart. In this way, we obtain shapes and lifetimes of SNRs as a function of the explosion distance from the SMBH, the gas density profile and the SMBH mass. As an application, we focus here exclusively on quiescent SMBHs, because their light may not hamper detections of SNRs and because we can take advantage of the unsurpassed detailed observations of our Galactic Centre. Assuming that properties such as gas and stellar content scale appropriately with the SMBH mass, we study SNR evolution around other quiescent SMBHs. We find that, for SMBH masses over ~10^7 Msun, tidal disruption of SNRs can occur at less than 10^4 yr, leading to a shortened X-ray emitting adiabatic phase, and to no radiative phase. On the other hand, only modest disruption is expected in our Galactic Centre for SNRs in their X-ray stage. This is in accordance with estimates of the lifetime of the Sgr A East SNR, which leads us to expect one supernova per 10^4 yr in the sphere of influence of Sgr A*.

We have simulated, for the first time, the long term evolution of the Milky Way Galaxy using 51 billion particles on the Swiss Piz Daint supercomputer with our $N$-body gravitational tree-code Bonsai. Herein, we describe the scientific motivation and numerical algorithms. The Milky Way model was simulated for 6 billion years, during which the bar structure and spiral arms were fully formed. This improves upon previous simulations by using 1000 times more particles, and provides a wealth of new data that can be directly compared with observations. We also report the scalability on both the Swiss Piz Daint and the US ORNL Titan. On Piz Daint the parallel efficiency of Bonsai was above 95%. The highest performance was achieved with a 242 billion particle Milky Way model using 18600 GPUs on Titan, thereby reaching a sustained GPU and application performance of 33.49 Pflops and 24.77 Pflops respectively.

The general consensus in the N-body community is that statistical results of an ensemble of collisional N-body simulations are accurate, even though individual simulations are not. A way to test this hypothesis is to make a direct comparison of an ensemble of solutions obtained by conventional methods with an ensemble of true solutions. In order to make this possible, we wrote an N-body code called Brutus, that uses arbitrary-precision arithmetic. In combination with the Bulirsch--Stoer method, Brutus is able to obtain converged solutions, which are true up to a specified number of digits. We perform simulations of democratic 3-body systems, where after a sequence of resonances and ejections, a final configuration is reached consisting of a permanent binary and an escaping star. We do this with conventional double-precision methods, and with Brutus; both have the same set of initial conditions and initial realisations. The ensemble of solutions from the conventional simulations is compared directly to that of the converged simulations, both as an ensemble and on an individual basis to determine the distribution of the errors. We find that on average at least half of the conventional simulations diverge from the converged solution, such that the two solutions are microscopically incomparable. For the solutions which have not diverged significantly, we observe that if the integrator has a bias in energy and angular momentum, this propagates to a bias in the statistical properties of the binaries. In the case when the conventional solution has diverged onto an entirely different trajectory in phase-space, we find that the errors are centred around zero and symmetric; the error due to divergence is unbiased, as long as the time-step parameter, eta <= 2^(-5) and when simulations which violate energy conservation by more than 10% are excluded.

The aim of this paper is to find lost siblings of the Sun by analyzing high resolution spectra. Finding solar siblings will enable us to constrain the parameters of the parental cluster and the birth place of the Sun in the Galaxy. The solar siblings can be identified by accurate measurements of metallicity, stellar age and elemental abundances for solar neighbourhood stars. The solar siblings candidates were kinematically selected based on their proper motions, parallaxes and colours. Stellar parameters were determined through a purely spectroscopic approach and partly physical method, respectively. Comparing synthetic with observed spectra, elemental abundances were computed based on the stellar parameters obtained using a partly physical method. A chemical tagging technique was used to identify the solar siblings. We present stellar parameters, stellar ages, and detailed elemental abundances for Na, Mg, Al, Si, Ca, Ti, Cr, Fe, and Ni for 32 solar sibling candidates. Our abundances analysis shows that four stars are chemically homogenous together with the Sun. Technique of chemical tagging gives us a high probability that they might be from the same open cluster. Only one candidate HIP 40317 which has solar metallicity and age could be a solar sibling. We performed simulations of the Sun's birth cluster in analytical Galactic model and found that most of the radial velocities of the solar siblings lie in the range $-10 \leq \mathrm{V_r}\leq 10$ $\mathrm{km~s^{-1}}$, which is smaller than the radial velocity of HIP 40317 $(\mathrm{V_r} = 34.2~\mathrm{km~s^{-1}})$, under different Galactic parameters and different initial conditions of the Sun's birth cluster. The sibling status for HIP 40317 is not directly supported by our dynamical analysis.

We study the formation and evolution of the cosmic web, using the high-resolution CosmoGrid \$\Lambda\$CDM simulation. In particular, we investigate the evolution of the large-scale structure around void halo groups, and compare this to observations of the VGS-31 galaxy group, which consists of three interacting galaxies inside a large void. The structure around such haloes shows a great deal of tenuous structure, with most of such systems being embedded in intra-void filaments and walls. We use the Nexus+ algorithm to detect walls and filaments in CosmoGrid, and find them to be present and detectable at every scale. The void regions embed tenuous walls, which in turn embed tenuous filaments. We hypothesize that the void galaxy group of VGS-31 formed in such an environment.

Flares in X-ray and near infrared are observed above the quiescent emission of the supermassive black hole (SBH) in the Galactic Centre (GC) at a rate of approximately once per day. One proposed energy source for these flares is the tidal disruption of planetesimals with radius $\gtrsim 10$ km passing within $\sim$ 1 AU of the SBH. Very little is known about the formation and evolution of planetesimals in galactic nuclei such as the GC, making predictions for flaring event rates uncertain. We explore two scenarios for the formation of planetesimals in the GC: (1) in a large-scale cloud bound to the SBH, and (2) in debris discs around stars. We model their orbital evolution around the SBH using the Fokker-Planck equation and investigate the effect of gravitational interactions with various relevant perturbers. Our predicted flaring rate, $\approx 0.6 \, \mathrm{day^{-1}}$, is nearly independent of the distribution of perturbers. Moreover, it is insensitive to scenarios (1) or (2). The assumed number of planetesimals per star is consistent with debris discs around stars in the Solar neighbourhood. In scenario (1) this implies that the number of planetesimals formed in the large-scale cloud is strongly correlated with the number of stars, and this requires finetuning for our results to be consistent with the observed flaring rate. We favour the alternative explanation that planetesimals in the GC are formed in debris discs around stars, similar to the Solar neighbourhood.

The determination of the birth radius of the Sun is important to understand the evolution and consequent disruption of the Sun's birth cluster in the Galaxy. Motivated by this fact, we study the motion of the Sun in the Milky Way during the last 4.6 Gyr in order to find its birth radius. We carried out orbit integrations backward in time using an analytical model of the Galaxy which includes the contribution of spiral arms and a central bar. We took into account the uncertainty in the parameters of the Milky Way potential as well as the uncertainty in the present day position and velocity of the Sun. We find that in general the Sun has not migrated from its birth place to its current position in the Galaxy (R_\odot). However, significant radial migration of the Sun is possible when: 1) The 2:1 Outer Lindblad resonance of the bar is separated from the corrotation resonance of spiral arms by a distance ~1 kpc. 2) When these two resonances are at the same Galactocentric position and further than the solar radius. In both cases the migration of the Sun is from outer regions of the Galactic disk to R_\odot, placing the Sun's birth radius at around 11 kpc. We find that in general it is unlikely that the Sun has migrated significantly from the inner regions of the Galactic disk to R_\odot.

We review the recent optimizations of gravitational $N$-body kernels for running them on graphics processing units (GPUs), on single hosts and massive parallel platforms. For each of the two main $N$-body techniques, direct summation and tree-codes, we discuss the optimization strategy, which is different for each algorithm. Because both the accuracy as well as the performance characteristics differ, hybridizing the two algorithms is essential when simulating a large $N$-body system with high-density structures containing few particles, and with low-density structures containing many particles. We demonstrate how this can be realized by splitting the underlying Hamiltonian, and we subsequently demonstrate the efficiency and accuracy of the hybrid code by simulating a group of 11 merging galaxies with massive black holes in the nuclei.

We investigate the effect of different initial virial temperatures, Q, on the dynamics of star clusters. We find that the virial temperature has a strong effect on many aspects of the resulting system, including among others: the fraction of bodies escaping from the system, the depth of the collapse of the system, and the strength of the mass segregation. These differences deem the practice of using "cold" initial conditions no longer a simple choice of convenience. The choice of initial virial temperature must be carefully considered as its impact on the remainder of the simulation can be profound. We discuss the pitfalls and aim to describe the general behavior of the collapse and the resultant system as a function of the virial temperature so that a well reasoned choice of initial virial temperature can be made. We make a correction to the previous theoretical estimate for the minimum radius, $R_{min}$, of the cluster at the deepest moment of collapse to include a Q dependency, $R_{min}\approx Q + N^{(-1/3)}$, where $N$ is the number of particles. We use our numerical results to infer more about the initial conditions of the young cluster R136. Based on our analysis, we find that R136 was likely formed with a rather cool, but not cold, initial virial temperature ($Q\approx 0.13$). Using the same analysis method, we examined 15 other young clusters and found the most common initial virial temperature to be between 0.18 and 0.25.

Aims: We study single and binary white dwarfs in the inner halo of the Milky Way in order to learn more about the conditions under which the population of halo stars was born, such as the initial mass function (IMF), the star formation history, or the binary fraction. Methods: We simulate the evolution of low-metallicity halo stars at distances up to ~ 3 kpc using the binary population synthesis code SeBa. We use two different white dwarf cooling models to predict the present-day luminosities of halo white dwarfs. We determine the white dwarf luminosity functions (WDLFs) for eight different halo models and compare these with the observed halo WDLF of white dwarfs in the SuperCOSMOS Sky Survey. Furthermore, we predict the properties of binary white dwarfs in the halo and determine the number of halo white dwarfs that is expected to be observed with the Gaia satellite. Results: By comparing the WDLFs, we find that a standard IMF matches the observations more accurately than a top-heavy one, but the difference with a bottom-heavy IMF is small. A burst of star formation 13 Gyr ago fits slightly better than a star formation burst 10 Gyr ago and also slightly better than continuous star formation $10-13$ Gyrs ago. Gaia will be the first instrument to constrain the bright end of the field halo WDLF, where contributions from binary WDs are considerable. Many of these will have He cores, of which a handful have atypical surface gravities ($\log g < 6$) and reach luminosities $\log(L/L_\odot) > 0$ in our standard model for WD cooling. These so called pre-WDs, if observed, can help us to constrain white dwarf cooling models and might teach us something about the fraction of halo stars that reside in binaries.

We show the effects of the perturbation caused by a passing by star on the Kuiper belt objects (KBOs) of our Solar System. The dynamics of the Kuiper belt (KB) is followed by direct $N$-body simulations. The sampling of the KB has been done with $N$ up to $131,062$, setting the KBOs on initially nearly circular orbits distributed in a ring of surface density $\Sigma \sim r^{-2}$. This modelization allowed us to investigate the secular evolution of the KB upon the encounter with the perturbing star. Actually, the encounter itself usually leads toward eccentricity and inclination distributions similar to observed ones, but tends also to excite the low-eccentricity population ($e < 0.1$ around $a\sim 40$\,$\mathrm{AU}$ from the Sun), depleting this region of low eccentricities. The following long-term evolution shows a "cooling" of the eccentricities repopulating the low-eccentricity area. In dependence on the assumed KBO mass spectrum and sampled number of bodies, this repopulation takes place in a time that goes from 0.5 Myr to 100 Myr. Due to the unavoidable limitation in the number of objects in our long-term simulations ($N \leq 16384$), we could not consider a detailed KBO mass spectrum, accounting for low mass objects, thus our present simulations are not reliable in constraining correlations among inclination distribution of the KBOs and other properties, such as their size distribution. However, our high precision long term simulations are a starting point for future larger studies on massively parallel computational platforms which will provide a deeper investigation of the secular evolution ($\sim 100\,$Myr) of the KB over its whole mass spectrum.

It is often assumed that the strong gravitational field of a super-massive black hole disrupts an adjacent molecular cloud preventing classical star formation in the deep potential well of the black hole. Yet, young stars have been observed across the entire nuclear star cluster of the Milky Way including the region close ($<$0.5~pc) to the central black hole, Sgr A*. Here, we focus particularly on small groups of young stars, such as IRS 13N located 0.1 pc away from Sgr A*, which is suggested to contain about five embedded massive young stellar objects ($<$1 Myr). We perform three dimensional hydrodynamical simulations to follow the evolution of molecular clumps orbiting about a $4\times10^6~M_{\odot}$ black hole, to constrain the formation and the physical conditions of such groups. The molecular clumps in our models assumed to be isothermal containing 100 $M_{\odot}$ in $<$0.2 pc radius. Such molecular clumps exist in the circumnuclear disk of the Galaxy. In our highly eccentrically orbiting clump, the strong orbital compression of the clump along the orbital radius vector and perpendicular to the orbital plane causes the gas densities to increase to values higher than the tidal density of Sgr A*, which are required for star formation. Additionally, we speculate that the infrared excess source G2/DSO approaching Sgr A* on a highly eccentric orbit could be associated with a dust enshrouded star that may have been formed recently through the mechanism supported by our models.

We present a novel method for efficient direct integration of gravitational N-body systems with a large variation in characteristic time scales. The method is based on a recursive and adaptive partitioning of the system based on the connected components of the graph generated by the particle distribution combined with an interaction-specific time step criterion. It uses an explicit and approximately time-symmetric time step criterion, and conserves linear and angular momentum to machine precision. In numerical tests on astrophysically relevant setups, the method compares favourably to both alternative Hamiltonian-splitting integrators as well as recently developed block time step-based GPU-accelerated Hermite codes. Our reference implementation is incorporated in the HUAYNO code, which is freely available as a part of the AMUSE framework.

General relativistic precession limits the ability of gravitational encounters to increase the eccentricity $e$ of orbits near a supermassive black hole (SBH). This "Schwarzschild barrier" (SB) has been shown to play an important role in the orbital evolution of stars like the galactic center S-stars. However, the evolution of orbits below the SB, $e>e_\mathrm{SB}$, is not well understood; the main current limitation is the computational complexity of detailed simulations. Here we present an $N$-body algorithm that allows us to efficiently integrate orbits of test stars around a SBH including general relativistic corrections to the equations of motion and interactions with a large ($\gtrsim 10^3$) number of field stars. We apply our algorithm to the S-stars and extract diffusion coefficients describing the evolution in angular momentum $L$. We identify three angular momentum regimes, in which the diffusion coefficients depend in functionally different ways on $L$. Regimes of lowest and highest $L$ are well-described in terms of non-resonant relaxation (NRR) and resonant relaxation (RR), respectively. In addition, we find a new regime of "anomalous relaxation" (AR). We present analytic expressions, in terms of physical parameters, that describe the diffusion coefficients in all three regimes, and propose a new, empirical criterion for the location of the SB in terms of the $L$-dependence of the diffusion coefficients. Subsequently we apply our results to obtain the steady-state distribution of angular momentum for orbits near a SBH.

Star clusters with multi-mass components dynamically evolve faster than those modeled with equal-mass components. Using a series of direct $N$-body simulations, we investigate the dynamical evolution of star clusters with mass functions, especially their core collapse time. Multi-mass clusters tend to behave like systems with a smaller number of particles, which we call the effective number of particles ($N_{eff}$) and for which $N_{eff} = M/m_{max}$ (here $M$ and $m_{max}$ are the total cluster mass and the mass of the most massive star in the cluster, respectively). We find that the time of core collapse is inversely proportional to the mass of the most massive star in the cluster and analytically confirm that this is because the core collapse of clusters with a mass function proceeds on the dynamical friction timescale of the most massive stars. As the mass of the most massive star increases, however, the core-collapse time, which is observed as a core bounce of the cluster core from the evolution of the core density or core radius, becomes ambiguous. We find that in that case the total binding energy of the hard binaries gives a good diagnosis for determining the moment of core collapses. Based on the results of our simulations, we argue that the core bounce becomes ambiguous when the mass of the most massive star exceeds 0.1\% of the total mass of the cluster.

Recent simulations and observations suggest that star clusters form via the assembling of smaller sub-clusters. Because of their short relaxation time, sub-clusters experience core collapse much earlier than virialized solo-clusters, which have similar properties of the merger remnant of the assembling clusters. As a consequence it seems that the assembling clusters result in efficient multiple collisions of stars in the cluster core. We performed a series of $N$-body simulations of ensemble and solitary clusters including stellar collisions and found that the efficiency of multiple collisions between stars are suppressed if sub-clusters assemble after they experience core collapse individually. In this case, sub-clusters form their own multiple collision stars which experienced a few collisions, but they fail to collide with each other after their host sub-clusters assemble. The multiple collision stars scatter each other and escape, and furthermore the central density of the remnant clusters had already been depleted for the stars to experience more collisions. On the other hand, if sub-clusters assemble before they experience core collapse, the multiple collisions of stars proceed efficiently in the remnant cluster, and the collision products are more massive than virialized solo-clusters and comparable in mass to cold solo-clusters.

The conservation of energy, linear momentum and angular momentum are important drivers for our physical understanding of the evolution of the Universe. These quantities are also conserved in Newton's laws of motion under gravity \citep{Newton:1687}. Numerical integration of the associated equations of motion is extremely challenging, in particular due to the steady growth of numerical errors (by round-off and discrete time-stepping, \cite{1981PAZh....7..752B,1993ApJ...415..715G,1993ApJ...402L..85H,1994LNP...430..131M}) and the exponential divergence \citep{1964ApJ...140..250M,2009MNRAS.392.1051U} between two nearby solution. As a result, numerical solutions to the general N-body problem are intrinsically questionable \citep{2003gmbp.book.....H,1994JAM....61..226L}. Using brute force integrations to arbitrary numerical precision we demonstrate empirically that ensembles of different realizations of resonant 3-body interactions produce statistically indistinguishable results. Although individual solutions using common integration methods are notoriously unreliable, we conjecture that an ensemble of approximate 3-body solutions accurately represents an ensemble of true solutions, so long as the energy during integration is conserved to better than 1/10. We therefore provide an independent confirmation that previous work on self-gravitating systems can actually be trusted, irrespective of the intrinsic chaotic nature of the N-body problem.

Astronomical phenomena are governed by processes on all spatial and temporal scales, ranging from days to the age of the Universe (13.8,Gyr) as well as from km size up to the size of the Universe. This enormous range in scales is contrived, but as long as there is a physical connection between the smallest and largest scales it is important to be able to resolve them all, and for the study of many astronomical phenomena this governance is present. Although covering all these scales is a challenge for numerical modelers, the most challenging aspect is the equally broad and complex range in physics, and the way in which these processes propagate through all scales. In our recent effort to cover all scales and all relevant physical processes on these scales we have designed the Astrophysics Multipurpose Software Environment (AMUSE). AMUSE is a Python-based framework with production quality community codes and provides a specialized environment to connect this plethora of solvers to a homogeneous problem solving environment.

We developed a Keplerian-based Hamiltonian splitting for solving the gravitational $N$-body problem. This splitting allows us to approximate the solution of a general $N$-body problem by a composition of multiple, independently evolved $2$-body problems. While the Hamiltonian splitting is exact, we show that the composition of independent $2$-body problems results in a non-symplectic non-time-symmetric first-order map. A time-symmetric second-order map is then constructed by composing this basic first-order map with its self-adjoint. The resulting method is precise for each individual $2$-body solution and produces quick and accurate results for near-Keplerian $N$-body systems, like planetary systems or a cluster of stars that orbit a supermassive black hole. The method is also suitable for integration of $N$-body systems with intrinsic hierarchies, like a star cluster with primordial binaries. The superposition of Kepler solutions for each pair of particles makes the method excellently suited for parallel computing; we achieve $\gtrsim 64\%$ efficiency for only $8$ particles per core, but close to perfect scaling for $16384$ particles on a $128$ core distributed-memory computer. We present several implementations in \texttt{Sakura}, one of which is publicly available via the AMUSE framework.

Black holes orbiting the Super Massive Black Hole (SMBH) Sgr A* in the Milky-way galaxy center (GC) generate gravitational waves. The spectrum, due to stars and black holes, is continuous below 40 nHz while individual BHs within about 200 AU of the central SMBH stick out in the spectrum at higher frequencies. The GWs can be detected by timing radio pulsars within a few parsecs of this region. Future observations with the Square Kilometer Array of such pulsars with sufficient timing accuracy may be sensitive to signals from intermediate mass BHs (IMBH) in a 3 year observation baseline. The recent detection of radio pulsations from the magnetar SGR J1745-29 very near the GC opens up the possibilities of detecting millisecond pulsars (which can be used as probes of the GWs) through lines of sight with only moderate pulse and angular broadening due to scattering.

We aim at modelling small groups of young stars such as IRS 13N, 0.1 pc away from Sgr A*, which is suggested to contain a few embedded massive young stellar objects. We perform hydrodynamical simulations to follow the evolution of molecular clumps orbiting about a $4\times10^6 ~ M_{\odot}$ black hole, to constrain the formation and the physical conditions of such groups. We find that, the strong compression due to the black hole along the orbital radius vector of clumps evolving on highly eccentric orbits causes the clumps densities to increase to higher than the tidal density of Sgr A*, and required for star formation. This suggests that the tidal compression from the black hole could support star formation. Additionally, we speculate that the infrared excess source G2/DSO approaching Sgr A* on a highly eccentric orbit could be associated with a dust enshrouded star that may have been formed recently through the mechanism supported by our models.

Super Luminous supernovae (SLSN) occur almost exclusively in small galaxies (SMC/LMC-like or smaller), and the few SLSN observed in larger star-forming galaxies always occur close to the nuclei of their hosts. Another type of peculiar and highly energetic supernovae are the broad-line type Ic SNe (SN Ic-BL) that are associated with long-duration gamma-ray bursts (LGRBs). Also these have a strong preference for occurring in small (SMC/LMC-like or smaller) star-forming galaxies, and in these galaxies LGRBs always occur in the brightest spots. Studies of nearby star-forming galaxies that are similar to the hosts of LGRBs show that these brightest spots are giant HII regions produced by massive dense young star clusters with many hundreds of O- and Wolf-Rayet-type stars. Such dense young clusters are also found in abundance within a few hundred parsecs from the nucleus of larger galaxies like our own. We argue that the SLSN and the SN Ic-BL/LGRBs are exclusive products of two types of dynamical interactions in dense young star clusters. In our model the high angular momentum of the collapsing stellar cores required for the "engines" of a SN Ic-BL results from the post-main sequence mergers of dynamically produced cluster binaries with almost equal-mass components. The merger produces a critically rotating single helium star with sufficent angular momentum to produce a LGRB; the observed "metal aversio" of LGRBs is a natural consequence of the model. We argue that, on the other hand, SLSN could be the products of runaway multiple collisions in dense clusters, and we present (and quantize) plausible scenarios of how the different types of SLSNs can be produced.

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.

The evolution of triples has not attracted much attention in the literature, although their evolution can be dramatically different from binaries and single stars. Triples are quite common, and we find that for about 1% of the triples in the Tokovinin catalogue of multiple stellar systems in the solar neighbourhood, the tertiary star will overflow its Roche lobe at some time in its evolution, before any of the inner stars leave the main sequence. For two of these systems, Xi Tauri and HD97131 we simulate in detail this phase of mass transfer, during which stellar evolution, gravitational dynamics and hydrodynamics all play an important role. We have used the Astrophysical Multi-purpose Software Environment (AMUSE) to solve these physical processes in a self-consistent way. The resulting evolution, mass transfer and the effects on the inner as well as on the outer orbit are profound, although it is not trivial to predict the eventual consequence of the phase of mass transfer and the appearance of the resulting system.

We study the formation and evolution of filamentary configurations of dark matter haloes in voids. Our investigation uses the high-resolution LambdaCDM simulation CosmoGrid to look for void systems resembling the VGS_31 elongated system of three interacting galaxies that was recently discovered by the Void Galaxy Survey (VGS) inside a large void in the SDSS galaxy redshift survey. HI data revealed these galaxies to be embedded in a common elongated envelope, possibly embedded in intravoid filament. In the CosmoGrid simulation we look for systems similar to VGS_31 in mass, size and environment. We find a total of eight such systems. For these systems, we study the distribution of neighbour haloes, the assembly and evolution of the main haloes and the dynamical evolution of the haloes, as well as the evolution of the large-scale structure in which the systems are embedded. The spatial distribution of the haloes follows that of the dark matter environment. We find that VGS_31-like systems have a large variation in formation time, having formed between 10 Gyr ago and the present epoch. However, the environments in which the systems are embedded evolved resemble each other substantially. Each of the VGS_31-like systems is embedded in an intra-void wall, that no later than z = 0.5 became the only prominent feature in its environment. While part of the void walls retain a rather featureless character, we find that around half of them are marked by a pronounced and rapidly evolving substructure. Five haloes find themselves in a tenuous filament of a few Mpc/h long inside the intra-void wall. Finally, we compare the results to observed data from VGS_31. Our study implies that the VGS_31 galaxies formed in the same (proto)filament, and did not meet just recently. The diversity amongst the simulated halo systems indicates that VGS_31 may not be typical for groups of galaxies in voids.

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.

In young dense clusters repeated collisions between massive stars may lead to the formation of a very massive star (above 100 Msun). In the past the study of the long-term evolution of merger remnants has mostly focussed on collisions between low-mass stars (up to about 2 Msun) in the context of blue-straggler formation. The evolution of collision products of more massive stars has not been as thoroughly investigated. In this paper we study the long-term evolution of a number of stellar mergers formed by the head-on collision of a primary star with a mass of 5-40 Msun with a lower mass star at three points in its evolution in order to better understand their evolution. We use smooth particle hydrodynamics (SPH) calculations to model the collision between the stars. The outcome of this calculation is reduced to one dimension and imported into a stellar evolution code. We follow the subsequent evolution of the collision product through the main sequence at least until the onset of helium burning. We find that little hydrogen is mixed into the core of the collision products, in agreement with previous studies of collisions between low-mass stars. For collisions involving evolved stars we find that during the merger the surface nitrogen abundance can be strongly enhanced. The evolution of most of the collision products proceeds analogously to that of normal stars with the same mass, but with a larger radius and luminosity. However, the evolution of collision products that form with a hydrogen depleted core is markedly different from that of normal stars with the same mass. They undergo a long-lived period of hydrogen shell burning close to the main-sequence band in the Hertzsprung-Russell diagram and spend the initial part of core helium burning as compact blue supergiants.

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.

[abridged] Aims. We model the number and properties of ultracompact X-ray binaries (UCXBs) in the Galactic Bulge. The objective is to learn about the formation of UCXBs and their evolution, such as the onset of mass transfer and late-time behavior. Methods. The binary population synthesis code SeBa and stellar evolutionary tracks are used to model the UCXB population. The luminosity behavior of UCXBs is predicted using long-term X-ray observations of the known UCXBs and the thermal-viscous disk instability model. Results. In our model, the majority of UCXBs initially have a helium burning star donor. In the absence of a mechanism that destroys old UCXBs, we predict (0.2 - 1.9) x 10^5 UCXBs in the Galactic Bulge, mostly at orbital periods longer than 60 min (a large number of long-period systems also follows from the observed short-period UCXB population). About 5 - 50 UCXBs should be brighter than 10^35 erg/s, mostly persistent sources with orbital periods shorter than 30 min and with degenerate helium and carbon-oxygen donors. This is one order of magnitude more than the observed number of (probably) three. Conclusions. This overprediction of short-period UCXBs implies that fewer systems are formed, or that a super-Eddington mass transfer rate is more difficult to survive than assumed. The very small number of observed long-period UCXBs with respect to short-period UCXBs, the surprisingly high luminosity of the observed UCXBs with orbital periods around 50 min, and the properties of the PSR J1719-1438 system all point to much faster UCXB evolution than expected from angular momentum loss via gravitational wave radiation alone. Old UCXBs, if they still exist, probably have orbital periods longer than 2 h and have become very faint due to either reduced accretion or quiescence, or have become detached. UCXBs are promising candidate progenitors of isolated millisecond radio pulsars.

Despite the significance of type Ia supernovae (SNeIa) in many fields in astrophysics, SNeIa lack a theoretical explanation. We investigate the potential contribution to the SNeIa rate from the most common progenitor channels using the binary population synthesis (BPS) code SeBa. Using SeBa, we aim constrain binary processes such as the common envelope phase and the efficiency of mass retention of white dwarf accretion. We find that the simulated rates are not sufficient to explain the observed rates. Further, we find that the mass retention efficiency of white dwarf accretion significantly influences the rates, but does not explain all the differences between simulated rates from different BPS codes.

Despite the significance of Type Ia supernovae (SNeIa) in many fields in astrophysics, SNeIa lack a theoretical explanation. The standard scenarios involve thermonuclear explosions of carbon/oxygen white dwarfs approaching the Chandrasekhar mass; either by accretion from a companion or by a merger of two white dwarfs. We investigate the contribution from both channels to the SNIa rate with the binary population synthesis (BPS) code SeBa in order to constrain binary processes such as the mass retention efficiency of WD accretion and common envelope evolution. We determine the theoretical rates and delay time distribution of SNIa progenitors and in particular study how assumptions affect the predicted rates.

Massive galaxies with a half-mass radius <~ 1kpc are observed in the early universe (z~>2), but not in the local universe. In the local universe similar-mass (within a factor of two) galaxies tend to be a factor of 4 to 5 larger. Dry minor mergers are known to drive the evolution of the size of a galaxy without much increasing the mass, but it is unclear if the growth in size is sufficient to explain the observations. We test the hypothesis that galaxies grow through dry minor mergers by simulating merging galaxies with mass ratios of q=1:1 (equal mass) to q=1:160. In our N-body simulations the total mass of the parent galaxy doubles. We confirm that major mergers do not cause a sufficient growth in size. The observation can be explained with mergers with a mass ratio of q=1:5--1:10. Smaller mass ratios cause a more dramatic growth in size, up to a factor of ~17 for mergers with a mass ratio of 1:80. For relatively massive minor mergers q ~> 1:20 the mass of the incoming child galaxies tend to settle in the halo of the parent galaxy. This is caused by the tidal stripping of the child galaxies by the time they enter the central portion of the parent. When the accretion of minor galaxies becomes more continuous, when q <~ 1:40, the foreign mass tends to concentrate more in the central region of the parent galaxy. We speculate that this is caused by dynamic interactions between the child galaxies inside the merger remnant and the longer merging times when the difference in mass is larger. These interactions cause dynamical heating which results in accretion of mass inside the galaxy core and a reduction of the parent's circular velocity and density.

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.

Cataclysmic variables (CVs) are binaries in which a compact white dwarf accretes material from a low-mass companion star. The discovery of two planets in orbit around the CV HU Aquarii opens unusual opportunities for understanding the formation and evolution of this system. In particular the orbital parameters of the planets constrains the past and enables us to reconstruct the evolution of the system through the common-envelope phase. During this dramatic event the entire hydrogen envelope of the primary star is ejected, passing the two planets on the way. The observed eccentricities and orbital separations of the planets in HU Aqr enable us to limit the common-envelope parameter $\alpha \lambda = 0.45\pm 0.17$ or $\gamma = 1.77\pm0.02$ and measure the rate at which the common envelope is ejected, which turns out to be copious. The mass in the common envelope is ejected from the binary system at a rate of ${\dot m} = 1.9\pm 0.3\,\MSun/yr$. The reconstruction of the initial conditions for HU Aqr indicates that the primary star had a mass of $M_{\rm ZAMS} = 1.6\pm0.2$\,\MSun\, and a $m_{\rm ZAMS} = 0.47\pm 0.04$\,\MSun\, companion in a $a=25$--160\,\RSun\, (best value $a=97$\,\RSun) binary. The two planets were born with an orbital separation of $a_a=541\pm44$\,\RSun\, and $a_b=750\pm72$\,\Rsun\, respectively. After the common envelope, the primary star turns into a $0.52\pm0.01$\,\MSun\, helium white dwarf, which subsequently accreted $\sim 0.30$\,\MSun\, from its Roche-lobe filling companion star, grinding it down to its current observed mass of $0.18\,\MSun$.

We examine the formation of planets around binary stars in light of the recently discovered systems Kepler 16, 34 and 35. We conduct hydrodynamical simulations of self gravitating disks around binary systems. The selected binary and disk parameters are chosen consistent with observed systems. The disks are evolved until they settle in a quasi-equilibrium and the resulting systems are compared with the parameters of Kepler 16, 34 and 35. We find a close correspondence of the peak density at the inner disk gap and the orbit of the observed planets. We conclude, based on our simulations, that the orbits of the observed Kepler planets are determined by the size of the inner disk gap which for these systems results from the binary driving. This mediates planet formation either through the density enhancement or through planetary trapping at the density gradient inversion in the inner disk. For all three systems the current eccentricity of the planetary orbit is less than the disk eccentricity in the simulations. This, together with the long term stability of the orbits argues against in situ formation (e.g. a direct collapse scenario of the material in the ring). Conducting additional simulations of systems with a wider range of parameters (taken from a survey of eclipsing binaries), we find that the planet semi-major axis and binary eccentricity in such a scenario should be tightly correlated providing an observational test of this formation mechanism.

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.

The study of Type Ia supernovae (SNIa) has lead to greatly improved insights into many fields in astrophysics, however a theoretical explanation of the origin of these events is still lacking. We investigate the potential contribution to the SNIa rate from the population of merging double carbon-oxygen white dwarfs. We aim to develope a model that fits the observed SNIa progenitors as well as the observed close double white dwarf population. We differentiate between two scenarios for the common envelope (CE) evolution; the alpha-formalism based on the energy equation and the gamma-formalism that is based on the angular momentum equation. In one model we apply the alpha-formalism always. In the second model the gamma-formalism is applied, unless the binary contains a compact object or the CE is triggered by a tidal instability for which the alpha-formalism is used. The binary population synthesis code SeBa was used to evolve binary systems from the zero-age main sequence to the formation of double white dwarfs and subsequent mergers. SeBa has been thoroughly updated since the last publication of the content of the code. The limited sample of observed double white dwarfs is better represented by the simulated population using the gamma-formalism than the alpha-formalism. For both CE formalisms, we find that although the morphology of the simulated delay time distribution matches that of the observations within the errors, the normalisation and time-integrated rate per stellar mass are a factor 7-12 lower than observed. Furthermore, the characteristics of the simulated populations of merging double carbon-oxygen white dwarfs are discussed and put in the context of alternative SNIa models for merging double white dwarfs.

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.

We develop a method to analyze the effect of an asymmetric supernova on hierarchical multiple star systems and we present analytical formulas to calculate orbital parameters for surviving binaries or hierarchical triples and runaway velocities for their dissociating equivalents. The effect of an asymmetric supernova on the orbital parameters of a binary system has been studied to great extent (e.g. Hills 1983; Kalogera 1996; Tauris & Takens 1998), but this effect on higher multiplicity hierarchical systems has not been explored before. With our method, the supernova effect can be computed by reducing the hierarchical multiple to an effective binary by means of recursively replacing the inner binary by an effective star at the center of mass of that binary. We apply our method to a hierarchical triple system similar to the progenitor of PSR J1903+0327 suggested by Portegies Zwart et al. (2011). We confirm their earlier finding that PSR J1903+0327 could have evolved from a hierarchical triple that became unstable and ejected the secondary star of the inner binary. Furthermore, if such as system did evolve via this mechanism the most probable configuration would be a small supernova kick velocity, an inner binary with a large semi-major axis, and the fraction of mass accreted onto the neutron star to the mass lost by the secondary would most likely be between 0.35 and 0.5

We review the implementation of individual particle time-stepping for N-body dynamics. We present a class of integrators derived from second order Hamiltonian splitting. In contrast to the usual implementation of individual time-stepping, these integrators are momentum conserving and show excellent energy conservation in conjunction with a symmetrized time step criterion. We use an explicit but approximate formula for the time symmetrization that is compatible with the use of individual time steps. No iterative scheme is necessary. We implement these ideas in the HUAYNO (available online at www.amusecode.org) code and present tests of the integrators and show that the presented integration schemes shows good energy conservation, with little or no systematic drift, while conserving momentum and angular momentum to machine precision for long term integrations.

Young star clusters like R136 in the Large Magellanic Cloud and NGC 3603, Westerlund 1, and 2 in the Milky Way are dynamically more evolved than expected based on their current relaxation times. In particular, the combination of a high degree of mass segregation, a relatively low central density, and the large number of massive runaway stars in their vicinity are hard to explain with the monolithic formation of these clusters. Young star clusters can achieve such a mature dynamical state if they formed through the mergers of a number of less massive clusters. The shorter relaxation times of less massive clusters cause them to dynamically evolve further by the time they merge, and the merger product preserves the memory of the dynamical evolution of its constituent clusters. With a series of $N$-body simulations, we study the dynamical evolution of single massive clusters and those that are assembled through merging smaller clusters together. We find that the formation of massive star clusters through the mergers of smaller clusters can reproduce the currently observed spatial distribution of massive stars, the density, and the characteristics (number and mass distribution) of the stars ejected as runaways from young dense clusters. We therefore conclude that these clusters and possibly other young massive star clusters formed through the mergers of smaller clusters.

In this short review we present the developments over the last 5 decades that have led to the use of Graphics Processing Units (GPUs) for astrophysical simulations. Since the introduction of NVIDIA's Compute Unified Device Architecture (CUDA) in 2007 the GPU has become a valuable tool for N-body simulations and is so popular these days that almost all papers about high precision N-body simulations use methods that are accelerated by GPUs. With the GPU hardware becoming more advanced and being used for more advanced algorithms like gravitational tree-codes we see a bright future for GPU like hardware in computational astrophysics.

We present parallel algorithms for constructing and traversing sparse octrees on graphics processing units (GPUs). The algorithms are based on parallel-scan and sort methods. To test the performance and feasibility, we implemented them in CUDA in the form of a gravitational tree-code which completely runs on the GPU.(The code is publicly available at: http://castle.strw.leidenuniv.nl/software.html) The tree construction and traverse algorithms are portable to many-core devices which have support for CUDA or OpenCL programming languages. The gravitational tree-code outperforms tuned CPU code during the tree-construction and shows a performance improvement of more than a factor 20 overall, resulting in a processing rate of more than 2.8 million particles per second.

We present a gravitational hierarchical N-body code that is designed to run efficiently on Graphics Processing Units (GPUs). All parts of the algorithm are executed on the GPU which eliminates the need for data transfer between the Central Processing Unit (CPU) and the GPU. Our tests indicate that the gravitational tree-code outperforms tuned CPU code for all parts of the algorithm and show an overall performance improvement of more than a factor 20, resulting in a processing rate of more than 2.8 million particles per second.

We examine the nHz gravitational wave (GW) foreground of stars and black holes (BHs) orbiting SgrA* in the Galactic Center. A cusp of stars and BHs generates a continuous GW spectrum below 40 nHz; individual BHs within 1 mpc to SgrA* stick out in the spectrum at higher GW frequencies. The GWs and gravitational near-field effects can be resolved by timing pulsars within a few pc of this region. Observations with the Square Kilometer Array (SKA) may be especially sensitive to intermediate mass black holes (IMBHs) in this region, if present. A 100ns-10mus timing accuracy is sufficient to detect BHs of mass 1000 Msun with pulsars at distance 0.1-1 pc in a 3 yr observation baseline. Unlike electromagnetic imaging techniques, the prospects for resolving individual objects through GW measurements improve closer to SgrA*, even if the number density of objects increases inwards steeply. Scattering by the interstellar medium will pose the biggest challenge for such observations.

High-performance scientific applications require more and more compute power. The concurrent use of multiple distributed compute resources is vital for making scientific progress. The resulting distributed system, a so-called Jungle Computing System, is both highly heterogeneous and hierarchical, potentially consisting of grids, clouds, stand-alone machines, clusters, desktop grids, mobile devices, and supercomputers, possibly with accelerators such as GPUs. One striking example of applications that can benefit greatly of Jungle Computing Systems are Multi-Model / Multi-Kernel simulations. In these simulations, multiple models, possibly implemented using different techniques and programming models, are coupled into a single simulation of a physical system. Examples include the domain of computational astrophysics and climate modeling. In this paper we investigate the use of Jungle Computing Systems for such Multi-Model / Multi-Kernel simulations. We make use of the software developed in the Ibis project, which addresses many of the problems faced when running applications on Jungle Computing Systems. We create a prototype Jungle-aware version of AMUSE, an astrophysical simulation framework. We show preliminary experiments with the resulting system, using clusters, grids, stand-alone machines, and GPUs.

Most recent progress in understanding the dynamical evolution of star clusters relies on direct N-body simulations. Owing to the computational demands, and the desire to model more complex and more massive star clusters, hardware calculational accelerators, such as GRAPE special-purpose hardware or, more recently, GPUs (i.e. graphics cards), are generally utilised. In addition, simulations can be accelerated by adjusting parameters determining the calculation accuracy (i.e. changing the internal simulation time step used for each star). We extend our previous thorough comparison (Anders et al. 2009) of basic quantities as derived from simulations performed either with STARLAB/KIRA or NBODY6. Here we focus on differences arising from using different hardware accelerations (including the increasingly popular graphic card accelerations/GPUs) and different calculation accuracy settings. We use the large number of star cluster models (for a fixed stellar mass function, without stellar/binary evolution, primordial binaries, external tidal fields etc) already used in the previous paper, evolve them with STARLAB/KIRA (and NBODY6, where required), analyse them in a consistent way and compare the averaged results quantitatively. For this quantitative comparison, we apply the bootstrap algorithm for functional dependencies developed in our previous study. In general we find very high comparability of the simulation results, independent of the used computer hardware (including the hardware accelerators) and the used N-body code. For the tested accuracy settings we find that for reduced accuracy (i.e. time step at least a factor 2.5 larger than the standard setting) most simulation results deviate significantly from the results using standard settings. The remaining deviations are comprehensible and explicable.

We present results of simulations aimed at tracing the formation of nuclear star clusters (NCs) and black hole (BH) seeds, in a cosmological context. We focus on two mechanisms for the formation of BHs at high redshifts: as end-products of (1) Population III stars in metal free halos, and of (2) runaway stellar collisions in metal poor NCs. Our model tracks the chemical, radiative and mechanical feedback of stars on the baryonic component of the evolving halos. This procedure allows us to evaluate when and where the conditions for BH formation are met, and to trace the emergence of BH seeds arising from the dynamical channel, in a cosmological context. BHs start to appear already at z~30 as remnants of Population III stars. The efficiency of this mechanism begins decreasing once feedbacks become increasingly important. Around redshift z~15, BHs mostly form in the centre of mildly metal enriched halos inside dense NCs. The seed BHs that form along the two pathways have at birth a mass around 100-1000M\odot. The occupation fraction of BHs is a function of both halo mass and mass growth rate: at a given z, heavier and faster growing halos have a higher chance to form a native BH, or to acquire an inherited BH via merging of another system. With decreasing z, the probability of finding a BH shifts toward progressively higher mass halo intervals. This is due to the fact that, at later cosmic times, low mass systems rarely form a seed, and already formed BHs are deposited into larger mass systems due to hierarchical mergers. Our model predict that at z=0, all halos above 10^11M\odot should host a BH (in agreement with observational results), most probably inherited during their lifetime. Halos less massive then 10^9M\odot have a higher probability to host a native BH, but their occupation fraction decreases below 10%.

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.

(Abridged) Interacting galaxies are well-known for their high star formation rates and rich star cluster populations, but the rapidly changing tidal field can also efficiently destroy clusters. We use numerical simulations of merging disc galaxies to investigate which mechanism dominates. The simulations include a model for the formation and dynamical disruption of the entire star cluster population. We find that the dynamical heating of clusters by tidal shocks is about an order of magnitude higher in interacting galaxies than in isolated galaxies. This is driven by the increased gas density, and is sufficient to destroy star clusters at a higher rate than new clusters are formed: the total number of clusters in the merger remnant is 2-50% of the amount in the progenitor discs, with low-mass clusters being disrupted preferentially. By adopting observationally motivated selection criteria, we find that the observed surplus of star clusters in nearby merging galaxies is caused by the bias to detect young, massive clusters. We provide a general expression for the survival fraction of clusters, which increases with the gas depletion time-scale. Due to the preferential disruption of low-mass clusters, the mass distribution of the surviving star clusters in a merger remnant develops a peak at a mass of about 10^3 Msun, which evolves to higher masses at a rate of 0.3-0.4 dex per Gyr. The peak mass initially depends weakly on the galactocentric radius, but this correlation disappears as the system ages. We discuss the similarities between the cluster populations of the simulated merger remnants and (young) globular cluster systems. Our results suggest that the combination of cluster formation and destruction should be widespread in the dense star-forming environments at high redshifts, which could provide a natural origin to present-day globular cluster systems.

We compute the gravitational wave background (GWB) generated by a cosmological population of (BH-BH) binaries using hybrid waveforms recently produced by numerical simulations of (BH-BH) coalescence, which include the inspiral, merger and ring-down contributions. A large sample of binary systems is simulated using the population synthesis code SeBa, and we extract fundamental statistical information on (BH-BH) physical parameters (primary and secondary BH masses, orbital separations and eccentricities, formation and merger timescales). We then derive the binary birth and merger rates using the theoretical cosmic star formation history obtained from a numerical study which reproduces the available observational data at redshifts $z < 8$. We evaluate the contributions of the inspiral, merger and ring-down signals to the GWB, and discuss how these depend on the parameters which critically affect the number of coalescing (BH-BH) systems. We find that Advanced LIGO/Virgo have a chance to detect the GWB signal from the inspiral phase with a $(S/N)=10$ only for the most optimistic model, which predicts the highest local merger rate of 0.85 Mpc$^{-3}$ Myr$^{-1}$. Third generation detectors, such as ET, could reveal the GWB from the inspiral phase predicted by any of the considered models. In addition, ET could sample the merger phase of the evolution at least for models which predict local merger rates between $[0.053 - 0.85]$ Mpc$^{-3}$ Myr$^{-1}$, which are more than a factor 2 lower the the upper limit inferred from the analysis of the LIGO S5 run\cite{Abadieetal2011}.

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.

About 20% of all massive stars in the Milky Way have unusually high velocities, the origin of which has puzzled astronomers for half a century. We argue that these velocities originate from strong gravitational interactions between single stars and binaries in the centers of star clusters. The ejecting binary forms naturally during the collapse of a young ($\aplt 1$\,Myr) star cluster. This model replicates the key characteristics of OB runaways in our galaxy and it explains the $\apgt 100$\,\Msun\, runaway stars around young star clusters, e.g. R136 and Westerlund~2. The high proportion and the distributions in mass and velocity of runaways in the Milky Way is reproduced if the majority of massive stars are born in dense and relatively low-mass (5000-10000 \Msun) clusters.

We study the evolution of embedded clusters. The equations of motion of the stars in the cluster are solved by direct N-body integration while taking the effects of stellar evolution and the hydrodynamics of the natal gas content into account. The gravity of the stars and the surrounding gas are coupled self consistently to allow the realistic dynamical evolution of the cluster. While the equations of motion are solved, a stellar evolution code keeps track of the changes in stellar mass, luminosity and radius. The gas liberated by the stellar winds and supernovae deposits mass and energy into the gas reservoir in which the cluster is embedded. We examine cluster models with 1000 stars, but we varied the star formation efficiency (between 0.05-0.5), cluster radius (0.1-1.0 parsec), the degree of virial support of the initial population of stars (0-100%) and the strength of the feedback. We find that an initial star fraction $M_\star/M_{\rm tot} > 0.05$ is necessary for cluster survival. Survival is more likely if gas is not blown out violently by a supernova and if the cluster has time to approach virial equilibrium during out-gassing. While the cluster is embedded, dynamical friction drives early and efficient mass segregation in the cluster. Stars of $m \gtrsim 2\,M_\odot$ are preferentially retained, at the cost of the loss of less massive stars. We conclude that the degree of mass segregation in open clusters such as the Pleiades is not the result of secular evolution but a remnant of its embedded stage.

This paper has been withdrawn by the authors because it had major revisions.

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.

We present results from our cosmological N-body simulation which consisted of 2048x2048x2048 particles and ran distributed across three supercomputers throughout Europe. The run, which was performed as the concluding phase of the Gravitational Billion Body Problem DEISA project, integrated a 30 Mpc box of dark matter using an optimized Tree/Particle Mesh N-body integrator. We ran the simulation up to the present day (z=0), and obtained an efficiency of about 0.93 over 2048 cores compared to a single supercomputer run. In addition, we share our experiences on using multiple supercomputers for high performance computing and provide several recommendations for future projects.

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.

The millisecond pulsar J1903+0327 is accompanied by an ordinary G-dwarf star in an unusually wide ($P_{\rm orb} \simeq 95.2$\,days) and eccentric ($e \simeq 0.44$) orbit. The standard model for producing MSPs fails to explain the orbital characteristics of this extraordinary binary, and alternative binary models are unable to explain the observables. We present a triple-star model for producing MSPs in relatively wide eccentric binaries with a normal (main-sequence) stellar companion. We start from a stable triple system consisting of a Low-Mass X-ray Binary (LMXB) with an orbital period of at least 1 day, accompanied by a G-dwarf in a wide and possibly eccentric orbit. Variations in the initial conditions naturally provide a satisfactory explanation for the unexplained triple component in the eclipsing soft X-ray transient 4U~2129+47 or the cataclysmic variable EC 19314-5915. The best explanation for J1903, however, results from the expansion of the orbit of the LMXB, driven by the mass transfer from the evolving donor star to its neutron star companion, which causes the triple eventually to becomes dynamically unstable. Using numerical computations we show that, depending on the precise system configuration at the moment the triple becomes dynamically unstable, the ejection of each of the three components is possible. If the donor star of the LMXB is ejected, a system resembling J1903, will result. If the neutron star is ejected, a single MSP results. This model therefore also provides a straightforward mechanism for forming single MSP in the Galactic disk. We conclude that the Galaxy contains some 30--300 binaries with characteristics similar to J1903, and about an order of magnitude fewer single millisecond pulsars produced with the proposed triple scenario.

An update is presented to the software package SeBa for simulating single star and binary evolution in which new stellar evolution tracks have been implemented. SeBa is applied to study the population of close double white dwarf and the delay time distribution of double white dwarf mergers that may lead to Supernovae Type Ia.

(Abridged) The formation and evolution of star cluster populations are related to the galactic environment. Cluster formation is governed by processes acting on galactic scales, and star cluster disruption is driven by the tidal field. In this paper, we present a self-consistent model for the formation and evolution of star cluster populations, for which we combine an N-body/SPH galaxy evolution code with semi-analytic models for star cluster evolution. The model includes star formation, feedback, stellar evolution, and star cluster disruption by two-body relaxation and tidal shocks. We apply the model by simulating a suite of 9 isolated disc galaxies and 24 galaxy mergers. The evolutionary histories of individual clusters in these simulations are discussed to illustrate how the environment of clusters changes in time and space. The resulting variability of the disruption rate with time and space affects the properties of star cluster populations. The combined effect of clusters escaping their dense formation sites (`cluster migration') and the preferential disruption of clusters residing in dense environments (`natural selection') implies that the mean disruption rate of the population decreases with cluster age. This affects the slope of the cluster age distribution, which becomes a function of the star formation rate density. The evolutionary histories of clusters in a galaxy merger vary widely and determine which clusters survive the merger. This impacts the age distributions and the locations of the surviving clusters at all times during a merger. We conclude that accounting for the interplay between the formation, disruption, and orbital histories of clusters enables a more sophisticated interpretation of observed properties of cluster populations, thereby extending the role of cluster populations as tracers of galaxy evolution.

We report on the performance of our cold-dark matter cosmological N-body simulation which was carried out concurrently using supercomputers across the globe. We ran simulations on 60 to 750 cores distributed over a variety of supercomputers in Amsterdam (the Netherlands, Europe), in Tokyo (Japan, Asia), Edinburgh (UK, Europe) and Espoo (Finland, Europe). Regardless the network latency of 0.32 seconds and the communication over 30.000 km of optical network cable we are able to achieve about 87% of the performance compared to an equal number of cores on a single supercomputer. We argue that using widely distributed supercomputers in order to acquire more compute power is technically feasible, and that the largest obstacle is introduced by local scheduling and reservation policies.

In Galactic studies a distinction is made between (open) star clusters and associations. For barely resolved objects at a distance of several Mpc this distinction is not trivial to make. Here we provide an objective definition by comparing the age of the stars to the crossing time of nearby stellar agglomerates. We find that a satisfactory separation can be made where this ratio equals unity. Stellar agglomerates for which the age of the stars exceeds the crossing time are bound, and are referred to as star clusters. Alternatively, those for which the crossing time exceeds the stellar age are unbound and are referred to as associations. This definition is useful whenever reliable measurements for the mass, radius and age are available.

We present a new very fast tree-code which runs on massively parallel Graphical Processing Units (GPU) with NVIDIA CUDA architecture. The tree-construction and calculation of multipole moments is carried out on the host CPU, while the force calculation which consists of tree walks and evaluation of interaction list is carried out on the GPU. In this way we achieve a sustained performance of about 100GFLOP/s and data transfer rates of about 50GB/s. It takes about a second to compute forces on a million particles with an opening angle of $\theta \approx 0.5$. The code has a convenient user interface and is freely available for use\footnote{{\tt http://castle.strw.leidenuniv.nl/software/octgrav.html}}.

We show that we can obtain a good fit to the present day stellar mass functions (MFs) of a large sample of young and old Galactic clusters in the range 0.1 - 10 Msolar with a tapered power law distribution function with an exponential truncation of the form dN/dm \propto m^alpha [1 - exp-(m/m_c)^beta]. The average value of the power-law index alpha is -2, that of beta is 2.5, whereas the characteristic mass m_c is in the range 0.1 - 0.8 Msolar and does not seem to vary in any systematic way with the present cluster parameters such as metal abundance, total cluster mass or central concentration. However, m_c shows a remarkable correlation with the dynamical age of the cluster, namely m_c/Msolar ~ 0.15 + 0.5 tau_dyn^0.75, where tau_dyn is the dynamical age taken as the ratio of cluster age and dissolution time. The small scatter seen around this correlation is consistent with the uncertainties on the estimated value of tau_dyn. We attribute the observed trend to the onset of mass segregation via two-body relaxation in a tidal environment, causing the preferential loss of low-mass stars from the cluster and hence a drift of the characteristic mass m_c towards higher values. If dynamical evolution is indeed at the origin of the observed trend, it would seem plausible that high-concentration globular clusters, now with median m_c ~ 0.33 Msolar, were born with a stellar MF very similar to that measured today in the youngest Galactic clusters and with a value of m_c ~ 0.15 Msolar. This hypothesis is consistent with the absence of a turn-over in the MF of the Galactic bulge down to the observational limit at ~0.2 Msolar and, if correct, it would carry the implication that the characteristic mass is not set by the thermal Jeans mass of the cloud.

Observations and theory both suggest that star clusters form sub-virial (cool) with highly sub-structured distributions. We perform a large ensemble of N-body simulations of moderate-sized (N=1000) cool, fractal clusters to investigate their early dynamical evolution. We find that cool, clumpy clusters dynamically mass segregate on a short timescale, that Trapezium-like massive higher-order multiples are commonly formed, and that massive stars are often ejected from clusters with velocities > 10 km/s (c.f. the average escape velocity of 2.5 km/s). The properties of clusters also change rapidly on very short timescales. Young clusters may also undergo core collapse events, in which a dense core containing massive stars is hardened due to energy losses to a halo of lower-mass stars. Such events can blow young clusters apart with no need for gas expulsion. The warmer and less substructured a cluster is initially, the less extreme its evolution.

We describe the results of a search for the remnants of the Sun's birth cluster among stars in the Hipparcos Catalogue. This search is based on the predicted phase space distribution of the Sun's siblings from simple simulations of the orbits of the cluster stars in a smooth Galactic potential. For stars within 100 pc the simulations show that it is interesting to examine those that have small space motions relative to the Sun. From amongst the candidate siblings thus selected there are six stars with ages consistent with that of the Sun. Considering their radial velocities and abundances only one potential candidate, HIP 21158, remains but essentially the result of the search is negative. This is consistent with predictions by Portegies Zwart (2009) on the number of siblings near the Sun. We discuss the steps that should be taken in anticipation of the data from the Gaia mission in order to conduct fruitful searches for the Sun's siblings in the future.

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.

Many young extra-galactic clusters have a measured velocity dispersion that is too high for the mass derived from their age and total luminosity, which has led to the suggestion that they are not in virial equilibrium. Most of these clusters are confined to a narrow age range centred around 10 Myr because of observational constraints. At this age the cluster light is dominated by luminous evolved stars, such as red supergiants, with initial masses of ~13-22 Msun for which (primordial) binarity is high. In this study we investigate to what extent the observed excess velocity dispersion is the result of the orbital motions of binaries. We demonstrate that estimates for the dynamical mass of young star clusters, derived from the observed velocity dispersion, exceed the photometric mass by up-to a factor of 10 and are consistent with a constant offset in the square of the velocity dispersion. This can be reproduced by models of virialised star clusters hosting a massive star population of which ~25 is in binaries, with typical mass ratios of ~0.6 and periods of ~1000 days. We conclude that binaries play a pivotal role in deriving the dynamical masses of young (~10 Myr) moderately massive and compact (<1e5 Msun; > 1 pc) star clusters.

The relation between the globular cluster luminosity function (GCLF, dN/dlogL) and globular cluster mass function (GCMF, dN/dlogM) is considered. Due to low-mass star depletion, dissolving GCs have mass-to-light (M/L) ratios that are lower than expected from their metallicities. This has been shown to lead to an M/L ratio that increases with GC mass and luminosity. We model the GCLF and GCMF and show that the power law slopes inherently differ (1.0 versus 0.7, respectively) when accounting for the variability of M/L. The observed GCLF is found to be consistent with a Schechter-type initial cluster mass function and a mass-dependent mass-loss rate.

We performed numerical simulations of dynamical encounters between hard massive binaries and a very massive star (VMS; formed through runaway mergers of ordinary stars in the dense core of a young massive star cluster), in order to explore the hypothesis that this dynamical process could be responsible for the origin of high-velocity (\geq 200-400 km/s) early or late B-type stars. We estimated the typical velocities produced in encounters between very tight massive binaries and VMSs (of mass of \geq 200 Msun) and found that about 3-4 per cent of all encounters produce velocities of \geq 400 km/s, while in about 2 per cent of encounters the escapers attain velocities exceeding the Milky Ways's escape velocity. We therefore argue that the origin of high-velocity (\geq 200-400 km/s) runaway stars and at least some so-called hypervelocity stars could be associated with dynamical encounters between the tightest massive binaries and VMSs formed in the cores of star clusters. We also simulated dynamical encounters between tight massive binaries and single ordinary 50-100 Msun stars. We found that from 1 to \simeq 4 per cent of these encounters can produce runaway stars with velocities of \geq 300-400 km/s (typical of the bound population of high-velocity halo B-type stars) and occasionally (in less than 1 per cent of encounters) produce hypervelocity (\geq 700 km/s) late B-type escapers.

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 discuss the observations and theory of star cluster formation to argue that clusters form dynamically cool (subvirial) and with substructure. We then perform an ensemble of simulations of cool, clumpy (fractal) clusters and show that they often dynamically mass segregate on timescales far shorter than expected from simple models. The mass segregation comes about through the production of a short-lived, but very dense core. This shows that in clusters like the Orion Nebula Cluster the stars >4 Msun can dynamically mass segregate within the current age of the cluster. Therefore, the observed mass segregation in apparently dynamically young clusters need not be primordial, but could be the result of rapid and violent early dynamical evolution.

The conversion of the globular cluster luminosity function (GCLF, dN/dlogL) to the globular cluster mass function (GCMF, dN/dlogM) is addressed. Dissolving globular clusters (GCs) become preferentially depleted in low-mass stars, which have a high mass-to-light ratio. This has been shown to result in a mass-to-light ratio (M/L) that increases with GC luminosity or mass, because more massive GCs have lost a smaller fraction of their stars than low-mass GCs. Using GC models, we study the influence of the luminosity dependency of M/L on the inferred GCMF. The observed GCLF is consistent with a powerlaw or Schechter type GC initial mass function in combination with a cluster mass-dependent mass loss rate. Below the peak, the logarithmic slope of the GCMF is shallower than that of the GCLF (0.7 versus 1.0), whereas the peak mass is 0.1-0.3 dex lower when accounting for the variability of M/L than in the case where a constant M/L is adopted.

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.

The onset of runaway stellar collisions in young star clusters is more likely to initiate with an encounter between a binary and a third star than between two single stars. Using the initial conditions of such three-star encounters from direct $N$-body simulations, we model the resulting interaction by means of Smoothed Particle Hydrodynamics (SPH). We find that, in the majority of the cases considered, all three stars merge together, and in such three star mergers, the hydrodynamic simulations reveal that: (1) mass lost as ejecta can be a considerable fraction of the total mass in the system (up to $\sim25$%); (2) due to asymmetric mass loss, the collision product can sometimes receive a kick velocity that exceeds 10 km/s, large enough to allow the collision product to escape the core of the cluster; and (3) the energy of the ejected matter can be large enough (up to $\sim 3\times 10^{50}$ erg) to remove or disturb the inter cluster gas appreciably.

We explore the hypothesis that some high-velocity runaway stars attain their peculiar velocities in the course of exchange encounters between hard massive binaries and a very massive star (either an ordinary 50-100 Msun star or a more massive one, formed through runaway mergers of ordinary stars in the core of a young massive star cluster). In this process, one of the binary components becomes gravitationally bound to the very massive star, while the second one is ejected, sometimes with a high speed. We performed three-body scattering experiments and found that early B-type stars (the progenitors of the majority of neutron stars) can be ejected with velocities of $\ga$ 200-400 km/s (typical of pulsars), while 3-4 Msun stars can attain velocities of $\ga$ 300-400 km/s (typical of the bound population of halo late B-type stars). We also found that the ejected stars can occasionally attain velocities exceeding the Milky Ways's escape velocity.

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.

The anomalous chemical abundances and the structure of the Edgewood-Kuiper belt observed in the solar system constrain the initial mass and radius of the star cluster in which the sun was born to $M\simeq500$ to 3000 \msun and $R\simeq 1$ to 3 pc. When the cluster dissolved the siblings of the sun dispersed through the galaxy, but they remained on a similar orbit around the Galactic center. Today these stars hide among the field stars, but 10 to 60 of them are still present within a distance of $\sim 100$ pc. These siblings of the sun can be identified by accurate measurements of their chemical abundances, positions and their velocities. Finding even a few will strongly constrain the parameters of the parental star cluster and the location in the Galaxy where we were born.

N-body simulations are widely used to simulate the dynamical evolution of a variety of systems, among them star clusters. Much of our understanding of their evolution rests on the results of such direct N-body simulations. They provide insight in the structural evolution of star clusters, as well as into the occurrence of stellar exotica. Although the major pure N-body codes STARLAB/KIRA and NBODY4 are widely used for a range of applications, there is no thorough comparison study yet. Here we thoroughly compare basic quantities as derived from simulations performed either with STARLAB/KIRA or NBODY4. We construct a large number of star cluster models for various stellar mass function settings (but without stellar/binary evolution, primordial binaries, external tidal fields etc), evolve them in parallel with STARLAB/KIRA and NBODY4, analyse them in a consistent way and compare the averaged results quantitatively. For this quantitative comparison we develop a bootstrap algorithm for functional dependencies. We find an overall excellent agreement between the codes, both for the clusters' structural and energy parameters as well as for the properties of the dynamically created binaries. However, we identify small differences, like in the energy conservation before core collapse and the energies of escaping stars, which deserve further studies. Our results reassure the comparability and the possibility to combine results from these two major N-body codes, at least for the purely dynamical models (i.e. without stellar/binary evolution) we performed. (abridged)

In the cores of young dense star clusters repeated stellar collisions involving the same object can occur, which has been suggested to lead to the formation of an intermediate-mass black hole. In order to verify this scenario we compute the detailed evolution of the merger remnant of three sequences. We follow the evolution until the onset of carbon burning and estimate the final remnant mass to determine the ultimate fate of a runaway merger sequence. We use a detailed stellar evolution code to follow the evolution of the collision product. At each collision, we mix the two colliding stars, taking account of mass loss during the collision. During the stellar evolution we apply mass loss rates from the literature, as appropriate for the evolutionary stage of the merger remnant. We compute models for high ($Z=0.02$) and low ($Z=0.001$) metallicity to quantify metallicity effects. We find that the merger remnant becomes a Wolf-Rayet star before the end of core hydrogen burning. Mass loss from stellar winds dominates over the mass increase due to repeated mergers for all three merger sequences that we consider. In none of our high metallicity models an intermediate-mass black hole is formed, instead our models have a mass of 10--14 \Msun{} at the onset of carbon burning. For low metallicity we expect the final remnant of the merger sequence to explode as a pair creation supernova. We find that our metal-rich models become inflated as a result of developing an extended low-density envelope. This may increase the probability of further collisions, but self-consistent $N$-body calculations with detailed evolution of runaway mergers are required to verify this.

We present a new method to detect and quantify mass segregation in star clusters. It compares the minimum spanning tree (MST) of massive stars with that of random stars. If mass segregation is present, the MST length of the most massive stars will be shorter than that of random stars. This difference can be quantified (with an associated significance) to measure the degree of mass segregation. We test the method on simulated clusters in both 2D and 3D and show that the method works as expected. We apply the method to the Orion Nebula Cluster (ONC) and show that the method is able to detect the mass segregation in the Trapezium with a `mass segregation ratio' \Lambda_{MSR}=8.0 \pm 3.5 (where \Lambda_{MSR}=1 is no mass segregation) down to 16 \Msun, and also that the ONC is mass segregated at a lower level (~2.0 \pm 0.5) down to 5 \Msun. Below 5 \Msun we find no evidence for any further mass segregation in the ONC.

Using new and archival radio data, we have measured the proper motion of the black hole X-ray binary V404 Cyg to be 9.2+/-0.3 mas/yr. Combined with the systemic radial velocity from the literature, we derive the full three-dimensional heliocentric space velocity of the system, which we use to calculate a peculiar velocity in the range 47-102 km/s, with a best fitting value of 64 km/s. We consider possible explanations for the observed peculiar velocity, and find that the black hole cannot have formed via direct collapse. A natal supernova is required, in which either significant mass (approximately 11 solar masses) was lost, giving rise to a symmetric Blaauw kick of up to 65 km/s, or, more probably, asymmetries in the supernova led to an additional kick out of the orbital plane of the binary system. In the case of a purely symmetric kick, the black hole must have been formed with a mass of approximately 9 solar masses, since when it has accreted 0.5-1.5 solar masses from its companion.

Knowledge of the binary population in stellar groupings provides important information about the outcome of the star forming process in different environments. Binarity is also a key ingredient in stellar population studies and is a prerequisite to calibrate the binary evolution channels. In these proceedings we present an overview of several commonly used methods to pair individual stars into binary systems, which we refer to as the pairing function. Many pairing functions are frequently used by observers and computational astronomers, either for the mathematical convenience, or because they roughly describe the expected outcome of the star forming process. We discuss the consequences of each pairing function for the interpretation of observations and numerical simulations. The binary fraction and mass ratio distribution generally depend strongly on the selection of the range in primary spectral type in a sample. These quantities, when derived from a binary survey with a mass-limited sample of target stars, are thus not representative for the population as a whole.

Knowledge of the binary population in stellar groupings provides important information about the outcome of the star forming process in different environments (see, e.g., Blaauw 1991, and references therein). Binarity is also a key ingredient in stellar population studies, and is a prerequisite to calibrate the binary evolution channels. In this paper we present an overview of several commonly used methods to pair individual stars into binary systems, which we refer to as pairing functions. These pairing functions are frequently used by observers and computational astronomers, either for their mathematical convenience, or because they roughly describe the expected outcome of the star forming process. We discuss the consequences of each pairing function for the interpretation of observations and numerical simulations. The binary fraction and mass ratio distribution generally depend strongly on the selection of the range in primary spectral type in a sample. The mass ratio distribution and binary fraction derived from a binarity survey among a mass-limited sample of targets is thus not representative for the population as a whole. Neither theory nor observations indicate that random pairing of binary components from the mass distribution, the simplest pairing function, is realistic. It is more likely that companion stars are formed in a disk around a star, or that a pre-binary core fragments into two binary components. The results of our analysis are important for (i) the interpretation of the observed mass ratio distribution and binary fraction for a sample of stars, (ii) a range of possible initial condition algorithms for star cluster simulations, and (iii) how to discriminate between the different star formation scenarios.

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

(Abridged) Aims: We obtained VLT/FLAMES+UVES high-resolution, fibre-fed spectroscopy (FFS) of five young massive clusters in M83 (NGC 5236). This forms the basis of a pilot study testing the feasibility of using FFS to measure the velocity dispersions of several clusters simultaneously, in order to determine their dynamical masses; Methods: We adopted two methods for determining the velocity dispersion of the star clusters: cross-correlating the cluster spectrum with the template spectra and minimising a chi^2 value between the cluster spectrum and the broadened template spectra. Cluster 805 in M83 was chosen as a control to test the reliability of the method, through a comparison with the results obtained from a standard echelle VLT/UVES spectrum obtained by Larsen & Richtler; Results: We find no dependence of the velocity dispersions measured for a cluster on the choice of red giant versus red supergiant templates, nor on the method adopted. We measure a velocity dispersion of sigma_los = 10.2+/-1.1 km/s for cluster 805 from our FFS. Our FLAMES+UVES velocity dispersion measurement gives M_vir = (6.6+/-1.7)e5 M_sun, consistent with previous results. This is a factor of ~3 greater than the cluster's photometric mass, indicating a lack of virial equilibrium. However, based on its effective star formation efficiency, the cluster is likely to virialise, and may survive for a Hubble time, in the absence of external disruptive forces; Conclusions: We find that reliable velocity dispersions can be determined from FFS. The advantages of observing several clusters simultaneously outweighs the difficulty of accurate galaxy background subtraction, providing that the targets are chosen to provide sufficient S/N ratios, and are much brighter than the galaxy background.

Recent proper motion and parallax measurements for the pulsar PSR B1508+55 indicate a transverse velocity of ~1100 km/s, which exceeds earlier measurements for any neutron star. The spin-down characteristics of PSR B1508+55 are typical for a non-recycled pulsar, which implies that the velocity of the pulsar cannot have originated from the second supernova disruption of a massive binary system. The high velocity of PSR B1508+55 can be accounted for by assuming that it received a kick at birth or that the neutron star was accelerated after its formation in the supernova explosion. We propose an explanation for the origin of hyperfast neutron stars based on the hypothesis that they could be the remnants of a symmetric supernova explosion of a high-velocity massive star which attained its peculiar velocity (similar to that of the pulsar) in the course of a strong dynamical three- or four-body encounter in the core of dense young star cluster. To check this hypothesis we investigated three dynamical processes involving close encounters between: (i) two hard massive binaries, (ii) a hard binary and an intermediate-mass black hole, and (iii) a single star and a hard binary intermediate-mass black hole. We find that main-sequence O-type stars cannot be ejected from young massive star clusters with peculiar velocities high enough to explain the origin of hyperfast neutron stars, but lower mass main-sequence stars or the stripped helium cores of massive stars could be accelerated to hypervelocities. Our explanation for the origin of hyperfast pulsars requires a very dense stellar environment of the order of 10^6 -10^7 stars pc^{-3}. Although such high densities may exist during the core collapse of young massive star clusters, we caution that they have never been observed.

We implement an efficient method to quantify time-dependent orbital complexity in gravitational $N$-body simulations. The technique, which we name DWaTIM, is based on a discrete wavelet transform of velocity orbital time series. The wavelet power-spectrum is used to measure trends in complexity continuously in time. We apply the method to the test cases N=3 Pythagorean- and a perturbed N=5 Caledonian configurations. The method recovers the well-known time-dependent complexity of the dynamics in these small-$N$ problems. We then apply the technique to an equal-mass collisional N=256 body simulation ran through core-collapse. We find that a majority of stars evolve on relatively complex orbits up to the time when the first hard binary forms, whereas after core-collapse, less complex orbits are found on the whole as a result of expanding mass shells.

We propose an explanation for the origin of hyperfast neutron stars (e.g. PSR B1508+55, PSR B2224+65, RX J0822-4300) based on the hypothesis that they could be the remnants of a symmetric supernova explosion of a high-velocity massive star (or its helium core) which attained its peculiar velocity (similar to that of the neutron star) in the course of a strong three- or four-body dynamical encounter in the core of a young massive star cluster. This hypothesis implies that the dense cores of star clusters (located either in the Galactic disk or near the Galactic centre) could also produce the so-called hypervelocity stars -- the ordinary stars moving with a speed of ~1000 km/s.

We study the circumstances under which first collisions occur in young and dense star clusters. The initial conditions for our direct $N$-body simulations are chosen such that the clusters experience core collapse within a few million years, before the most massive stars have left the main-sequence. It turns out that the first collision is typically driven by the most massive stars in the cluster. Upon arrival in the cluster core, by dynamical friction, massive stars tend to form binaries. The enhanced cross section of the binary compared to a single star causes other stars to engage the binary. A collision between one of the binary components and the incoming third star is then mediated by the encounters between the binary and other cluster members. Due to the geometry of the binary-single star engagement the relative velocity at the moment of impact is substantially different than in a two-body encounter. This may have profound consequences for the further evolution of the collision product.

Supernova 2006gy in the galaxy NGC 1260 is the most luminous one recorded \cite{2006CBET..644....1Q, 2006CBET..647....1H, 2006CBET..648....1P, 2006CBET..695....1F}. Its progenitor might have been a very massive ($>100$ \msun) star \cite{2006astro.ph.12617S}, but that is incompatible with hydrogen in the spectrum of the supernova, because stars $>40$ \msun are believed to have shed their hydrogen envelopes several hundred thousand years before the explosion \cite{2005A&A...429..581M}. Alternatively, the progenitor might have arisen from the merger of two massive stars \cite{2007ApJ...659L..13O}. Here we show that the collision frequency of massive stars in a dense and young cluster (of the kind to be expected near the center of a galaxy) is sufficient to provide a reasonable chance that SN 2006gy resulted from such a bombardment. If this is the correct explanation, then we predict that when the supernova fades (in a year or so) a dense cluster of massive stars becomes visible at the site of the explosion.

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 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 study the evolution and fate of solar composition supermassive stars in the mass range 60 - 1000 Msun. Our study is relevant both for very massive objects observed in young stellar complexes as well as supermassive stars that may potentially form through runaway stellar collisions. We predict the outcomes of stellar evolution employing a mass-loss prescription that is consistent with the observed Hertzsprung-Russell Diagram location of the most massive stars. We compute a series of stellar models with an appropriately modified version of the Eggleton evolutionary code. We find that super-massive stars with initial masses up to 1000 Msun end their lives as objects less massive than ~150 Msun. These objects are expected to collapse into black holes (with M < 70 Msun) or explode as pair-instability supernovae. We argue that if ultraluminous X-ray sources (ULXs) contain intermediate-mass black holes, these are unlikely to be the result of runaway stellar collisions in the cores of young clusters.

The early evolution of dense star clusters is possibly dominated by close interactions between stars, and physical collisions between stars may occur quite frequently. Simulating a stellar collision event can be an intensive numerical task, as detailed calculations of this process require hydrodynamic simulations in three dimensions. We present a computationally inexpensive method in which we approximate the merger process, including shock heating, hydrodynamic mixing and mass loss, with a simple algorithm based on conservation laws and a basic qualitative understanding of the hydrodynamics of stellar mergers. The algorithm relies on Archimedes' principle to dictate the distribution of the fluid in the stable equilibrium situation. We calibrate and apply the method to mergers of massive stars, as these are expected to occur in young and dense star clusters. We find that without the effects of microscopic mixing, the temperature and chemical composition profiles in a collision product can become double-valued functions of enclosed mass. Such an unphysical situation is mended by simulating microscopic mixing as a post-collision effect. In this way we find that head-on collisions between stars of the same spectral type result in substantial mixing, while mergers between stars of different spectral type, such as type B and O stars ($\sim$10 and $\sim$40\msun respectively), are subject to relatively little hydrodynamic mixing.

We characterize the binary population in the young and nearby OB association Scorpius OB2 using available observations of visual, spectroscopic, and astrometric binaries with intermediate-mass primaries. We take into account observational biases by comparing the observations with simulated observations of model associations. The available data indicate a large binary fraction (> 70% with 3sigma confidence), with a large probability that all intermediate mass stars in Sco OB2 are part of a binary system. The binary systems have a mass ratio distribution of the form f(q) ~ q^-0.4. Sco OB2 has a semi-major axis distribution of the form f(log a) ~ constant (Opik's law), in the range 5-5e6 Rsun. The log-normal period distribution of Duquennoy & Mayor results in too few spectroscopic binaries, even if the model binary fraction is 100%. Sco OB2 is a young association with a low stellar density; its current population is expected to be very similar to the primordial population. The fact that practically all stars in Sco OB2 are part of a binary (or multiple) system demonstrates that multiplicity is a fundamental factor in the star formation process, at least for intermediate mass stars.

We present the results of gravitational direct $N$-body simulations using the Graphics Processing Unit (GPU) on a commercial NVIDIA GeForce 8800GTX designed for gaming computers. The force evaluation of the $N$-body problem is implemented in ``Compute Unified Device Architecture'' (CUDA) using the GPU to speed-up the calculations. We tested the implementation on three different $N$-body codes: two direct $N$-body integration codes, using the 4th order predictor-corrector Hermite integrator with block time-steps, and one Barnes-Hut treecode, which uses a 2nd order leapfrog integration scheme. The integration of the equations of motions for all codes is performed on the host CPU. We find that for $N > 512$ particles the GPU outperforms the GRAPE-6Af, if some softening in the force calculation is accepted. Without softening and for very small integration time steps the GRAPE still outperforms the GPU. We conclude that modern GPUs offer an attractive alternative to GRAPE-6Af special purpose hardware. Using the same time-step criterion, the total energy of the $N$-body system was conserved better than to one in $10^6$ on the GPU, only about an order of magnitude worse than obtained with GRAPE-6Af. For $N \apgt 10^5$ the 8800GTX outperforms the host CPU by a factor of about 100 and runs at about the same speed as the GRAPE-6Af.

We re-analyze the age distribution (dN/dt) of star clusters in the Small Magellanic Cloud (SMC) using age determinations based on the Magellanic Cloud Photometric Survey. For ages younger than 3x10^9 yr the dN/dt distribution can be approximated by a power-law distribution, dN/dt propto t^-beta, with -beta=-0.70+/-0.05 or -beta=-0.84+/-0.04, depending on the model used to derive the ages. Predictions for a cluster population without dissolution limited by a V-band detection result in a power-law dN/dt distribution with an index of ~-0.7. This is because the limiting cluster mass increases with age, due to evolutionary fading of clusters, reducing the number of observed clusters at old ages. When a mass cut well above the limiting cluster mass is applied, the dN/dt distribution is flat up to 1 Gyr. We conclude that cluster dissolution is of small importance in shaping the dN/dt distribution and incompleteness causes dN/dt to decline. The reason that no (mass independent) infant mortality of star clusters in the first ~10-20 Myr is found is explained by a detection bias towards clusters without nebular emission, i.e. cluster that have survived the infant mortality phase. The reason we find no evidence for tidal (mass dependent) cluster dissolution in the first Gyr is explained by the weak tidal field of the SMC. Our results are in sharp contrast to the interpretation of Chandar et al. (2006), who interpret the declining dN/dt distribution as rapid cluster dissolution. This is due to their erroneous assumption that the sample is limited by cluster mass, rather than luminosity.

We study the effect of spiral arm passages on the evolution of star clusters on planar and circular orbits around the centres of galaxies. Individual passages with different relative velocity (V_drift) and arm width are studied using N-body simulations. When the ratio of the time it takes the cluster to cross the density wave to the crossing time of stars in the cluster is much smaller than one, the energy gain of stars can be predicted accurately in the impulsive approximation. When this ratio is much larger than one, the cluster is heated adiabatically and the net effect of heating is largely damped. For a given duration of the perturbation, this ratio is smaller for stars in the outer parts of the cluster compared to stars in the inner part. The cluster energy gain due to perturbations of various duration as obtained from our N-body simulations is in good agreement with theoretical predictions taking into account the effect of adiabatic damping. Perturbations by the broad stellar component of the spiral arms on a cluster are in the adiabatic regime and, therefore, hardly contribute to the energy gain and mass loss of the cluster. We consider the effect of crossings through the high density shocked gas in the spiral arms, which result in a more impulsive compression of the cluster. The time scale of disruption is shortest at ~0.8-0.9 R_CR since there V_drift is low. This location can be applicable to the solar neighbourhood. In addition, the four-armed spiral pattern of the Milky Way makes spiral arms contribute more to the disruption of clusters than in a similar but two-armed galaxy. Still, the disruption time due to spiral arm perturbations there is about an order of magnitude higher than what is observed for the solar neighbourhood.[ABRIDGED]

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 present the results of gravitational direct $N$-body simulations using the commercial graphics processing units (GPU) NVIDIA Quadro FX1400 and GeForce 8800GTX, and compare the results with GRAPE-6Af special purpose hardware. The force evaluation of the $N$-body problem was implemented in Cg using the GPU directly to speed-up the calculations. The integration of the equations of motions were, running on the host computer, implemented in C using the 4th order predictor-corrector Hermite integrator with block time steps. We find that for a large number of particles ($N \apgt 10^4$) modern graphics processing units offer an attractive low cost alternative to GRAPE special purpose hardware. A modern GPU continues to give a relatively flat scaling with the number of particles, comparable to that of the GRAPE. The GRAPE is designed to reach double precision, whereas the GPU is intrinsically single-precision. For relatively large time steps, the total energy of the N-body system was conserved better than to one in $10^6$ on the GPU, which is impressive given the single-precision nature of the GPU. For the same time steps, the GRAPE gave somewhat more accurate results, by about an order of magnitude. However, smaller time steps allowed more energy accuracy on the grape, around $10^{-11}$, whereas for the GPU machine precision saturates around $10^{-6}$ For $N\apgt 10^6$ the GeForce 8800GTX was about 20 times faster than the host computer. Though still about a factor of a few slower than GRAPE, modern GPUs outperform GRAPE in their low cost, long mean time between failure and the much larger onboard memory; the GRAPE-6Af holds at most 256k particles whereas the GeForce 8800GTX can hold 9 million particles in memory.

We describe XMM-Newton observations taken 4.3 days prior to and 1.5 days subsequent to two remarkable events that were detected with Swift on 2006 September 21 from the candidate magnetar CXOU J164710.2-455216: (1) a 20 ms burst with an energy of 1e37 erg (15-150 keV), and (2) a rapid spin-down (glitch) with a fractionap period change of 1e-4. We find that the luminosity of the pulsar increased by a factor of 100 in the interval between observations, from 1e33 to 1e35 erg/s (0.5-8.0 keV), and that its spectrum hardened. The pulsed count rate increased by a factor of 10 (0.5-8.0 keV), but the fractional rms amplitude of the pulses decreased from 65 to 11 per cent, and their profile changed from being single-peaked to exhibiting three peaks. Similar changes have been observed from other magnetars in response to outbursts, such as that of 1E 2259+586 in 2002 June. We suggest that a plastic deformation of the neutron star's crust induced a very slight twist in the external magnetic field, which in turn generated currents in the magnetosphere that were the direct cause of the X-ray outburst.

In this chapter we will argue that studying such multi-scale multi-science systems gives rise to inherently hybrid models containing many different algorithms best serviced by different types of computing environments (ranging from massively parallel computers, via large-scale special purpose machines to clusters of PC's) whose total integrated computing capacity can easily reach the PFlop/s scale. Such hybrid models, in combination with the by now inherently distributed nature of the data on which the models `feed' suggest a distributed computing model, where parts of the multi-scale multi-science model are executed on the most suitable computing environment, and/or where the computations are carried out close to the required data (i.e. bring the computations to the data instead of the other way around). We presents an estimate for the compute requirements to simulate the Galaxy as a typical example of a multi-scale multi-physics application, requiring distributed Petaflop/s computational power.

We study the evolution of the mass function in young and dense star clusters by means of direct N-body simulations. Our main aim is to explain the recent observations of the relatively flat mass function observed near the centre of the Arches star cluster. In this region, the power law index of the mass function for stars more massive than about 5-6 solar mass, is larger than the Salpeter value by about unity; whereas further out, and for the lower mass stars, the mass function resembles the Salpeter distribution. We show that the peculiarities in the Arches mass function can be explained satisfactorily without primordial mass segregation. We draw two conclusions from our simulations: 1) The Arches initial mass function is consistent with a Salpeter slope down to ~1 solar mass, 2) The cluster is about half way towards core collapse. The cores of other star clusters with characteristics similar to those of the Arches are expected to show similar flattening in the mass functions for the high mass (>5 solar mass) stars.

We study the acceleration of the star HE0437-5439, to hypervelocity and discuss its possible origin in the Large Magellanic Cloud (LMC). The star has a radial velocity of 723 km/s and is located at a distance of 61 kpc from the Sun. With a mass of about 8 Msun, the travel time from the Galactic centre is of about 100 Myr, much longer than its main sequence lifetime. Given the relatively small distance to the LMC (18 kpc), we consider it likely that HE0437-5439 originated in the cloud rather than in the Galactic centre, like the other hypervelocity stars. The minimum ejection velocity required to travel from the LMC to its current location within its lifetime is of about 500 km/s. Such a high velocity can only be obtained in a dynamical encounter with a massive black hole. We perform 3-body scattering simulations in which a stellar binary encounters a massive black hole and find that a black hole more massive than 1000 Msun is necessary to explain the high velocity of HE0437-5439. We look for possible parent clusters for HE0437-5439 and find that NGC 2100 and NGC 2004 are young enough to host stars coeval to HE0437-5439 and dense enough to produce an intermediate mass black hole able to eject an 8 Msun star with hypervelocity.

We discuss the performance of direct summation codes used in the simulation of astrophysical stellar systems on highly distributed architectures. These codes compute the gravitational interaction among stars in an exact way and have an O(N^2) scaling with the number of particles. They can be applied to a variety of astrophysical problems, like the evolution of star clusters, the dynamics of black holes, the formation of planetary systems, and cosmological simulations. The simulation of realistic star clusters with sufficiently high accuracy cannot be performed on a single workstation but may be possible on parallel computers or grids. We have implemented two parallel schemes for a direct N-body code and we study their performance on general purpose parallel computers and large computational grids. We present the results of timing analyzes conducted on the different architectures and compare them with the predictions from theoretical models. We conclude that the simulation of star clusters with up to a million particles will be possible on large distributed computers in the next decade. Simulating entire galaxies however will in addition require new hybrid methods to speedup the calculation.

We reconstruct the initial two-body relaxation time at the half mass radius for a sample of young $\aplt 300$ Myr star clusters in the large Magellanic cloud. We achieve this by simulating star clusters with 12288 to 131072 stars using direct $N$-body integration. The equations of motion of all stars are calculated with high precision direct $N$-body simulations which include the effects of the evolution of single stars and binaries. We find that the initial relaxation times of the sample of observed clusters in the large Magellanic cloud ranges from about 200 Myr to about 2 Gyr. The reconstructed initial half-mass relaxation times for these clusters has a much narrower distribution than the currently observed distribution, which ranges over more than two orders of magnitude.

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 study the evolution of bound pairs of star clusters by means of direct N-body simulations. Our simulations include mass loss by stellar evolution. The initial conditions are selected to mimic the observed binary star cluster NGC 2136 and NGC 2137 in the Large Magellanic Cloud. Based on the rather old ages ($\sim 100$ Myr), masses, sizes of the two clusters and their projected separation, we conclude that the cluster pair must have been born with an initial separation of 15--20 pc. Clusters with a smaller initial separation tend to merge in $\aplt 60$ Myr due to loss of angular momentum from escaping stars. Clusters with a larger initial separation tend to become even more widely separated due to mass loss from the evolving stellar populations. The early orbital evolution of a binary cluster is governed by mass loss from the evolving stellar population and by loss of angular momentum from escaping stars. Mass loss by stellar winds and supernovae explosions in the first $\sim 30$ Myr causes the binary to expand and the orbit to become eccentric. The initially less massive cluster expands more quickly than the binary separation increases, and is therefore bound to initiate mass transfer to the more massive cluster. This process is quite contrary to stellar binaries in which the more massive star tends to initiate mass transfer. Since mass transfer proceeds on a thermal timescale from the less massive to the more massive cluster, this semi-detached phase is quite stable, even in an eccentric orbit until the orbital separation reaches the gyration radius of the two clusters, at which point both clusters merge to one.

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)

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 have performed N-body simulations of the formation of hyper-velocity stars (HVS) in the centre of the Milky Way due to inspiralling intermediate-mass black holes (IMBHs). We considered IMBHs of different masses, all starting from circular orbits at an initial distance of 0.1 pc. We find that the IMBHs sink to the centre of the Galaxy due to dynamical friction, where they deplete the central cusp of stars. Some of these stars become HVS and are ejected with velocities sufficiently high to escape the Galaxy. Since the HVS carry with them information about their origin, in particular in the moment of ejection, the velocity distribution and the direction in which they escape the Galaxy, detecting a population of HVS will provide insight in the ejection processes and could therefore provide indirect evidence for the existence of IMBHs. Our simulations show that HVS are generated in short bursts which last only a few Myrs until the IMBH is swallowed by the supermassive black hole (SMBH). HVS are ejected almost isotropically, which makes IMBH induced ejections hard to distinguish from ejections due to encounters of stellar binaries with a SMBH. After the HVS have reached the galactic halo, their escape velocities correlate with the distance from the Galactic centre in the sense that the fastest HVS can be found furthest away from the centre. The velocity distribution of HVS generated by inspiralling IMBHs is also nearly independent of the mass of the IMBH and can be quite distinct from one generated by binary encounters. Finally, our simulations show that the presence of an IMBH in the Galactic centre changes the stellar density distribution inside r<0.02 pc into a core profile, which takes at least 100 Myrs to replenish.

We present a qualitative diagnostic based on the continuous wavelet transform, for the detection of irregular behavior in time series of particle simulations. We apply the method to three qualitatively different gravitational 3-body encounters. The intrinsic irregular behavior of these encounters is well reproduced by the presented method, and we show that the method accurately identifies the irregular regime in these encounters. We also provide an instantaneous quantification for the degree of irregularity in these simulations. Furthermore we demonstrate how the method can be used to analyze larger systems by applying it to simulations with 100-particles. It turns out that the number of stars on irregular orbits is systematically larger for clusters in which all stars have the same mass compared to a multimass system. The proposed method provides a quick and sufficiently accurate diagnostic for identifying stars on irregular orbits in large scale N-body simulations.

We present the diffuse X-ray emission identified in Chandra observations of the young, massive Galactic star cluster Westerlund 1. After removing point-like X-ray sources down to a completeness limit of 2e31 erg/s, we identify 3e34 erg/s (2--8 keV) of diffuse emission. The spatial distribution of the emission can be described as a slightly-elliptical Lorentzian core with a half-width half-maximum along the major axis of 25+/-1", similar to the distribution of point sources in the cluster, plus a 5' halo of extended emission. The spectrum of the diffuse emission is dominated by a hard continuum component that can be described as a kT>3 keV thermal plasma that has a low iron abundance (<0.3 solar), or as non-thermal emission that could be stellar light that is inverse-Compton scattered by MeV electrons. Only 5% of the flux is produced by a kT=0.7 keV plasma. The low luminosity of the thermal emission and the lack of a 6.7 keV iron line suggests that <40,000 unresolved stars with masses between 0.3 and 2 Msun are present in the cluster. Moreover, the flux in the diffuse emission is a factor of two lower than would be expected from a supersonically-expanding cluster wind, and there is no evidence for thermal remnants produced by supernovae. Less than 1e-5 of the mechanical luminosity of the cluster is dissipated as 2--8 keV X-rays, leaving a large amount of energy that either is radiated at other wavelengths, is dissipated beyond the bounds of our image, or escapes into the intergalactic medium.

We investigate encounters between giant molecular clouds (GMCs) and star clusters. We propose a single expression for the energy gain of a cluster due to an encounter with a GMC, valid for all encounter distances and GMC properties. This relation is verified with N-body simulations of cluster-GMC encounters and excellent agreement is found. The fractional mass loss from the cluster is 0.25 times the fractional energy gain. This is because 75% of the injected energy goes to the velocities of escaping stars, that are higher than the escape velocity. We derive an expression for the cluster disruption time (t_dis) based on the mass loss from the simulations, taking into account the effect of gravitational focusing by the GMC. The disruption time depends on the cluster mass (M_c) and half-mass radius (r_h) as t_dis=2.0 S (M_c/10^4 M_sun)(3.75 pc/r_h)^3 Gyr, with S=1 for the solar neighbourhood and inversely proportional to the GMC density. The observed shallow relation between cluster radius and mass gives t_dis a power-law dependence on the mass with index 0.7, similar to that found from observations and from simulations of clusters dissolving in tidal fields (0.62). The constant of 2.0 Gyr is about a factor of 3.5 shorter than found from earlier simulations of clusters dissolving under the combined effect of galactic tidal field and stellar evolution. It is somewhat higher than the observationally determined value of 1.3 Gyr. It suggests, however, that the combined effect of tidal field and encounters with GMCs can explain the lack of old open clusters in the solar neighbourhood. GMC encounters can also explain the (very) short disruption time that was observed for star clusters in the central region of M51, since there rho_n is an order of magnitude higher than in the solar neighbourhood.

We study the encounters between stars clusters and giant molecular clouds (GMCs). The effect of these encounters has previously been studied analytically for two cases: 1) head-on encounters, for which the cluster moves through the centre of the GMC and 2) distant encounters, where the encounter distance p > 3*R_n, with p the encounter parameter and R_n the radius of the GMC. We introduce an expression for the energy gain of the cluster due to GMC encounters valid for all values of p and R_n. This analytical result is confronted with results from N-body simulations and excellent agreement is found. From the simulations we find that the fractional mass loss is only 25% of the fractional energy gain. This is because stars escape with velocities much higher than the escape velocity. Based on the mass loss, we derive a disruption time for star clusters due to encounters with GMCs of the form t_dis [Gyr] = 2.0*S*(M_c/10^4 M_sun)^gamma, with S=1 for the solar neighbourhood and inversely proportional with the global GMC density and gamma=1-3lambda, with lambda the index that relates the cluster half-mass radius to the cluster mass (r_h ~ M_c^lambda). The observed shallow relation between cluster radius and mass (e.g. lambda=0.1), makes the index (gamma=0.7) similar to the index found both from observations and from simulations of clusters dissolving in tidal fields (gamma=0.62). The constant of 2.0 Gyr, which is the disruption time of a 10^4 M_sun cluster in the solar neighbourhood, is close to the value of 1.3 Gyr which was empirically determined from the age distribution of open clusters. This suggests that the combined effect of GMC encounters, stellar evolution and galactic tidal field can explain the lack of old open clusters in the solar neighbourhood.

Short-hard gamma-ray bursts (SHBs) may arise from gravitational wave (GW) driven mergers of double neutron star (DNS) systems. DNSs may be "primordial" or can form dynamically by binary exchange interactions in globular clusters during core-collapse. For primordial binaries, the time delay between formation and merger is expected to be short, tau~0.1 Gyr, implying that the redshift distribution of merger events should follow that of star-formation. We point out here that for dynamically formed DNSs, the time delay between star-formation and merger is dominated by the cluster core-collapse time, rather than by the GW inspiral time, yielding delays comparable to the Hubble time. We derive the redshift distribution of merger events of dynamically formed DNSs, and find it to differ significantly from that typically expected for primordial binaries. The observed redshift distribution of SHBs favors dynamical formation, although a primordial origin cannot be ruled out due to possible detection biases. Future red-shift observations of SHBs may allow to determine whether they are dominated by primordial or dynamically formed DNSs.

We present results of a population synthesis study for semidetached short orbital period binaries which contain low-mass(<1.5 Msun) donors and black hole (>4 Msun) accretors. Evolution of these binaries is determined by nuclear evolution of the donors and/or orbital angular momentum loss due to magnetic braking by the stellar wind of the donors and gravitational wave radiation. According to our model, the estimated total number of this type of black-hole binaries in the Galaxy is about 10000. If the magnetic braking is described by the Verbunt & Zwaan formula, the model predicts around 3000 transient systems with periods >2 hours and around 300 luminous stable systems with periods between 3 and 8 hours. Several dozens of these bright systems should be above the RXTE ASM sensitivity limit. The absence of such systems implies that angular momentum losses are reduced by a factor more than 2 with respect to the Verbunt & Zwaan prescription. We show that it is unlikely that the transient behaviour of black-hole short-period X-ray binaries is explained by the evolved nature of the stellar companion. A substantial fraction of black-hole binaries with periods >3 hours could be faint with truncated, stable cold accretion discs as proposed by Menou et al. Most of the semidetached black-hole binaries are expected to have periods shorter than ~2 hours. Properties of such, still to be observed, very small mass-ratio (q<0.02) binaries are different from those of their longer period cousins.

We describe a simple computer package which illustrates a method of combining stellar dynamics with stellar evolution. Though the method is intended for elaborate applications (especially the dynamical evolution of rich star clusters) it is illustrated here in the context of three-body scattering, i.e. interactions between a binary star and a field of single stars. We describe the interface between the dynamics and the two independent packages which describe the internal evolution of single stars and binaries. We also give an example application, and introduce a stand alone utility for the visual presentation of simulation results.

The starburst galaxy M82 at a distance of 12 million light years, is the host of an unusually bright $2.4-16\times 10^{40}$ erg/s X-ray point source, which is best explained by an accreting black hole $10^{2}$ to $10^4$ times more massive than the Sun. Though the strongest candidate for a so called intermediate mass black hole, the only support stems from the observed luminosity and the 0.05-0.1 Hz quasi periodicity in its signal. Interestingly, the $7-12 $Myr old star cluster MGG-11 which has been associated with the X-ray source is sufficiently dense that an intermediate mass black hole could have been produced in the cluster core via collision runaway. The recently discovered $62.0\pm 2.5$ day periodicity in the X-ray source X-1 further supports the hypothesis that this source is powered by a black hole several hundred times more massive than the Sun. We perform detailed binary evolution simulations with an accreting compact object of $10-5000 \msun $ and find that the X-ray luminosity, the age of the cluster, the observed quasi periodic oscillations and the now observed orbital period are explained best by a black hole of $200-5000 \msun $ that accretes material from a $22-25\msun$ giant companion in a state of Roche-lobe contact. Interestingly such a companion star is consistent with the expectation based on the tidal capture in a young and dense star cluster like MGG-11, making the picture self consistent.

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.

Recent X-ray observations and theoretical modelling have made it plausible that some ultraluminous X-ray sources (ULX) are powered by intermediate-mass black holes (IMBHs). N-body simulations have also shown that runaway merging of stars in dense star clusters is a way to form IMBHs. In the present paper we have performed N-body simulations of young clusters such as MGG-11 of M82 in which IMBHs form through runaway merging. We took into account the effect of tidal heating of stars by the IMBH to study the tidal capture and disruption of stars by IMBHs. Our results show that the IMBHs have a high chance of capturing stars through tidal heating within a few core relaxation times and we find that 1/3 of all runs contain a ULX within the age limits of MGG-11, a result consistent with the fact that a ULX is found in this galaxy. Our results strengthen the case for some ULX being powered by intermediate-mass black holes.

Ultraluminous X-ray sources (ULXs) with X-ray luminosities larger than the Eddington luminosity of stellar mass objects may be powered by intermediate mass black holes (IBHs) of masses Mbh~10^3Msun. If IBHs form in young dense stellar clusters, they can be fed by Roche lobe overflow from a tidally captured massive (Ms>10Msun) stellar companion. After the donor leaves the main sequence it forms a compact remnant, which spirals in due to gravitational wave (GW) emission. We show that space based detectors such as the Laser Interferometer Space Antenna are likely to detect several of these sources. GW sources stemming from this scenario have small eccentricities which give distinct GW signals. Detection of such a GW signal will unambiguously prove the existence of IBHs, and support the hypothesis that some ULXs are powered by IBHs with captured companions.

We investigate the effects of the addition of pre-main sequence evolution to star cluster simulations. We allowed stars to follow pre-main sequence tracks that begin at the deuterium burning birthline and end at the zero age main sequence. We compared our simulations to ones in which the stars began their lives at the zero age main sequence, and also investigated the effects of particular choices for initial binary orbital parameters. We find that the inclusion of the pre-main sequence phase results in a slightly higher core concentration, lower binary fraction, and fewer hard binary systems. In general, the global properties of star clusters remain almost unchanged, but the properties of the binary star population in the cluster can be dramatically modified by the correct treatment of the pre-main sequence stage.

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.

Hills (1988) predicted that runaway stars could be accelerated to velocities larger than 1000 km/s by dynamical encounters with the supermassive black hole (SMBH) in the Galactic center. The recently discovered hypervelocity star SDSS J090745.0+024507 (hereafter HVS) is escaping the Galaxy at high speed and could be the first object in this class. With the measured radial velocity and the estimated distance to the HVS, we trace back its trajectory in the Galactic potential. Assuming it was ejected from the center, we find that a $\sim$ 2 mas/yr proper motion is necessary for the star to have come within a few parsecs of the SMBH. We perform three-body scattering experiments to constrain the progenitor encounter which accelerated the HVS. As proposed by Yu & Tremaine (2003), we consider the tidal disruption of binary systems by the SMBH and the encounter between a star and a binary black hole, as well as an alternative scenario involving intermediate mass black holes. We find that the tidal disruption of a stellar binary ejects stars with a larger velocity compared to the encounter between a single star and a binary black hole, but has a somewhat smaller ejection rate due to the greater availability of single stars.

Cygnus X-3 is a strong X-ray source (L_X about 10^38 erg/s) which is thought to consist of a compact object, accreting matter from a helium star. We find analytically that the estimated ranges of mass-loss rate and orbital-period derivative for Cyg X-3 are consistent with two models: i) the system is detached and the mass loss from the system comes from the stellar wind of a massive helium star, of which only a fraction that allows for the observed X-ray luminosity is accreted, or ii) the system is semidetached and a Roche-lobe-overflowing low- or moderate-mass helium donor transfers mass to the compact object, followed by ejection of its excess over the Eddington rate from the system. These analytical results appear to be consistent with evolutionary calculations. By means of population synthesis we find that currently in the Galaxy there may exist ~1 X-ray binary with a black hole that accretes from a >~ 7 MSun Wolf-Rayet star and ~1 X-ray binary in which a neutron star accretes matter from a Roche-lobe-overflowing helium star with mass <~ 1.5 MSun. Cyg X-3 is probably one of these systems.

We present a simple analytical description of the disruption of star clusters in a tidal field, which agrees excellently with detailed N-body simulations. The analytic expression can be used to predict the mass and age histograms of surviving clusters for any cluster initial mass function and any cluster formation history. The method is applied to open clusters in the solar neighbourhood, based on the new cluster sample of Kharchenko et al. From a comparison between the observed and predicted age distributions in the age range between 10 Myr to 3 Gyr we find the following results: (1) The disruption time of a 10^4 M_sun cluster in the solar neighbourhood is about 1.3+/-0.5 Gyr. This is a factor 5 shorter than derived from N-body simulations of clusters in the tidal field of the galaxy. (2) The present starformation rate in bound clusters within 600 pc from the Sun is 5.9+/-0.8 * 10^2 M_sun / Myr, which corresponds to a surface star formation rate in bound clusters of 5.2+/-0.7 10^(-10) M_sun/yr/pc^2. (3) The age distribution of open clusters shows a bump between 0.26 and 0.6 Gyr when the cluster formation rate was 2.5 times higher than before and after. (4) The present star formation rate in bound clusters is half as small as that derived from the study of embedded clusters. The difference suggests that half of the clusters in the solar neighbourhood become unbound within 10 Myr. (5) The most massive clusters within 600 pc had an initial mass of 3*10^4 M_sun. This is in agreement with the statistically expected value based on a cluster initial mass function with a slope of -2, even if the physical upper mass limit is as high as 10^6 M_sun.

We discuss the signature of a peculiar constellation of very massive stars at a projected distance of 2-3 pc around R136a. We discuss various scenarios for its possible origin, such as independent clusters, triggered star formation, supernovae, and ejections via dynamical interactions. If the latter scenario were the correct one this would have significant implications on the way to probe the conditions in dense stellar cores, and on the evolution of massive clusters in general.

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 study the origin of the ultraluminous X-ray source M82-X1 in the nearby starburst galaxy M82. This X-ray source is of particular interest as it is currently the best candidate intermediate mass black hole; it is associated with a 54mHz quasi periodic oscillations with a relatively low (~1keV) black-body temperature. We perform detailed binary evolution calculations of 2--15Msun stars which transfer mass to a 100--2000Msun black hole and present an empirical model for the X-ray characteristics expected for such binaries. Based on the binary evolution calculations and the assumption in our simulations we conclude that the most likely candidate for the bright X-ray source M82-X1 is a 10-15Msun star near the end of its main-sequence or slightly evolved, which transfers mass to a ~1000Msun black hole. We expect the system to be in the high/soft state. In that case the binary will not be visible as a source of gravitational wave radiation, but other transient X-ray binaries with lower mass donors way be rather bright sources of gravitational wave radiation.

Young star clusters are often born with such high stellar densities that stellar collisions play an important role in their further evolution. In such environments the same star may participate in several tens to hundreds of collisions ultimately leading to the collapse of the star to a black hole of intermediate mass. At later time the black hole may acquire a companion star by tidal capture or by dynamical --3-body-- capture. When this companion star evolves it will fill its Roche-lobe and transfer mass to its accompanying black hole. This then leads to a bright phase of X-ray emission but the binary may also be visible as a source of gravitational wave radiation in the LISA band. If the star captured by the intermediate mass black hole is relatively low mass ($\aplt 2$ \msun) the binary will be visible as a bright source in gravitational waves and X-rays for the entire lifetime of the binary. The majority of compact binaries which formed from the population of primordial binaries in young and dense star clusters do not lead to detectable gravitational wave sources.

We explore the origin of the Galactic high latitude black hole X-ray binary XTE J1118+480, and in particular its birth location and the magnitude of the kick received by the black hole upon formation in the supernova explosion. We constrain the age of the companion to the black hole using stellar evolution calculations between 2 Gyr and 5 Gyr, making an origin in a globular cluster unlikely. We therefore argue that the system was born in the Galactic disk and the supernova propelled it in its current high latitude orbit. Given the current estimates on its distance, proper motion and radial velocity, we back-trace the orbit of XTE J1118+480 in the Galactic potential to infer the peculiar velocity of the system at different disk crossings over the last 5 Gyr. Taking into account the uncertainties on the velocity components, we infer an average peculiar velocity of 183 \pm 31 km/s. The maximum velocity which the binary can acquire by symmetric supernova mass loss is about 100 km/s, which is 2.7 sigma away from the mean of the peculiar velocity distribution. We therefore argue that an additional asymmetric kick velocity is required. By considering the orientation of the system relative to the plane of the sky, we derive a 95% probability for a non null component of the kick perpendicular to the orbital plane of the binary. The distribution of perpendicular velocities is skewed to lower velocities with an average of 93^{+55}_{-60} km/s.

We present the results of a near-infrared adaptive optics survey with the aim to detect close companions to Hipparcos members in the three subgroups of the nearby OB association Sco OB2: Upper Scorpius (US), Upper Centaurus Lupus (UCL) and Lower Centaurus Crux (LCC). We have targeted 199 A-type and late B-type stars in the Ks band, and a subset also in the J and H band. We find 151 stellar components other than the target stars. A brightness criterion is used to separate these components into 77 background stars and 74 candidate physical companion stars. Out of these 74 candidate companions, 41 have not been reported before (14 in US; 13 in UCL; 14 in LCC). Companion star masses range from 0.1 to 3 Msun. The mass ratio distribution follows f(q) = q^-0.33, which excludes random pairing. No close (rho < 3.75'') companion stars or background stars are found in the magnitude range 12 < Ks < 14. The lack of stars with these properties cannot be explained by low-number statistics, and may imply a lower limit on the companion mass of ~ 0.1 Msun. Close stellar components with Ks > 14 are observed. If these components are very low-mass companion stars, a gap in the companion mass distribution might be present. The small number of close low-mass companion stars could support the embryo-ejection formation scenario for brown dwarfs. Our findings are compared with and complementary to visual, spectroscopic, and astrometric data on binarity in Sco OB2. We find an overall companion star fraction of 0.52 in this association. This paper is the first step toward our goal to derive the primordial binary population in Sco OB2.

We show that one can obtain a good fit to the measured main sequence mass function (MF) of a large sample of Galactic clusters (young and old) with a tapered Salpeter power law distribution function with an exponential truncation. The average value of the power law index is very close to Salpeter (~2.3), whereas the characteristic mass is in the range 0.1 - 0.5 Msolar and does not seem to vary in a systematic way with the present cluster parameters such as metal abundance and central concentration. However, a remarkable correlation with age is seen, in that the peak mass of young clusters increases with it. This trend does not extend to globular clusters, whose peak mass is firmly at ~0.35 Msolar. This correlation is due to the onset of mass segregation following early dynamical interactions in the loose cluster cores. Differences between globular and younger clusters may depend on the initial environment of star formation, which in turn affects their total mass.

The observed average lifetime of the population of star clusters in the Solar Neighbourhood, the Small Magellanic Cloud and in selected regions of M51 and M33 is compared with simple theoretical predictions and with the results of N-body simulations. The empirically derived lifetimes (or disruption times) of star clusters depend on their initial mass as t_dis ~ Mcl^0.60 in all four galaxies. N-body simulations have shown that the predicted disruption time of clusters in a tidal field scales as t_dis^pred ~ t_rh^0.75 t_cr^0.25, where t_rh is the initial half-mass relaxation time and t_cr is the crossing time for a cluster in equilibrium. We show that this can be approximated accurately by t_dis^pred ~ M_cl^0.62 for clusters in the mass range of about 10^3 to 10^6 M_sun, in excellent agreement with the observations. Observations of clusters in different extragalactic environments show that t_dis also depends on the ambient density in the galaxies where the clusters reside. Linear analysis predicts that the disruption time will depend on the ambient density of the cluster environment as t_dis ~ rho_amb^-0.5. This relation is consistent with N-body simulations.

We investigate the possibility that mass transfer early in the evolution of a massive binary can effect a reversal of the end states of the two components, resulting in a neutron star which forms before a black hole. In this sense, such systems would comprise the high-mass analogues of white dwarf-neutron star systems such as PSR B2303+46. One consequence of this reversal is that a second episode of mass transfer from the black hole progenitor star can recycle the nascent neutron star, extending the life of the pulsar. An estimate of the formation rate through this channel is first performed via a simple analytical approach, and then refined using the results of the SeBa binary evolution package. The central role of kicks in determining the survival rate of these binaries is clearly demonstrated. The final result is expressed in terms of the number of field pulsars one can expect for every single neutron star-black hole (ns,bh) binary. We also calculate this figure for black hole-neutron star (bh,ns) systems formed through the usual channel. Assuming kicks drawn from the distribution of Cordes & Chernoff (1998), we find an expectation value of one (ns,bh) binary per 4 x 10^4 pulsars, and one (bh,ns) system for every 1500 pulsars. This helps to explain why neither system has been seen to date, though it suggests that detection of a (bh,ns) binary is imminent.

In this lecture we investigate the formation and evolution of black holes in star clusters. The star clusters under consideration are generally rich, containing more than 10^4 stars, and with a density exceeding 10^4 stars/pc^3. Among these are young dense clusters (YoDeCs), globular cluster and the nuclei of galaxies. We will also address the the possible evolutionary link between stellar mass black holes, via intermediate mass black holes to supermassive black holes, mainly focus on the ecology of star clusters.

We propose a new non-parametric algorithm that can be implemented to study and characterize stellar clusters. The scheme attempts to simultaneously recover the stellar distribution function and the cluster potential by using projected radii and velocity information about the cluster members. The pair of these functions that is most consistent with the input data is detected by the Metropolis algorithm. In this work, the cluster characteristics recovered by CHASSIS are calibrated against the N-body realizations of two clusters, namely Hyades and Arches. The cluster mass and line-of-sight projected velocity dispersion profiles are correctly reproduced by the algorithm when the cluster obeys the assumption used in the code, namely isotropy in phase space. The results recovered by the code are shown to be insensitive to the choice of the initial parameters. The results are also not influenced by increasing the number of input data points as long as this number exceeds a minimum value which is moderately low for an input data set that obeys the assumptions of isotropy and sphericity.

Physical collisions between stars occur frequently in dense star clusters, either via close encounters between two single stars, or during strong dynamical interactions involving binary stars. Here we study stellar collisions that occur during binary-single and binary-binary interactions, by performing numerical scattering experiments. Our results include cross sections, branching ratios, and sample distributions of parameters for various outcomes. For interactions of hard binaries containing main-sequence stars, we find that the normalized cross section for at least one collision to occur (between any two of the four stars involved) is essentially unity, and that the probability of collisions involving more than two stars is significant. Hydrodynamic calculations have shown that the effective radius of a collision product can be 2-30 times larger than the normal main-sequence radius for a star of the same total mass. We study the effect of this expansion, and find that it increases the probability of further collisions considerably. We discuss these results in the context of recent observations of blue stragglers in globular clusters with masses exceeding twice the main-sequence turnoff mass. We also present Fewbody, a new, freely available numerical toolkit for simulating small-N gravitational dynamics that is particularly suited to performing scattering experiments.

Core radii of globular clusters in the Large and Small Magellanic Clouds show an increasing trend with age. We propose that this trend is a dynamical effect resulting from the accumulation of massive stars and stellar-mass black holes at the cluster centers. The black holes are remnants of stars with initial masses exceeding 20-25 solar masses; as their orbits decay by dynamical friction, they heat the stellar background and create a core. Using analytical estimates and N-body experiments, we show that the sizes of the cores so produced and their growth rates are consistent with what is observed. We propose that this mechanism is responsible for the formation of cores in all globular clusters and possibly in other systems as well.

The nature of ultraluminous X-ray sources (ULXs) is presently unknown. A possible explanation is that they are accreting intermediate mass black holes (IBHs) that are fed by Roche lobe overflow from a tidally captured stellar companion. We show that a star can circularize around an IBH without being destroyed by tidal heating (in contrast to the case of M_bh> 10^6 M_sun massive black holes in galactic centers, where survival is unlikely). We find that the capture and circularization rate is of the order of 5 \times 10^-8 yr^-1, almost independently of the cluster's relaxation time. We follow the luminosity evolution of the binary system during the main sequence Roche lobe overflow phase and show it can maintain ULX-like luminosities for >10 Myr. In particular, we show that the ULX in the young cluster MGG-11 in star-burst galaxy M82, which possibly harbors an IBH, is well explained by this mechanism, and we predict that \gtrsim 10% of similar clusters with IBHs have a tidally captured circularized star. The cluster can evaporate on a time-scale shorter than the lifetime of the binary. This raises the possibility of a ULX that outlives its host cluster, or even lights up only after the cluster has evaporated, in agreement with observations of host-less ULXs.

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 study the dynamical interaction in which the two single runaway stars AE Aurigae and mu Columbae and the binary iota Orionis acquired their unusually high space velocity. The two single runaways move in almost opposite directions with a velocity greater than 100 km/s away from the Trapezium cluster. The star iota Ori is an eccentric (e=0.8) binary moving with a velocity of about 10 km/s at almost right angles with respect to the two single stars. The kinematic properties of the system suggest that a strong dynamical encounter occurred in the Trapezium cluster about 2.5 Myr ago. Curiously enough, the two binary components have similar spectral type but very different masses, indicating that their ages must be quite different. This observation leads to the hypothesis that an exchange interaction occurred in which an older star was swapped into the original iota Orionis binary. We test this hypothesis by a combination of numerical and theoretical techniques, using N-body simulations to constrain the dynamical encounter, binary evolution calculations to constrain the high orbital eccentricity of iota Orionis and stellar evolution calculations to constrain the age discrepancy of the two binary components. We find that an encounter between two low eccentricity (0.4<e<0.6) binaries with comparable binding energy, leading to an exchange and the ionization of the wider binary, provides a reasonable solution to this problem.

We discuss the model for the Galactic sample of the AM CVn systems with orbital periods P less than 1500 s, that can be detected in the optical and/or X-ray bands and may be resolved by the gravitational waves detector LISA. At P=3 to 10 min all detectable systems are X-ray sources. At P greater than 10 min most systems are only detectable in the V-band. About 30% of the X-ray sources is also detectable in the V-band. About 10,000 AM CVn systems might be resolved by LISA; this is comparable to the number of detached double white dwarfs that can be resolved. Several hundreds of AM CVn LISA sources might be also detectable in the V- and/or X-ray bands.

We model the population of AM CVn systems in the Galaxy and discuss the detectability of these systems with optical, X-ray and gravitational wave detectors. We concentrate on the short period (P < 1500 s) systems, some of which are expected to be in a phase of direct impact accretion. Using a self-consistent model for the star formation history and radial distribution of stars in the Galaxy plus simple models for the emission of optical and X-ray radiation from the AM CVn systems and interstellar absorption, we derive the sample of short-period AM CVn systems that can be detected in the optical and/or X-ray bands. At the shortest periods the detectable systems are all X-ray sources, some with periods as short as three minutes. At periods above 10 minutes most detectable systems are optical sources. About one third of the X-ray sources are also detectable in the optical band. We also calculate the gravitational wave signal of the short-period AM CVn systems. We find that potentially several thousand AM CVn systems can be resolved by the gravitational wave detector LISA, comparable to the expected number of detached double white dwarfs that can be resolved. We estimate that several hundreds of the AM CVn systems resolvable by LISA, are also detectable in the optical and/or X-ray bands.

Using adaptive optics we study the binary population in the nearby OB association Scorpius OB2. We present the first results of our near-infrared adaptive optics survey among 199 (mainly) A- and B-type stars in Sco OB2. In total 151 components other than the target stars are found, out of which 77 are probably background stars. Our findings are compared with data collected from literature. Out of the remaining 74 candidate physical companions 42 are new, demonstrating that many stars A/B stars have faint, close companions.

We study the efficiency at which a black hole or dense star cluster spirals-in to the Galactic centre. This process takes place on a dynamical friction time scale, which depends on the value of the Coulomb logarithm ln(L). We determine the accurate value of this parameter using the direct N-body method, a tree algorithm and a particle-mesh technique with up to 2 million plus one particles. The three different techniques are in excellent agreement. Our result for the Coulomb logarithm appears to be independent of the number of particles. We conclude that ln(L) = 6.6 +/- 0.6 for a massive point particle in the inner few parsec of the Galactic bulge. For an extended object, like a dense star cluster, ln(L) is smaller, with a value of the logarithm argument L inversely proportional to the object size.

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.

For understanding the process of star formation it is essential to know how many stars are formed as singles or in multiple systems, as a function of environment and binary parameters. This requires a characterization of the primordial binary population, which we define as the population of binaries that is present just after star formation has ceased, but before dynamical and stellar evolution have significantly altered its characteristics. In this article we present the first results of our adaptive optics survey of 200 (mainly) A-type stars in the nearby OB association Sco OB2. We report the discovery of 47 new candidate companions of Sco OB2 members. The next step will be to combine these observations with detailed simulations of young star clusters, in order to find the primordial binary population.

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.

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.

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 investigate the choice of stellar population for use as the Astrometric Grid for the Space Interferometry Mission (SIM). SIM depends on the astrometric stability of about 2000 stars, the so called Grid, against which the science measures are referenced. Low metallicity, and thus relatively high luminosity K giants are shown to be the population of choice, when available. The alternative, nearby G dwarfs, are shown to be suseptable to unmodeled motions induced by gas-giant planetary companions, should there be a significant population of such companions. Radial velocity filtering is quite efficient in selecting Grid members from the K giants with yields exceeding 50% if filtering at 30m/s (1-sigma) is available. However if the binary fraction of the G dwarfs approaches 100% as some studies suggest, the yield of stable systems would be in the range of 15% at best (with 10m/s filtering). Use of the initial SIM measurement as a final filter is shown not to be critical in either case, although it could improve the yield of stable grid members. For a Grid composed of weak-lined K giants, the residual contamination by large unmodeled motions will amount to about 3% (and rises to about 6% if a 60m/s radial velocity criterion is used). The selective introduction of quadratic terms in the proper motion solutions during the post-mission phase of data reduction can reduce contamination to a remarkable 1% or better in either case. Analytic estimates based on circular orbits are developed which show how these results come about.

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.

We analyzed archival Chandra X-ray observations of the central portion of the 30 Doradus region in the Large Magellanic Cloud. The image contains 20 X-ray point sources with luminosities between $5 \times 10^{32}$ and $2 \times 10^{35}$ erg s$^{-1}$ (0.2 -- 3.5 keV). A dozen sources have bright WN Wolf-Rayet or spectral type O stars as optical counterparts. Nine of these are within $\sim 3.4$pc of R136, the central star cluster of NGC2070. We derive an empirical relation between the X-ray luminosity and the parameters for the stellar wind of the optical counterpart. The relation gives good agreement for known colliding wind binaries in the Milky Way Galaxy and for the identified X-ray sources in NGC2070. We conclude that probably all identified X-ray sources in NGC2070 are colliding wind binaries and that they are not associated with compact objects. This conclusion contradicts Wang (1995) who argued, using ROSAT data, that two earlier discovered X-ray sources are accreting black-hole binaries. Five of the eighteen brightest stars in R136 are not visible in our X-ray observations. These stars are either single, have low mass companions or very wide orbits. The resulting binary fraction among early type stars is then unusually high (at least 70%).

We present the results of a new study of mass segregation in two-component star clusters, based on a large number of numerical N-body simulations using our recently developed dynamical Monte Carlo code. Specifically, we follow the dynamical evolution of clusters containing stars with individual masses m_1 as well as a tracer population of objects with individual masses m_2=\mu m_1, using N=10^5 total stars. For heavy tracers, which could represent stellar remnants such as neutron stars or black holes in a globular cluster, we characterize in a variety of ways the tendency for these objects to concentrate in or near the cluster core. In agreement with simple theoretical arguments, we find that the characteristic time for this mass segregation process varies as 1/\mu. For models with very light tracers (\mu <~ 10^-2), which could represent free-floating planets or brown dwarfs, we find the expected depletion of light objects in the cluster core, but also sometimes a significant enhancement in the halo. Using these results we estimate the optical depth to gravitational microlensing by planetary mass objects or brown dwarfs in typical globular clusters. For some initial conditions, the optical depth in the halo due to very low-mass objects could be much greater than that of luminous stars. If we apply our results to M22, using the recent null detection of Sahu, Anderson, & King (2001), we find an upper limit of ~25% at the 63% confidence level for the current mass fraction of M22 in the form of very low-mass objects.

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.

We present an internal shock model with external characteristics for explaining the complicated light curves of gamma-ray bursts. Shocks produce gamma-rays in the interaction between a precessing beam of relativistic particles and the interstellar medium. Each time the particle beam passes the same line of sight with the observer the inter stellar medium is pushed outward. Subsequent interactions between the medium and the beam are delayed by the extra distance to be traveled for the particles before the shock can form. This results in a natural retardation and leads to an intrinsic asymmetry in the produced light curves for gamma-ray bursts. In addition we account for the cooling of the electron-proton plasma in the shocked region, which gives rise to an exponential decay in the gamma-ray flux. The combination of these effects and the precessing jet of ultra relativistic particles produces light curves which can be directly compared with observed gamma-ray burst light curves. We illustrate the model by fitting a number of observed gamma-ray bursts which are hard to explain with only a precessing jet. With a genetic algorithm we are able to fit several observed gamma-ray bursts with remarkable accuracy. We find that for different bursts the observed fluence, assuming isotropic emission, easily varies over four orders of magnitude from the energy generated intrinsically.

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.

We review the properties of Galactic binaries containing two compact objects, as derived by means of population synthesis. Using this information we calculate the gravitational wave signal of these binaries. At frequencies below f ~< 2 mHz the double white dwarf population forms an unresolved background for the low-frequency gravitational wave detector LISA. Above this limit some few thousand double white dwarfs and few tens of binaries containing neutron stars will be resolved. Of the resolved double white dwarfs ~500 have a total mass above the Chandrasekhar limit. About ~95 of these have a measurable frequency change allowing a determination of their chirp mass. We discuss the properties of the resolved systems.

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 two models for AM CVn stars: white dwarfs accreting (i) from a helium white dwarf companion and (ii) from a helium-star donor. We show that in the first model possibly no accretion disk forms at the onset of mass transfer. The stability and the rate of mass transfer then depend on the tidal coupling between the accretor and the orbital motion. In the second model the formation of AM CVn stars may be prevented by detonation of the CO white dwarf accretor and the disruption of the system. With the most favourable conditions for the formation of AM CVn stars we find a current Galactic birth rate of 6.8 10^{-3} yr^{-1}. Unfavourable conditions give 1.1 10^{-3} yr^{-1}. The expected total number of the systems in the Galaxy is 9.4 10^{7} and 1.6 10^{7}, respectively. We model very simple selection effects to get some idea about the currently expected observable population and discuss the (quite good) agreement with the observed systems.

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 present a new model for the Galactic population of close double white dwarfs. The model accounts for the suggestion of the avoidance of a substantial spiral-in during mass transfer between a giant and a main-sequence star of comparable mass and for detailed cooling models. It agrees well with the observations of the local sample of white dwarfs if the initial binary fraction is close to 50% and an ad hoc assumption is made that white dwarfs with mass less than about 0.3 solar mass cool faster than the models suggest. About 1000 white dwarfs brighter than V=15 have to be surveyed for detection of a pair which has total mass greater than the Chandrasekhar mass and will merge within 10 Gyr.

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.

We model the population of double white dwarfs in the Galaxy and find a better agreement with observations compared to earlier studies, due to two modifications. The first is the treatment of the first phase of unstable mass transfer and the second the modelling of the cooling of the white dwarfs. A satisfactory agreement with observations of the local sample of white dwarfs is achieved if we assume that the initial binary fraction is ~ 50% and that the lowest mass white dwarfs (M < 0.3 Msun) cool faster than the most recently published cooling models predict. With this model we find a Galactic birth rate of close double white dwarfs of 0.05 yr^{-1}, a birth rate of AM CVn systems of 0.005 yr^{-1}, a merger rate of pairs with a combined mass exceeding the Chandrasekhar limit (which may be progenitors of SNe Ia) of 0.003 yr^{-1} and a formation rate of planetary nebulae of 1 yr^{-1}. We estimate the total number of double white dwarfs in the Galaxy as 2.5 10^8. In an observable sample with a limiting magnitude V_lim = 15 we predict the presence of ~855 white dwarfs of which ~220 are close pairs. Of these 10 are double CO white dwarfs of which one has a combined mass exceeding the Chandrasekhar limit and will merge within a Hubble time.

Globular clusters house a population of compact binaries that will be interesting gravitational wave sources for LISA. We provide estimates for the numbers of sources of several categories and discuss the sensitivity of LISA to detecting these sources. The estimated total number of detectable sources ranges from about 10 to about 1000 with gravitational wave frequencies above 1 mHz. These sources are typically undetectable by any other means and thus offer an opportunity for doing true gravitational-wave astronomy. The detection of these sources would provide information about both binary star evolution and the dynamics of globular clusters.

At least one, but more likely two or more, eccentric neutron-star, carbon-oxygen white-dwarf binaries with an unrecycled pulsar have been observed. According to the standard scenario for evolving neutron stars which are recycled in common envelope evolution we expect to observe $\gsim 50$ such circular neutron star-carbon oxygen white dwarf binaries, since their formation rate is roughly equal to that of the eccentric binaries and the time over which they can be observed is two orders of magnitude longer, as we shall outline. We observe at most one or two such circular binaries and from that we conclude that the standard scenario must be revised. Introducing hypercritical accretion into common envelope evolution removes the discrepancy by converting the neutron star into a black hole which does not emit radio waves, and therefore would not be observed.

We performed populations synthesis calculations of single stars and binaries and show that binary evolution is extremely important for Galactic astronomy. We review several binary evolution models and conclude that they give quite different results. These differences can be understood from the assumptions related to how mass is transfered in the binary systems. Most important are 1) the fraction of mass that is accreted by the companion star during mass transfer, 2) the amount of specific angular momentum which is carried away with the mass that leaves the binary system.

We study the evolution and observability of young and compact star clusters near the Galactic center, such as the Arches cluster and the Quintuplet. The star clusters are modeled with a combination of techniques; using direct N-body, integration to calculate the motions of all stars and detailed stellar and binary evolution to follow the evolution of the stars. The modeled star clusters dissolve within 10 to 60 million years in the tidal field of the Galaxy. The projected stellar density in the modeled clusters drops within 5% to 70% of the lifetime to a level comparable to the projected background density towards the Galactic center. And it will be very hard to distinguished these clusters at later age among the background stars. This effect is more severe for clusters at larger distance from the Galactic center but in projection at the same distance. Based on these arguments we conclude that the Galactic center easily hides 10 to 40 clusters with characteristics similar to the Arches and the Quintuplet cluster.

The unique core-mass - radius relation for giants with degenerate helium cores enables us to reconstruct the evolution of three observed double helium white dwarfs with known masses of both components. The last mass transfer phase in their evolution must have been a spiral-in. In the formalism proposed by Webbink (1984), we can constrain the efficiency of the deposition of orbital energy into the envelope to be $1 \la \alpha \la 6$, for an envelope structure parameter $\lambda=0.5$. We find that the two standard mass transfer types (stable mass transfer and spiral-in) are both unable to explain the first phase of mass transfer for these three binaries. We use a parametric approach to describe mass transfer in low-mass binaries, where both stars are of comparable mass and find that the orbital characteristics of the observed double helium white dwarfs can be well reproduced if the envelope of the primary is lost with ~1.5 times the specific angular momentum of the initial binary. In this case no substantial spiral-in occurs.

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.

The recent finding by Chevalier & Ilovaisky (1998) that OB-supergiant X-ray binaries have relatively large runaway velocities whereas Be/X-ray binaries have low runaway velocities, provides confirmation of the current models for the formation of these two types of systems. These predict a difference in runaway velocity of an order of magnitude. This difference basically results from the variation of the fractional helium core mass as a function of stellar mass, in combination with the conservation of orbital angular momentum during the mass transfer phase that preceded the formation of the compact object in the system. This combination results into: (i) Systematically narrower pre-supernova orbits in the OB-supergiant systems than in the Be-systems, and (ii) A larger fractional amount of mass ejected in the supernovae in high-mass systems relative to systems of lower mass. Regardless of possible kick velocities imparted to neutron stars at birth, this combination leads to a considerable difference in average runaway velocity between these two groups. The observed low runaway velocities of the Be/X-ray binaries confirm that in most cases not more than 1 to 2Msun was ejected in the supernovae that produced their neutron stars. This, in combination with the --on average-- large orbital eccentricities of these systems, indicates that their neutron stars must have received a velocity kick in the range 60 - 250 km/s at birth.

We perform binary population synthesis calculations to study the origin and the characteristics of runaway O and B stars which are ejected by the supernova explosion of the companion star in a binary system. The number of OB runaways can be explained from supernova ejections only if: high mass stars are preferentially formed in binaries, the initial mass ratio distribution is strongly peaked to unity and stars are rejuvenated to zero age upon accretion of mass from a companion star. Taking these requirements into consideration we conclude that at most 30% of the runaway O stars but possibly all runaway B stars obtain high velocities due to supernovae in evolving binaries. Stars which obtain high velocities via supernova ejections have the following characteristics: 1) at least 10% of the high velocity B stars and half the O stars have a mass greater than the turn off mass of the cluster in which they are born and would be observed as blue stragglers in the parent cluster, 2) their equatorial rotational velocities are proportional to their space velocity and 3) between 20% and 40% of the runaways have neutron star companions but less than 1% are visible as radio pulsars in part of the orbit.

We consider gravitational waves emitted by various populations of compact binaries at cosmological distances. We use population synthesis models to characterize the properties of double neutron stars, double black holes and double white dwarf binaries as well as white dwarf-neutron star, white dwarf-black hole and black hole-neutron star systems. We use the observationally determined cosmic star formation history to reconstruct the redshift distribution of these sources and their merging rate evolution. The gravitational signals emitted by each source during its early-inspiral phase add randomly to produce a stochastic background in the low frequency band with spectral strain amplitude between 10^{-18} Hz^{-1/2} and 5 10^{-17} Hz^{-1/2} at frequencies in the interval [5 10^{-6}-5 10^{-5}] Hz. The overall signal which, at frequencies above 10^{-4}Hz, is largely dominated by double white dwarf systems, might be detectable with LISA in the frequency range [1-10] mHz and acts like a confusion limited noise component which might limit the LISA sensitivity at frequencies above 1 mHz.

We investigate the evolution of globular clusters using N-body calculations and anisotropic Fokker-Planck (FP) calculations. The models include a mass spectrum, mass loss due to stellar evolution, and the tidal field of the parent galaxy. Recent N-body calculations have revealed a serious discrepancy between the results of N-body calculations and isotropic FP calculations. The main reason for the discrepancy is an oversimplified treatment of the tidal field employed in the isotropic FP models. In this paper we perform a series of calculations with anisotropic FP models with a better treatment of the tidal boundary and compare these with N-body calculations. The new tidal boundary condition in our FP model includes one free parameter. We find that a single value of this parameter gives satisfactory agreement between the N-body and FP models over a wide range of initial conditions. Using the improved FP model, we carry out an extensive survey of the evolution of globular clusters over a wide range of initial conditions varying the slope of the mass function, the central concentration, and the relaxation time. The evolution of clusters is followed up to the moment of core collapse or the disruption of the clusters in the tidal field of the parent galaxy. In general, our model clusters, calculated with the anisotropic FP model with the improved treatment for the tidal boundary, live longer than isotropic models. The difference in the lifetime between the isotropic and anisotropic models is particularly large when the effect of mass loss via stellar evolution is rather significant. On the other hand the difference is small for relaxation- dominated clusters which initially have steep mass functions and high central concentrations.

We present a new parallel supercomputer implementation of the Monte-Carlo method for simulating the dynamical evolution of globular star clusters. Our method is based on a modified version of Henon's Monte-Carlo algorithm for solving the Fokker-Planck equation. Our code allows us to follow the evolution of a cluster containing up to 5x10^5 stars to core collapse in < 40 hours of computing time. In this paper we present the results of test calculations for clusters with equal-mass stars, starting from both Plummer and King model initial conditions. We consider isolated as well as tidally truncated clusters. Our results are compared to those obtained from approximate, self-similar analytic solutions, from direct numerical integrations of the Fokker-Planck equation, and from direct N-body integrations performed on a GRAPE-4 special-purpose computer with N=16384. In all cases we find excellent agreement with other methods, establishing our new code as a robust tool for the numerical study of globular cluster dynamics using a realistic number of stars.

The high stellar densities in young compact star clusters, such as the star cluster R 136 in the 30 Doradus region, may lead to a large number of stellar collisions. Such collisions were recently found to be much more frequent than previous estimates. The number of collisions scales with the number of stars for clusters with the same initial relaxation time. These collisions take place in a few million years. The collision products may finally collapse into massive black holes. The fraction of the total mass in the star cluster which ends up in a single massive object scales with the total mass of the cluster and its relaxation time. This mass fraction is rather constant, within a factor two or so. Wild extrapolation from the relatively small masses of the studied systems to the cores of galactic nuclei may indicate that the massive black holes in these systems have formed in a similar way.

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.

We discuss the formation of pulsars with massive companions in eccentric orbits. We demonstrate that the probability for a non-recycled radio pulsar to have a white dwarf as a companion is comparable to that of having an old neutron star as a companion. Special emphasis is given to PSR B1820-11 and PSR B2303+46. Based on population synthesis calculations we argue that PSR B1820-11 and PSR B2303+46 could very well be accompanied by white dwarfs with mass > 1.1Msun. For PSR B1820-11, however, we can not exclude the possibility that its companion is a main-sequence star with a mass between about 0.7Msun and 5Msun.

Double white dwarfs could be important sources for space based gravitational wave detectors like OMEGA and LISA. We use population synthesis to predict the current population of double white dwarfs in the Galaxy and the gravitational waves produced by this population. We simulate a detailed power spectrum for an observation with an integration time of 10^6 s. At frequencies below ~3 mHz confusion limited noise dominates. At higher frequencies a few thousand double white dwarfs are resolved individually. Including compact binaries containing neutron stars and black holes in our calculations yields a further few hundred resolved binaries and some tens which can be detected above the double white dwarf noise at low frequencies. We find that binaries in which one white dwarf transfers matter to another white dwarf are rare, and thus unimportant for gravitational wave detectors. We discuss the uncertainties and compare our results with other authors.

A phenomenological model is presented to explain the light curves of gamma-ray bursts. Gamma-rays are produced in a narrow beam which sweeps through space due to the precession of a slaved accretion disc. The light curve expected from such a precessing luminosity cone can explain the complex temporal behavior of bright gamma-ray bursts.

We present a phenomenological model to explain the light curves of gamma-ray bursts. In the model a black hole is orbited by a precessing accretion disc which is fed by a neutron star. Gamma-rays are produced in a highly collimated beam via the Blandford-Znajek mechanism. The gamma-ray beam sweeps through space due to the precession of the slaved accretion disc. The light curve expected from such a precessing luminosity cone can explain the complex temporal behavior of observed bright gamma-ray bursts.

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.

By comparing relative rates of supernovae versus formation rates of single radio pulsars, recycled pulsars, binary pulsars and X-ray binaries we put strong limits on the progenitors of radio pulsars and on the requirement of an asymmetry in the supernova. The assumption that radio pulsars are only formed in type Ib and type Ic supernovae from interacting binaries (Iben & Tutukov 1996) breaks down on the implication that in that case either the formation rate of binary pulsars (double neutron stars) should be of the order of 20%, of the single pulsar birthrate or, alternatively, almost all single pulsars (85% to 98%) should originate from Thorne-Zytkow stars. In the latter case the pulsar velocity distribution is inconsistent with observations. Also, in that case the difference between the supernova rate and the pulsar formation rate would be about one order of magnitude, i.e.: much larger than observed. Allowing type II supernovae from single stars and non-interacting binaries to form radio pulsars solves this conundrum, but then a kick is required in order to explain the high velocities of single radio pulsars. A kick is also required to understand the small birthrate, relative to the supernova rate, of binary pulsars consisting of two neutron stars.

Gamma-ray bursts are characterized by a duration of milliseconds to several minutes in which an enormous amount of radiation is emitted. The origin of these phenomena is still unknown because proposed models fail to explain all the observed features. Our proposed solution to this conundrum is a new class of mass-exchanging binaries in which a neutron star transfers mass to a black hole. According to recent studies binaries which contain a neutron star and a black hole are much more frequent than was previously believed. Mass exchange is driven by the emission of gravitational waves but the redistribution of mass in the binary system prevents coalescence. The phase of mass transfer is surprisingly stable and lasts for several thousands of orbital revolutions (about a minute). With a simple analytic model we demonstrated that this new class of binaries could provide an excellent candidate for the observed phenomena known as gamma-ray bursts.

Recent N-body simulations have shown that there is a serious discrepancy between the results of the N-body simulations and the results of Fokker-Planck simulations for the evolution of globular and rich open clusters under the influence of the galactic tidal field. In some cases, the lifetime obtained by Fokker-Planck calculations is more than an order of magnitude smaller than those by N-body simulations. In this letter we show that the principal cause for this discrepancy is an over-simplified treatment of the tidal field used in previous Fokker-Planck simulations. We performed new Fokker-Planck calculations using a more appropriate implementation of the boundary condition of the tidal field. The implementation is only possible with anisotropic Fokker-Planck models, while all previous Fokker-Planck calculations rely on the assumption of isotropy. Our new Fokker-Planck results agree well with N-body results. Comparison of the two types of simulations gives a better understanding of the cluster evolution.

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.

The evolution of the comoving cosmic merger-rate density of neutron star binaries n_c(z) is calculated using a distribution of their merging times provided by population-synthesis computations of binary stars. We adopt an exponential law for the star formation rate with various timescales for different morphological types of galaxies. For elliptical galaxies also an initial burst of star formation, lasting one Gyr, is considered. The resulting n_c(z) of most models agree with the form n_c(z) \propto (1+z)^{1.5-2} for z < 2, which has been proposed for the source population of gamma-ray bursts. Assuming a standard candle luminosity, the computed brightness distribution is consistent with the BATSE results if bursts at the peak flux threshold, P = 0.4 photons cm^{-2}/s, are located at a limiting redshift of 1.9 to 3.3. Progenitors of the systems producing gamma-ray bursts at small redshift (bright) are likely to host in spiral galaxies and star forming regions whereas these at high redshift (dim) reside mainly in elliptical galaxies. The location of a burst may be up to a Mpc away from the host galaxy.

We discuss the formation and evolution of binaries which contain neutron stars or black holes. It is shown that in a stellar system which for $10^{10}$ yr had star formation rate similar to the current one in the Galactic disc, the rate of neutron star mergers is $\sim 2 \times 10^{-5} yr{-1}$, consistent with the observational estimates. The rate of black hole and neutron star mergers is by an order of magnitude lower. The pairs of black holes form but do not merge because heavy mass loss by their progenitors efficiently moves the binary components apart.

The formation and evolution of binaries which contain two neutron stars or a neutron star with a black hole are discussed in detail. The evolution of the distributions in orbital period and eccentricity for neutron star binaries are studied as a function of time. In the model which fits the observations of high mass binary pulsars best the deposition of orbital energy into common envelopes has to be very efficient and a kick velocity distribution has to contain a significant contribution of low velocity kicks. The estimated age of the population has to be between several 100 Myr and 1 Gyr. The birthrate of binary neutron stars is about 3.4 10^{-5} per year (assuming 100% binarity) and their merger rate is about 2 x 10^{-5} per year. The merger rate of neutron star binaries is consistent with the estimated rate of gamma-ray bursts, if the latter are beamed into an opening angle of a few degrees. We argue that PSR B2303+46 is possibly formed in a scenario in which the common envelope is avoided while for the other three known high-mass binary pulsars a common envelope is required to explain their orbital period.

Dynamical models of star clusters are maturing in the sense that effects other than simple point particle dynamics are taken into account. We summarize the relevance of and prospects for this new generation of N-body models.

We study the formation of low-mass X-ray binaries with a black hole as accreting object. The performed semi-analytic analysis reveals that the formation rate of black holes in low-mass X-ray binaries is about two orders of magnitude smaller than that of systems with a neutron star as accretor. This is contradicted by the six observed systems, which are all transients, which suggest that the majority of these systems has not been seen jet. The birthrate for both type of objects are expected to be similar (for reviews see Cowley 1992, Tanaka & Lewin 1995). Solution of this conundrum requires either that stars with masses less than 20~\msun can evolve into black holes, or that stellar wind from a member of a binary is accompanied by a much larger loss of angular momentum than hitherto assumed.

The absence or presence of extremely wide binaries with a radio pulsar and an optical counterpart imposes a strong constraint on the existence and magnitude of kicks in supernova explosions. We search for such systems by comparing the positions of radio pulsars which are not known to be in binaries with the positions of visible stars, and find that the number of associations is negligible. According to the performed population synthesis, this implies that kicks must occur, with a lower limit of at least 10 to 20 km/s.

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.

Star clusters with a high central density contain an ecological network of evolving binaries, affected by interactions with passing stars, while in turn affecting the energy budget of the cluster as a whole by giving off binding energy. This is the first paper in a series aimed at providing the tools for increasingly realistic simulations of these ecological networks.

We calculate the formation rates of low-mass X-ray binaries with a black hole. Both a semi-analytic and a more detailed model predict formation rates two orders of magnitude lower than derived from the observations. Solution of this conundrum requires either that stars with masses less than 20M$_\odot$ can evolve into a black hole, or that stellar wind from a member of a binary is accompanied by a much larger loss of angular momentum than hitherto assumed.