Hierarchical star formation: stars and stellar clusters in the Gould Belt

We perform a study of the spatial and kinematical distribution of young open clusters in the solar neighbourhood, discerning between bound clusters and transient stellar condensations within our sample. Then, we discriminate between Gould Belt (GB) and local Galactic disc (LGD) members, using our previous estimate of the structural parameters of both systems obtained from a sample of O-B6 Hipparcos stars. Single membership probabilities of the clusters are also calculated in the separation process. Using this classified sample, we analyse the spatial structure and the kinematic behaviour of the cluster system in the GB. The two star formation regions that dominate and give the GB its characteristic-inclined shape show a striking difference in their content of star clusters: while Ori OB1 is richly populated by open clusters, not a single one can be found within the boundaries of Sco OB2. This is mirrored in the velocity space, translating again into an abundance of clusters in the region of the kinematic space populated by the members of Ori OB1, and a marginal number of them associated with Sco OB2. We interpret all these differences by characterizing the Orion region as a cluster complex typically surrounded by a stellar halo, and the Sco-Cen region as an OB association in the outskirts of the complex. In the light of these results, we study the nature of the GB with respect to the optical segment of the Orion Arm, and we propose that the different content of star clusters, the different heights over the Galactic plane and the different residual velocities of Ori OB1 and Sco OB2 can be explained in terms of their relative position to the density maximum of the Local Arm in the solar neighbourhood. Although morphologically intriguing, the GB appears to be the result of our local and biased view of a larger star cluster complex in the Local Arm, that could be explained by the internal dynamics of the Galactic disc.     

Ratios of star cluster core and half-mass radii: a cautionary note on intermediate-mass black holes in star clusters

There is currently much interest in the possible presence of intermediate-mass black holes (IMBHs) in the cores of globular clusters (GCs). Based on theoretical arguments and simulation results it has previously been suggested that a large core radius – or particularly a large ratio of the core radius to half-mass radius – is a promising indicator for finding such a black hole (BH) in a star cluster. In this study N -body models of 100 000 stars with and without primordial binaries are used to investigate the long-term structural evolution of star clusters. Importantly, the simulation data are analysed using the same processes by which structural parameters are extracted from observed star clusters. This gives a ratio of the core and half-mass (or half-light) radii that are directly comparable to the Galactic GC sample. As a result, it is shown that the ratios observed for the bulk of this sample can be explained without the need for an IMBH. Furthermore, it is possible that clusters with large core to half-light radius ratios harbour a BH binary (comprising stellar mass BHs) rather than a single massive BH. This work does not rule out the existence of IMBHs in the cores of at least some star clusters.     

Quantitative analysis of clumps in the tidal tails of star clusters

Tidal tails of star clusters are not homogeneous but show well-defined clumps in observations as well as in numerical simulations. Recently, an epicyclic theory for the formation of these clumps was presented. A quantitative analysis was still missing. We present a quantitative derivation of the angular momentum and energy distribution of escaping stars from a star cluster in the tidal field of the Milky Way and derive the connection to the position and width of the clumps. For the numerical realization we use star-by-star N -body simulations. We find a very good agreement of theory and models. We show that the radial offset of the tidal arms scales with the tidal radius, which is a function of cluster mass and the rotation curve at the cluster orbit. The mean radial offset is 2.77 times the tidal radius in the outer disc. Near the Galactic Centre the circumstances are more complicated, but to lowest order the theory still applies. We have also measured the Jacobi energy distribution of bound stars and showed that there is a large fraction of stars (about 35 per cent) above the critical Jacobi energy at all times, which can potentially leave the cluster. This is a hint that the mass loss is dominated by a self-regulating process of increasing Jacobi energy due to the weakening of the potential well of the star cluster, which is induced by the mass loss itself.     

Captured older stars as the reason for apparently prolonged star formation in young star clusters

The existence of older stars within a young star cluster can be interpreted to imply that star formation occurs on time-scales longer than a free-fall time of a pre-cluster cloud core. Here, the idea is explored that these older stars are not related to the star formation process forming the young star cluster but rather that the orbits of older field stars are focused by the collapsing pre-cluster cloud core. Two effects appear: the focusing of stellar orbits leads to an enhancement of the density of field stars in the vicinity of the centre of the young star cluster; and due to the time-dependent potential of the forming cluster some of these stars can get bound gravitationally to the cluster. These stars exhibit similar kinematical properties to the newly formed stars and cannot be distinguished from them on the basis of radial velocity or proper motion surveys. Such contaminations may lead to a wrong apparent star formation history of a young cluster. In the case of the ONC, the theoretical number of gravitationally bound older low-mass field stars agrees with the number of observed older low-mass stars.     

The effect of gas loss on the formation of bound stellar clusters

The effect of gas ejection on the structure and binding energy of newly formed stellar clusters is investigated. The star formation efficiency (SFE), necessary for forming a gravitationally bound stellar cluster, is determined.
Two sets of numerical N -body simulations are presented. As a first simplified approach we treat the residual gas as an external potential. The gas expulsion is approximated by reducing the gas mass to zero on a given time-scale, which is treated as a free parameter. In a second set of simulations we use smoothed particle hydrodynamics (SPH) to follow the dynamics of the outflowing residual gas self-consistently. We investigate cases where gas outflow is induced by an outwards propagating shock front and where the whole gas cloud is heated homogeneously, leading to ejection.
If the stars are in virial equilibrium with the gaseous environment initially, bound clusters only form in regions where the local SFE is larger than 50 per cent or where the gas expulsion time-scale is long compared with the dynamical time-scale. A small initial velocity dispersion of the stars leads to a compaction of the cluster during the expulsion phase and reduces the SFE needed to form bound clusters to less than 10 per cent.     

