The kinematics and dynamics of young stellar populations enable us to test theories of star formation. With this aim, we continue our analysis of the SDSS-III/APOGEE IN-SYNC survey, a high resolution near infrared spectroscopic survey of young clusters. We focus on the Orion A star-forming region, for which IN-SYNC obtained spectra of $\sim2700$ stars. In Paper IV we used these data to study the young stellar population. Here we study the kinematic properties through radial velocities ($v_r$). The young stellar population remains kinematically associated with the molecular gas, following a $\sim10\:{\rm{km\:s}}^{-1}$ gradient along filament. However, near the center of the region, the $v_r$ distribution is slightly blueshifted and asymmetric; we suggest that this population, which is older, is slightly in foreground. We find evidence for kinematic subclustering, detecting statistically significant groupings of co-located stars with coherent motions. These are mostly in the lower-density regions of the cloud, while the ONC radial velocities are smoothly distributed, consistent with it being an older, more dynamically evolved cluster. The velocity dispersion $\sigma_v$ varies along the filament. The ONC appears virialized, or just slightly supervirial, consistent with an old dynamical age. Here there is also some evidence for on-going expansion, from a $v_r$--extinction correlation. In the southern filament, $\sigma_v$ is $\sim2$--$3$ times larger than virial in the L1641N region, where we infer a superposition along the line of sight of stellar sub-populations, detached from the gas. On the contrary, $\sigma_v$ decreases towards L1641S, where the population is again in agreement with a virial state.

We present radial velocity measurements for 70 high confidence, and 34 potential binary systems in fields containing the Perseus Molecular Cloud, Pleiades, NGC 2264, and the Orion A star forming region. 18 of these systems have been previously identified as binaries in the literature. Candidate double-lined spectroscopic binaries (SB2s) are identified by analyzing the cross-correlation functions (CCFs) computed during the reduction of each APOGEE spectrum. We identify sources whose CCFs are well fit as the sum of two Lorentzians as likely binaries, and provide an initial characterization of the system based on the radial velocities indicated by that dual fit. For systems observed over several epochs, we present mass ratios and systemic velocities; for two systems with observations on eight or more epochs, and which meet our criteria for robust orbital coverage, we derive initial orbital parameters. The distribution of mass ratios for multi-epoch sources in our sample peaks at q=1, but with a significant tail toward lower q values. Tables reporting radial velocities, systemic velocities, and mass ratios are provided online. We discuss future improvements to the radial velocity extraction method we employ, as well as limitations imposed by the number of epochs currently available in the APOGEE database. The Appendix contains brief notes from the literature on each system in the sample, and more extensive notes for select sources of interest.

The Gaia ESO Public Spectroscopic Survey (GES) is providing the astronomical community with high-precision measurements of many stellar parameters including radial velocities (RVs) of stars belonging to several young clusters and star-forming regions. One of the main goals of the young cluster observations is to study of their dynamical evolution and provide insight into their future, revealing if they will eventually disperse to populate the field, rather than evolve into bound open clusters. In this paper we report the analysis of the dynamical state of L1688 in the $\rho$~Ophiuchi molecular cloud using the dataset provided by the GES consortium. We performed the membership selection of the more than 300 objects observed. Using the presence of the lithium absorption and the location in the Hertzspung-Russell diagram, we identify 45 already known members and two new association members. We provide accurate RVs for all 47 confirmed members.A dynamical analysis, after accounting for unresolved binaries and errors, shows that the stellar surface population of L1688 has a velocity dispersion $\sigma \sim$1.14$\pm$0.35 km s$^{-1}$ that is consistent with being in virial equilibrium and is bound with a $\sim$80% probability. We also find a velocity gradient in the stellar surface population of $\sim$1.0 km s$^{-1}$pc$^{-1}$ in the northwest/southeast direction, which is consistent with that found for the pre-stellar dense cores, and we discuss the possibility of sequential and triggered star formation in L1688.

