The velocity distribution of stars is a sensitive probe of the gravitational potential of the Galaxy, and hence of its dark matter distribution. In particular, the shape of the dark halo (e.g. spherical, oblate, or prolate) determines velocity correlations, and different halo geometries are expected to result in measurable differences. Here we explore and interpret the correlations in the $(v_R, v_z)$-velocity distribution as a function of position in the Milky Way. We selected a high-quality sample of stars from the Gaia DR2 catalogue and characterised the orientation of the velocity distribution or tilt angle over a radial distance range of $[4-13]~$kpc and up to $3.5~$kpc away from the Galactic plane while taking into account the effects of the measurement errors. We find that the tilt angles change from spherical alignment in the inner Galaxy ($R\sim4~$kpc) towards more cylindrical alignments in the outer Galaxy ($R\sim11~$kpc) when using distances that take a global zero-point offset in the parallax of $-29~\mu$as. However, if the amplitude of this offset is underestimated, then the inferred tilt angles in the outer Galaxy only appear shallower and are intrinsically more consistent with spherical alignment for an offset as large as $-54~\mu$as. We further find that the tilt angles do not seem to strongly vary with Galactic azimuth and that different stellar populations depict similar tilt angles. Therefore we introduce a simple analytic function that describes the trends found over the full radial range. Since the systematic parallax errors in Gaia DR2 depend on celestial position, magnitude, and colour in complex ways, it is not possible to fully correct for them. Therefore it will be particularly important for dynamical modelling of the Milky Way to thoroughly characterise the systematics in astrometry in future Gaia data releases.

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 investigate the kinematics of red clump stars in the Solar Neighbourhood by combining data from TGAS and RAVE to constrain the local dark matter density. After calibrating the absolute magnitude of red clump stars, we characterize their velocity distribution over a radial distance range of $6$-$10$ kpc and up to $1.5$ kpc away from the Galactic plane. We then apply the axisymmetric Jeans equations on subsets representing the thin and thick disks to determine the (local) distribution of mass near the disk of our Galaxy. Our kinematic maps are well-behaved permitting a straightforward local determination of the vertical force, which we find to be $K_z^{\rm thin} = -2454 \pm 619$ and $K_z^{\rm thick} = -2141 \pm 774$ $(\mathrm{km/s})^2 \! /\mathrm{kpc}$ at $1.5$ kpc away from the Galactic plane for the thin and thick disk samples and for thin and thick disk scale heights of 0.28 kpc and 1.12 kpc respectively. These measurements can be translated into a local dark matter density $\rho_\textrm{DM} \sim 0.018 \pm 0.002 M_\odot/\textrm{pc}^3$. The systematic error on this estimate is much larger than the quoted statistical error, since even a 10% difference in the scale height of the thin disk leads to a 30% change in the value of $\rho_\textrm{DM}$, and a nearly equally good fit to the data.