Collaborative research of open star clusters

Preliminary results on observations of open clusters are presented. The project has been initiated in the framework of the Uzbek-Taiwan and Taiwan-Baltic collaboration, mainly to upgrade and make use of facilities at Maidanak Observatory. We present detailed, multiwavelength studies of the young cluster NGC 6823 and the associated complex nebulosity, to diagnose the young stellar population and star formation history in the region. In addition, 7 compact open clusters have been monitored for stellar variability. We show how observations like these could feasibly be used to look for exoplanet transit events. We also expect to join the Whole-Earth Telescope effort in future campaigns for asteroseismology.     

The time-scale of escape from star clusters

In this paper a cluster is modelled as a smooth potential (due to the cluster stars) plus the steady tidal field of the Galaxy. In this model there is a minimum energy below which stars cannot escape. Above this energy, however, the time-scale on which a star escapes varies with the orbital parameters of the star (mainly its energy) in a way which we attempt to quantify, with both theoretical arguments and computer simulations. Within the limitations of the model we show that the time-scale is long enough to complicate the interpretation of full N -body simulations of clusters, and that stars above the escape energy may remain bound to the cluster for about a Hubble time.     

Analysing the structure of elongated star clusters

We have previously reported a measure     which both quantifies and distinguishes between a (relatively smooth) large-scale radial density gradient and multiscale (fractal) subclustering. Here, we extend the applicability of     to clusters which deviate significantly from an overall circular shape.
    varies systematically as clusters assume a more elongated shape, and it is therefore possible to correct for the effect, if the elongation of the cluster is also quantified.     therefore remains a useful and robust analytical technique for classifying and quantifying the internal structure of star clusters, even when their overall shape is far from circular.
The corrections required are small for individual clusters which are not extremely elongated (not more than three times longer than they are wide) of the same order as the uncertainty in the value of     for a particular cluster type. We therefore recommend that no correction be applied to the calculation of     for individual clusters, unless they are more than three times longer than their width, but that correction for elongation be applied when     is used for statistical analyses of large numbers of observed or simulated clusters.     

A Kinematical Calibration of the Galactocentric Distance

1 INTRODUCTION The distance to the Galactic center R0 is a fundamental constant for astronomy and astrophysics. Most determinations of astronomical quantities are directly connected with the Galactic distance scale, e.g., the rotational speed of our Galax…     

Measurement of the epicycle frequency in the Galactic disc and initial velocities of open clusters

A new method to measure the epicycle frequency κ in the Galactic disc is presented. We make use of the large data base on open clusters completed by our group to derive the observed velocity vector (amplitude and direction) of the clusters in the Galactic plane. In the epicycle approximation, this velocity is equal to the circular velocity given by the rotation curve, plus a residual or perturbation velocity, of which the direction rotates as a function of time with the frequency κ. Due to the non-random direction of the perturbation velocity at the birth time of the clusters, a plot of the present-day direction angle of this velocity as a function of the age of the clusters reveals systematic trends from which the epicycle frequency can be obtained. Our analysis considers that the Galactic potential is mainly axis-symmetric, or in other words, that the effect of the spiral arms on the Galactic orbits is small; in this sense, our results do not depend on any specific model of the spiral structure. The values of κ that we obtain provide constraints on the rotation velocity of the disc; in particular, V ₀ is found to be  230 ± 15 km s⁻¹  even if the short scale  ( R ₀= 7.5 kpc)  of the Galaxy is adopted. The measured κ at the solar radius is  43 ± 5 km s⁻¹ kpc ⁻¹  . The distribution of initial velocities of open clusters is discussed.     

Roche volume filling and the dissolution of open star clusters

From direct N‐body simulations we find that the dynamical evolution of star clusters is strongly influenced by the Roche volume filling factor. We present a parameter study of the dissolution of open star clusters with different Roche volume filling factors and different particle numbers. We study both Roche volume underfilling and overfilling models and compare with the Roche volume filling case. We find that in the Roche volume overfilling limit of our simulations two‐body relaxation is no longer the dominant dissolution mechanism but the changing cluster potential. We call this mechanism “mass‐loss driven dissolution” in contrast to “two‐body relaxation driven dissolution” which occurs in the Roche volume underfilling regime. We have measured scaling exponents of the dissolution time with the two‐body relaxation time. In this experimental study we find a decreasing scaling exponent with increasing Roche volume filling factor. The evolution of the escaper number in the Roche volume overfilling limit can be described by a log‐logistic differential equation. We report the finding of a resonance condition which may play a role for the evolution of star clusters and may be calibrated by the main periodic orbit in the large island of retrograde quasiperiodic orbits in the Poincaré surfaces of section. We also report on the existence of a stability curve which may be of relevance with respect to the structure of star clusters. (© 2015 WILEY‐VCH Verlag GmbH & Co. KGaA, Weinheim)