Source author record

Simon Portegies Zwart

Simon Portegies Zwart appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

Catalog footprint

What is connected

81works

21topics

4close collaborators

Actions

Connect this record

Open graph Browse works

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

preprint2022arXiv

Birth cluster simulations of planetary systems with multiple super-Earths: initial conditions for white dwarf pollution drivers

Previous investigations have revealed that eccentric super-Earths represent a class of planets which are particularly effective at transporting minor bodies towards white dwarfs and subsequently polluting their atmospheres with observable chemical signatures. However, the lack of discoveries of these planets beyond a few astronomical units from their host stars prompts a better understanding of their orbital architectures from their nascent birth cluster. Here, we perform stellar cluster simulations of 3-planet and 7-planet systems containing super-Earths on initially circular, coplanar orbits. We adopt the typical stellar masses of main-sequence progenitors of white dwarfs ($1.5\,\mathrm{M}_{\odot}$-$2.5\,\mathrm{M}_{\odot}$) as host stars and include 8,000 main-sequence stars following a Kroupa initial mass function in our clusters. Our results reveal that about 30 per cent of the simulated planets generate eccentricities of at least 0.1 by the time of cluster dissolution, which would aid white dwarf pollution. We provide our output parameters to the community for potential use as initial conditions for subsequent evolution simulations.

preprint2022arXiv

Impact of bar resonances in the velocity-space distribution of the solar neighbourhood stars in a self-consistent $N$-body Galactic disc simulation

The velocity-space distribution of the solar neighbourhood stars shows complex substructures. Most of the previous studies use static potentials to investigate their origins. Instead we use a self-consistent $N$-body model of the Milky Way, whose potential is asymmetric and evolves with time. In this paper, we quantitatively evaluate the similarities of the velocity-space distributions in the $N$-body model and that of the solar neighbourhood, using Kullback-Leibler divergence (KLD). The KLD analysis shows the time evolution and spatial variation of the velocity-space distribution. The KLD fluctuates with time, which indicates the velocity-space distribution at a fixed position is not always similar to that of the solar neighbourhood. Some positions show velocity-space distributions with small KLDs (high similarities) more frequently than others. One of them locates at $(R,ϕ)=(8.2\;\mathrm{kpc}, 30^{\circ})$, where $R$ and $ϕ$ are the distance from the galactic centre and the angle with respect to the bar's major axis, respectively. The detection frequency is higher in the inter-arm regions than in the arm regions. In the velocity maps with small KLDs, we identify the velocity-space substructures, which consist of particles trapped in bar resonances. The bar resonances have significant impact on the stellar velocity-space distribution even though the galactic potential is not static.

preprint2022arXiv

Lucky planets: how circum-binary planets survive the supernova in one of the inner-binary components

A planet hardly ever survives the supernova of the host star in a bound orbit, because mass loss in the supernova and the natal kick imparted to the newly formed compact object cause the planet to be ejected. A planet in orbit around a binary has a considerably higher probability to survive the supernova explosion of one of the inner binary stars. In those cases, the planet most likely remains bound to the companion of the exploding star, whereas the compact object is ejected. We estimate this to happen to $\sim 1/33$ the circum-binary planetary systems. These planetary orbits tend to be highly eccentric ($e \apgt 0.9$), and $\sim 20$\,\% of these planets have retrograde orbits compared to their former binary. The probability that the planet as well as the binary (now with a compact object) remains bound is about ten times smaller ($\sim 3\cdot 10^{-3}$). We then expect the Milky way Galaxy to host $\aplt 10$ x-ray binaries that are still orbited by a planet, and $\aplt 150$ planets that survived in orbit around the compact object's companion. These numbers should be convolved with the fraction of massive binaries that is orbited by a planet.

preprint2022arXiv

Revisiting Common Envelope Evolution -- A New Semi-Analytic Model for N-body and Population Synthesis Codes

