We test the ability of Schwarzschild's orbit superposition method in measuring the mass content, scale radius and shape of a flattened dwarf spheroidal galaxy. Until now, most dynamical model efforts have assumed that dwarf spheroidal galaxies and their host halos are spherical. We use an Evans model (1993) to construct an isothermal mock galaxy whose properties somewhat resemble those of the Sculptor dwarf spheroidal galaxy. This mock galaxy contains flattened luminous and dark matter components, resulting in a logarithmic profile for the global potential. We have tested how well our Schwarzschild method could constrain the characteristic parameters of the system for different sample sizes, and also if the functional form of the potential was unknown. When assuming the true functional form of the potential, the Schwarzschild modelling technique is able to provide an accurate and precise measurement of the characteristic mass parameter of the system and reproduces well the light distribution and the stellar kinematics of our mock galaxy. When assuming a different functional form for the potential, such as a flattened NFW profile, we also constrain the mass and scale radius to their expected values. However in both cases, we find that the flattening parameter remains largely unconstrained. This is likely because the information content of the velocity dispersion on the geometric shape of the potential is too small, since $\sigma$ is constant across our mock dSph. Our results using Schwarzschild's method indicate that the mass enclosed can be derived reliably, even if the flattening parameter is unknown, and already for samples containing 2000 line-of-sight radial velocities, such as those currently available. Further applications of the method to more general distribution functions of flattened systems are needed to establish how well the flattening of dSph dark halos can be determined.

We highlight the power of the Gaia DR2 in studying many fine structures of the Hertzsprung-Russell diagram (HRD). Gaia allows us to present many different HRDs, depending in particular on stellar population selections. We do not aim here for completeness in terms of types of stars or stellar evolutionary aspects. Instead, we have chosen several illustrative examples. We describe some of the selections that can be made in Gaia DR2 to highlight the main structures of the Gaia HRDs. We select both field and cluster (open and globular) stars, compare the observations with previous classifications and with stellar evolutionary tracks, and we present variations of the Gaia HRD with age, metallicity, and kinematics. Late stages of stellar evolution such as hot subdwarfs, post-AGB stars, planetary nebulae, and white dwarfs are also analysed, as well as low-mass brown dwarf objects. The Gaia HRDs are unprecedented in both precision and coverage of the various Milky Way stellar populations and stellar evolutionary phases. Many fine structures of the HRDs are presented. The clear split of the white dwarf sequence into hydrogen and helium white dwarfs is presented for the first time in an HRD. The relation between kinematics and the HRD is nicely illustrated. Two different populations in a classical kinematic selection of the halo are unambiguously identified in the HRD. Membership and mean parameters for a selected list of open clusters are provided. They allow drawing very detailed cluster sequences, highlighting fine structures, and providing extremely precise empirical isochrones that will lead to more insight in stellar physics. Gaia DR2 demonstrates the potential of combining precise astrometry and photometry for large samples for studies in stellar evolution and stellar population and opens an entire new area for HRD-based studies.

The second Gaia data release (DR2), contains very precise astrometric and photometric properties for more than one billion sources, astrophysical parameters for dozens of millions, radial velocities for millions, variability information for half a million of stellar sources and orbits for thousands of solar system objects. Before the Catalogue publication, these data have undergone dedicated validation processes. The goal of this paper is to describe the validation results in terms of completeness, accuracy and precision of the various Gaia DR2 data. The validation processes include a systematic analysis of the Catalogue content to detect anomalies, either individual errors or statistical properties, using statistical analysis, and comparisons to external data or to models. Although the astrometric, photometric and spectroscopic data are of unprecedented quality and quantity, it is shown that the data cannot be used without a dedicated attention to the limitations described here, in the Catalogue documentation and in accompanying papers. A particular emphasis is put on the caveats for the statistical use of the data in scientific exploitation.

