sort results by

Use logical operators AND, OR, NOT and round brackets to construct complex queries. Whitespace-separated words are treated as ANDed.

Show articles per page in mode

Kouwenhoven, M. B. N.

Normalized to: Kouwenhoven, M.

56 article(s) in total. 116 co-authors, from 1 to 16 common article(s). Median position in authors list is 2,0.

[1]  [pdf] - 2061750
Linking the formation and fate of exo-Kuiper belts within solar system analogues
Comments: Accepted for publication in MNRAS
Submitted: 2020-02-19
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.
[2]  [pdf] - 1958096
Planets In Young Massive Clusters: On the survivability of planets in young massive clusters and its implication of planet orbital architectures in globular clusters
Comments: 12 pages, 11 figures, 1 table, accepted by MNRAS
Submitted: 2019-03-06, last modified: 2019-09-09
As of August 2019, among the more than 4000 confirmed exoplanets, only one has been detected in a globular cluster (GC) M4. The scarce of exoplanet detections motivates us to employ direct $N$-body simulations to investigate the dynamical stability of planets in young massive clusters (YMCs), which are potentially the progenitors of GCs. In an $N=128{\rm k}$ cluster of virial radius 1.7 pc (comparable to Westerlund-1), our simulations show that most wide-orbit planets ($a\geq 20$~au) will be ejected within a timescale of 10 Myr. Interestingly, more than $70\%$ of planets with $a<5$~au survive in the 100 Myr simulations. Ignoring planet-planet scattering and tidal damping, the survivability at $t$ Myr as a function of initial semi-major axis $a_0$ in au in such a YMC can be described as $f_{\rm surv}(a_0, t)=-0.33 \log_{10}(a_0) \left(1 - e^{-0.0482t} \right) + 1$. Upon ejection, about $28.8\%$ of free-floating planets (FFPs) have sufficient speeds to escape from the host cluster at a crossing timescale. The other FFPs will remain bound to the cluster potential, but the subsequent dynamical evolution of the stellar system can result in the delayed ejection of FFPs from the host cluster. Although a full investigation of planet population in GCs requires extending the simulations to multi-Gyr, our results suggest that wide-orbit planets and free-floating planets are unlikely to be found in GCs.
[3]  [pdf] - 1948854
Planetary systems in a star cluster I: the Solar system scenario
Comments: 20 pages, 46 figures, accepted for pubblication on Monthly Notices of the Royal Astronomical Society on date 13 August 2019
Submitted: 2019-08-21
Young stars are mostly found in dense stellar environments, and even our own Solar system may have formed in a star cluster. Here, we numerically explore the evolution of planetary systems similar to our own Solar system in star clusters. We investigate the evolution of planetary systems in star clusters. Most stellar encounters are tidal, hyperbolic, and adiabatic. A small fraction of the planetary systems escape from the star cluster within 50 Myr; those with low escape speeds often remain intact during and after the escape process. While most planetary systems inside the star cluster remain intact, a subset is strongly perturbed during the first 50 Myr. Over the course of time, 0.3 % - 5.3 % of the planets escape, sometimes up to tens of millions of years after a stellar encounter occurred. Survival rates are highest for Jupiter, while Uranus and Neptune have the highest escape rates. Unless directly affected by a stellar encounter itself, Jupiter frequently serves as a barrier that protects the terrestrial planets from perturbations in the outer planetary system. In low-density environments, Jupiter provides protection from perturbations in the outer planetary system, while in high-density environments, direct perturbations of Jupiter by neighbouring stars is disruptive to habitable-zone planets. The diversity amongst planetary systems that is present in the star clusters at 50 Myr, and amongst the escaping planetary systems, is high, which contributes to explaining the high diversity of observed exoplanet systems in star clusters and in the Galactic field
[4]  [pdf] - 1860049
Kepler-411: a four-planet system with an active host star
Comments: 11 pages, 9 figures
Submitted: 2019-02-25
We present a detailed characterization of the Kepler-411 system (KOI 1781). This system was previously known to host two transiting planets: one with a period of 3 days ($R=2.4R_\oplus$; Kepler-411b) and one with a period of 7.8 days ($R=4.4R_\oplus$; Kepler-411c), as well as a transiting planetary candidate with a 58-day period ($R=3.3R_\oplus$; KOI 1781.03) from Kepler photometry. Here, we combine Kepler photometry data and new transit timing variation (TTV) measurements from all the Kepler quarters with previous adaptive-optics imaging results, and dynamical simulations, in order to constrain the properties of the Kepler-411 system. From our analysis, we obtain masses of 25.6$\pm$2.6$M_\oplus$ for Kepler-411b and 26.4$\pm$5.9 $M_\oplus$ for Kepler-411c, and we confirm the planetary nature of KOI 1781.03 with a mass of 15.2$\pm$5.1$M_\oplus$, hence the name Kepler-411d. Furthermore, by assuming near-coplanarity of the system (mutual inclination below $30^\circ$), we discover a nontransiting planet, Kepler-411e, with a mass of 10.8$\pm$1.1$M_\oplus$ on a 31.5-day orbit, which has a strong dynamical interaction with Kepler-411d. With densities of $1.71\pm0.39$~g\,cm$^{-3}$ and $2.32\pm0.83$~g\,cm$^{-3}$, both Kepler-411c and Kepler-411d belong to the group of planets with a massive core and a significant fraction of volatiles. Although Kepler-411b has a sub-Neptune size, it belongs to the group of rocky planets.
[5]  [pdf] - 1860030
The substructure and halo population of the Double Cluster $h$ and $\chi$ Persei
Comments: 9 pagges, 9 figures, Accecpted to A&A
Submitted: 2019-02-18
In order to study the stellar population and possible substructures in the outskirts of Double Cluster $h$ and $\chi$ Persei, we investigate using the GAIA DR2 data a sky area of about 7.5 degrees in radius around the Double Cluster cores. We identify member stars using various criteria, including their kinematics (viz, proper motion), individual parallaxes, as well as photometric properties. A total of 2186 member stars in the parameter space were identified as members. Based on the spatial distribution of the member stars, we find an extended halo structure of $h$ and $\chi$ Persei, about 6 - 8 times larger than their core radii. We report the discovery of filamentary substructures extending to about 200 pc away from the Double Cluster. The tangential velocities of these distant substructures suggest that they are more likely to be the remnants of primordial structures, instead of a tidally disrupted stream from the cluster cores. Moreover, the internal kinematic analysis indicates that halo stars seems to be experiencing a dynamic stretching in the RA direction, while the impact of the core components is relatively negligible. This work also suggests that the physical scale and internal motions of young massive star clusters may be more complex than previously thought.
[6]  [pdf] - 1791836
Planetary Systems in Star Clusters: the dynamical evolution and survival
Comments: 2 pages, 2 figures, Proceedings for Symposia IAU 345 "Origins: from the Protosun to the First Steps of Life" at the XXX general assembly in Vienna
Submitted: 2018-11-30
Most stars, perhaps even all stars, form in crowded stellar environments. Such star forming regions typically dissolve within ten million years, while others remain bound as stellar groupings for hundreds of millions to billions of years, and then become the open clusters or globular clusters that are present in our Milky Way galaxy today. A large fraction of stars in the Galaxy hosts planetary companions. To understand the origin and dynamical evolution of such exoplanet systems, it is necessary to carefully study the effect of their environments. Here, we combine theoretical estimates with state-of-the-art numerical simulations of evolving planetary systems similar to our own solar system in different star cluster environments. We combine the REBOUND planetary system evolution code, and the NBODY6++GPU star cluster evolution code, integrated in the AMUSE multi-physics environment. With our study we can constrain the effect of external perturbations of different environments on the planets and debris structures of a wide variety of planetary systems, which may play a key role for the habitability of exoplanets in the Universe.
[7]  [pdf] - 1588682
Clearing residual planetesimals by sweeping secular resonances in transitional disks: a lone-planet scenario for the wide gaps in debris disks around Vega and Fomalhaut
Comments: 20 pages, 12 figures. Accepted for publication in ApJ
Submitted: 2017-09-21
Extended gaps in the debris disks of both Vega and Fomalhaut have been observed. These structures have been attributed to tidal perturbations by multiple super-Jupiter gas giant planets. Within the current observational limits, however, no such massive planets have been detected. Here we propose a less stringent `lone-planet' scenario to account for the observed structure with a single eccentric gas giant and suggest that clearing of these wide gaps is induced by its sweeping secular resonance. During the depletion of the disk gas, the planet's secular resonance propagates inward and clears a wide gap over an extended region of the disk. Although some residual intermediate-size planetesimals may remain in the gap, their surface density is too low to either produce super-Earths or lead to sufficiently frequent disruptive collisions to generate any observable dusty signatures. The main advantage of this lone-planet sweeping-secular-resonance model over the previous multiple gas giant tidal truncation scenario is the relaxed requirement on the number of gas giants. The observationally inferred upper mass limit can also be satisfied provided the hypothetical planet has a significant eccentricity. A significant fraction of solar or more massive stars bear gas giant planets with significant eccentricities. If these planets acquired their present-day kinematic properties prior to the depletion of their natal disks, their sweeping secular resonance would effectively impede the retention of neighboring planets and planetesimals over a wide range of orbital semi-major axes.
[8]  [pdf] - 1586502
The dynamical origin of multiple populations in intermediate-age clusters in the Magellanic Clouds
Comments: 11 pages, 8 figures, 1 table; accepted for publication in MNRAS
Submitted: 2017-07-28
Numerical simulations were carried out to study the origin of multiple stellar populations in the intermediate-age clusters NGC 411 and NGC 1806 in the Magellanic Clouds. We performed NBODY6++ simulations based on two different formation scenarios, an ad hoc formation model where second-generation (SG) stars are formed inside a cluster of first-generation (FG) stars using the gas accumulated from the external intergalactic medium and a minor merger model of unequal mass ($M_{SG}$/$M_{FG}$ ~5-10%) clusters with an age difference of a few hundred million years. We compared our results such as the radial profile of the SG-to-FG number ratio with observations on the assumption that the SG stars in the observations are composed of cluster members, and confirmed that both the ad hoc formation and merger scenarios reproduce the observed radial trend of the SG-to-FG number ratio which shows less centrally concentrated SG than FG stars. It is difficult to constrain the formation scenario for the multiple populations by only using the spatial distribution of the SG stars. SG stars originating from the merger scenario show a significant velocity anisotropy and rotational features compared to those from the ad hoc formation scenario. Thus, observations aimed at kinematic properties like velocity anisotropy or rotational velocities for SG stars should be obtained to better understand the formation of the multiple populations in these clusters. This is, however, beyond current instrumentation capabilities.
[9]  [pdf] - 1584577
Stability of Multiplanetary Systems in Star Clusters
Comments: 19 pages, 13 figures, 4 tables, accepted for publication in MNRAS
Submitted: 2017-06-12
Most stars form in star clusters and stellar associated. To understand the roles of star cluster environments in shaping the dynamical evolution of planetary systems, we carry out direct $N$-body simulations of four planetary systems models in three different star cluster environments with respectively N=2k, 8k and 32k stars. In each cluster, an ensemble of initially identical planetary systems are assigned to solar-type stars with $\sim 1 M_{\odot}$ and evolved for 50~Myr. We found that following the depletion of protoplanetary disks, external perturbations and planet-planet interactions are two driving mechanisms responsible for the destabilization of planetary systems. The planet survival rate varies from $\sim 95\%$ in the N=2k cluster to $\sim 60\%$ in the N=32k cluster, which suggests that most planetary systems can indeed survive in low-mass clusters, except in the central regions. We also find that planet ejections through stellar encounters are cumulative processes, as only $\sim 3\%$ of encounters are strong enough to excite the eccentricity by $\Delta e \geq 0.5$. Short-period planets can be perturbed through orbit crossings with long-period planets. When taking into account planet-planet interactions, the planet ejection rate nearly doubles, and therefore multiplicity contributes to the vulnerability of planetary systems. In each ensemble, $\sim 0.2\%$ of planetary orbits become retrograde due to random directions of stellar encounters. Our results predict that young low-mass star clusters are promising sites for next-generation planet surveys, yet low planet detection rates are expected in dense globular clusters such as 47 Tuc. Nevertheless, planets in denser stellar environments are likely to have shorter orbital periods, which enhances their detectability.
[10]  [pdf] - 1542773
Planetesimal clearing and size-dependent asteroid retention by secular resonance sweeping during the depletion of the solar nebula
Comments: 26 pages, 18figures, accepted for publication in the ApJ
Submitted: 2016-10-30
The distribution of heavy elements is anomalously low in the asteroid main belt region compared with elsewhere in the solar system. Observational surveys also indicate a deficit in the number of small ($ \le 50$~km size) asteroids that is two orders of magnitude lower than what is expected from the single power-law distribution that results from a collisional coagulation and fragmentation equilibrium. Here, we consider the possibility that a major fraction of the original asteroid population may have been cleared out by Jupiter's secular resonance, as it swept through the main asteroid belt during the depletion of the solar nebula. This effect leads to the excitation of the asteroids' orbital eccentricities. Concurrently, hydrodynamic drag and planet-disk tidal interaction effectively damp the eccentricities of sub-100 km-size and of super-lunar-size planetesimals, respectively. These combined effects lead to the asteroids' orbital decay and clearing from the present-day main belt region ($\sim 2.1-3.3$~AU). The intermediate-size (50 to several hundreds of kilometers) planetesimals therefore preferentially remain as main belt asteroids near their birthplaces, with modest asymptotic eccentricities. The smaller asteroids are the fragments of subsequent disruptive collisions at later times as suggested by the present-day asteroid families. This scenario provides a natural explanation for both the observed low surface density and the size distribution of asteroids in the main belt. It also offers an explanation for the confined spatial extent of the terrestrial planet building blocks without the requirement of extensive migration of Jupiter.
[11]  [pdf] - 1474412
Planetary systems in star clusters
Comments: Proceedings of Cosmic-Lab: Star Clusters as Cosmic Laboratories for Astrophysics, Dynamics, and Fundamental Physics - MODEST16
Submitted: 2016-09-04
Thousands of confirmed and candidate exoplanets have been identified in recent years. Consequently, theoretical research on the formation and dynamical evolution of planetary systems has seen a boost, and the processes of planet-planet scattering, secular evolution, and interaction between planets and gas/debris disks have been well-studied. Almost all of this work has focused on the formation and evolution of isolated planetary systems, and neglect the effect of external influences, such as the gravitational interaction with neighbouring stars. Most stars, however, form in clustered environments that either quickly disperse, or evolve into open clusters. Under these conditions, young planetary systems experience frequent close encounters with other stars, at least during the first 1-10 Myr, which affects planets orbiting at any period range, as well as their debris structures.
[12]  [pdf] - 1580195
The long-term dynamical evolution of disc-fragmented multiple systems in the Solar Neighborhood
Submitted: 2016-09-01
The origin of very low-mass hydrogen-burning stars, brown dwarfs, and planetary-mass objects at the low-mass end of the initial mass function is not yet fully understood. Gravitational fragmentation of circumstellar discs provides a possible mechanism for the formation of such low-mass objects. The kinematic and binary properties of very low-mass objects formed through disc fragmentation at early times (< 10 Myr) were discussed in Li et al. (2015). In this paper we extend the analysis by following the long-term evolution of disc-fragmented systems, up to an age of 10 Gyr, covering the ages of the stellar and substellar population in the Galactic field. We find that the systems continue to decay, although the rates at which companions escape or collide with each other are substantially lower than during the first 10 Myr, and that dynamical evolution is limited beyond 1 Gyr. By t = 10 Gyr, about one third of the host stars is single, and more than half have only one companion left. Most of the other systems have two companions left that orbit their host star in widely separated orbits. A small fraction of companions have formed binaries that orbit the host star in a hierarchical triple configuration. The majority of such double companion systems have internal orbits that are retrograde with respect to their orbits around their host stars. Our simulations allow a comparison between the predicted outcomes of disc-fragmentation with the observed low-mass hydrogen-burning stars, brown dwarfs, and planetary-mass objects in the Solar neighborhood. Imaging and radial velocity surveys for faint binary companions among nearby stars are necessary for verification or rejection for the formation mechanism proposed in this paper.
[13]  [pdf] - 1381572
The DRAGON simulations: globular cluster evolution with a million stars
Comments: 18 pages, 8 figures, MNRAS accepted
Submitted: 2016-02-01
Introducing the DRAGON simulation project, we present direct $N$-body simulations of four massive globular clusters (GCs) with $10^6$ stars and 5$\%$ primordial binaries at a high level of accuracy and realism. The GC evolution is computed with NBODY6++GPU and follows the dynamical and stellar evolution of individual stars and binaries, kicks of neutron stars and black holes, and the effect of a tidal field. We investigate the evolution of the luminous (stellar) and dark (faint stars and stellar remnants) GC components and create mock observations of the simulations (i.e. photometry, color-magnitude diagrams, surface brightness and velocity dispersion profiles). By connecting internal processes to observable features we highlight the formation of a long-lived 'dark' nuclear subsystem made of black holes (BHs), which results in a two-component structure. The inner core is dominated by the BH subsystem and experiences a core collapse phase within the first Gyr. It can be detected in the stellar (luminous) line-of-sight velocity dispersion profiles. The outer extended core - commonly observed in the (luminous) surface brightness profiles - shows no collapse features and is continuously expanding. We demonstrate how a King (1966) model fit to observed clusters might help identify the presence of post core-collapse BH subsystems. For global observables like core and half-mass radii the direct simulations agree well with Monte-Carlo models. Variations in the initial mass function can result in significantly different GC properties (e.g. density distributions) driven by varying amounts of early mass loss and the number of forming BHs.
[14]  [pdf] - 1359170
The dynamical fate of binary star clusters in the Galactic tidal field
Comments: 14 pages, 7 figures, accepted for publication in MNRAS
Submitted: 2016-01-07
Fragmentation and fission of giant molecular clouds occasionally results in a pair of gravitationally bound star clusters that orbit their mutual centre of mass for some time, under the influence of internal and external perturbations. We investigate the evolution of binary star clusters with different orbital configurations, with a particular focus on the Galactic tidal field. We carry out $N$-body simulations of evolving binary star clusters and compare our results with estimates from our semi-analytic model. The latter accounts for mass loss due to stellar evolution and two-body relaxation, and for evolution due to external tides. Using the semi-analytic model we predict the long-term evolution for a wide range of initial conditions. It accurately describes the global evolution of such systems, until the moment when a cluster merger is imminent. $N$-body simulations are used to test our semi-analytic model and also to study additional features of evolving binary clusters, such as the kinematics of stars, global cluster rotation, evaporation rates, and the cluster merger process. We find that the initial orientation of a binary star cluster with respect to the Galactic field, and also the initial orbital phase, are crucial for its fate. Depending on these properties, the binaries may experience orbital reversal, spiral-in, or vertical oscillation about the Galactic plane before they actually merge at $t\approx100$~Myr, and produce rotating star clusters with slightly higher evaporation rates. The merger process of a binary cluster induces an outburst that ejects $\sim10\%$ of the stellar members into the Galactic field.
[15]  [pdf] - 1277167
The dynamical fate of planetary systems in young star clusters
Comments: 14 pages, 9 figures, accepted for publication in MNRAS
Submitted: 2015-08-06
We carry out N-body simulations to examine the effects of dynamical interactions on planetary systems in young open star clusters. We explore how the planetary populations in these star clusters evolve, and how this evolution depends on the initial amount of substructure, the virial ratio, the cluster mass and density, and the initial semi-major axis of the planetary systems. The fraction of planetary systems that remains intact as a cluster member, fbps, is generally well-described by the functional form fbps=f0(1+[a/a0]^c)^-1, where (1-f0) is the fraction of stars that escapes from the cluster, a0 the critical semi-major axis for survival, and c a measure for the width of the transition region. The effect of the initial amount of substructure over time can be quantified as fbps=A(t)+B(D), where A(t) decreases nearly linearly with time, and B(D) decreases when the clusters are initially more substructured. Provided that the orbital separation of planetary systems is smaller than the critical value a0, those in clusters with a higher initial stellar density (but identical mass) have a larger probability of escaping the cluster intact. These results help us to obtain a better understanding of the difference between the observed fractions of exoplanets-hosting stars in star clusters and in the Galactic field. It also allows us to make predictions about the free-floating planet population over time in different stellar environments.
[16]  [pdf] - 1288404
Block Time Step Storage Scheme for Astrophysical N-body Simulations
Comments: 13 pages, 10 figures, 5 tables, Accepted for publication in The Astrophysical Journal Supplement Series
Submitted: 2015-06-24
Astrophysical research in recent decades has made significant progress thanks to the availability of various $N$-body simulation techniques. With the rapid development of high-performance computing technologies, modern simulations have been able to take the computing power of massively parallel clusters with more than $10^5$ GPU cores. While unprecedented accuracy and dynamical scales have been achieved, the enormous amount of data being generated continuously poses great challenges for the subsequent procedures of data analysis and archiving. As an urgent response to these challenges, in this paper we propose an adaptive storage scheme for simulation data, inspired by the block time step integration scheme found in a number of direct $N$-body integrators available nowadays. The proposed scheme, namely the block time step storage scheme, works by minimizing the data redundancy with assignments of data with individual output frequencies as required by the researcher. As demonstrated by benchmarks, the proposed scheme is applicable to a wide variety of simulations. Despite the main focus of developing a solution for direct $N$-body simulation data, the methodology is transferable for grid-based or tree-based simulations where hierarchical time stepping is used.
[17]  [pdf] - 1155002
The dynamical evolution of low-mass hydrogen-burning stars, brown dwarfs and planetary-mass objects formed through disc fragmentation
Comments: 38 pages, 15 figures
Submitted: 2015-06-10
Theory and simulations suggest that it is possible to form low-mass hydrogen-burning stars, brown dwarfs and planetary-mass objects via disc fragmentation. As disc fragmentation results in the formation of several bodies at comparable distances to the host star, their orbits are generally unstable. Here, we study the dynamical evolution of these objects. We set up the initial conditions based on the outcomes of the SPH simulations of Stamatellos & Whitworth, and for comparison we also study the evolution of systems resulting from lower-mass fragmenting discs. We refer to these two sets of simulations as set 1 and set 2. At 10 Myr, approximately half of the host stars have one companion left, and approximately 22% (set 1) to 9.8% (set 2) of the host stars are single. Systems with multiple secondaries in relatively stable configurations are common (about 30% and 44%, respectively). The majority of the companions are ejected within 1 Myr with velocities mostly below 5 km/s, with some runaway escapers with velocities over 30 km/s. About 6% (set 1) and 2% (set 2) of the companions pair up into very low-mass binary systems. The majority of these pairs escape as very low-mass binaries, while others remain bound to the host star in hierarchical configurations (often with retrograde inner orbits). Physical collisions with the host star (0.43 and 0.18 events per host star for set 1 and set 2) and between companions (0.08 and 0.04 events per host star for set 1 and set 2) are relatively common and their frequency increases with increasing disc mass. Our study predicts observable properties of very low-mass binaries, low-mass hierarchical systems, the brown dwarf desert, and free-floating brown dwarfs and planetary-mass objects in and near young stellar groupings, which can be used to distinguish between different formation scenarios of very low-mass stars, brown dwarfs and planetary-mass objects.
[18]  [pdf] - 1280830
NBODY6++GPU: Ready for the gravitational million-body problem
Comments: 13 pages, 9 figures, 3 tables
Submitted: 2015-04-14, last modified: 2015-05-21
Accurate direct $N$-body simulations help to obtain detailed information about the dynamical evolution of star clusters. They also enable comparisons with analytical models and Fokker-Planck or Monte-Carlo methods. NBODY6 is a well-known direct $N$-body code for star clusters, and NBODY6++ is the extended version designed for large particle number simulations by supercomputers. We present NBODY6++GPU, an optimized version of NBODY6++ with hybrid parallelization methods (MPI, GPU, OpenMP, and AVX/SSE) to accelerate large direct $N$-body simulations, and in particular to solve the million-body problem. We discuss the new features of the NBODY6++GPU code, benchmarks, as well as the first results from a simulation of a realistic globular cluster initially containing a million particles. For million-body simulations, NBODY6++GPU is $400-2000$ times faster than NBODY6 with 320 CPU cores and 32 NVIDIA K20X GPUs. With this computing cluster specification, the simulations of million-body globular clusters including $5\%$ primordial binaries require about an hour per half-mass crossing time.
[19]  [pdf] - 1045014
Close encounters involving free-floating planets in star clusters
Comments: Accepted for publication in MNRAS; 18 pages, 12 figures
Submitted: 2015-03-10
Instabilities in planetary systems can result in the ejection of planets from their host system, resulting in free-floating planets (FFPs). If this occurs in a star cluster, the FFP may remain bound to the star cluster for some time and interact with the other cluster members until it is ejected. Here, we use $N$-body simulations to characterise close star-planet and planet-planet encounters and the dynamical fate of the FFP population in star clusters containing $500-2000$ single or binary star members. We find that FFPs ejected from their planetary system at low velocities typically leave the star cluster 40% earlier than their host stars, and experience tens of close ($<1000$ AU) encounters with other stars and planets before they escape. The fraction of FFPs that experiences a close encounter depends on both the stellar density and the initial velocity distribution of the FFPs. Approximately half of the close encounters occur within the first 30 Myr, and only 10% occur after 100 Myr. The periastron velocity distribution for all encounters is well-described by a modified Maxwell-Bolzmann distribution, and the periastron distance distribution is linear over almost the entire range of distances considered, and flattens off for very close encounters due to strong gravitational focusing. Close encounters with FFPs can perturb existing planetary systems and their debris structures, and they can result in re-capture of FFPs. In addition, these FFP populations may be observed in young star clusters in imaging surveys; a comparison between observations and dynamical predictions may provide clues to the early phases of stellar and planetary dynamics in star clusters.
[20]  [pdf] - 1579584
Planetary Systems in Star Clusters
Comments: Paper submitted for proceedings of IAU Symp. No. 312, 2 pages, 3 figures
Submitted: 2015-01-07, last modified: 2015-01-27
In the solar neighborhood, where the typical relaxation timescale is larger than the cosmic age, at least 10\% to 15\% of Sun-like stars have planetary systems with Jupiter-mass planets. In contrast, dense star clusters, charactered by frequent close encounters, have been found to host very few planets. We carry out numerical simulations with different initial conditions to investigate the dynamical stability of planetary systems in star cluster environments.
[21]  [pdf] - 1216821
How does a low-mass cut-off in the stellar IMF affect the evolution of young star clusters?
Comments: 13 pages, 7 figures, accepted for publication in MNRAS
Submitted: 2014-09-09
We investigate how different stellar initial mass functions (IMFs) can affect the mass loss and survival of star clusters. We find that IMFs with radically different low-mass cut-offs (between 0.1 and 2 Msun) do not change cluster destruction time-scales as much as might be expected. Unsurprisingly, we find that clusters with more high-mass stars lose relatively more mass through stellar evolution, but the response to this mass loss is to expand and hence significantly slow their dynamical evolution. We also argue that it is very difficult, if not impossible, to have clusters with different IMFs that are initially "the same", since the mass, radius and relaxation times depend on each other and on the IMF in a complex way. We conclude that changing the IMF to be biased towards more massive stars does speed up mass loss and dissolution, but that it is not as dramatic as might be thought.
[22]  [pdf] - 1209329
Spitzer Mid-IR Spectroscopy of Powerful 2Jy and 3CRR Radio Galaxies. II. AGN Power Indicators and Unification
Comments: Matches the ApJ accepted version
Submitted: 2014-05-04, last modified: 2014-08-26
It remains uncertain which continuum and emission line diagnostics best indicate the bolometric powers of active galactic nuclei (AGNs), especially given the attenuation caused by the circumnuclear material and the possible contamination by components related to star formation. Here we use mid-IR spectra along with multiwavelength data to investigate the merit of various diagnostics of AGN radiative power, including the mid-IR [Ne III] lambda25.89 mum and [O IV] lambda25.89 mum fine-structure lines, the optical [O III] lambda5007 forbidden line, and mid-IR 24 mum, 5 GHz radio, and X-ray continuum emission, for complete samples of 46 2Jy radio galaxies (0.05 < z < 0.7) and 17 3CRR FRII radio galaxies (z < 0.1). We find that the mid-IR [O IV] line is the most reliable indicator of AGN power for powerful radio-loud AGNs. By assuming that the [O IV] is emitted isotropically, and comparing the [O III] and 24 mum luminosities of the broad- and narrow-line AGNs in our samples at fixed [O IV] luminosity, we show that the [O III] and 24 mum emission are both mildly attenuated in the narrow-line compared to the broad-line objects by a factor of #2. However, despite this attenuation, the [O III] and 24 mum luminosities are better AGN power indicators for our sample than either the 5 GHz radio or the X-ray continuum luminosities. We also detect the mid-IR 9.7 mum silicate feature in the spectra of many objects but not ubiquitously: at least 40% of the sample shows no clear evidence for these features. We conclude that, for the majority of powerful radio galaxies, the mid-IR lines are powered by AGN photoionization.
[23]  [pdf] - 977824
The Link Between Ejected Stars, Hardening and Eccentricity Growth of Super Massive Black Holes in Galactic Nuclei
Comments: 22 pages, 8 figures, accepted for publication in ApJ
Submitted: 2013-11-18
The hierarchical galaxy formation picture suggests that super massive black holes (MBHs) observed in galactic nuclei today have grown from coalescence of massive black hole binaries (MBHB) after galaxy merging. Once the components of a MBHB become gravitationally bound, strong three-body encounters between the MBHB and stars dominate its evolution in a "dry" gas free environment, and change the MBHB's energy and angular momentum (semi-major axis, eccentricity and orientation). Here we present high accuracy direct N-body simulations of spherical and axisymmetric (rotating) galactic nuclei with order a million stars and two massive black holes that are initially unbound. We analyze the properties of the ejected stars due to slingshot effects from three-body encounters with the MBHB in detail. Previous studies have investigated the eccentricity and energy changes of MBHs using approximate models or Monte-Carlo three body scatterings. We find general agreement with the average results of previous semi-analytic models for spherical galactic nuclei, but our results show a large statistical variation. Our new results show many more phase space details of how the process works, and also show the influence of stellar system rotation on the process. We detect that the angle between the orbital plane of the MBHBs and that of the stellar system (when it rotates) influences the phase-space properties of the ejected stars. We also find that massive MBHB tend to switch stars with counter-rotating orbits into co-rotating orbits during their interactions.
[24]  [pdf] - 1166453
The dynamical evolution of multi-planet systems in open clusters
Comments: 12 pages, 7 figures, accepted for publication in MNRAS
Submitted: 2013-05-07
The majority of stars form in star clusters and many are thought to have planetary companions. We demonstrate that multi-planet systems are prone to instabilities as a result of frequent stellar encounters in these star clusters much more than single-planet systems. The cumulative effect of close and distant encounters on these planetary systems are investigated using Monte Carlo scattering experiments. We consider two types of planetary configurations orbiting Sun-like stars: (i) five Jupiter-mass planets in the semi-major axis range 1-42 AU orbiting a Solar mass star, with orbits that are initially co-planar, circular, and separated by 10 mutual Hill radii, and (ii) the four gas giants of our Solar system. Planets with short orbital periods are not directly affected by encountering stars. However, secular evolution of perturbed systems may result in the ejection of the innermost planets or in physical collisions of the innermost planets with the host star, up to many thousands of years after a stellar encounter. The simulations of the Solar system-like systems indicate that Saturn, Uranus and Neptune are affected by both direct interactions with encountering stars, as well as planet-planet scattering. Jupiter, on the other hand, is almost only affected by direct encounters with neighbouring stars, as its mass is too large to be substantially perturbed by the other three planets. Our results indicate that stellar encounters can account for the apparent scarcity of exoplanets in star clusters, not only for those on wide-orbit that are directly affected by stellar encounters, but also planets close to the star which can disappear long after a stellar encounter has perturbed the planetary system.
[25]  [pdf] - 1158923
Gravitational conundrum? Dynamical mass segregation versus disruption of binary stars in dense stellar systems
Comments: Accepted for publication in The Astrophysical Journal; 19 pages in AASTeX format; 3 figures
Submitted: 2013-01-09
Upon their formation, dynamically cool (collapsing) star clusters will, within only a few million years, achieve stellar mass segregation for stars down to a few solar masses, simply because of gravitational two-body encounters. Since binary systems are, on average, more massive than single stars, one would expect them to also rapidly mass segregate dynamically. Contrary to these expectations and based on high-resolution Hubble Space Telescope observations, we show that the compact, 15-30 Myr-old Large Magellanic Cloud cluster NGC 1818 exhibits tantalizing hints at the >= 2 sigma level of significance (> 3 sigma if we assume a power-law secondary-to-primary mass-ratio distribution) of an increasing fraction of F-star binary systems (with combined masses of 1.3-1.6 Msun) with increasing distance from the cluster center, specifically between the inner 10 to 20" (approximately equivalent to the cluster's core and half-mass radii) and the outer 60 to 80". If confirmed, this will offer support of the theoretically predicted but thus far unobserved dynamical disruption processes of the significant population of 'soft' binary systems---with relatively low binding energies compared to the kinetic energy of their stellar members---in star clusters, which we have access to here by virtue of the cluster's unique combination of youth and high stellar density.
[26]  [pdf] - 1116562
On the origin of planets at very wide orbits from the re-capture of free floating planets
Comments: ApJ, in press. Added two figures
Submitted: 2012-02-10, last modified: 2012-03-13
In recent years several planets have been discovered at wide orbits (>100 AU) around their host stars. Theoretical studies encounter difficulties in explaining their formation and origin. Here we propose a novel scenario for the production of planetary systems at such orbits, through the dynamical recapture of free floating planets (FFPs) in dispersing stellar clusters. This process is a natural extension of the recently suggested scenario for the formation of wide stellar binaries. We use N-body simulations of dispersing clusters with 10-1000 stars and comparable numbers of FFPs to study this process. We find that planets are captured into wide orbits in the typical range ~100-10^6 AU, and have a wide range of eccentricities (thermal distribution). Typically, 3-6 x (f_FFP/1) % of all stars capture a planetary companion with such properties (where f_FFP is the number of FFP per star). The planetary capture efficiency is comparable to that of capture-formed stellar-binaries, and shows a similar dependence on the cluster size and structure. It is almost independent of the specific planetary mass; planets as well as sub-stellar companions of any mass can be captured. The capture efficiency decreases with increasing cluster size, and for a given cluster size the it increases with the host/primary mass. More than one planet can be captured around the same host and planets can be captured into binary systems. Planets can also be captured into pre-existing planetary and into orbits around black holes and massive white dwarfs, if these formed early enough before the cluster dispersal. In particular, stellar black holes have a high capture efficiency (>50 % and 5-10 x (f_FFP/1) % for capture of stars and planetary companions, respectively) due to their large mass. Finally, although rare, two FFPs or brown dwarfs can become bound and form a FFP-binary system with no stellar host.
[27]  [pdf] - 1091785
Spitzer mid-IR spectroscopy of powerful 2Jy and 3CRR radio galaxies. I. Evidence against a strong starburst-AGN connection in radio-loud AGN
Comments: 33 pages, 14 Figures, Accepted for publication in ApJ
Submitted: 2011-11-18, last modified: 2011-11-22
We present deep Spitzer/IRS spectra for complete samples of 46 2Jy radio galaxies (0.05<z<0.7) and 19 3CRR FRII radio galaxies (z<0.1), and use the detection of polycyclic aromatic hydrocarbon (PAH) features to examine the incidence of contemporaneous star formation and radio-loud AGN activity. Our analysis reveals PAH features in only a minority (30%) of the objects with good IRS spectra. Using the wealth of complementary data available for the 2Jy and 3CRR samples we make detailed comparisons between a range of star formation diagnostics: optical continuum spectroscopy, mid- to far-IR (MFIR) color, far-IR excess and PAH detection. There is good agreement between the various diagnostic techniques: most candidates identified to have star formation activity on the basis of PAH detection are also identified using at least two of the other techniques. We find that only 35% of the combined 2Jy and 3CRR sample show evidence for recent star formation activity (RSFA) at optical and/or MFIR wavelengths. This result argues strongly against the idea of a close link between starburst and powerful radio-loud AGN activity, reinforcing the view that, although a large fraction of powerful radio galaxies may be triggered in galaxy interactions, only a minority are triggered at the peaks of star formation activity in major, gas-rich mergers. However, we find that compact radio sources (D < 15 kpc) show a significantly higher incidence of RSFA (>75%) than their more extended counterparts (=15 -- 25%). We discuss this result in the context of a possible bias towards the selection of compact radio sources triggered in gas-rich environments.
[28]  [pdf] - 1083998
The Helium Abundance in the Ejecta of U Scorpii
Comments: Accepted for publication in MNRAS
Submitted: 2011-09-12, last modified: 2011-09-13
U Scorpii is a recurrent nova which has been observed in outburst on 10 occasions, most recently in 2010. We present near-infrared and optical spectroscopy of the 2010 outburst of U Sco. The reddening of U Sco is found to be $E(B-V) = 0.14\pm0.12$, consistent with previous determinations, from simultaneous optical and near-IR observations. The spectra show the evolution of the line widths and profiles to be consistent with previous outbursts. Velocities are found to be up to 14000\,kms$^{-1}$ in broad components and up to 1800\,kms$^{-1}$ in narrow line components, which become visible around day 8 due to changes in the optical depth. From the spectra we derive a helium abundance of $N$(He)/$N$(H)$ = 0.073\pm0.031$ from the most reliable lines available; this is lower than most other estimates and indicates that the secondary is not helium-rich, as previous studies have suggested.
[29]  [pdf] - 401456
The origin of very wide binary stars
Comments: 5 pages, 3 figures, to be published in the proceedings of the 9th Pacific-Rim Conference on Stellar Astrophysics
Submitted: 2011-08-23
A large population of fragile, wide (> 1000 AU) binary systems exists in the Galactic field and halo. These wide binary stars cannot be primordial because of the high stellar density in star forming regions, while formation by capture in the Galactic field is highly improbable. We propose that these binary systems were formed during the dissolution phase of star clusters (see Kouwenhoven et al. 2010, for details). Stars escaping from a dissolving star cluster can have very similar velocities, which can lead to the formation of a wide binary systems. We carry out N-body simulations to test this hypothesis. The results indicate that this mechanism explains the origin of wide binary systems in the Galaxy. The resulting wide binary fraction and semi-major axis distribution depend on the initial conditions of the dissolving star cluster, while the distributions in eccentricity and mass ratio are universal. Finally, since most stars are formed in (relatively tight) primordial binaries, we predict that a large fraction of the wide "binary stars" are in fact higher-order multiple systems.
[30]  [pdf] - 1025171
Sequential Star Formation in RCW 34: A Spectroscopic Census of the Stellar Content of High-mass Star-forming Regions
Comments: 19 pages, 11 figures, accepted by ApJ
Submitted: 2010-02-12
We present VLT/SINFONI integral field spectroscopy of RCW 34 along with Spitzer/IRAC photometry of the surroundings. RCW 34 consists of three different regions. A large bubble has been detected on the IRAC images in which a cluster of intermediate- and low-mass class II objects is found. At the northern edge of this bubble, an HII region is located, ionized by 3 OB stars. Intermediate mass stars (2 - 3 Msun) are detected of G- and K- spectral type. These stars are still in the pre-main sequence (PMS) phase. North of the HII region, a photon-dominated region is present, marking the edge of a dense molecular cloud traced by H2 emission. Several class 0/I objects are associated with this cloud, indicating that star formation is still taking place. The distance to RCW 34 is revised to 2.5 +- 0.2 kpc and an age estimate of 2 - 1 Myrs is derived from the properties of the PMS stars inside the HII region. The most likely scenario for the formation of the three regions is that star formation propagates from South to North. First the bubble is formed, produced by intermediate- and low-mass stars only, after that, the HII region is formed from a dense core at the edge of the molecular cloud, resulting in the expansion as a champagne flow. More recently, star formation occurred in the rest of the molecular cloud. Two different formation scenarios are possible: (a) The bubble with the cluster of low- and intermediate mass stars triggered the formation of the O star at the edge of the molecular cloud which in turn induces the current star-formation in the molecular cloud. (b) An external triggering is responsible for the star-formation propagating from South to North. [abridged]
[31]  [pdf] - 1024761
The formation of very wide binaries during the star cluster dissolution phase
Comments: 14 pages, 12 figures, accepted by MNRAS
Submitted: 2010-01-22
Over the past few decades, numerous wide (>1000 au) binaries in the Galactic field and halo have been discovered. Their existence cannot be explained by the process of star formation or by dynamical interactions in the field, and their origin has long been a mystery. We explain the origin of these wide binaries by formation during the dissolution phase of young star clusters: an initially unbound pair of stars may form a binary when their distance in phase-space is small. Using N-body simulations, we find that the resulting wide binary fraction in the semi-major axis range 1000 au - 0.1 pc for individual clusters is 1-30%, depending on the initial conditions. The existence of numerous wide binaries in the field is consistent with observational evidence that most clusters start out with a large degree of substructure. The wide binary fraction decreases strongly with increasing cluster mass, and the semi-major axis of the newly formed binaries is determined by the initial cluster size. The resulting eccentricity distribution is thermal, and the mass ratio distribution is consistent with gravitationally-focused random pairing. As a large fraction of the stars form in primordial binaries, we predict that a large number of the observed 'wide binaries' are in fact triple or quadruple systems. By integrating over the initial cluster mass distribution, we predict a binary fraction of a few per cent in the semi-major axis range 1000 au - 0.1 pc in the Galactic field, which is smaller than the observed wide binary fraction. However, this discrepancy may be solved when we consider a broad range of cluster morphologies.
[32]  [pdf] - 1018153
The origin of very wide binary systems
Comments: 4 pages, 1 figure, to appear in the proceedings of IAU Symposium 266, eds. R. de Grijs & J.R.D. Lepine
Submitted: 2009-10-30
The majority of stars in the Galactic field and halo are part of binary or multiple systems. A significant fraction of these systems have orbital separations in excess of thousands of astronomical units, and systems wider than a parsec have been identified in the Galactic halo. These binary systems cannot have formed through the 'normal' star-formation process, nor by capture processes in the Galactic field. We propose that these wide systems were formed during the dissolution phase of young star clusters. We test this hypothesis using N-body simulations of evolving star clusters and find wide binary fractions of 1-30%, depending on initial conditions. Moreover, given that most stars form as part of a binary system, our theory predicts that a large fraction of the known wide 'binaries' are, in fact, multiple systems.
[33]  [pdf] - 1018063
What does the IMF really tell us about star formation?
Comments: To be published in the proceedings of IAU Symposium 262 "Stellar Populations - Planning for the Next Decade"
Submitted: 2009-10-22
Obtaining accurate measurements of the initial mass function (IMF) is often considered to be the key to understanding star formation, and a universal IMF is often assumed to imply a universal star formation process. Here, we illustrate that different modes of star formation can result in the same IMF, and that, in order to truly understand star formation, a deeper understanding of the primordial binary population is necessary. Detailed knowledge on the binary fraction, mass ratio distribution, and other binary parameters, as a function of mass, is a requirement for recovering the star formation process from stellar population measurements.
[34]  [pdf] - 28072
The formation of very wide binaries
Comments: 4 pages, to appear in the proceedings of "BINKEY: Binaries - a Key to Comprehension of the Universe", held 8-12 June 2009 in Brno, Czech Republic
Submitted: 2009-09-07
Over the last decades, numerous wide (>1000 AU) binaries have been discovered in the Galactic field and halo. The origin of these wide binaries cannot be explained by star formation or by dynamical interactions in the Galactic field. We explain their existence by wide binary formation during the dissolution phase of young star clusters. In this scenario, two single stars that leave the dissolving cluster at the same time, in the same direction, and with similar velocities, form a new, very wide binary. Using N-body simulations we study how frequently this occurs, and how the orbital parameters of such binaries depend on the properties of the cluster from which they originate. The resulting wide binary fraction for individual star clusters is 1-30%, depending on the initial conditions. As most stars form as part of a binary or multiple system, we predict that a large fraction of these wide binaries are in fact wide triple and quadruple systems.
[35]  [pdf] - 25683
Dynamical mass segregation on a very short timescale
Comments: ApJ Letter, accepted
Submitted: 2009-06-26
We discuss the observations and theory of star cluster formation to argue that clusters form dynamically cool (subvirial) and with substructure. We then perform an ensemble of simulations of cool, clumpy (fractal) clusters and show that they often dynamically mass segregate on timescales far shorter than expected from simple models. The mass segregation comes about through the production of a short-lived, but very dense core. This shows that in clusters like the Orion Nebula Cluster the stars >4 Msun can dynamically mass segregate within the current age of the cluster. Therefore, the observed mass segregation in apparently dynamically young clusters need not be primordial, but could be the result of rapid and violent early dynamical evolution.
[36]  [pdf] - 1002217
Do binaries in clusters form in the same way as in the field?
Comments: 12 pages, 6 figures, accepted for publication in MNRAS
Submitted: 2009-05-13
We examine the dynamical destruction of binary systems in star clusters of different densities. We find that at high densities (10^4 - 10^5 Msun pc^-3) almost all binaries with separations > 10^3 AU are destroyed after a few crossing times. At low densities (order(10^2) Msun pc^-3) many binaries with separations > 10^3 AU are destroyed, and no binaries with separations > 10^4 AU survive after a few crossing times. Therefore the binary separations in clusters can be used as a tracer of the dynamical age and past density of a cluster. We argue that the central region of the Orion Nebula Cluster was around 100 times denser in the past with a half-mass radius of only 0.1 - 0.2 pc as (a) it is expanding, (b) it has very few binaries with separations > 10^3 AU, and (c) it is well-mixed and therefore dynamically old. We also examine the origin of the field binary population. Binaries with separations < 10^2 AU are not significantly modified in any cluster, therefore at these separations the field reflects the sum of all star formation. Binaries with separations in the range 10^2 - 10^4 AU are progressively more and more heavily affected by dynamical disruption in increasingly dense clusters. If most star formation is clustered, these binaries must be over-produced relative to the field. Finally, no binary with a separation > 10^4 AU can survive in any cluster and so must be produced by isolated star formation, but only if all isolated star formation produces extremely wide binaries.
[37]  [pdf] - 1002066
What does a universal IMF imply about star formation?
Comments: 6 pages, 4 figures. MNRAS, in press
Submitted: 2009-04-29
We show that the same initial mass function (IMF) can result from very different modes of star formation from very similar underlying core and/or system mass functions. In particular, we show that the canonical IMF can be recovered from very similar system mass functions, but with very different mass ratio distributions within those systems. This is a consequence of the basically log-normal shapes of all of the distributions. We also show that the relationships between the shapes of the core, system, and stellar mass functions may not be trivial. Therefore, different star formation in different regions could still result in the same IMF.
[38]  [pdf] - 20319
Using the minimum spanning tree to trace mass segregation
Comments: Accepted in MNRAS
Submitted: 2009-01-14
We present a new method to detect and quantify mass segregation in star clusters. It compares the minimum spanning tree (MST) of massive stars with that of random stars. If mass segregation is present, the MST length of the most massive stars will be shorter than that of random stars. This difference can be quantified (with an associated significance) to measure the degree of mass segregation. We test the method on simulated clusters in both 2D and 3D and show that the method works as expected. We apply the method to the Orion Nebula Cluster (ONC) and show that the method is able to detect the mass segregation in the Trapezium with a `mass segregation ratio' \Lambda_{MSR}=8.0 \pm 3.5 (where \Lambda_{MSR}=1 is no mass segregation) down to 16 \Msun, and also that the ONC is mass segregated at a lower level (~2.0 \pm 0.5) down to 5 \Msun. Below 5 \Msun we find no evidence for any further mass segregation in the ONC.
[39]  [pdf] - 18687
Pairing mechanisms for binary stars
Comments: 4 pages, 2 figures, in: Galactic and stellar dynamics in the era of high-resolution surveys, Boily C., Combes F., Hensler G., eds., Strasbourg (France), March 2008, in press (Astron. Nachr.)
Submitted: 2008-11-19
Knowledge of the binary population in stellar groupings provides important information about the outcome of the star forming process in different environments. Binarity is also a key ingredient in stellar population studies and is a prerequisite to calibrate the binary evolution channels. In these proceedings we present an overview of several commonly used methods to pair individual stars into binary systems, which we refer to as the pairing function. Many pairing functions are frequently used by observers and computational astronomers, either for the mathematical convenience, or because they roughly describe the expected outcome of the star forming process. We discuss the consequences of each pairing function for the interpretation of observations and numerical simulations. The binary fraction and mass ratio distribution generally depend strongly on the selection of the range in primary spectral type in a sample. These quantities, when derived from a binary survey with a mass-limited sample of target stars, are thus not representative for the population as a whole.
[40]  [pdf] - 18639
Exploring the consequences of pairing algorithms for binary stars
Comments: 43 pages, 18 figures, accepted for publication in A&A
Submitted: 2008-11-18
Knowledge of the binary population in stellar groupings provides important information about the outcome of the star forming process in different environments (see, e.g., Blaauw 1991, and references therein). Binarity is also a key ingredient in stellar population studies, and is a prerequisite to calibrate the binary evolution channels. In this paper we present an overview of several commonly used methods to pair individual stars into binary systems, which we refer to as pairing functions. These pairing functions are frequently used by observers and computational astronomers, either for their mathematical convenience, or because they roughly describe the expected outcome of the star forming process. We discuss the consequences of each pairing function for the interpretation of observations and numerical simulations. The binary fraction and mass ratio distribution generally depend strongly on the selection of the range in primary spectral type in a sample. The mass ratio distribution and binary fraction derived from a binarity survey among a mass-limited sample of targets is thus not representative for the population as a whole. Neither theory nor observations indicate that random pairing of binary components from the mass distribution, the simplest pairing function, is realistic. It is more likely that companion stars are formed in a disk around a star, or that a pre-binary core fragments into two binary components. The results of our analysis are important for (i) the interpretation of the observed mass ratio distribution and binary fraction for a sample of stars, (ii) a range of possible initial condition algorithms for star cluster simulations, and (iii) how to discriminate between the different star formation scenarios.
[41]  [pdf] - 10932
Open cluster stability and the effects of binary stars
Comments: 9 pages, 3 figures, accepted for publication in A&A (replacement of v1 in which our results were based on a circular argument)
Submitted: 2008-03-13, last modified: 2008-10-17
The diagnostic age versus mass-to-light ratio diagram is often used in attempts to constrain the shape of the stellar initial mass function (IMF), and the potential longevity of extragalactic young to intermediate-age massive star clusters. Here, we explore its potential for Galactic open clusters. On the basis of a small, homogenised cluster sample we provide useful constraints on the presence of significant binary fractions. Using the massive young Galactic cluster Westerlund 1 as a key example, we caution that stochasticity in the IMF introduces significant additional uncertainties. We conclude that for an open cluster to survive for any significant length of time, and in the absence of substantial external perturbations, it is a necessary but not a sufficient condition to be located close to or (in the presence of a significant binary population) somewhat BELOW the predicted photometric evolutionary sequences for "normal" simple stellar populations (although such a location may be dominated by a remaining "bound" cluster core and thus not adequately reflect the overall cluster dynamics).
[42]  [pdf] - 17062
What do dynamical cluster masses really tell us about dynamics?
Comments: 4 pages (incl. 3 postscript figures), in: Galactic and stellar dynamics in the era of high-resolution surveys, Boily C., Combes F., Hensler G., eds., Strasbourg (France), March 2008, in press (Astron. Nachr.)
Submitted: 2008-10-03
The diagnostic age versus mass-to-light ratio diagram is often used in attempts to constrain the shape of the stellar initial mass function, and the stability and the potential longevity of extragalactic young to intermediate-age massive star clusters. Here, we explore the pitfalls associated with this approach and its potential for use with Galactic open clusters. We conclude that for an open cluster to survive for any significant fraction of a Hubble time (in the absence of substantial external perturbations), it is a necessary but not a sufficient condition to be located close to the predicted photometric evolutionary sequences for "normal" simple stellar populations.
[43]  [pdf] - 11258
Probing the Early Evolution of Young High-Mass Stars
Comments: To appear in ASP Conf. Proceedings of "Massive Star Formation: Observations confront Theory", H. Beuther et al. (eds.), held in Heidelberg, September 2007
Submitted: 2008-03-27
Near-infrared imaging surveys of high-mass star-forming regions reveal an amazingly complex interplay between star formation and the environment (Churchwell et al. 2006; Alvarez et al. 2004). By means of near-IR spectroscopy the embedded massive young stars can be characterized and placed in the context of their birth site. However, so far spectroscopic surveys have been hopelessly incomplete, hampering any systematic study of these very young massive stars. New integral field instrumentation available at ESO has opened the possibility to take a huge step forward by obtaining a full spectral inventory of the youngest massive stellar populations in star-forming regions currently accessible. Simultaneously, the analysis of the extended emission allows the characterization of the environmental conditions. The Formation and Early Evolution of Massive Stars (FEMS) collaboration aims at setting up a large observing campaign to obtain a full census of the stellar content, ionized material, outflows and PDR's over a sample of regions that covers a large parameter space. Complementary radio, mm and infrared observations will be used for the characterization of the deeply embedded population. For the first eight regions we have obtained 40 hours of SINFONI observations. In this contribution, we present the first results on three regions that illustrate the potential of this strategy.
[44]  [pdf] - 9394
How do binaries affect the derived dynamical mass of a star cluster?
Comments: 6 pages, 5 figures, to appear in ``Young Massive Star Clusters - Initial Conditions and Environments'', 2008, Astrophysics & Space Science, eds. E. Perez, R. de Grijs, R. M. Gonzalez Delgado
Submitted: 2008-01-25
The dynamical mass of a star cluster can be derived from the virial theorem, using the measured half-mass radius and line-of-sight velocity dispersion of the cluster. However, this dynamical mass may be a significant overestimation of the cluster mass if the contribution of the binary orbital motion is not taken into account. In these proceedings we describe the mass overestimation as a function of cluster properties and binary population properties, and briefly touch the issue of selection effects. We find that for clusters with a measured velocity dispersion of sigma > 10 km/s the presence of binaries does not affect the dynamical mass significantly. For clusters with sigma < 1 km/s (i.e., low-density clusters), the contribution of binaries to sigma is significant, and may result in a major dynamical mass overestimation. The presence of binaries may introduce a downward shift of Delta log(L/Mdyn) = 0.05-0.4 in the log(L/Mdyn) vs. age diagram.
[45]  [pdf] - 8003
The effect of binaries on the dynamical mass determination of star clusters
Comments: 12 pages, 10 figures, accepted by A&A
Submitted: 2007-12-11
The total mass of distant star clusters is often derived from the virial theorem, using line-of-sight velocity dispersion measurements and half-light radii. Although most stars form in binary systems, this is mostly ignored when interpreting the observations. The components of binary stars exhibit orbital motion, which may increase the measured velocity dispersion, and may therefore result in a dynamical mass overestimation. In this paper we quantify the effect of neglecting the binary population on the derivation of the dynamical mass of a star cluster. We simulate star clusters numerically, and study the dependence of the derived dynamical mass on the properties of the binary population. We find that the presence of binaries plays a crucial role for very sparse clusters with a stellar density comparable to that of the field star population (~0.1 stars/pc3), as the velocity dispersion is fully dominated by the binary orbital motion. For such clusters, the dynamical mass may overestimate the true mass by up to an order of magnitude. For very dense clusters (>10^7 stars/pc3), binaries do not affect the dynamical mass estimation significantly. For clusters of intermediate density (0.1-10^7 stars/pc3), the dynamical mass can be overestimated by 10-100%, depending on the properties of the binary population.
[46]  [pdf] - 5729
Binaries and the dynamical mass of star clusters
Comments: Comments: 2 pages. Conference proceedings for IAUS246 'Dynamical Evolution of Dense Stellar Systems', ed. E. Vesperini (Chief Editor), M. Giersz, A. Sills, Capri, Sept. 2007
Submitted: 2007-10-05
The total mass of a distant star cluster is often derived from the virial theorem, using line-of-sight velocity dispersion measurements and half-light radii, under the implicit assumption that all stars are single (although it is known that most stars form part of binary systems). The components of binary stars exhibit orbital motion, which increases the measured velocity dispersion, resulting in a dynamical mass overestimation. In this article we quantify the effect of neglecting the binary population on the derivation of the dynamical mass of a star cluster. We find that the presence of binaries plays an important role for clusters with total mass M < 10^5 Msun; the dynamical mass can be significantly overestimated (by a factor of two or more). For the more massive clusters, with Mcl > 10^5 Msun, binaries do not affect the dynamical mass estimation significantly, provided that the cluster is significantly compact (half-mass radius < 5 pc).
[47]  [pdf] - 3221
The primordial binary population II: Recovering the binary population for intermediate mass stars in Sco OB2
Comments: 36 pages, 11 figures, accepted by A&A
Submitted: 2007-07-18
We characterize the binary population in the young and nearby OB association Scorpius OB2 using available observations of visual, spectroscopic, and astrometric binaries with intermediate-mass primaries. We take into account observational biases by comparing the observations with simulated observations of model associations. The available data indicate a large binary fraction (> 70% with 3sigma confidence), with a large probability that all intermediate mass stars in Sco OB2 are part of a binary system. The binary systems have a mass ratio distribution of the form f(q) ~ q^-0.4. Sco OB2 has a semi-major axis distribution of the form f(log a) ~ constant (Opik's law), in the range 5-5e6 Rsun. The log-normal period distribution of Duquennoy & Mayor results in too few spectroscopic binaries, even if the model binary fraction is 100%. Sco OB2 is a young association with a low stellar density; its current population is expected to be very similar to the primordial population. The fact that practically all stars in Sco OB2 are part of a binary (or multiple) system demonstrates that multiplicity is a fundamental factor in the star formation process, at least for intermediate mass stars.
[48]  [pdf] - 675
Optical BVI Imaging and HI Synthesis Observations of the Dwarf Irregular Galaxy ESO 364-G 029
Comments: 10 pages, 11 figures, accepted by A&A
Submitted: 2007-04-23
As part of an effort to enlarge the number of well-studied Magellanic-type galaxies, we obtained broadband optical imaging and neutral hydrogen radio synthesis observations of the dwarf irregular galaxy ESO 364-G 029. The optical morphology characteristically shows a bar-like main body with a one-sided spiral arm, an approximately exponential light distribution, and offset photometric and kinematic centers. The HI distribution is mildly asymmetric and, although slightly offset from the photometric center, roughly follows the optical brightness distribution, extending to over 1.2 Holmberg radii (where mu_B = 26.5 mag/arcsec^2). In particular, the highest HI column densities closely follow the bar, one-arm spiral, and a third optical extension. The rotation is solid-body in the inner parts but flattens outside of the optical extent. The total HI flux F_HI = 23.1 pm 1.2 Jy km/s, yielding a total HI mass M_HI= (6.4 pm 1.7) x 10^8 Msun (for a distance D = 10.8 pm 1.4 Mpc) and a total HI mass-to-blue-luminosity ratio M_HI/L_B = (0.96 pm 0.14) Msun / Lsun,B (distance independent). The HI data suggest a very complex small-scale HI structure, with evidence of large shells and/or holes, but deeper observations are required for a detailed study. Follow-up observations are also desirable for a proper comparison with the Large Magellanic Cloud, where despite an optical morphology very similar to ESO 364-G 029 the HI bears little resemblance to the optical.
[49]  [pdf] - 1234014
A brown dwarf desert for intermediate mass stars in Sco OB2?
Comments: 29 pages, 9 figures, accepted by A&A
Submitted: 2006-11-29
We present JHK observations of 22 intermediate-mass stars in Sco OB2, obtained with VLT/NACO. The survey was performed to determine the status of (sub)stellar candidate companions of A and late-B members. The distinction between companions and background stars is by a comparison with isochrones and statistical arguments. We are sensitive to companions in the separation range 0.1''-11'' (13-1430 AU) and K<17. We detect 62 secondaries of which 18 are physical companions (3 new), 11 candidates, and 33 background stars. The companion masses are in the range 0.03<M<1.19 Msun, with mass ratios 0.06<q<0.55. We include in our sample a subset of 9 targets with multi-color ADONIS observations from Kouwenhoven et al. (2005). In the ADONIS survey secondaries with K<12 were classified as companions; those with K>12 as background stars. Our multi-color analysis demonstrates that the simple K=12 criterion correctly classifies the secondaries in ~80% of the cases. We reanalyse the total ADONIS/NACO sample and conclude that of the 176 secondaries, 25 are physical companions, 55 are candidates, and 96 are background stars. Although we are sensitive and complete to brown dwarfs as faint as K=14 in the separation range 130-520 AU, we detect only one, giving a brown dwarf companion fraction of 0.5% (M>30 MJ). However, the number of brown dwarfs is consistent with an extrapolation of the stellar companion mass distribution. This indicates that the physical mechanism for the formation of brown dwarfs around intermediate mass stars is similar to that of stellar companions, and that the embryo ejection mechanism does not need to be invoked in order to explain the small number of brown dwarf companions among these stars.
[50]  [pdf] - 86216
The primordial binary population in the association Sco OB2
Comments: PhD thesis, University of Amsterdam, 2006 (very low resolution). Full-resolution PDF available at
Submitted: 2006-10-26
Observations over the last decade have indicated that a large fraction of the stars are part of a binary or multiple system. For our understanding of star formation it is therefore of crucial importance to characterise the outcome of the star forming process in terms of binary parameters. This thesis describes the recovery of the primordial binary population in the nearby OB association Sco OB2. The current binary population in Sco OB2 is first recovered using two adaptive optics surveys (ADONIS, NAOS-CONICA), a literature study, and a detailed analysis of the selection effects of visual, spectroscopic, and astrometric binary surveys. Our results indicate a binary fraction close to 100%. The mass ratio distribution (among A/B primaries) has the form f(q) = q^-0.4, while random pairing is excluded. The semi-major axis distribution has the form f(log a) = constant. Due to its youth and low stellar density, the current binary population of Sco OB2 is very similar to its primordial binary population. Our study further indicates a small brown dwarf companion frequency and a small substellar-to-stellar companion frequency among A and B type stars. These properties, often referred to as the brown dwarf desert, are a natural result of the mass ratio distribution in Sco OB2. The embryo ejection scenario is not necessary to explain observations. The brown dwarf desert may be ascribed to an excess of planetary companions, rather than by a lack of brown dwarf companions.
[51]  [pdf] - 67956
The Primordial Binary Population - I: A near-infrared adaptive optics search for close visual companions to A star members of Scorpius OB2
Comments: 27 pages, to accepted by A&A
Submitted: 2004-10-05
We present the results of a near-infrared adaptive optics survey with the aim to detect close companions to Hipparcos members in the three subgroups of the nearby OB association Sco OB2: Upper Scorpius (US), Upper Centaurus Lupus (UCL) and Lower Centaurus Crux (LCC). We have targeted 199 A-type and late B-type stars in the Ks band, and a subset also in the J and H band. We find 151 stellar components other than the target stars. A brightness criterion is used to separate these components into 77 background stars and 74 candidate physical companion stars. Out of these 74 candidate companions, 41 have not been reported before (14 in US; 13 in UCL; 14 in LCC). Companion star masses range from 0.1 to 3 Msun. The mass ratio distribution follows f(q) = q^-0.33, which excludes random pairing. No close (rho < 3.75'') companion stars or background stars are found in the magnitude range 12 < Ks < 14. The lack of stars with these properties cannot be explained by low-number statistics, and may imply a lower limit on the companion mass of ~ 0.1 Msun. Close stellar components with Ks > 14 are observed. If these components are very low-mass companion stars, a gap in the companion mass distribution might be present. The small number of close low-mass companion stars could support the embryo-ejection formation scenario for brown dwarfs. Our findings are compared with and complementary to visual, spectroscopic, and astrometric data on binarity in Sco OB2. We find an overall companion star fraction of 0.52 in this association. This paper is the first step toward our goal to derive the primordial binary population in Sco OB2.
[52]  [pdf] - 56117
The Primordial Binary Population in OB Associations
Comments: 2 pages, 1 figure, poster paper to appear in proceedings of IAU Coll. 191 "The environments and evolution of binary and multiple stars"
Submitted: 2003-04-11
For understanding the process of star formation it is essential to know how many stars are formed as singles or in multiple systems, as a function of environment and binary parameters. This requires a characterization of the primordial binary population, which we define as the population of binaries that is present just after star formation has ceased, but before dynamical and stellar evolution have significantly altered its characteristics. In this article we present the first results of our adaptive optics survey of 200 (mainly) A-type stars in the nearby OB association Sco OB2. We report the discovery of 47 new candidate companions of Sco OB2 members. The next step will be to combine these observations with detailed simulations of young star clusters, in order to find the primordial binary population.
[53]  [pdf] - 47971
Null-induced mode changes in PSR B0809+74
Comments: 10 pages, 10 low-res figures, published in Astronomy & Astrophysics. High-res versions available at [pdf (2.7M), ps (4.2M)]. (v2: removed typos)
Submitted: 2002-02-26, last modified: 2002-05-02
We have found that there are two distinct emission modes in PSR B0809+74. Beside its normal and most common mode, the pulsar emits in a significantly different quasi-stable mode after most or possibly all nulls, occasionally for over 100 pulses. In this mode the pulsar is brighter, the subpulse separation is less, the subpulses drift more slowly and the pulse window is shifted towards earlier longitudes. We can now account for several previously unexplained phenomena associated with the nulling-drifting interaction: the unexpected brightness of the first active pulse and the low post-null driftrate. We put forward a new interpretation of the subpulse-position jump over the null, indicating that the speedup time scale of the post-null drifting is much shorter than previously thought. The speedup time scale we find is no longer discrepant with the time scales found for the subpulse-drift slowdown and the emission decay around the null.
[54]  [pdf] - 45789
Pulsar "Drifting"-Subpulse Polarization: No Evidence for Systematic Polarization-Angle Rotations
Comments: Accepted for publication in Astronomy & Astrophysics
Submitted: 2001-11-01
Polarization-angle density displays are given for pulsars B0809+74 and B2303+30, which exhibit no evidence of the systematic polarization-angle rotation within individual subpulses previously reported for these two stars. The ``drifting'' subpulses of both pulsars exhibit strikingly linear and circular polarization which appears to reflect the characteristics of two nearly orthogonally polarized emission ``modes''--along which the severe average-profile depolarization that is characteristic of their admixture at comparable overall intensities.
[55]  [pdf] - 1232482
Pulsars in the Westerbork Northern Sky Survey
Comments: 12 pages, 9 figures. Accepted for publication in A&A
Submitted: 2000-05-26
I have searched the Westerbork Northern Sky Survey (WENSS) source list for detections of known radio pulsars. A source with a flux density greater than five times the local noise level is found near the positions of 25 pulsars. The probability that one out of these 25 sources is a chance coincidence is about 10%. I have looked at the WENSS maps of the non-detected pulsars. A flux density between three and five times the local noise level is found near the positions of 14 of these non-detected pulsars. There is a 50 percent probability that (at least) one of these marginal detections is just a noise fluctuation. Fourteen radio pulsars, which according to earlier flux measurements have flux densities above three times the WENSS noise level, are not detected. Of the 39 pulsars detected in the WENSS 19 are also detected in the NRAO VLA Sky Survey (NVSS). By combining the WENSS and NVSS flux densities for these 19 pulsars spectral indices are obtained that differ by up to 50% from the long term averaged values reported in the literature. This affects the reliability of pulsar candidates that are selected on the basis of their WENSS-NVSS spectral index.
[56]  [pdf] - 98660
Observational Constraints on Kicks in Supernovae
Comments: 4 pages including 2 figs., LaTeX. To appear in A&A letters
Submitted: 1997-09-22
The absence or presence of extremely wide binaries with a radio pulsar and an optical counterpart imposes a strong constraint on the existence and magnitude of kicks in supernova explosions. We search for such systems by comparing the positions of radio pulsars which are not known to be in binaries with the positions of visible stars, and find that the number of associations is negligible. According to the performed population synthesis, this implies that kicks must occur, with a lower limit of at least 10 to 20 km/s.