We present a novel way of modeling common envelope evolution in binary and few-body systems. We consider the common envelope inspiral as driven by a drag force with a power-law dependence in relative distance and velocity. The orbital motion is resolved either by direct N-body integration or by solving the set of differential equations for the orbital elements as derived using perturbation theory. Our formalism can model the eccentricity during the common envelope inspiral, and it gives results consistent with smoothed particles hydrodynamical simulations. We apply our formalism to common envelope events from binary population synthesis models and find that the final eccentricity distribution resembles the observed distribution of post-common-envelope binaries. Our model can be used for time-resolved common-envelope evolution in population synthesis calculations or as part of binary interactions in direct N-body simulations of star clusters.

preprint2022arXiv

Statistical features of multiple Compton scattering in a strong magnetic field

Compton scattering is a key process shaping spectra formation and accretion flow dynamics in accreting strongly magnetized neutron stars. A strong magnetic field affects the scattering cross section and makes it dependent on photon energy, momentum, and polarization state. Using Monte Carlo simulations, we investigate statistical features of Compton scattering of polarized X-ray radiation in a strong magnetic field. Our analysis is focused on photon gas behaviour well inside the scattering region. We take into account the resonant scattering at the fundamental cyclotron frequency, thermal distribution of electrons at the ground Landau level, and bulk velocity of the electron gas. We show that (i) the photons scattered around the cyclotron energy by the electron gas at rest tend to acquire the final energy close to the cyclotron one with a very small dispersion measure; (ii) the redistribution of photons within the Doppler core of cyclotron resonance differs significantly from the complete redistribution; (iii) the efficiency of momentum transfer from photons to the electron gas is affected by the temperature of electron gas both for photons at cyclotron energy and below it; (iv) the momentum transfer from photons to the electron gas of non-zero bulk velocity is more efficient in the case of magnetic scattering.

preprint2022arXiv

The paradox of infinitesimal granularity: Chaos and the reversibility of time in Newton's theory of gravity

The fundamental laws of physics are time-symmetric, but our macroscopic experience contradicts this. The time reversibility paradox is partly a consequence of the unpredictability of Newton's equations of motion. We measure the dependence of the fraction of irreversible, gravitational N-body systems on numerical precision and find that it scales as a power law. The stochastic wave packet reduction postulate then introduces fundamental uncertainties in the Cartesian phase space coordinates that propagate through classical three-body dynamics to macroscopic scales within the triple's lifetime. The spontaneous collapse of the wave function then drives the global chaotic behavior of the Universe through the superposition of triple systems (and probably multi-body systems). The paradox of infinitesimal granularity then arises from the superposition principle, which states that any multi-body system is composed of an ensemble of three-body problems.

preprint2020arXiv

Bonsai-SPH: A GPU accelerated astrophysical Smoothed Particle Hydrodynamics code

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.

preprint2020arXiv

Deep-learning enhancement of large scale numerical simulations

Traditional simulations on High-Performance Computing (HPC) systems typically involve modeling very large domains and/or very complex equations. HPC systems allow running large models, but limits in performance increase that have become more prominent in the last 5-10 years will likely be experienced. Therefore new approaches are needed to increase application performance. Deep learning appears to be a promising way to achieve this. Recently deep learning has been employed to enhance solving problems that traditionally are solved with large-scale numerical simulations using HPC. This type of application, deep learning for high-performance computing, is the theme of this whitepaper. Our goal is to provide concrete guidelines to scientists and others that would like to explore opportunities for applying deep learning approaches in their own large-scale numerical simulations. These guidelines have been extracted from a number of experiments that have been undertaken in various scientific domains over the last two years, and which are described in more detail in the Appendix. Additionally, we share the most important lessons that we have learned.

preprint2020arXiv

Gargantuan chaotic gravitational three-body systems and their irreversibility to the Planck length

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.

preprint2020arXiv

Identifying Earth-impacting asteroids using an artificial neural network

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.

preprint2020arXiv

Implementing Primordial Binaries in Simulations of Star Cluster Formation with a Hybrid MHD and Direct N-Body Method