Extended stellar haloes are a natural by-product of the hierarchical formation of massive galaxies. If merging is a non-negligible factor in the growth of our Galaxy, evidence of such events should be encoded in its stellar halo. Reliable identification of genuine halo stars is a challenging task however. The 1st Gaia data release contains the positions, parallaxes and proper motions for over 2 million stars, mostly in the Solar neighbourhood. Gaia DR2 will enlarge this sample to over 1.5 billion stars, the brightest ~5 million of which will have a full phase-space information. Our aim is to develop a machine learning model to reliably identify halo stars, even when their full phase-space information is not available. We use the Gradient Boosted Trees algorithm to build a supervised halo star classifier. The classifier is trained on a sample extracted from the Gaia Universe Model Snapshot, convolved with the errors of TGAS, as well as with the expected uncertainties of the upcoming Gaia DR2. We also trained our classifier on the cross-match between the TGAS and RAVE catalogues, where the halo stars are labelled in an entirely model independent way. We then use this model to identify halo stars in TGAS. When full phase- space information is available and for Gaia DR2-like uncertainties, our classifier is able to recover 90% of the halo stars with at most 30% distance errors, in a completely unseen test set, and with negligible levels of contamination. When line-of-sight velocity is not available, we recover ~60% of such halo stars, with less than 10% contamination. When applied to the TGAS data, our classifier detects 337 high confidence RGB halo stars. Although small, this number is consistent with the expectation from models given the data uncertainties. The large parallax errors are the biggest limitation to identify a larger number of halo stars in all the cases studied.

We study the dynamical properties of halo stars located in the Solar Neighbourhood. Our goal is to explore how the properties of the halo depend on the selection criteria used for defining a sample of halo stars. Once this is understood we proceed to measure the shape and orientation of the halo's velocity ellipsoid and we use this information to put constraints on the gravitational potential of the Galaxy. We use the recently released Gaia DR1 catalogue cross-matched to the RAVE dataset for our analysis. We develop a dynamical criterion based on the distribution function of stars in various Galactic components, using action integrals to identify halo members, and compare this to metallicity and to kinematically selected samples. With this new method, we find 1156 stars in the Solar Neighbourhood to be likely members of the stellar halo. Our dynamically selected sample consists mainly of distant giants on elongated orbits. Their metallicity distribution is rather broad, with roughly half of the stars having [M/H] $\ge -1$ dex. The use of different selection criteria has an important impact on the characteristics of the velocity distributions obtained. Nonetheless, for our dynamically selected and for the metallicity selected samples, we find the local velocity ellipsoid to be aligned in spherical coordinates in a Galactocentric reference frame: this suggests that the total gravitational potential is rather spherical in the region spanned by the orbits of the halo stars in these samples.

We present a new Python library called vaex, to handle extremely large tabular datasets, such as astronomical catalogues like the Gaia catalogue, N-body simulations or any other regular datasets which can be structured in rows and columns. Fast computations of statistics on regular N-dimensional grids allows analysis and visualization in the order of a billion rows per second. We use streaming algorithms, memory mapped files and a zero memory copy policy to allow exploration of datasets larger than memory, e.g. out-of-core algorithms. Vaex allows arbitrary (mathematical) transformations using normal Python expressions and (a subset of) numpy functions which are lazily evaluated and computed when needed in small chunks, which avoids wasting of RAM. Boolean expressions (which are also lazily evaluated) can be used to explore subsets of the data, which we call selections. Vaex uses a similar DataFrame API as Pandas, a very popular library, which helps migration from Pandas. Visualization is one of the key points of vaex, and is done using binned statistics in 1d (e.g. histogram), in 2d (e.g. 2d histograms with colormapping) and 3d (using volume rendering). Vaex is split in in several packages: vaex-core for the computational part, vaex-viz for visualization mostly based on matplotlib, vaex-jupyter for visualization in the Jupyter notebook/lab based in IPyWidgets, vaex-server for the (optional) client-server communication, vaex-ui for the Qt based interface, vaex-hdf5 for hdf5 based memory mapped storage, vaex-astro for astronomy related selections, transformations and memory mapped (column based) fits storage. Vaex is open source and available under MIT license on github, documentation and other information can be found on the main website: https://vaex.io, https://docs.vaex.io or https://github.com/maartenbreddels/vaex

