The density of dark matter near the Sun, ρ DM, , is important for experiments hunting for dark matter particles in the laboratory, and for constraining the local shape of the Milky Way's dark matter halo. Estimates to date have typically assumed that the Milky Way's stellar disc is axisymmetric and in a steady-state. Yet the Milky Way disc is neither, exhibiting prominent spiral arms and a bar, and vertical and radial oscillations. We assess the impact of these assumptions on determinations of ρ DM, by applying a free-form, steady-state, Jeans method to two different N -body simulations of Milky Way-like galaxies. In one, the galaxy has experienced an ancient major merger, similar to the hypothesized Gaia-Sausage-Enceladus; in the other, the galaxy is perturbed more recently by the repeated passage and slow merger of a Sagittarius-like dwarf galaxy. We assess the impact of each of the terms in the Jeans-Poisson equations on our ability to correctly extract ρ DM, from the simulated data. We find that common approximations employed in the literature -axisymmetry and a locally flat rotation curve -can lead to significant systematic errors of up to a factor ∼ 1.5 in the recovered surface mass density ∼ 2 kpc above the disc plane, implying a fractional error on ρ DM, of order unity. However, once we add in the tilt term and the rotation curve term in our models, we obtain an unbiased estimate of ρ DM, , consistent with the true value within our 95% confidence intervals for realistic 20% uncertainties on the baryonic surface density of the disc. Other terms -the axial tilt, 2:nd Poisson and time dependent terms -contribute less than 10% to ρ DM, (given current data) and can be safely neglected for now. In the future, as more data become available, these terms will need to be included in the analysis. Recent estimates of ρDM, make use of a variety of meth-ods (e.g. de Salas & Widmark 2020). 'Local' methods use Jeans and/or distribution function models over small volumes near the Sun (e.g. Xia et al. 2016;Sivertsson et al. 2018;Guo et al. 2020; Salomon et al. 2020). 'Global' methods use stellar tracers over larger volumes (e.g. Cole & Binney 2017;Wegg et al. 2019), or use the Milky Way's circular velocity, assuming that its dark matter halo is spherically symmetric (e.g. Benito et al. 2019;Karukes et al. 2019;de Salas et al. 2019). More recently, newer methods have been introduced, like 'made to measure ' (e.g. Syer & Tremaine 1996;Bovy et al. 2018), and a recent attempt to harness the Milky Way disc's disequilibrium to estimate ρDM, (Widmark et al. 2021). Some studies exploit the vertical kinematics of stars up to a height z ∼ 1 − 3 kpc above the Galactic disc (e.g. Xia et al. 2016;Sivertsson et al. 2018;Guo et al. 2020; Salomon et al. 2020), while others use a similar methodology much closer to the plane (< 200 pc;