No Arabic abstract
We have coupled a fast, parametrized star cluster evolution code to a Markov Chain Monte Carlo code to determine the distribution of probable initial conditions of observed star clusters, which may serve as a starting point for future $N$-body calculations. In this paper we validate our method by applying it to a set of star clusters which have been studied in detail numerically with $N$-body simulations and Monte Carlo methods: the Galactic globular clusters M4, 47 Tucanae, NGC 6397, M22, $omega$ Centauri, Palomar 14 and Palomar 4, the Galactic open cluster M67, and the M31 globular cluster G1. For each cluster we derive a distribution of initial conditions that, after evolution up to the clusters current age, evolves to the currently observed conditions. We find that there is a connection between the morphology of the distribution of initial conditions and the dynamical age of a cluster and that a degeneracy in the initial half-mass radius towards small radii is present for clusters which have undergone a core collapse during their evolution. We find that the results of our method are in agreement with $N$-body and Monte Carlo studies for the majority of clusters. We conclude that our method is able to find reliable posteriors for the determined initial mass and half-mass radius for observed star clusters, and thus forms an suitable starting point for modeling an observed clusterrq{}s evolution.
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. Aims: The VVV observations will foster the construction of a sample of Galactic star clusters with reliable and homogeneously derived physical parameters (e.g., age, distance, and mass, etc.). In this first paper in a series, the methodology employed to establish cluster parameters for the envisioned database are elaborated upon by analyzing a subsample of 4 known young open clusters: Danks 1, Danks 2, RCW 79, and DBS 132. The analysis offers a first glimpse of the information that can be gleaned for the final cluster database from the VVV observations. Methods: Wide-field, deep JHKs VVV observations, combined with new infrared spectroscopy, are employed to constrain fundamental parameters for a subset of clusters. Results: Results inferred from the deep near-infrared photometry which features mitigated uncertainties (e.g. the accuracy of the photometry is better than 0.1mag for Ks<18mag), the wide field-of-view of the VVV survey, and numerous high quality low resolution spectra (typically more than 10 per cluster), are used to establish independent cluster parameters which enable existing determinations to be assessed. An anomalous reddening law in the direction toward the Danks clusters is found, i.e. E(J-H)/E(H-Ks)=2.20+/-0.06, which exceeds published values for the inner Galaxy. The G305 star forming complex, which includes the Danks clusters, lies beyond the Sagittarius-Carina spiral arm and occupies the Centaurus arm. Finally, the first deep infrared color-magnitude diagram of RCW 79 is presented which reveal a sizable pre-main sequence population. A list of candidate variable stars in G305 region is reported.
Cool giant and supergiant star atmospheres are characterized by complex velocity fields originating from convection and pulsation processes which are not fully understood yet. The velocity fields impact the formation of spectral lines, which thus contain information on the dynamics of stellar atmospheres. The tomographic method allows to recover the distribution of the component of the velocity field projected on the line of sight at different optical depths in the stellar atmosphere. The computation of the contribution function to the line depression aims at correctly identifying the depth of formation of spectral lines in order to construct numerical masks probing spectral lines forming at different optical depths. The tomographic method is applied to 1D model atmospheres and to a realistic 3D radiative hydrodynamics simulation performed with CO5BOLD in order to compare their spectral line formation depths and velocity fields. In 1D model atmospheres, each spectral line forms in a restricted range of optical depths. On the other hand, in 3D simulations, the line formation depths are spread in the atmosphere mainly because of temperature and density inhomogeneities. Comparison of CCF profiles obtained from 3D synthetic spectra with velocities from the 3D simulation shows that the tomographic method correctly recovers the distribution of the velocity component projected on the line of sight in the atmosphere.
We describe the development and implementation of the SEGUE (Sloan Extension for Galactic Exploration and Understanding) Stellar Parameter Pipeline (SSPP). The SSPP derives, using multiple techniques, radial velocities and the fundamental stellar atmospheric parameters (effective temperature, surface gravity, and metallicity) for AFGK-type stars, based on medium-resolution spectroscopy and $ugriz$ photometry obtained during the course of the original Sloan Digital Sky Survey (SDSS-I) and its Galactic extension (SDSS-II/SEGUE). The SSPP also provides spectral classification for a much wider range of stars, including stars with temperatures outside of the window where atmospheric parameters can be estimated with the current approaches. This is Paper I in a series of papers on the SSPP; it provides an overview of the SSPP, and initial tests of its performance using multiple data sets. Random and systematic errors are critically examined for the current version of the SSPP, which has been used for the sixth public data release of the SDSS (DR-6).
We present a new technique to quantify cluster-to-cluster variations in the observed present-day stellar mass functions of a large sample of star clusters. Our method quantifies these differences as a function of both the stellar mass and the total cluster mass, and offers the advantage that it is insensitive to the precise functional form of the mass function. We applied our technique to data taken from the ACS Survey for Globular Clusters, from which we obtained completeness-corrected stellar mass functions in the mass range 0.25-0.75 M$_{odot}$ for a sample of 27 clusters. The results of our observational analysis were then compared to Monte Carlo simulations for globular cluster evolution spanning a range of initial mass functions, total numbers of stars, concentrations, and virial radii. We show that the present-day mass functions of the clusters in our sample can be reproduced by assuming an universal initial mass function for all clusters, and that the cluster-to-cluster differences are consistent with what is expected from two-body relaxation. A more complete exploration of the initial cluster conditions will be needed in future studies to better constrain the precise functional form of the initial mass function. This study is a first step toward using our technique to constrain the dynamical histories of a large sample of old Galactic star clusters and, by extension, star formation in the early Universe.
We present combined interferometer and single dish telescope data of NH3 (J,K) = (1,1) and (2,2) emission towards the clustered star forming Ophiuchus B, C and F Cores at high spatial resolution (~1200 AU) using the Australia Telescope Compact Array, the Very Large Array, and the Green Bank Telescope. While the large scale features of the NH3 (1,1) integrated intensity appear similar to 850 micron continuum emission maps of the Cores, on 15 (1800 AU) scales we find significant discrepancies between the dense gas tracers in Oph B, but good correspondence in Oph C and F. Using the Clumpfind structure identifying algorithm, we identify 15 NH3 clumps in Oph B, and 3 each in Oph C and F. Only five of the Oph B NH3 clumps are coincident within 30 (3600 AU) of a submillimeter clump. We find v_LSR varies little across any of the Cores, and additionally varies by only ~1.5 km/s between them. The observed NH3 line widths within the Oph B and F Cores are generally large and often mildly supersonic, while Oph C is characterized by narrow line widths which decrease to nearly thermal values. We find several regions of localized narrow line emission (Delta v < 0.4 km/s), some of which are associated with NH3 clumps. We derive the kinetic temperatures of the gas, and find they are remarkably constant across Oph B and F, with a warmer mean value (T_K = 15 K) than typically found in isolated regions and consistent with previous results in clustered regions. Oph C, however, has a mean T_K = 12 K, decreasing to a minimum T_K = 9.4 K towards the submillimeter continuum peak, similar to previous studies of isolated starless cores. There is no significant difference in temperature towards protostars embedded in the Cores. [Abridged]