The 3D motions of stars in small galaxies beyond our own are minute and yet they are crucial for our understanding of the nature of gravity and dark matter. Even for the dwarf galaxy Sculptor which is one of the best studied systems and inferred to be strongly dark matter dominated, there are conflicting reports on its mean motion around the Milky Way and the 3D internal motions of its stars have never been measured. Here we report, based on data from the Gaia space mission and the Hubble Space Telescope, a new precise measurement of Sculptor's mean proper motion. From this we deduce that Sculptor is currently at its closest approach to the Milky Way and moving on an elongated high-inclination orbit that takes it much farther away than previously thought. For the first time we are also able to measure the internal motions of stars in Sculptor. We find $\sigma_{R}=11.5 \pm 4.3$ km/s and $\sigma_{T}=8.5\pm3.2$ km/s along the projected radial and tangential directions, implying that the stars in our sample move preferentially on radial orbits as quantified by the anisotropy parameter, which we find to be $\beta\sim 0.86^{+0.12}_{-0.83}$ at a location beyond the core radius. Taken at face value such a high radial anisotropy requires abandoning conventional models for the mass distribution in Sculptor. Our sample is dominated by metal-rich stars and for these we find $\beta^{MR} \sim 0.95^{+0.04}_{-0.27}$, a value consistent with multi-component models where Sculptor is embedded in a cuspy dark halo as expected for cold dark matter.

Parallaxes for 331 classical Cepheids, 31 Type II Cepheids and 364 RR Lyrae stars in common between Gaia and the Hipparcos and Tycho-2 catalogues are published in Gaia Data Release 1 (DR1) as part of the Tycho-Gaia Astrometric Solution (TGAS). In order to test these first parallax measurements of the primary standard candles of the cosmological distance ladder, that involve astrometry collected by Gaia during the initial 14 months of science operation, we compared them with literature estimates and derived new period-luminosity ($PL$), period-Wesenheit ($PW$) relations for classical and Type II Cepheids and infrared $PL$, $PL$-metallicity ($PLZ$) and optical luminosity-metallicity ($M_V$-[Fe/H]) relations for the RR Lyrae stars, with zero points based on TGAS. The new relations were computed using multi-band ($V,I,J,K_{\mathrm{s}},W_{1}$) photometry and spectroscopic metal abundances available in the literature, and applying three alternative approaches: (i) by linear least squares fitting the absolute magnitudes inferred from direct transformation of the TGAS parallaxes, (ii) by adopting astrometric-based luminosities, and (iii) using a Bayesian fitting approach. TGAS parallaxes bring a significant added value to the previous Hipparcos estimates. The relations presented in this paper represent first Gaia-calibrated relations and form a "work-in-progress" milestone report in the wait for Gaia-only parallaxes of which a first solution will become available with Gaia's Data Release 2 (DR2) in 2018.

Context. The first Gaia Data Release contains the Tycho-Gaia Astrometric Solution (TGAS). This is a subset of about 2 million stars for which, besides the position and photometry, the proper motion and parallax are calculated using Hipparcos and Tycho-2 positions in 1991.25 as prior information. Aims. We investigate the scientific potential and limitations of the TGAS component by means of the astrometric data for open clusters. Methods. Mean cluster parallax and proper motion values are derived taking into account the error correlations within the astrometric solutions for individual stars, an estimate of the internal velocity dispersion in the cluster, and, where relevant, the effects of the depth of the cluster along the line of sight. Internal consistency of the TGAS data is assessed. Results. Values given for standard uncertainties are still inaccurate and may lead to unrealistic unit-weight standard deviations of least squares solutions for cluster parameters. Reconstructed mean cluster parallax and proper motion values are generally in very good agreement with earlier Hipparcos-based determination, although the Gaia mean parallax for the Pleiades is a significant exception. We have no current explanation for that discrepancy. Most clusters are observed to extend to nearly 15 pc from the cluster centre, and it will be up to future Gaia releases to establish whether those potential cluster-member stars are still dynamically bound to the clusters. Conclusions. The Gaia DR1 provides the means to examine open clusters far beyond their more easily visible cores, and can provide membership assessments based on proper motions and parallaxes. A combined HR diagram shows the same features as observed before using the Hipparcos data, with clearly increased luminosities for older A and F dwarfs.

