BZPEERPapers, people, evidence and research paths
Menu

Discover

SearchFind papers, people, topics, institutions and opportunities from one place.GraphSee how papers, authors, topics and evidence connect.

Research tools

HomeContinue from the papers, people and topics that matter now.CollectionsTurn saved papers into reading lists, evidence maps and shareable context.ResearchReconstruct the path that led you from search to paper to question.CollaborationMove from a paper, topic or expert request into a focused thread.PublishCreate full-text research articles with metadata, formulas and versions.

Network

ExpertsFind specialists whose work and reviews explain why they can help.CommunitiesJoin research groups organized around papers, topics and shared evidence.

Opportunities

ListingsBrowse roles, calls and collaborations connected to real research context.

Account

Log inReturn to your saved papers, graph views and active threads.Create accountStart saving papers, building a research trail and joining teams.
Log inCreate account

Source author record

S. Portegies Zwart

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

ResearcherUnclaimed source record
astro-ph.GAastro-ph.SRastro-ph.EPastro-ph.IMastro-ph.COastro-ph.HEastro-phDistributed, Parallel, and Cluster Computinggr-qcNetworking and Internet Architecture

Catalog footprint

What is connected

35works
10topics
4close collaborators

Actions

Connect this record

Log in to claim
Open graphBrowse works

Research graph

See the researcher in context

Open full explorer

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.

Published work

35 published item(s)

preprint2022arXiv

A Monte Carlo study of early gas expulsion and evolution of star clusters: new simulations with the MOCCA code in the AMUSE framework

We introduce a new prescription for the evolution of globular clusters (GCs) during the initial embedded gas phase into a Monte Carlo method. With a simplified version of the Monte Carlo MOCCA code embedded in the AMUSE framework, we study the survival of GCs after the removal of primordial gas. We first test our code and show that our results for the evolution of mass and Lagrangian radii are in good agreement with those obtained with N-body simulations. The Monte Carlo code enables a more rapid exploration of the evolution of systems with a larger number of stars than N-body simulations. We have carried out a new survey of simulations to explore the evolution of globular clusters with up to $N = 500000$ stars for a range of different star formation efficiencies and half-mass radii. Our study shows the range of initial conditions leading to the clusters' dissolution and those for which the clusters can survive this early evolutionary phase.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Stellar triples on the edge; Comprehensive overview of the evolution of destabilised triples leading to stellar and binary exotica

Hierarchical triple stars are ideal laboratories for studying the interplay between orbital dynamics and stellar evolution. Both stellar wind mass loss and three-body dynamics cooperate to destabilise triples, which can lead to a variety of astrophysical exotica. So far our understanding of their evolution was mainly built upon results from extensive binary-single scattering experiments. Starting from generic initial conditions, we evolve an extensive set of hierarchical triples using a combination of the triple evolution code TRES and an N-body code. We find that the majority of triples preserve their hierarchy throughout their evolution, which is in contradiction with the commonly adopted picture that unstable triples always experience a chaotic, democratic resonant interaction. The duration of the unstable phase is much longer than expected, so that stellar evolution cannot be neglected. Typically an unstable triple dissolve into a single star and a binary; sometimes democratically (the initial hierarchy is lost and the lightest body usually escapes), but also in a hierarchical way (the tertiary is ejected in a slingshot, independent of its mass). Collisions are common, and mostly involve the two original inner binary components still on the main-sequence. This contradicts the idea that collisions with a giant during democratic encounters dominate. Together with collisions in stable triples, we find that triple evolution is the dominant mechanism for stellar collisions in the Milky Way. Furthermore, our simulations produce runaway and walk-away stars with speeds up to several tens km/s, with a maximum of a few 100km/s. We suggest that destabilised triples can alleviate the tension behind the origin of the observed run-away stars. Lastly, we present a promising indicator to make general predictions for the fate of a specific triple, based on the initial inclination of the system.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Survey of Orion Disks with ALMA (SODA) I: Cloud-level demographics of 873 protoplanetary disks