The fraction of stars in binary systems within star clusters is important for their evolution, but what proportion of binaries form by dynamical processes after initial stellar accretion remains unknown. In previous work, we showed that dynamical interactions alone produced too few low-mass binaries compared to observations. We therefore implement an initial population of binaries in the coupled MHD and direct N-body star cluster formation code Torch. We compare simulations with, and without, initial binary populations and follow the dynamical evolution of the binary population in both sets of simulations, finding that both dynamical formation and destruction of binaries take place. Even in the first few million years of star formation, we find that an initial population of binaries is needed at all masses to reproduce observed binary fractions for binaries with mass ratios above the $q \geq 0.1$ detection limit. Our simulations also indicate that dynamical interactions in the presence of gas during cluster formation modify the initial distributions towards binaries with smaller primary masses, larger mass ratios, smaller semi-major axes and larger eccentricities. Systems formed dynamically do not have the same properties as the initial systems, and systems formed dynamically in the presence of an initial population of binaries differ from those formed in simulations with single stars only. Dynamical interactions during the earliest stages of star cluster formation are important for determining the properties of binary star systems.

preprint2020arXiv

Linking the formation and fate of exo-Kuiper belts within solar system analogues

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.

preprint2020arXiv

Non-intrusive hierarchical coupling strategies for multi-scale simulations in gravitational dynamics

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

preprint2020arXiv

On the survival of resonant and non-resonant planetary systems in star clusters

Despite the discovery of thousands of exoplanets in recent years, the number of known exoplanets in star clusters remains tiny. This may be a consequence of close stellar encounters perturbing the dynamical evolution of planetary systems in these clusters. Here, we present the results from direct $N$-body simulations of multiplanetary systems embedded in star clusters containing $N = 8k, 16k, 32k$, and $64k$ stars. The planetary systems, which consist of the four Solar system giant planets Jupiter, Saturn, Uranus, and Neptune, are initialized in different orbital configurations, to study the effect of the system architecture on the dynamical evolution of the entire planetary system, and on the escape rate of the individual planets. We find that the current orbital parameters of the Solar system giants (with initially circular orbits, as well as with present-day eccentricities) and a slightly more compact configuration, have a high resilience against stellar perturbations. A configuration with initial mean-motion resonances of 3:2, 3:2, and 5:4 between the planets, which is inspired by the Nice model, and for which the two outermost planets are usually ejected within the first $10^5$ yr, is in many cases stabilized due to the removal of the resonances by external stellar perturbation and by the rapid ejection of at least one planet. Assigning all planets the same mass of 1 Jovian mass almost equalizes the survival fractions. Our simulations reproduce the broad diversity amongst observed exoplanet systems. We find not only many very wide and/or eccentric orbits, but also a significant number of (stable) retrograde orbits.

preprint2020arXiv

Searching for Solar Siblings in APOGEE and $Gaia$ DR2 with N-body Simulations

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.

preprint2019arXiv

Collisional N-Body Dynamics Coupled to Self-Gravitating Magnetohydrodynamics Reveals Dynamical Binary Formation

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.

preprint2019arXiv

Newton vs the machine: solving the chaotic three-body problem using deep neural networks

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.

preprint2016arXiv

Mass transfer between debris discs during close stellar encounters

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.

preprint2016arXiv

The contribution of young core-collapse supernova remnants to the X-ray emission near quiescent supermassive black holes

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.

preprint2015arXiv

A debris disk under the influence of a wide planetary mass companion: The system of HD106906

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.

preprint2015arXiv

From Thread to Transcontinental Computer: Disturbing Lessons in Distributed Supercomputing

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.

preprint2015arXiv

GPU-Enabled Particle-Particle Particle-Tree Scheme for Simulating Dense Stellar Cluster System

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.

preprint2015arXiv

How Sedna and family were captured in a close encounter with a solar sibling

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.

preprint2015arXiv

Sapporo2: A versatile direct $N$-body library

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.

preprint2015arXiv

Simulations of stripped core-collapse supernovae in close binaries

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.

preprint2015arXiv

Simulations of the tidal interaction and mass transfer of a star in an eccentric orbit around an intermediate-mass black hole: the case of HLX-1

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.

preprint2015arXiv

Stellar disk destruction by dynamical interactions in the Orion Trapezium star cluster

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.

preprint2015arXiv

Stellar winds near massive black holes: The case of the S-stars

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.

preprint2015arXiv

