We apply a new method to determine the local disc matter and dark halo matter density to kinematic and position data for \sim2000 K dwarf stars taken from the literature. Our method assumes only that the disc is locally in dynamical equilibrium, and that the 'tilt' term in the Jeans equations is small up to \sim1 kpc above the plane. We present a new calculation of the photometric distances to the K dwarf stars, and use a Monte Carlo Markov Chain to marginalise over uncertainties in both the baryonic mass distribution, and the velocity and distance errors for each individual star. We perform a series of tests to demonstrate that our results are insensitive to plausible systematic errors in our distance calibration, and we show that our method recovers the correct answer from a dynamically evolved N-body simulation of the Milky Way. We find a local dark matter density of {\rho}dm = 0.025+0.014-0.013 M\odotpc^{-3} (0.95+0.53-0.49 GeV cm^{-3}) at 90% confidence assuming no correction for the non-flatness of the local rotation curve, and {\rho}dm = 0.022+0.015-0.013 M\odotpc^-3 (0.85+0.57-0.50 GeV cm^{-3}) if the correction is included. Our 90% lower bound on {\rho}dm is larger than the canonical value typically assumed in the literature, and is at mild tension with extrapolations from the rotation curve that assume a spherical halo. Our result can be explained by a larger normalisation for the local Milky Way rotation curve, an oblate dark matter halo, a local disc of dark matter, or some combination of these.

We study the systematic problems in determining the local dark matter density $\rho_{dm}(R_\odot)$ from kinematics of stars in the Solar Neighbourhood, using a simulated Milky Way-like galaxy. We introduce a new unbiased method for recovering $\rho_{dm}(R_\odot)$ based on the moments of the Jeans equations, combined with a Monte Carlo Markov Chain technique and apply it to real data. ov

We revisit systematics in determining the local dark matter density (rho_dm) from the vertical motion of stars in the Solar Neighbourhood. Using a simulation of a Milky Way-like galaxy, we determine the data-quality required to detect the dark matter density at its expected local value. We introduce a new method for recovering rho_dm that uses moments of the Jeans equations, combined with a Monte Carlo Markov Chain technique to marginalise over the unknown parameters. Given sufficiently good data, we show that our method can recover the correct local dark matter density even in the face of disc inhomogeneities, non-isothermal tracers and a non-separable distribution function. We illustrate the power of our technique by applying it to Hipparcos data [Holmberg & Flynn 2000,2004]. We first make the assumption that the A and F star tracer populations are isothermal. This recovers rho_dm=0.003^{+0.009}_{-0.007}Msun/pc^3 (with 90 per cent confidence), consistent with previous determinations. However, the vertical dispersion profile of these tracers is poorly known. If we assume instead a non-isothermal profile similar to the blue disc stars from SDSS DR-7 [Abazajian et al. 2009] measured by Bond et al. (2009), we obtain a fit with a very similar chi^2 value, but with rho_dm=0.033^{+0.008}_{-0.009}Msun/pc^3 (with 90 per cent confidence). This highlights that it is vital to measure the vertical dispersion profile of the tracers to recover an unbiased estimate of the local dark matter density.

We examine systematic problems in determining the local matter density from the vertical motion of stars, i.e. the 'Oort limit'. Using collisionless simulations and a Monte Carlo Markov Chain technique, we determine the data quality required to detect local dark matter at its expected density. We find that systematic errors are more important than observational errors and apply our technique to Hipparcos data to reassign realistic error bars to the local dark matter density.