We present the results of the SDSS APOGEE INfrared Spectroscopy of Young Nebulous Clusters program (IN-SYNC) survey of the Orion A molecular cloud. This survey obtained high resolution near infrared (NIR) spectroscopy of about 2700 young pre-main sequence stars throughout the region, acquired across five distinct fields spanning 6deg field of view (FOV). With these spectra, we have measured accurate stellar parameters (T_eff, log g, v sin i) and extinctions, and placed the sources in the Hertzsprung-Russel Diagram (HRD). We have also extracted radial velocities for the kinematic characterization of the population. We compare our measurements with literature results for a sub-sample of targets in order to assess the performances and accuracy of the survey. Source extinction shows evidence for dust grains that are larger than those in the diffuse interstellar medium (ISM): we estimate an average R_V=5.5 in the region. Importantly, we find a clear correlation between HRD inferred ages and spectroscopic surface-gravity inferred ages. This clearly indicates a real spread of stellar radii at fixed temperature, and together with additional correlations with extinction and with disk presence, strongly suggests a real spread of ages large than a few Myr. Focussing on the young population around NGC1980 iota Ori, which has previously been suggested to be a separate, foreground, older cluster, we confirm its older (5Myr) age and low A_V, but considering that its radial velocity distribution is indistinguishable from the Orion A's population, we suggest that NGC1980 is part of Orion A's star formation activity. Based on their stellar parameters and kinematic properties, we identify 383 new candidate members of Orion A, most of which are diskless sources in areas of the region poorly studied by previous works.

Most field stars will have encountered the highest stellar density and hence the largest number of interactions in their birth environment. Yet the stellar dynamics during this crucial phase are poorly understood. Here we analyze the radial velocities measured for 152 out of 380 observed stars in the 2-6 Myr old star cluster IC 348 as part of the SDSS-III APOGEE. The radial velocity distribution of these stars is fitted with one or two Gaussians, convolved with the measurement uncertainties including binary orbital motions. Including a second Gaussian improves the fit; the high-velocity outliers that are best fit by this second component may either (1) be contaminants from the nearby Perseus OB2 association, (2) be a halo of ejected or dispersing stars from IC 348, or (3) reflect that IC 348 has not relaxed to a Gaussian velocity distribution. We measure a velocity dispersion for IC 348 of $0.72 \pm 0.07$ km s$^{-1}$ (or $0.64 \pm 0.08$ km s$^{-1}$ if two Gaussians are fitted), which implies a supervirial state, unless the gas contributes more to the gravitational potential than expected. No evidence is found for a dependence of this velocity dispersion on distance from the cluster center or stellar mass. We also find that stars with lower extinction (in the front of the cloud) tend to be redshifted compared with stars with somewhat higher extinction (towards the back of the cloud). This data suggests that the stars in IC 348 are converging along the line of sight. We show that this correlation between radial velocity and extinction is unlikely to be spuriously caused by the small cluster rotation of $0.024 \pm 0.013$ km s$^{-1}$ arcmin$^{-1}$ or by correlations between the radial velocities of neighboring stars. This signature, if confirmed, will be the first detection of line-of-sight convergence in a star cluster(...)

This paper describes the analysis of UVES and GIRAFFE spectra acquired by the Gaia-ESO Public Spectroscopic Survey in the fields of young clusters whose population includes pre-main sequence (PMS) stars. Both methods that have been extensively used in the past and new ones developed in the contest of the Gaia-ESO survey enterprise are available and used. The internal precision of these quantities is estimated by inter-comparing the results obtained by such different methods, while the accuracy is estimated by comparison with independent external data, like effective temperature and surface gravity derived from angular diameter measurements, on a sample of benchmarks stars. Specific strategies are implemented to deal with fast rotation, accretion signatures, chromospheric activity, and veiling. The analysis carried out on spectra acquired in young clusters' fields during the first 18 months of observations, up to June 2013, is presented in preparation of the first release of advanced data products. Stellar parameters obtained with the higher resolution and larger wavelength coverage from UVES are reproduced with comparable accuracy and precision using the smaller wavelength range and lower resolution of the GIRAFFE setup adopted for young stars, which allows us to provide with confidence stellar parameters for the much larger GIRAFFE sample. Precisions are estimated to be $\approx$ 120 K r.m.s. in Teff, $\approx$0.3 dex r.m.s. in logg, and $\approx$0.15 dex r.m.s. in [Fe/H], for both the UVES and GIRAFFE setups.

