Source author record

S. F. Portegies Zwart

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

12works
7topics
4close collaborators

Actions

Connect this record

Log in to claim

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

12 published item(s)

preprint2020arXiv

Mergers of Equal-Mass Binaries with Compact Object Companions from Mass Transfer in Triple Star Systems

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

preprint2015arXiv

The initial conditions of observed star clusters - I. Method description and validation

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

preprint2014arXiv

A Keplerian-based Hamiltonian Splitting for Gravitational $N$-body Simulations

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

preprint2014arXiv

Quest for finding the lost siblings of the Sun

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

preprint2013arXiv

Population synthesis of ultracompact X-ray binaries in the Galactic Bulge

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

preprint2013arXiv

The Astrophysical Multipurpose Software Environment

We present the open source Astrophysical Multi-purpose Software Environment (AMUSE, www.amusecode.org), a component library for performing astrophysical simulations involving different physical domains and scales. It couples existing codes within a Python framework based on a communication layer using MPI. The interfaces are standardized for each domain and their implementation based on MPI guarantees that the whole framework is well-suited for distributed computation. It includes facilities for unit handling and data storage. Currently it includes codes for gravitational dynamics, stellar evolution, hydrodynamics and radiative transfer. Within each domain the interfaces to the codes are as similar as possible. We describe the design and implementation of AMUSE, as well as the main components and community codes currently supported and we discuss the code interactions facilitated by the framework. Additionally, we demonstrate how AMUSE can be used to resolve complex astrophysical problems by presenting example applications.

preprint2012arXiv

Asymmetric supernova in hierarchical multiple star systems and application to J1903+0327

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

preprint2009arXiv

On the velocity dispersion of young star clusters: super-virial or binaries?

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

preprint2006arXiv

The effect of giant molecular clouds on star clusters

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

preprint2001arXiv

Population synthesis for double white dwarfs. II. Semi-detached systems: AM CVn stars

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

preprint2001arXiv

The gravitational wave signal from the Galactic disk population of binaries containing two compact objects

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

preprint2000arXiv

Evolution of Neutron-Star, Carbon-Oxygen White-Dwarf Binaries

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