The fate of supernova remnants near quiescent supermassive black holes

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

preprint2015arXiv

The formation and dynamical evolution of young star clusters

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.

preprint2015arXiv

The fragility of planetary systems

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.

preprint2015arXiv

The influence of dense gas rings on the dynamics of a stellar disk in the Galactic center

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.

preprint2015arXiv

Was the nineteenth century giant eruption of Eta Carinae a merger event in a triple system?

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

preprint2014arXiv

24.77 Pflops on a Gravitational Tree-Code to Simulate the Milky Way Galaxy with 18600 GPUs

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.

preprint2014arXiv

Binary white dwarfs in the halo of the Milky Way

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.

preprint2014arXiv

Computational Gravitational Dynamics with Modern Numerical Accelerators

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.

preprint2014arXiv

Detecting gravitational waves from the galactic center with Pulsar Timing

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.

preprint2014arXiv

On the Effects of Subvirial Initial Conditions and the Birth Temperature of R136

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.

preprint2014arXiv

On the minimal accuracy required for simulating self-gravitating systems by means of direct N-body methods

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.

preprint2014arXiv

On the Reliability of N-body Simulations

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.

preprint2014arXiv

The cosmic web in CosmoGrid void regions

We study the formation and evolution of the cosmic web, using the high-resolution CosmoGrid \$Λ\$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.

preprint2013arXiv

Are Super-Luminous supernovae and Long GRBs produced exclusively in young dense star clusters?

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.

preprint2013arXiv

Assembly of filamentary void galaxy configurations

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.

preprint2013arXiv

Evolution of star clusters in a cosmological tidal field

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.

preprint2013arXiv

MPWide: a light-weight library for efficient message passing over wide area networks

We present MPWide, a light weight communication library which allows efficient message passing over a distributed network. MPWide has been designed to connect application running on distributed (super)computing resources, and to maximize the communication performance on wide area networks for those without administrative privileges. It can be used to provide message-passing between application, move files, and make very fast connections in client-server environments. MPWide has already been applied to enable distributed cosmological simulations across up to four supercomputers on two continents, and to couple two different bloodflow simulations to form a multiscale simulation.

preprint2013arXiv

Simulating Star Clusters with the AMUSE Software Framework: I. Dependence of Cluster Lifetimes on Model Assumptions and Cluster Dissolution Modes

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.

preprint2013arXiv

Structure and evolution of high-mass stellar mergers

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.

preprint2013arXiv

The Cosmogrid Simulation: Statistical Properties of Small Dark Matter Halos

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.

preprint2013arXiv

The Effect of Many Minor Mergers on the Size Growth of Compact Quiescent Galaxies

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.

preprint2013arXiv

The evolution of triples with a Roche-lobe filling outer star

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.

preprint2012arXiv

A pilgrimage to gravity on GPUs

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.

preprint2012arXiv

A sparse octree gravitational N-body code that runs entirely on the GPU processor

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.

preprint2012arXiv

Bonsai: A GPU Tree-Code

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.

preprint2012arXiv

High-Performance Distributed Multi-Model / Multi-Kernel Simulations: A Case-Study in Jungle Computing

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.

preprint2012arXiv

Mapping the Galactic Center with Gravitational Wave Measurements using Pulsar Timing

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.

preprint2012arXiv

Multi-physics simulations using a hierarchical interchangeable software interface

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.

preprint2012arXiv

Planet-Mediated Precision-Reconstruction of the Evolution of the Cataclysmic Variable HU Aquarius

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 $αλ= 0.45\pm 0.17$ or $γ= 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$.

preprint2012arXiv

PSDF: Particle Stream Data Format for N-Body Simulations

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.

preprint2012arXiv

Supernova Type Ia progenitors from merging double white dwarfs: Using a new population synthesis model

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.

preprint2011arXiv

High performance cosmological simulations on a grid of supercomputers

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.

preprint2011arXiv

High Performance Gravitational N-body Simulations on a Planet-wide Distributed Supercomputer

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.

preprint2011arXiv

Multi-physics simulations using a hierarchical interchangeable software interface

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.

preprint2011arXiv

New population synthesis model: Preliminary results for close double white dwarf populations

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.