Before the publication of the Gaia Catalogue, the contents of the first data release have undergone multiple dedicated validation tests. These tests aim at analysing in-depth the Catalogue content to detect anomalies, individual problems in specific objects or in overall statistical properties, either to filter them before the public release, or to describe the different caveats of the release for an optimal exploitation of the data. Dedicated methods using either Gaia internal data, external catalogues or models have been developed for the validation processes. They are testing normal stars as well as various populations like open or globular clusters, double stars, variable stars, quasars. Properties of coverage, accuracy and precision of the data are provided by the numerous tests presented here and jointly analysed to assess the data release content. This independent validation confirms the quality of the published data, Gaia DR1 being the most precise all-sky astrometric and photometric catalogue to-date. However, several limitations in terms of completeness, astrometric and photometric quality are identified and described. Figures describing the relevant properties of the release are shown and the testing activities carried out validating the user interfaces are also described. A particular emphasis is made on the statistical use of the data in scientific exploitation.

The hierarchical structure formation model predicts that stellar halos should form, at least partly, via mergers. If this was a predominant formation channel for the Milky Way's halo, imprints of this merger history in the form of moving groups or streams should exist also in the vicinity of the Sun. Here we study the kinematics of halo stars in the Solar neighbourhood using the very recent first data release from the Gaia mission, and in particular the TGAS dataset, in combination with data from the RAVE survey. Our aim is to determine the amount of substructure present in the phase-space distribution of halo stars that could be linked to merger debris. To characterise kinematic substructure, we measure the velocity correlation function in our sample of halo (low metallicity) stars. We also study the distribution of these stars in the space of energy and two components of the angular momentum, in what we call "Integrals of Motion" space. The velocity correlation function reveals substructure in the form of an excess of pairs of stars with similar velocities, well above that expected for a smooth distribution. Comparison to cosmological simulations of the formation of stellar halos indicate that the levels found are consistent with the Galactic halo having been built fully via accretion. Similarly, the distribution of stars in the space of "Integrals of motion" is highly complex. A strikingly high fraction (between 58% and upto 73%) of the stars that are somewhat less bound than the Sun are on (highly) retrograde orbits. A simple comparison to Milky Way-mass galaxies in cosmological hydrodynamical simulations suggests that less than 1% have such prominently retrograde outer halos. We also identify several other statistically significant structures in "Integrals of Motion" space that could potentially be related to merger events.

With new catalogues arriving such as the Gaia DR1, containing more than a billion objects, new methods of handling and visualizing these data volumes are needed. In visualization, one problem is that the number of datapoints can become so large, that a scatter plot becomes cluttered. Another problem is that with over a billion objects, only a few cpu cycles are available per object if one wants to process them within a second, making traditional methods by rendering glyphs not viable. Instead, we show that by calculating statistics on a regular (N-dimensional) grid, visualizations of a billion objects can be done within a second on a modern desktop computer. This is achieved using memory mapping of hdf5 files together with a simple binning algorithm, which are part of a Python library called vaex. This enables efficient exploration or large datasets interactively, making science exploration of large catalogues feasible. Vaex is a Python library, which also integrates well in the Jupyter/Numpy/Astropy/matplotlib stack. Build on top of this is the vaex application, which allows for interactive exploration and visualization. The motivation for developing vaex is the catalogue of the Gaia satellite, however, vaex can also be used on SPH or N-body simulations, any other (future) catalogues such as SDSS, Pan-STARRS, LSST, WISE, 2MASS, etc. or other tabular data. The homepage for vaex is http://vaex.astro.rug.nl.

