M. B. N. Kouwenhoven

The primordial binary population. II. Recovering the binary population for intermediate mass stars in Scorpius OB2

We characterize the binary population in the young and nearby OB association Scorpius OB2 (Sco 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 3σ 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 fq(q) ∝ q γq , with γq ≈− 0.4. Sco OB2 has a semi-major axis distribution of the form fa(a) ∝ a γ a with γa ≈− 1.0 (Opik’s law), in the range 5R� < a < 5 × 10 6 R� . The log-normal period distribution of Duquennoy & Mayor (1991) 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.

Dynamical Mass Segregation on a Very Short Timescale

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

The primordial binary population. I: A near-infrared adaptive optics search for close visual companions to A star members of Scorpius OB2

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 K S 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). The angular separation between primaries and observed companion stars ranges from 0.22 to 12.4. At the mean distance of Sco OB2 (130 pc) this corresponds to a projected separation of 28.6 AU to 1612 AU. Absolute magnitudes are derived for all primaries and observed companions using the parallax and interstellar extinction for each star individually. For each object we derive the mass from K S , assuming an age of 5 Myr for the US subgroup, and 20 Myr for the UCL and LCC subgroups. Companion star masses range from 0.10 M ○. to 3.0 M ○. . The mass ratio distribution follows f(q) = q -Γ with r = 0.33, which excludes random pairing. No close (ρ ≤ 3.75) companion stars or background stars are found in the magnitude range 12 mag 14 mag 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 is a lower limit since the data from the observations and from literature are hampered by observational biases and selection effects. This paper is the first step toward our goal to derive the primordial binary population in Sco OB2.

Do binaries in clusters form in the same way as in the field

We examine the dynamical destruction of binary systems in star clusters of different densities. We find that at high densities (10 4 -10 5 M ⊙ pc -3 ) almost all binaries with separations > 10 3 au are destroyed after a few crossing times. At low densities [O(10 2 )M ⊙ )pc -3 ], many binaries with separations > 1 0 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 ∼100 times denser in the past with a half-mass radius of only 0.1-0.2 pc as (i) it is expanding, (ii) it has very few binaries with separations > 1 0 3 au and (iii) it is well mixed and therefore dynamically old. We also examine the origin of the field binary population. Binaries with separations 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.

The dragon simulations: globular cluster evolution with a million stars

Introducing the dragon simulation project, we present direct N-body simulations of four massive globular clusters (GCs) with 106 stars and 5 per cent 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 (BHs), 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, colour–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 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 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. (Less)

Spitzer Mid-IR Spectroscopy of Powerful 2?Jy and 3CRR Radio Galaxies. I. Evidence against a Strong Starburst-AGN Connection in Radio-loud AGN

We present deep Spitzer/Infrared Spectrograph (IRS) spectra for complete samples of 46 2 Jy radio galaxies (0.05 75%) than their more extended counterparts (≈15%-25%). We discuss this result in the context of a possible bias toward the selection of compact radio sources triggered in gas-rich environments.

Sequential star formation in rcw 34: a spectroscopic census of the stellar content of high-mass star-forming regions

In this paper, 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 in 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, of which the most massive star has spectral type O8.5V. Intermediate-mass stars (2–3 M� ) are detected of G- and K-spectral type. These stars are still in the premain-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 Myr is derived from the properties of the PMS stars inside the Hii region. Between the class II sources in the bubble and the PMS stars in the Hii region, no age difference could be detected with the present data. The presence of the class 0/I sources in the molecular cloud, however, suggests that the objects inside the molecular cloud are significantly younger. The most likely scenario for the formation of the three regions is that star formation propagated 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 similar to a champagne flow. More recently, star formation occurred in the rest of the molecular cloud. Two different formation scenarios are possible. (1) 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 its turn induces the current star formation in the molecular cloud. (2) An external triggering is responsible for the star formation propagating from south to north.

nbody6++gpu: ready for the gravitational million-body problem

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

The dynamical evolution of multiplanet systems in open clusters

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.

Spitzer Mid-IR Spectroscopy of Powerful 2Jy and 3CRR Radio Galaxies. II. AGN Power Indicators and Unification

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] λ25.89 μm and [O IV] λ25.89 μm fine-structure lines, the optical [O III] λ5007 forbidden line, and mid-IR 24 μm, 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 μm luminosities of the broad- and narrow-line AGNs in our samples at fixed [O IV] luminosity, we show that the [O III] and 24 μm 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 μm 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 μm 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.