preprint2011arXiv

Simulations of Dense Stellar Systems with the AMUSE Software Toolkit

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.

preprint2011arXiv

The formation of the eccentric-orbit millisecond pulsar J1903+0327 and the origin of single millisecond pulsars

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.

preprint2011arXiv

The Origin of OB Runaway Stars

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.

preprint2010arXiv

A Light-Weight Communication Library for Distributed Computing

We present MPWide, a platform independent communication library for performing message passing between computers. Our library allows coupling of several local MPI applications through a long distance network and is specifically optimized for such communications. The implementation is deliberately kept light-weight, platform independent and the library can be installed and used without administrative privileges. The only requirements are a C++ compiler and at least one open port to a wide area network on each site. In this paper we present the library, describe the user interface, present performance tests and apply MPWide in a large scale cosmological N-body simulation on a network of two computers, one in Amsterdam and the other in Tokyo.

preprint2010arXiv

Gravitational tree-code on graphics processing units: implementation in CUDA

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 $θ\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}.

preprint2010arXiv

On the temporal evolution of the stellar mass function in Galactic clusters

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.

preprint2010arXiv

Simulating the universe on an intercontinental grid of supercomputers

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.

preprint2010arXiv

The distinction between star clusters and associations

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.

preprint2010arXiv

Young massive star clusters

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.

preprint2009arXiv

Effects of Primordial Mass Segregation on the Dynamical Evolution of Star Clusters

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.

preprint2009arXiv

On the onset of runaway stellar collisions in dense star clusters - II. Hydrodynamics of three-body interactions

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.

preprint2009arXiv

SAPPORO: A way to turn your graphics cards into a GRAPE-6

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.

preprint2007arXiv

High Performance Direct Gravitational N-body Simulations on Graphics Processing Units

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. 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, which is only about an order of magnitude worse than obtained with GRAPE. For $N\apgt 10^6$ the GeForce 8800GTX was about 20 times faster than the host computer. Though still about an order of magnitude slower than GRAPE, modern GPU's 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 8800GTF can hold 9 million particles in memory.

preprint2006arXiv

The redshift distribution of short gamma-ray bursts from dynamically formed neutron star binaries

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.

preprint2005arXiv

Short gamma-ray bursts from binary neutron star mergers in globular clusters

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.

preprint2004arXiv

The formation of massive black holes through collision runaway in dense young star clusters

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.

preprint2004arXiv

The stellar IMF of Galactic clusters and its evolution

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.

preprint1999arXiv

The possible companions of young radio pulsars

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.

Simon Portegies Zwart

What is connected

Connect this record

See the researcher in context

Building this map preview

81 published item(s)

Birth cluster simulations of planetary systems with multiple super-Earths: initial conditions for white dwarf pollution drivers

Impact of bar resonances in the velocity-space distribution of the solar neighbourhood stars in a self-consistent $N$-body Galactic disc simulation

Lucky planets: how circum-binary planets survive the supernova in one of the inner-binary components

Revisiting Common Envelope Evolution -- A New Semi-Analytic Model for N-body and Population Synthesis Codes

Statistical features of multiple Compton scattering in a strong magnetic field

The paradox of infinitesimal granularity: Chaos and the reversibility of time in Newton's theory of gravity

Bonsai-SPH: A GPU accelerated astrophysical Smoothed Particle Hydrodynamics code

Deep-learning enhancement of large scale numerical simulations

Gargantuan chaotic gravitational three-body systems and their irreversibility to the Planck length

Identifying Earth-impacting asteroids using an artificial neural network

Implementing Primordial Binaries in Simulations of Star Cluster Formation with a Hybrid MHD and Direct N-Body Method

Linking the formation and fate of exo-Kuiper belts within solar system analogues

Non-intrusive hierarchical coupling strategies for multi-scale simulations in gravitational dynamics

On the survival of resonant and non-resonant planetary systems in star clusters

Searching for Solar Siblings in APOGEE and $Gaia$ DR2 with N-body Simulations

Collisional N-Body Dynamics Coupled to Self-Gravitating Magnetohydrodynamics Reveals Dynamical Binary Formation