The nearby (distance~350-400 pc), rich Vela OB2 association, includes $\gamma^2$ Velorum, one of the most massive binaries in the solar neighbourhood and an excellent laboratory for investigating the formation and early evolution of young clusters. Recent Gaia-ESO survey observations have led to the discovery of two kinematically distinct populations in the young (10-15 Myr) cluster immediately surrounding $\gamma^2$ Velorum. Here we analyse the results of Gaia-ESO survey observations of NGC 2547, a 35 Myr cluster located two degrees south of $\gamma^2$ Velorum. The radial velocity distribution of lithium-rich pre-main sequence stars shows a secondary population that is kinematically distinct from and younger than NGC 2547. The radial velocities, lithium absorption lines, and the positions in a colour-magnitude diagram of this secondary population are consistent with those of one of the components discovered around $\gamma^2$ Velorum. This result shows that there is a young, low-mass stellar population spread over at least several square degrees in the Vela OB2 association. This population could have originally been part of a cluster around $\gamma^2$ Velorum that expanded after gas expulsion or formed in a less dense environment that is spread over the whole Vela OB2 region.

The initial velocity dispersion of newborn stars is a major unconstrained aspect of star formation theory. Using near-infrared spectra obtained with the APOGEE spectrograph, we show that the velocity dispersion of young (1-2 Myr) stars in NGC 1333 is 0.92+/-0.12 km/s after correcting for measurement uncertainties and the effect of binaries. This velocity dispersion is consistent with the virial velocity of the region and the diffuse gas velocity dispersion, but significantly larger than the velocity dispersion of the dense, star-forming cores, which have a sub-virial velocity dispersion of 0.5 km/s. Since the NGC 1333 cluster is dynamically young and deeply embedded, this measurement provides a strong constraint on the initial velocity dispersion of newly-formed stars. We propose that the difference in velocity dispersion between stars and dense cores may be due to the influence of a 70 micro-Gauss magnetic field acting on the dense cores, or be the signature of a cluster with initial sub-structure undergoing global collapse.

Over two years 8,859 high-resolution H-band spectra of 3,493 young (1 - 10 Myr) stars were gathered by the multi-object spectrograph of the APOGEE project as part of the IN-SYNC ancillary program of that SDSS-III survey. Here we present the forward modeling approach used to derive effective temperatures, surface gravities, radial velocities, rotational velocities, and H-band veiling from these near-infrared spectra. We discuss in detail the statistical and systematic uncertainties in these stellar parameters. In addition we present accurate extinctions by measuring the E(J-H) of these young stars with respect to the single-star photometric locus in the Pleiades. Finally we identify an intrinsic stellar radius spread of about 25% for late-type stars in IC 348 using three (nearly) independent measures of stellar radius, namely the extinction-corrected J-band magnitude, the surface gravity and the $R \sin i$ from the rotational velocities and literature rotation periods. We exclude that this spread is caused by uncertainties in the stellar parameters by showing that the three estimators of stellar radius are correlated, so that brighter stars tend to have lower surface gravities and larger $R \sin i$ than fainter stars at the same effective temperature. Tables providing the spectral and photometric parameters for the Pleiades and IC 348 have been provided online.

