BZPEERReading, trust and collaboration for research
Menu

Discover

GraphFollow connections around a paper, topic or saved view.

Workspaces

HomePick up where your research flow left off.InboxHandle requests, replies and notifications from one place.CollectionsCurate reading lists, evidence sets and shareable dossiers.ResearchSee your private provenance graph, trail and imports.CollaborationMove from discovery into focused discussion rooms.PublishDraft, preview and publish full-text research articles.

Network

ExpertsFind specialists worth following or asking for help.CommunitiesJoin research circles organized around shared work.PartnersSee external organizations active around the same topics.

Opportunities

ListingsBrowse roles, calls and collaborations with clearer fit signals.

Account

Log inReturn to your reading and collaboration workspace.Create accountStart saving work, following people and joining teams.
Log inCreate account

Researcher profile

Simon Portegies Zwart

Simon Portegies Zwart contributes to research discovery and scholarly infrastructure.

ResearcherAffiliation not importedOpen to collaborate
astro-ph.GAastro-ph.SRDistributed, Parallel, and Cluster Computingastro-ph.IMastro-ph.EPastro-ph.HEnlin.CDphysics.comp-phastro-ph.COComputational Engineering, Finance, and Sciencegr-qchep-phMachine Learningphysics.space-ph

Trust snapshot

Quick read

Trust 21 - EmergingVerification L1Unclaimed author
27works
0followers
14topics
4close collaborators

Actions

Decide how to stay connected

Follow researcher0

Identity and collaboration

How to connect with this researcher

Claiming links this public author record to a researcher profile and unlocks direct collaboration workflows.

Log in to claim

Direct collaboration

Open a focused conversation when the fit is right

Claim this author entity first to unlock direct invitations.

Research graph

See the researcher in context

Open full explorer

Inspect adjacent work, topics, institutions and collaborators without jumping out to a separate graph page.

Building this graph slice

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

Published work

27 published item(s)

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.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
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 ...

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
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&#39;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.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
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}.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access