BZPEERPapers, people, evidence and research paths
Menu

Discover

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

Research tools

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

Network

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

Opportunities

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

Account

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

Source author record

M. B. N. Kouwenhoven

M. B. N. Kouwenhoven appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record
astro-ph.GAastro-ph.SRastro-ph.EPastro-ph.COastro-ph.IMastro-ph

Catalog footprint

What is connected

31works
6topics
4close collaborators

Actions

Connect this record

Log in to claim
Open graphBrowse works

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

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

Published work

31 published item(s)

preprint2022arXiv

3D Morphology of Open Clusters in the Solar Neighborhood with Gaia EDR3 II: Hierarchical Star Formation Revealed by Spatial and Kinematic Substructures

We identify members of 65 open clusters in the solar neighborhood using the machine-learning algorithm StarGO based on Gaia EDR3 data. After adding members of twenty clusters from previous studies (Pang et al. 2021a,b; Li et al. 2021) we obtain 85 clusters, and study their morphology and kinematics. We classify the substructures outside the tidal radius into four categories: filamentary (f1) and fractal (f2) for clusters $<100$ Myr, and halo (h) and tidal-tail (t) for clusters $>100$ Myr. The kinematical substructures of f1-type clusters are elongated; these resemble the disrupted cluster Group X. Kinematic tails are distinct in t-type clusters, especially Pleiades. We identify 29 hierarchical groups in four young regions (Alessi 20, IC 348, LP 2373, LP 2442); ten among these are new. The hierarchical groups form filament networks. Two regions (Alessi 20, LP 2373) exhibit global "orthogonal" expansion (stellar motion perpendicular to the filament), which might cause complete dispersal. Infalling-like flows (stellar motion along the filament) are found in UBC 31 and related hierarchical groups in the IC 348 region. Stellar groups in the LP 2442 region (LP 2442 gp 1-5) are spatially well-mixed but kinematically coherent. A merging process might be ongoing in the LP 2442 subgroups. For younger systems ($\lesssim30$ Myr), the mean axis ratio, cluster mass and half-mass radius tend to increase with age values. These correlations between structural parameters may imply two dynamical processes occurring in the hierarchical formation scenario in young stellar groups: (1) filament dissolution and (2) sub-group mergers.

0 citations0 reviewsNo code linkedOpen access
preprint2022arXiv

Investigating the UV-excess in star clusters with $N$-body simulations: predictions for future CSST observations

We study the origin of the UV-excess in star clusters by performing N-body simulations of six clusters with N=10k and N=100k (single stars & binary systems) and metallicities of Z=0.01, 0.001, and 0.0001, using PETAR. All models initially have a 50 percent primordial binary fraction. Using GalevNB we convert the simulated data into synthetic spectra and photometry for the China Space Station Telescope (CSST) and Hubble Space Telescope (HST). From the spectral energy distributions we identify three stellar populations that contribute to the UV-excess: (1) second asymptotic giant branch stars, which contribute to the UV flux at early times; (2) naked helium stars, and (3) white dwarfs, which are long-term contributors to the FUV spectra. Binary stars consisting of a white dwarf and a main-sequence star are cataclysmic variable (CV) candidates. The magnitude distribution of CV candidates is bimodal up to 2 Gyr. The bright CV population is particularly bright in FUV-NUV. The FUV-NUV color of our model clusters is 1-2 mag redder than the UV-excess globular clusters in M 87 and in the Milky Way. This discrepancy may be induced by helium enrichment in observed clusters. Our simulations are based on simple stellar evolution; we do not include the effects of variations in helium and light elements or multiple stellar populations. A positive radial color gradient is present in CSST NUV-y for main-sequence stars of all models with a color difference of 0.2-0.5 mag, up to 4 half-mass radii. The CSST NUV-g color correlates strongly with HST FUV-NUV for NUV-g>1 mag, with the linear relation $FUV-NUV=(1.09\pm0.12)\times(NUV-g)+(-1.01\pm0.22)$. This allows for conversion of future CSST NUV-g colors into HST FUV-NUV colors, which are sensitive to UV-excess features. We find that CSST will be able to detect UV-excess in galactic/extra-galactic star clusters with ages >200 Myr.

0 citations0 reviewsNo code linkedOpen access
preprint2021arXiv

A highly mutually-inclined, compact warm-Jupiter system KOI-984 ?