Context: A key science goal of the Gaia-ESO survey (GES) is to use the kinematics of low-mass stars in young clusters to probe their dynamical histories and how they populate the field as they become unbound. The clustering of low-mass stars around the massive W-R binary gamma^2 Velorum was one of the first GES targets. Aims: To empirically determine the radial velocity (RV) precision of GES data, construct a kinematically unbiased sample of cluster members and characterise their dynamical state. Methods: Targets were selected from colour-magnitude diagrams and intermediate resolution spectroscopy used to derive RVs and assess membership from the strength of the Li6708A line. The RV distribution was analysed using a maximum likelihood technique that accounts for unresolved binaries. Results: The GES RV precision is about 0.25km/s and sufficient to resolve velocity structure in the low-mass population around gamma^2 Vel. The structure is well fitted by two kinematic components with roughly equal numbers of stars; the first has an intrinsic dispersion of 0.34+/-0.16km/s, consistent with virial equilibrium. The second has a broader dispersion of 1.60+/-0.37km/s and is offset from the first by ~2km/s. The first population is older by 1-2Myr based on a greater level of Li depletion seen among its M-stars and is probably more centrally concentrated around gamma^2 Vel. Conclusions: We consider several formation scenarios, concluding that the two kinematic components are a bound remnant of the original, denser cluster that formed gamma^2 Vel, and a dispersed population from the wider Vela OB2 association, of which gamma^2 Vel is the most massive member. The apparent youth of gamma^2 Vel compared to the older (>=10Myr) low-mass population surrounding it suggests a scenario where the massive binary formed in a clustered environment after the formation of the bulk of the low-mass stars.[abridged]

We have used Sloan Digital Sky Survey-III (SDSS-III) Apache Point Observatory Galactic Evolution Experiment (APOGEE) radial velocity observations in the near-infrared $H$-band to explore the membership of the nearby ($86.7 \pm 0.9$ pc) open cluster Coma Berenices (Melotte 111), concentrating on the poorly-populated low-mass end of the main sequence. Using SDSS-III APOGEE radial velocity measurements, we confirm the membership of eight K/M dwarf members, providing the first confirmed low-mass members of the Coma Berenices cluster. Using $R\sim2000$ spectra from IRTF-SpeX, we confirm the independently luminosity classes of these targets, and find their metallicities to be consistent with the known solar mean metallicity of Coma Berenices and of M dwarfs in the Solar neighborhood. In addition, the APOGEE spectra have enabled measurement of $v\sin i$ for each target and detection for the first time of the low-mass secondary components of the known binary systems Melotte 111 102 and Melotte 111 120, as well as identification of the previously unknown binary system 2MASS J12214070+2707510. Finally, we use Kilodegree Extremely Little Telescope (KELT) photometry to measure photometric variability and rotation periods for a subset of the Coma Berenices members.

Orbital motions from binary stars can broaden the observed line-of-sight velocity distribution of a stellar system, artificially inflating the measured line-of-sight velocity dispersion, which can in turn lead to erroneous conclusions about the dynamical state of the system. Cottaar et al. (2012b) proposed a maximum likelihood procedure to recover the intrinsic velocity dispersion of a resolved star cluster from a single epoch of radial velocity data of individual stars, which they achieved by simultaneously fitting the intrinsic velocity distribution of the single stars and the centres of mass of the binaries along with the velocity shifts caused by binary orbital motions. Assuming well-characterized binary properties, they showed that this procedure can accurately reproduce intrinsic velocity dispersions below 1 km s$^{-1}$ for solar-type stars. Here we investigate the systematic offsets induced in cases where the binary properties are uncertain, and we show how two epochs of radial velocity data with an appropriate baseline can help to mitigate these systematic effects. We first test the method above using Monte Carlo simulations, taking into account the large uncertainties in the binary properties of OB stars. We then apply it to radial velocity data in the young massive cluster R136, an example for which the intrinsic velocity dispersion of O-type stars is known from an intensive multi-epoch approach. For typical velocity dispersions of young massive clusters ($\gtrsim 4$ km s$^{-1}$) and with a single epoch of data, we demonstrate that the method can just about distinguish between a cluster in virial equilibrium and an unbound cluster. This is due to the higher spectroscopic binary fraction and more loosely constrained distributions of orbital parameters of OB stars compared to solar-type stars. By extending the maximum likelihood method to multi-epoch data, <abridged>

