No Arabic abstract
Until now it has been impossible to observationally measure how star cluster scale height evolves beyond 1Gyr as only small samples have been available. Here we establish a novel method to determine the scale height of a cluster sample using modelled distributions and Kolmogorov-Smirnov tests. This allows us to determine the scale height with a 25% accuracy for samples of 38 clusters or more. We apply our method to investigate the temporal evolution of cluster scale height, using homogeneously selected sub-samples of Kharchenko et al. (MWSC), Dias et al. (DAML02), WEBDA, and Froebrich et al. (FSR). We identify a linear relationship between scale height and log(age/yr) of clusters, considerably different from field stars. The scale height increases from about 40pc at 1Myr to 75pc at 1Gyr, most likely due to internal evolution and external scattering events. After 1Gyr, there is a marked change of the behaviour, with the scale height linearly increasing with log(age/yr) to about 550pc at 3.5Gyr. The most likely interpretation is that the surviving clusters are only observable because they have been scattered away from the mid-plane in their past. A detailed understanding of this observational evidence can only be achieved with numerical simulations of the evolution of cluster samples in the Galactic Disk. Furthermore, we find a weak trend of an age-independent increase in scale height with galactocentric distance. There are no significant temporal or spatial variations of the cluster distribution zero point. We determine the Suns vertical displacement from the Galactic Plane as $Z_odot=18.5pm1.2$pc.
Context: The ESO Public Survey VISTA Variables in the Via Lactea (VVV) provides deep multi-epoch infrared observations for an unprecedented 562 sq. degrees of the Galactic bulge, and adjacent regions of the disk. In this survey nearly 150 new open clusters and cluster candidates have been discovered. Aims: This is the second in a series of papers about young, massive open clusters observed using the VVV survey. We present the first study of six recently discovered clusters. These clusters contain at least one newly discovered Wolf-Rayet (WR) star. Methods: Following the methodology presented in the first paper of the series, wide-field, deep JHKs VVV observations, combined with new infrared spectroscopy, are employed to constrain fundamental parameters for a subset of clusters. Results: We affirm that the six studied stellar groups are real young (2-7 Myr) and massive (between 0.8 and 2.2 10^3 Msol) clusters. They are highly obscured (Av ~ 5-24 mag) and compact (1-2 pc). In addition to WR stars, two of the six clusters also contain at least one red supergiant star. We claim the discovery of 8 new WR stars, and 3 stars showing WR-like emission lines which could be classified WR or OIf. Preliminary analysis provides initial masses of ~30-50 Msol for the WR stars. Finally,we discuss the spiral structure of the Galaxy using as tracers the six new clusters together with the previously studied VVV clusters.
Based on recent findings of a formation mechanism of substructure in tidal tails by Kuepper, Macleod & Heggie (2008) we investigate a more comprehensive set of N-body models of star clusters on orbits about a Milky-Way-like potential. We find that the predicted epicyclic overdensities arise in any tidal tail no matter which orbit the cluster follows as long as the cluster lives long enough for the overdensities to build up. The distance of the overdensities along the tidal tail from the cluster centre depends for circular orbits only on the mass of the cluster and the strength of the tidal field, and therefore decreases monotonically with time, while for eccentric orbits the orbital motion influences the distance, causing a periodic compression and stretching of the tails and making the distance oscillate with time. We provide an approximation for estimating the distance of the overdensities in this case. We describe an additional type of overdensity which arises in extended tidal tails of clusters on eccentric orbits, when the acceleration of the tidal field on the stellar stream is no longer homogeneous. Moreover, we conclude that a pericentre passage or a disk shock is not the direct origin of an overdensity within a tidal tail. Escape due to such tidal perturbations does not take place immediately after the perturbation but is rather delayed and spread over the orbit of the cluster. All observable overdensities are therefore of the mentioned two types. In particular, we note that substructured tidal tails do not imply the existence of dark-matter sub-structures in the haloes of galaxies.
We study the evolution of embedded clusters. The equations of motion of the stars in the cluster are solved by direct N-body integration while taking the effects of stellar evolution and the hydrodynamics of the natal gas content into account. The gravity of the stars and the surrounding gas are coupled self consistently to allow the realistic dynamical evolution of the cluster. While the equations of motion are solved, a stellar evolution code keeps track of the changes in stellar mass, luminosity and radius. The gas liberated by the stellar winds and supernovae deposits mass and energy into the gas reservoir in which the cluster is embedded. We examine cluster models with 1000 stars, but we varied the star formation efficiency (between 0.05-0.5), cluster radius (0.1-1.0 parsec), the degree of virial support of the initial population of stars (0-100%) and the strength of the feedback. We find that an initial star fraction $M_star/M_{rm tot} > 0.05$ is necessary for cluster survival. Survival is more likely if gas is not blown out violently by a supernova and if the cluster has time to approach virial equilibrium during out-gassing. While the cluster is embedded, dynamical friction drives early and efficient mass segregation in the cluster. Stars of $m gtrsim 2,M_odot$ are preferentially retained, at the cost of the loss of less massive stars. We conclude that the degree of mass segregation in open clusters such as the Pleiades is not the result of secular evolution but a remnant of its embedded stage.
We present a novel and flexible tensor approach to computing the effect of a time-dependent tidal field acting on a stellar system. The tidal forces are recovered from the tensor by polynomial interpolation in time. The method has been implemented in a direct-summation stellar dynamics integrator (NBODY6) and test-proved through a set of reference calculations: heating, dissolution time and structural evolution of model star clusters are all recovered accurately. The tensor method is applicable to arbitrary configurations, including the important situation where the background potential is a strong function of time. This opens up new perspectives in stellar population studies reaching to the formation epoch of the host galaxy or galaxy cluster, as well as for star-burst events taking place during the merger of large galaxies. A pilot application to a star cluster in the merging galaxies NGC 4038/39 (the Antennae) is presented.
Dense stellar clusters are natural sites for the origin and evolution of exotic objects such as relativistic binaries (potential gravitational wave sources), blue stragglers, etc. We investigate the secular dynamics of a binary system driven by the global tidal field of an axisymmetric stellar cluster in which the binary orbits. In a companion paper (Hamilton & Rafikov 2019a) we developed a general Hamiltonian framework describing such systems. The effective (doubly-averaged) Hamiltonian derived there encapsulates all information about the tidal potential experienced by the binary in its orbit around the cluster in a single parameter $Gamma$. Here we provide a thorough exploration of the phase-space of the corresponding secular problem as $Gamma$ is varied. We find that for $Gamma > 1/5$ the phase-space structure and the evolution of binary orbital element are qualitatively similar to the Lidov-Kozai problem. However, this is only one of four possible regimes, because the dynamics are qualitatively changed by bifurcations at $Gamma = 1/5,0,-1/5$. We show how the dynamics are altered in each regime and calculate characteristics such as secular evolution timescale, maximum possible eccentricity, etc. We verify the predictions of our doubly-averaged formalism numerically and find it to be very accurate when its underlying assumptions are fulfilled, typically meaning that the secular timescale should exceed the period of the binary around the cluster by $gtrsim 10-10^2$ (depending on the cluster potential and binary orbit). Our results may be relevant for understanding the nature of a variety of exotic systems harboured by stellar clusters.