Newton vs the machine: solving the chaotic three-body problem using deep neural networks

Mass transfer between debris discs during close stellar encounters

The contribution of young core-collapse supernova remnants to the X-ray emission near quiescent supermassive black holes

A debris disk under the influence of a wide planetary mass companion: The system of HD106906

From Thread to Transcontinental Computer: Disturbing Lessons in Distributed Supercomputing

GPU-Enabled Particle-Particle Particle-Tree Scheme for Simulating Dense Stellar Cluster System

How Sedna and family were captured in a close encounter with a solar sibling

Sapporo2: A versatile direct $N$-body library

Simulations of stripped core-collapse supernovae in close binaries

Simulations of the tidal interaction and mass transfer of a star in an eccentric orbit around an intermediate-mass black hole: the case of HLX-1

Stellar disk destruction by dynamical interactions in the Orion Trapezium star cluster

Stellar winds near massive black holes: The case of the S-stars

The fate of supernova remnants near quiescent supermassive black holes

The formation and dynamical evolution of young star clusters

The fragility of planetary systems

The influence of dense gas rings on the dynamics of a stellar disk in the Galactic center

Was the nineteenth century giant eruption of Eta Carinae a merger event in a triple system?

24.77 Pflops on a Gravitational Tree-Code to Simulate the Milky Way Galaxy with 18600 GPUs

Binary white dwarfs in the halo of the Milky Way

Computational Gravitational Dynamics with Modern Numerical Accelerators

Detecting gravitational waves from the galactic center with Pulsar Timing

On the Effects of Subvirial Initial Conditions and the Birth Temperature of R136

On the minimal accuracy required for simulating self-gravitating systems by means of direct N-body methods

On the Reliability of N-body Simulations

The cosmic web in CosmoGrid void regions

Are Super-Luminous supernovae and Long GRBs produced exclusively in young dense star clusters?

Assembly of filamentary void galaxy configurations

Evolution of star clusters in a cosmological tidal field

MPWide: a light-weight library for efficient message passing over wide area networks

Simulating Star Clusters with the AMUSE Software Framework: I. Dependence of Cluster Lifetimes on Model Assumptions and Cluster Dissolution Modes

Structure and evolution of high-mass stellar mergers

The Cosmogrid Simulation: Statistical Properties of Small Dark Matter Halos

The Effect of Many Minor Mergers on the Size Growth of Compact Quiescent Galaxies

The evolution of triples with a Roche-lobe filling outer star

A pilgrimage to gravity on GPUs

A sparse octree gravitational N-body code that runs entirely on the GPU processor

Bonsai: A GPU Tree-Code

High-Performance Distributed Multi-Model / Multi-Kernel Simulations: A Case-Study in Jungle Computing

Mapping the Galactic Center with Gravitational Wave Measurements using Pulsar Timing

Multi-physics simulations using a hierarchical interchangeable software interface

Planet-Mediated Precision-Reconstruction of the Evolution of the Cataclysmic Variable HU Aquarius

PSDF: Particle Stream Data Format for N-Body Simulations

Supernova Type Ia progenitors from merging double white dwarfs: Using a new population synthesis model

High performance cosmological simulations on a grid of supercomputers

High Performance Gravitational N-body Simulations on a Planet-wide Distributed Supercomputer

Multi-physics simulations using a hierarchical interchangeable software interface

New population synthesis model: Preliminary results for close double white dwarf populations

Simulations of Dense Stellar Systems with the AMUSE Software Toolkit

The formation of the eccentric-orbit millisecond pulsar J1903+0327 and the origin of single millisecond pulsars

The Origin of OB Runaway Stars

A Light-Weight Communication Library for Distributed Computing

Gravitational tree-code on graphics processing units: implementation in CUDA

On the temporal evolution of the stellar mass function in Galactic clusters

Simulating the universe on an intercontinental grid of supercomputers

The distinction between star clusters and associations

Young massive star clusters

Effects of Primordial Mass Segregation on the Dynamical Evolution of Star Clusters

On the onset of runaway stellar collisions in dense star clusters - II. Hydrodynamics of three-body interactions