We compare the dark matter halos' structural parameters derived for four Milky Way dwarf spheroidal galaxies to those of subhalos found in cosmological $N$-body simulations. We confirm that estimates of the mass at a single fixed radius are fully consistent with the observations. However, when a second structural parameter such as the logarithmic slope of the dark halo density profile measured close to the half-light radius is included in the comparison, we find little to no overlap between the satellites and the subhalos. Typically the right mass subhalos have steeper profiles at these radii than measurements of the dSph suggest. Using energy arguments we explore if it is possible to solve this discrepancy by invoking baryonic effects. Assuming that feedback from supernovae can lead to a reshaping of the halos, we compute the required efficiency and find entirely plausible values for a significant fraction of the subhalos and even as low as 0.1%. This implies that care must be taken not to exaggerate the effect of supernovae feedback as this could make the halos too shallow. These results could be used to calibrate and possibly constrain feedback recipes in hydrodynamical simulations.

We present a proper-motion study on models of the dwarf spheroidal galaxy Sculptor, based on the predicted proper-motion accuracy of Gaia measurements. Gaia will measure proper motions of several hundreds of stars for a Sculptor-like system. Even with an uncertainty on the proper motion of order 1.5 times the size of an individual proper-motion value of ~10 mas/century, we find that it is possible to recover Sculptor's systemic proper motion at its distance of 79 kpc.

In previous work we have presented Schwarzschild models of the Sculptor dSph, demonstrating that this system could be embedded in dark matter halos that are either cusped or cored. Here we show that the non-parametric distribution function recovered through Schwarschild's method is bimodal in energy and angular momentum space for all best fitting mass models explored. We demonstrate that this bimodality is directly related to the two components known to be present in Sculptor through stellar populations analysis, although our method is purely dynamical in nature and does not use this prior information. It therefore constitutes independent confirmation of the existence of two physically distinct dynamical components in Sculptor and suggests a rather complex assembly history for this dwarf galaxy.

We demonstrate the existence of distribution functions that can be used to represent spherical massless cored stellar systems embedded in cuspy dark matter halos with constant mildly tangential velocity anisotropy. In particular, we derive analytically the functional form of the distribution function for a Plummer stellar sphere in a Hernquist dark halo, for \beta_0 = -0.5 and for different degrees of embedding. This particular example satisfies the condition that the central logarithmic slope of the light profile \gamma_0 > 2 \beta_0. Our models have velocity dispersion profiles similar to those observed in nearby dwarf spheroidal galaxies. Hence they can be used to generate initial conditions for a variety of problems, including N-body simulations that may represent dwarf galaxies in the Local Group.

We have developed spherically symmetric dynamical models of dwarf spheroidal galaxies using Schwarzschild's orbit superposition method. This type of modelling yields constraints both on the total mass distribution (e.g. enclosed mass and scale radius) as well as on the orbital structure of the system (e.g. velocity anisotropy). This method is thus less prone to biases introduced by assumptions in comparison to the more commonly used Jeans modelling, and it allows us to derive the dark matter content in a robust way. Here we present our results for the Sculptor dwarf spheroidal galaxy, after testing our methods on mock data sets. We fit both the second and fourth velocity moment profile to break the mass-anisotropy degeneracy. For an NFW dark matter halo profile, we find that the mass of Sculptor within 1 kpc is M_1kpc = 1.03 \pm 0.07 x 10^8 Msol, and that its velocity anisotropy profile is tangentially biased and nearly constant with radius. The preferred concentration (c ~ 15) is low for its dark matter mass but consistent within the scatter found in N-body cosmological simulations. When we let the value of the central logarithmic slope \alpha{} vary, we find that the best-fit model has \alpha{} = 0, although an NFW cusp or shallower is consistent at 1 \sigma{} confidence level. On the other hand, very cuspy density profiles with logarithmic central slopes \alpha{} < -1.5 are strongly disfavoured for Sculptor.