(Abridged) Surveys of protoplanetary disks in nearby star-forming regions (SFRs) have provided important information on their demographics, but due to their sample sizes, they cannot be used to study how disk properties vary with the environment. We conduct a survey of the unresolved millimeter continuum emission of 873 protoplanetary disks identified by Spitzer in the L1641 and L1647 regions of the Orion A cloud. This is the largest such survey yet, allowing us to identify even weak trends in the median disk mass as a function of position in the cloud and cluster membership. The sample detection rates and median masses are also compared to those of nearby (<300 pc) SFRs. The sample was observed with ALMA at 225 GHz, with a median rms of 0.08 mJy/beam (1.5 M$_{\oplus}$). The data were reduced and imaged using an innovative parallel data processing approach. We detect 58% (502/873) of the observed disks. This includes 20 disks with dust masses >100 M$_{\oplus}$, and two objects associated with extended dust emission. We infer a median disk dust mass in the full sample of $2.2^{+0.2}_{-0.2}$ M$_{\oplus}$. In L1641 and L1647 median dust masses are $2.1^{+0.2}_{-0.2}$ M$_{\oplus}$ and $2.6^{+0.4}_{-0.5}$ M$_{\oplus}$, respectively. The disk mass distribution of the full sample is similar to that of nearby low-mass SFRs at similar ages of 1-3 Myr. We find only weak trends in disk (dust) masses with galactic longitude and between the YSO clusters identified in the sample, with median masses varying by $\lesssim$ 50%. Age differences may explain these median disk mass variations. Apart from this, disk masses are essentially constant at scales of ~100 pc. This also suggests that the majority of disks, even in different SFRs, are formed with similar initial masses and evolve at similar rates, assuming no external irradiation, with disk mass loss rates of $\sim 10^{-8}$ M$_{\odot}$/yr.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

APEX-SEPIA660 Early Science: Gas at densities above $10^7$ cm$^{-3}$ towards OMC-1

Context. The star formation rates and stellar densities found in young massive clusters suggest that these stellar systems originate from gas at densities n(H$_2$) $>10^7$ cm$^{-3}$. Until today, however, the physical characterization of this ultra high density material remains largely unconstrained in observations. Aims. We investigated the density properties of the star-forming gas in the OMC-1 region located in the vicinity of the Orion Nebula Cluster (ONC). Methods. We mapped the molecular emission at 652 GHz in OMC-1 as part of the APEX-SEPIA660 Early Science. Results. We detect bright and extended N$_2$H$^+$ (J=7-6) line emission along the entire OMC-1 region. Comparisons with previous ALMA data of the (J=1-0) transition and radiative transfer models indicate that the line intensities observed in this N$_2$H$^+$ (7-6) line are produced by large mass reservoirs of gas at densities n(H$_2$) $>10^7$ cm$^{-3}$. Conclusions. The first detection of this N$_2$H$^+$ (7-6) line at parsec-scales demonstrates the extreme density conditions of the star-forming gas in young massive clusters such as the ONC. Our results highlight the unique combination of sensitivity and mapping capabilities of the new SEPIA660 receiver for the study of the ISM properties at high frequencies.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

The evolution of stellar triples: The most common evolutionary pathways

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

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Face-on accretion onto a protoplanetary disc

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

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

The evolution of hierarchical triple star-systems

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

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

The evolution of the Sun's birth cluster and the search for the solar siblings with Gaia

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

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

The rate of stellar encounters along a migrating orbit of the Sun

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

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

The formation of the solar system

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

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

The initial mass function of star clusters that form in turbulent molecular clouds

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

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

The solar siblings in the Gaia era

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

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Multi-scale and multi-domain computational astrophysics

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

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Radial Migration of the Sun in the Milky Way: a Statistical Study

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

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Star Formation in the vicinity of Nuclear Black Holes: Young Stellar Objects close to Sgr A*

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

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

The growth of massive stars via stellar collisions in ensemble star clusters

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

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

The moment of core collapse in star clusters with a mass function

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

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

The secular evolution of the Kuiper belt after a close stellar encounter

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

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Binary population synthesis and SNIa rates

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

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Progenitors of Supernovae Type Ia

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

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Young Stellar Objects close to Sgr A*

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

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

High-redshift formation and evolution of central massive objects II: The census of BH seeds

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

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

How well do STARLAB and NBODY compare? II: Hardware and accuracy

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

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Reconstructing the Arches I: Constraining the Initial Conditions

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

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

The Evolution of Protoplanetary Disks in the Arches Cluster

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

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

The formation of planets in circumbinary disks

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

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Dynamics in Young Star Clusters: From Planets to Massive Stars

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

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Imprint of the merger and ring-down on the gravitational wave background from black hole binaries coalescence

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

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

The evolution of embedded star clusters

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

0 citations0 reviewsNo code linkedOpen access
preprint2009arXiv

High-velocity runaway stars from three-body encounters

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

0 citations0 reviewsNo code linkedOpen access
preprint2009arXiv

How well do STARLAB and NBODY4 compare? I: Simple models

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

0 citations0 reviewsNo code linkedOpen access
preprint2009arXiv

On the origin of high-velocity runaway stars

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

0 citations0 reviewsNo code linkedOpen access
preprint2009arXiv

The Living Application: a Self-Organising System for Complex Grid Tasks

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

0 citations0 reviewsNo code linkedOpen access
preprint2009arXiv

The lost siblings of the Sun

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

0 citations0 reviewsNo code linkedOpen access
preprint2002arXiv

Selection of the SIM Astrometric Grid

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

0 citations0 reviewsNo code linkedOpen access