Radial velocity measurements can be used to constrain the dynamical state of a stellar cluster. However, for clusters with velocity dispersions smaller than a few km/s the observed radial velocity distribution tends to be dominated by the orbital motions of binaries rather than the stellar motions through the potential well of the cluster. Our goal is to characterize the intrinsic velocity distribution of a cluster from a single epoch of radial velocity data, even for a cluster with a velocity dispersion of a fraction of a km/s, using a maximum likelihood procedure. Assuming a period, mass ratio, and eccentricity distribution for the binaries in the observed cluster this procedure fits a dynamical model describing the velocity distribution for the single stars and center of masses of the binaries, simultaneously with the radial velocities caused by binary orbital motions, using all the information available in the observed velocity distribution. We find that the fits to the intrinsic velocity distribution depend only weakly on the binary properties assumed, so the uncertainty in the fitted parameters tends to be dominated by statistical uncertainties. Based on Monte Carlo simulations we provide an estimate of how these statistical uncertainties vary with the velocity dispersion, binary fraction, and the number of observed stars, which can be used to estimate the sample size needed to reach a specific accuracy. Finally we test the method on the well-studied open cluster NGC 188, showing that it can reproduce a velocity dispersion of only 0.5 km/s using a single epoch of the multi-epoch radial velocity data. If the binary period, mass ratio, and eccentricity distribution of the observed stars are roughly known, this procedure can be used to correct for the effect of binary orbital motions on an observed velocity distribution. [Abridged]

Context. Westerlund I is the richest young cluster currently known in our Galaxy, making it one of the most massive clusters for which we can resolve the individual stars even in the crowded centre. This makes it an ideal target to assess whether massive clusters formed currently will remain bound or will disperse and contribute significantly to the stellar field population. Aims. Here we report a measurement of the radial velocity dispersion of Westerlund I to explore whether the cluster is currently in virial equilibrium, if it is in the process of collapse or if it is expanding and dispersing into the field. Methods. We obtained MIKE/Magellan high resolution optical spectra of 22 post main-sequence stars jn Westerlund I for 2 or 3 epochs with a maximum baseline of about one year. Radial velocities variations between these spectra have been measured through cross correlation. Results. We calculate the velocity dispersion from the cross correlation of five yellow hypergiants and one luminous blue variable, that show little radial velocity variations between epochs and have many spectral features in common. After taking into account the effect of small number statistics and undetected binaries, we estimate the velocity dispersion for the massive stars in Westerlund I to be 2.1 (+3.3, -2.1) km s-1. For several different assumptions concerning possible mass segregation and the elongation of the cluster, we find that Westerlund I is subvirial at the 90% confidence level. Conclusions. We can rule out that the cluster is significantly supervirial at the 97% confidence level, indicating that Westerlund I is currently bound. This implies that Westerlund I has survived past the point where any gas expulsion has taken place and is expected to survive for billions of years.

Rapid mass transfer in a binary system can drive the accreting star out of thermal equilibrium, causing it to expand. This can lead to a contact system, strong mass loss from the system and possibly merging of the two stars. In low metallicity stars the timescale for heat transport is shorter due to the lower opacity. The accreting star can therefore restore thermal equilibrium more quickly and possibly avoid contact. We investigate the effect of accretion onto main sequence stars with radiative envelopes with different metallicities. We find that a low metallicity (Z<0.001), 4 solar mass star can endure a 10 to 30 times higher accretion rate before it reaches a certain radius than a star at solar metallicity. This could imply that up to two times fewer systems come into contact during rapid mass transfer when we compare low metallicity. This factor is uncertain due to the unknown distribution of binary parameters and the dependence of the mass transfer timescale on metallicity. In a forthcoming paper we will present analytic fits to models of accreting stars at various metallicities intended for the use in population synthesis models.