The discovery of a population of close-orbiting giant planets ($\le$ 1 au) has raised a number of questions about their origins and dynamical histories. These issues have still not yet been fully resolved, despite over 20 years of exoplanet detections and a large number of discovered exoplanets. In particular, it is unclear whether warm Jupiters (WJs) form in situ, or whether they migrate from further outside and are even currently migrating to form hot Jupiters (HJs). Here, we report the possible discovery and characterization of the planets in a highly mutually-inclined ($I_{\rm mut}\simeq 45^\circ$), compact two-planet system (KOI-984), in which the newly discovered warm Jupiter KOI-984$c$ is on a 21.5-day, moderately eccentric ($e\simeq 0.4$) orbit, in addition to a previously known 4.3-day planet candidate KOI-984$b$. Meanwhile, the orbital configuration of a moderately inclined ($I_{\rm mut}\simeq 15^\circ$), low-mass ($m_{c}\simeq 24 M_{\oplus}$;$P_b\simeq 8.6$ days) perturbing planet near 1:2 mean motion resonace with KOI-984$b$ could also well reproduce observed transit timing variations and transit duration variations of KOI-984$b$. Such an eccentric WJ with a close-in sibling would pose a challenge to the proposed formation and migration mechanisms of WJs, if the first scenario is supported with more evidences in near future; this system with several other well-measured inclined WJ systems (e.g., Kepler-419 and Kepler-108) may provide additional clues for the origin and dynamical histories of WJs.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Different Fates of Young Star Clusters After Gas Expulsion

We identify structures of the young star cluster NGC 2232 in the solar neighborhood (323.0 pc), and a newly discovered star cluster LP 2439 (289.1 pc). Member candidates are identified using the Gaia DR2 sky position, parallax and proper motion data, by an unsupervised machine learning method, \textsc{StarGO}. Member contamination from the Galactic disk is further removed using the color magnitude diagram. The four identified groups (NGC 2232, LP 2439 and two filamentary structures) of stars are coeval with an age of 25 Myr and were likely formed in the same giant molecular cloud. We correct the distance asymmetry from the parallax error with a Bayesian method. The 3D morphology shows the two spherical distributions of clusters NGC 2232 and LP 2439. Two filamentary structures are spatially and kinematically connected to NGC 2232. Both NGC 2232 and LP 2439 are expanding. The expansion is more significant in LP 2439, generating a loose spatial distribution with shallow volume number and mass density profiles. The expansion is suggested to be mainly driven by gas expulsion. NGC 2232, with 73~percent of the cluster mass bound, is currently experiencing a process of re-virialization, However, LP 2439, with 52 percent cluster mass being unbound, may fully dissolve in the near future. The different survivability traces different dynamical states of NGC 2232 and LP 2439 prior to the onset of gas expulsion. NGC 2232 may have been substructured and subvirial, while LP 2439 may either have been virial/supervirial, or it has experienced a much faster rate of gas removal.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Linking the formation and fate of exo-Kuiper belts within solar system analogues

Escalating observations of exo-minor planets and their destroyed remnants both passing through the solar system and within white dwarf planetary systems motivate an understanding of the orbital history and fate of exo-Kuiper belts and planetesimal discs. Here we explore how the structure of a 40-1000 au annulus of planetesimals orbiting inside of a solar system analogue that is itself initially embedded within a stellar cluster environment varies as the star evolves through all of its stellar phases. We attempt this computationally challenging link in four parts: (1) by performing stellar cluster simulations lasting 100 Myr, (2) by making assumptions about the subsequent quiescent 11 Gyr main-sequence evolution, (3) by performing simulations throughout the giant branch phases of evolution, and (4) by making assumptions about the belt's evolution during the white dwarf phase. Throughout these stages, we estimate the planetesimals' gravitational responses to analogues of the four solar system giant planets, as well as to collisional grinding, Galactic tides, stellar flybys, and stellar radiation. We find that the imprint of stellar cluster dynamics on the architecture of $\gtrsim 100$ km-sized exo-Kuiper belt planetesimals is retained throughout all phases of stellar evolution unless violent gravitational instabilities are triggered either (1) amongst the giant planets, or (2) due to a close ($\ll 10^3$ au) stellar flyby. In the absence of these instabilities, these minor planets simply double their semimajor axis while retaining their primordial post-cluster eccentricity and inclination distributions, with implications for the free-floating planetesimal population and metal-polluted white dwarfs.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

On the survival of resonant and non-resonant planetary systems in star clusters

Despite the discovery of thousands of exoplanets in recent years, the number of known exoplanets in star clusters remains tiny. This may be a consequence of close stellar encounters perturbing the dynamical evolution of planetary systems in these clusters. Here, we present the results from direct $N$-body simulations of multiplanetary systems embedded in star clusters containing $N = 8k, 16k, 32k$, and $64k$ stars. The planetary systems, which consist of the four Solar system giant planets Jupiter, Saturn, Uranus, and Neptune, are initialized in different orbital configurations, to study the effect of the system architecture on the dynamical evolution of the entire planetary system, and on the escape rate of the individual planets. We find that the current orbital parameters of the Solar system giants (with initially circular orbits, as well as with present-day eccentricities) and a slightly more compact configuration, have a high resilience against stellar perturbations. A configuration with initial mean-motion resonances of 3:2, 3:2, and 5:4 between the planets, which is inspired by the Nice model, and for which the two outermost planets are usually ejected within the first $10^5$ yr, is in many cases stabilized due to the removal of the resonances by external stellar perturbation and by the rapid ejection of at least one planet. Assigning all planets the same mass of 1 Jovian mass almost equalizes the survival fractions. Our simulations reproduce the broad diversity amongst observed exoplanet systems. We find not only many very wide and/or eccentric orbits, but also a significant number of (stable) retrograde orbits.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Planetary Systems in a Star Cluster II: intermediate-mass black holes and planetary systems