We review our current understanding of the internal dynamical properties of the dwarf spheroidal galaxies surrounding the Milky Way. These are the most dark matter dominated galaxies, and as such may be considered ideal laboratories to test the current concordance cosmological model, and in particular provide constraints on the nature of the dominant form of dark matter. We discuss the latest observations of the kinematics of stars in these systems, and how these may be used to derive their mass distribution. We tour through the various dynamical techniques used, with emphasis on the complementarity and limitations, and discuss what the results imply also in the context of cosmological models. Finally we provide an outlook on exciting developments in this field.

We use orbit based dynamical models to fit the 2nd and 4th moments of the line of sight velocity distributions of the Fornax, Sculptor, Carina and Sextans dwarf spheroidal galaxies. Our goal is to compare dark matter profile models of these four systems using Bayesian evidence. We consider NFW, Einasto and several cored profiles for their dark halos and present the probability distribution functions of the model parameters. When considering each system separately, we find there is no preference for one of these specific parametric density profiles. However, the combined evidence shows that it is unlikely that all galaxies are embedded in the same type of cored profiles of the form \rho_{DM} \propto 1/(1 + r^2)^{\beta/2}, where \beta = 3, 4. For each galaxy, we also obtain an almost model independent, and therefore accurate, measurement of the logarithmic slope of the dark matter density distribution at a radius ~ r_{-3}, i.e. where the logarithmic slope of the stellar density profile is -3. This slope ranges from ~ -1.4 for Fornax to ~ -1.1 for Sextans, both at ~ 1 kpc. All our best fit models essentially have the same mass distribution over a large range in radius (from just below r_{-3} to the last measured data point). This remarkable finding likely implies much stronger constraints on the characteristics that subhalos extracted from cosmological simulations should have in order to host the dSph galaxies around the Milky Way.

We present a new analysis of the Aquarius simulations done in combination with a semi-analytic galaxy formation model. Our goal is to establish whether the subhalos present in LCDM simulations of Milky Way-like systems could host the dwarf spheroidal (dSph) satellites of our Galaxy. Our analysis shows that, contrary to what has been assumed in most previous work, the mass profiles of subhalos are generally not well fit by NFW models but that Einasto profiles are preferred. We find that for shape parameters alpha = 0.2 - 0.5 and Vmax = 10 - 30 km/s there is very good correspondence with the observational constraints obtained for the nine brightest dSph of the Milky Way. However, to explain the internal dynamics of these systems as well as the number of objects of a given circular velocity the total mass of the Milky Way should be ~ 8x10^11 Msun, a value that is in agreement with many recent determinations, and at the low mass end of the range explored by the Aquarius simulations. Our simulations show important scatter in the number of bright satellites, even when the Aquarius Milky Way-like hosts are scaled to a common mass, and we find no evidence for a missing population of massive subhalos in the Galaxy. This conclusion is also supported when we examine the dynamics of the satellites of M31.

We have developed a spherically symmetric dynamical model of a dwarf spheroidal galaxy using the Schwarzschild method. This type of modelling yields constraints both on the total mass distribution (e.g. enclosed mass and scale radius) as well as on the orbital structure of the system modelled (e.g. velocity anisotropy). Therefore not only can we derive the dark matter content of these systems, but also explore possible formation scenarios. Here we present preliminary results for the Sculptor dSph. We find that the mass of Sculptor within 1kpc is 8.5\times10^(7\pm0.05) M\odot, its anisotropy profile is tangentially biased and slightly more isotropic near the center. For an NFW profile, the preferred concentration (~15) is compatible with cosmological models. Very cuspy density profiles (steeper than NFW) are strongly disfavoured for Sculptor.

