Context. A key science goal of the Gaia-ESO survey (GES) at the VLT is to use the kinematics of low-mass stars in young clusters and star forming regions to probe their dynamical histories and how they populate the field as they become unbound. The clustering of low-mass stars around the massive Wolf-Rayet binary system γ 2 Velorum was one of the first GES targets. Aims. We empirically determine the radial velocity 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 was used to derive radial velocities and assess membership from the strength of the Li i 6708 Å line. The radial velocity distribution was analysed using a maximum likelihood technique that accounts for unresolved binaries. Results. The GES radial velocity precision is about 0.25 km s −1 and sufficient to resolve velocity structure in the low-mass population around γ 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.16 km s −1 , consistent with virial equilibrium. The second has a broader dispersion of 1.60 ± 0.37 km s −1 and is offset from the first by 2 km s −1 . The first population is older by 1-2 Myr based on a greater level of Li depletion seen among its M-type stars and is probably more centrally concentrated around γ 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 γ 2 Vel, and a dispersed population from the wider Vela OB2 association, of which γ 2 Vel is the most massive member. The apparent youth of γ 2 Vel compared to the older (≥10 Myr) low-mass population surrounding it suggests a scenario in which the massive binary formed in a clustered environment after the formation of the bulk of the low-mass stars.

We present new Spitzer Space Telescope observations of stars in the young ($5 Myr) Velorum stellar cluster. Combining optical and 2MASS photometry, we have selected 579 stars as candidate members of the cluster. With the addition of the Spitzer mid-infrared data, we have identified five debris disks around A-type stars and five to six debris disks around solar-type stars, indicating that the strong radiation field in the cluster does not completely suppress the production of planetesimals in the disks of cluster members. However, we find some evidence that the frequency of circumstellar primordial disks is lower, and the infrared flux excesses are smaller than for disks around stellar populations with similar ages. This could be evidence for a relatively fast dissipation of circumstellar dust by the strong radiation field from the highest mass star(s) in the cluster. Another possibility is that Velorum stellar cluster is slightly older than reported ages and the low frequency of primordial disks reflects the fast disk dissipation observed at $5 Myr.

We have derived ages for 13 young (< 30 Myr) star-forming regions and find they are up to a factor two older than the ages typically adopted in the literature. This result has wide-ranging implications, including that circumstellar discs survive longer ( 10 − 12 Myr) and that the average Class I lifetime is greater ( 1 Myr) than currently believed.For each star-forming region we derived two ages from colour-magnitude diagrams. First we fitted models of the evolution between the zero-age main-sequence and terminal-age main-sequence to derive a homogeneous set of main-sequence ages, distances and reddenings with statistically meaningful uncertainties. Our second age for each star-forming region was derived by fitting pre-main-sequence stars to new semi-empirical model isochrones. For the first time (for a set of clusters younger than 50 Myr) we find broad agreement between these two ages, and since these are derived from two distinct mass regimes that rely on different aspects of stellar physics, it gives us confidence in the new age scale. This agreement is largely due to our adoption of empirical colour-T eff relations and bolometric corrections for pre-main-sequence stars cooler than 4000 K.The revised ages for the star-forming regions in our sample are -∼ 2 Myr for NGC 6611 (Eagle Nebula; M 16), IC 5146 (Cocoon Nebula), NGC 6530 (Lagoon Nebula; M 8), and NGC 2244 (Rosette Nebula); ∼ 6 Myr for σ Ori, Cep OB3b, and IC 348; 10 Myr for λ Ori (Collinder 69); 11 Myr for NGC 2169; 12 Myr for NGC 2362; 13 Myr for NGC 7160; 14 Myr for χ Per (NGC 884); and 20 Myr for NGC 1960 (M 36).

We have selected pre-main-sequence stars in 12 groups of notional ages ranging from 1 Myr to 35 Myrs, using heterogeneous membership criteria. Using these members we have constructed empirical isochrones in V, V-I colour magnitude diagrams (CMDs). This allows us to identify clearly the gap between the radiative main sequence and the convective pre-main-sequence (the R-C gap). We follow the evolution of this gap with age and show that it can be a useful age indicator for groups less than ≃ 15 Myrs old. We also observe a reduction in absolute spreads about the sequences with age. Finally the empirical isochrones allow us to place the groups in order of age, independently of theory. The youngest groups can be collated into three sets of similar ages. The youngest set is the ONC, NGC6530 and IC5146 (nominally 1 Myrs); next Cep OB3b, NGC2362, λ Ori and NGC2264 (nominally 3 Myrs); and finally σ Ori and IC348 (nominally 4-5 Myrs). This suggests Cep OB3b is younger than previously thought, and IC348 older. For IC348 the stellar rotation rate distribution and fraction of stars with discs imply a younger age than we derive. We suggest this is because of the absence of O-stars in this cluster, whose winds and/or ionising radiation may be an important factor in the removal of discs in other clusters.

We present a maximum-likelihood method for fitting two-dimensional model distributions to stellar data in colour-magnitude space. This allows one to include (for example) binary stars in an isochronal population. The method also allows one to derive formal uncertainties for fitted parameters, and assess the likelihood that a good fit has been found. We use the method to derive an age of 38.5 +3.5 −6.5 Myr and a true distance modulus of 7.79 +0.11 −0.05 mag from the V versus V − I diagram of NGC 2547 (the uncertainties are 67 per cent confidence limits, and the parameters are insensitive to the assumed binary fraction). These values are consistent with those previously determined from low-mass isochronal fitting, and are the first measurements to have statistically meaningful uncertainties. The age is also consistent with the lithium depletion age of NGC 2547, and the Hipparcos distance to the cluster is consistent with our value.The method appears to be quite general and could be applied to any N-dimensional data set, with uncertainties in each dimension. However, it is particularly useful when the data are sparse, in the sense that both the typical uncertainties for a data point and the size of structure in the function being fitted are small compared with the typical distance between data points. In this case binning the data will lose resolution, whilst the method presented here preserves it.Software implementing the methods described in this paper is available from