Most stars form in dense stellar environments. It is speculated that some dense star clusters may host intermediate-mass black holes (IMBHs), which may have formed from runaway collisions between high-mass stars, or from the mergers of less massive black holes. Here, we numerically explore the evolution of populations of planets in star clusters with an IMBH. We study the dynamical evolution of single-planet systems and free-floating planets, over a period of 100~Myr, in star clusters without an IMBH, and in clusters with a central IMBH of mass $100~M_\odot$ or $200~M_\odot$. In the central region ($r\lesssim 0.2$~pc), the IMBH's tidal influence on planetary systems is typically 10~times stronger than the average neighbour star. For a star cluster with a $200M_\odot$ IMBH, the region in which the IMBH's influence is stronger within the virial radius ($\sim 1$~pc). The IMBH quenches mass segregation, and the stars in the core tend to move towards intermediate regions. The ejection rate of both stars and planets is higher when an IMBH is present. The rate at which planets are expelled from their host star rate is higher for clusters with higher IMBH masses, for $t<0.5 t_{rh}$, while remains mostly constant while the star cluster fills its Roche lobe, similar to a star cluster without an IMBH. The disruption rate of planetary systems is higher in initially denser clusters, and for wider planetary orbits, but this rate is substantially enhanced by the presence of a central IMBH.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

The long-term evolution of main-sequence binaries in DRAGON simulations

We present a comprehensive investigation of main-sequence (MS) binaries in the DRAGON simulations, which are the first one-million particles direct $N$-body simulations of globular clusters. We analyse the orbital parameters of the binary samples in two of the DRAGON simulations, D1-R7- IMF93 and D2-R7-IMF01, focusing on their secular evolution and correlations up to 12 Gyr. These two models have different initial stellar mass functions: Kroupa 1993 (D1-R7-IMF93) and Kroupa 2001 (D2-R7-IMF01); and different initial mass ratio distributions: random paring (D1-R7-IMF93) and a power-law (D1-R7-IMF93). In general, the mass ratio of a population of binaries increases over time due to stellar evolution, which is less significant in D2-R7-IMF01. In D1-R7-IMF93, primordial binaries with mass ratio $q \approx$ 0.2 are most common, and the frequency linearly declines with increasing $q$ at all times. Dynamical binaries of both models have higher eccentricities and larger semi-major axes than primordial binaries. They are preferentially located in the inner part of the star cluster. Secular evolution of binary orbital parameters does not depend on the initial mass-ratio distribution, but is sensitive to the initial binary distribution of the system. At t = 12 Gyr, the binary fraction decreases radially outwards, and mass segregation is present. A color difference of 0.1 mag in $F330W-F814W$ and 0.2 mag in $NUV-y$ between the core and the outskirts of both clusters is seen, which is a reflection of the binary radial distribution and the mass segregation in the cluster. The complete set of data for primordial and dynamical binary systems at all snapshot intervals is made publicly available.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Planetary systems in star clusters

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.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

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 $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.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

The dynamical fate of binary star clusters in the Galactic tidal field

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.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Block Time Step Storage Scheme for Astrophysical N-body Simulations

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.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Close encounters involving free-floating planets in star clusters

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.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

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

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

The dynamical evolution of low-mass hydrogen-burning stars, brown dwarfs and planetary-mass objects formed through disc fragmentation

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.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

The dynamical fate of planetary systems in young star clusters

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.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

How does a low-mass cut-off in the stellar IMF affect the evolution of young star clusters?

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.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

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] 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.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Gravitational conundrum? Dynamical mass segregation versus disruption of binary stars in dense stellar systems

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.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

The dynamical evolution of multi-planet 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.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

The Link Between Ejected Stars, Hardening and Eccentricity Growth of Super Massive Black Holes in Galactic Nuclei

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.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

On the origin of planets at very wide orbits from the re-capture of free floating planets

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.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Spitzer mid-IR spectroscopy of powerful 2Jy and 3CRR radio galaxies. I. Evidence against a strong starburst-AGN connection in radio-loud AGN

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.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

The Helium Abundance in the Ejecta of U Scorpii

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.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

The origin of very wide binary stars

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.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Sequential Star Formation in RCW 34: A Spectroscopic Census of the Stellar Content of High-mass Star-forming Regions

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]

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

The formation of very wide binaries during the star cluster dissolution phase

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.

0 citations0 reviewsNo code linkedOpen access
preprint2009arXiv

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

0 citations0 reviewsNo code linkedOpen access
preprint2009arXiv

The origin of very wide binary systems

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.

0 citations0 reviewsNo code linkedOpen access
preprint2009arXiv

What does the IMF really tell us about star formation?

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.

0 citations0 reviewsNo code linkedOpen access
preprint2006arXiv

A brown dwarf desert for intermediate mass stars in Sco OB2?

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.

0 citations0 reviewsNo code linkedOpen access