The RAdial Velocity Experiment (RAVE) is a spectroscopic survey of the Milky Way. We use the subsample of spectra with spectroscopically determined values of stellar parameters to determine the distances to these stars. The list currently contains 235,064 high quality spectra which show no peculiarities and belong to 210,872 different stars. The numbers will grow as the RAVE survey progresses. The public version of the catalog will be made available through the CDS services along with the ongoing RAVE public data releases. The distances are determined with a method based on the work by Breddels et al.~(2010). Here we assume that the star undergoes a standard stellar evolution and that its spectrum shows no peculiarities. The refinements include: the use of either of the three isochrone sets, a better account of the stellar ages and masses, use of more realistic errors of stellar parameter values, and application to a larger dataset. The derived distances of both dwarfs and giants match within ~21% to the astrometric distances of Hipparcos stars and to the distances of observed members of open and globular clusters. Multiple observations of a fraction of RAVE stars show that repeatability of the derived distances is even better, with half of the objects showing a distance scatter of \simlt 11%. RAVE dwarfs are ~300 pc from the Sun, and giants are at distances of 1 to 2 kpc, and up to 10 kpc. This places the RAVE dataset between the more local Geneva-Copenhagen survey and the more distant and fainter SDSS sample. As such it is ideal to address some of the fundamental questions of Galactic structure and evolution in the pre-Gaia era. Individual applications are left to separate papers, here we show that the full 6-dimensional information on position and velocity is accurate enough to discuss the vertical structure and kinematic properties of the thin and thick disks.

We study the eccentricity distribution of a thick disc sample of stars observed in the Radial Velocity Experiment (RAVE) and compare it to that expected in four simulations of thick disc formation in the literature (accretion of satellites, heating of a primordial thin disc during a merger, radial migration, and gas-rich mergers), as compiled by Sales et al. (2009). We find that the distribution of our sample is peaked at low eccentricities and falls off smoothly and rather steeply to high eccentricities. This distribution is fairly robust to changes in distances, thin disc contamination, and the particular thick disc sample used. Our results are inconsistent with what is expected for the pure accretion simulation, since we find that the dynamics of local thick disc stars implies that the majority must have formed "in situ". Of the remaining models explored, the eccentricity distribution of our stars appears to be most consistent with the gas-rich merger case.

(Abridged) Aims:We develop a method for deriving distances from spectroscopic data and obtaining full 6D phase-space coordinates for the RAVE survey's second data release. Methods: We used stellar models combined with atmospheric properties from RAVE (Teff, logg and [Fe/H]) and (J-Ks) photometry from archival sources to derive absolute magnitudes. We are able to derive the full 6D phase-space coordinates for a large sample of RAVE stars. This method is tested with artificial data, Hipparcos trigonometric parallaxes and observations of the open cluster M67. Results: When we applied our method to a set of 16 146 stars, we found that 25% (4 037) of the stars have relative (statistical) distance errors of < 35%, while 50% (8 073) and 75% (12 110) have relative (statistical) errors smaller than 45% and 50%, respectively. Our various tests show that we can reliably estimate distances for main-sequence stars, but there is an indication of potential systematic problems with giant stars. For the main-sequence star sample (defined as those with log(g) > 4), 25% (1 744) have relative distance errors < 31%, while 50% (3 488) and 75% (5 231) have relative errors smaller than 36% and 42%, respectively. Our full dataset shows the expected decrease in the metallicity of stars as a function of distance from the Galactic plane. The known kinematic substructures in the U and V velocity components of nearby dwarf stars are apparent in our dataset, confirming the accuracy of our data and the reliability of our technique. We provide independent measurements of the orientation of the UV velocity ellipsoid and of the solar motion, and they are in very good agreement with previous work. Conclusions: The distance catalogue for the RAVE second data release is available at http://www.astro.rug.nl/~rave