Adiabatic Contraction Research Articles

Direct Collapse Supermassive Black Holes from Relic Particle Decay.

We investigate the formation of high-redshift supermassive black holes (SMBHs) via the direct collapse of baryonic clouds, where the unwanted formation of molecular hydrogen is successfully suppressed by a Lyman-Werner (LW) photon background from relic particle decay. We improve on existing studies by dynamically simulating the collapse, accounting for the adiabatic contraction of the DM halo, as well as the insitu production of the LW photons within the cloud which reduce the impact of the cloud's shielding. We find a viable parameter space where the decay of either some of the dark matter or all of a subdominant decaying species successfully allows direct collapse of the cloud to a SMBH.

Cold Dark Matter and Self-interacting Dark Matter Interpretations of the Strong Gravitational Lensing Object JWST-ER1

van Dokkum et al. reported the discovery of JWST-ER1, a strong lensing object at redshift z ≈ 2, using data from the James Webb Space Telescope. The lens mass within the Einstein ring is 5.9 times higher than the expected stellar mass from a Chabrier initial mass function, indicating a high dark matter density. In this work, we show that a cold dark matter halo, influenced by gas-driven adiabatic contraction, can account for the observed lens mass. We interpret the measurement of JWST-ER1 in the self-interacting dark matter scenario and show that the cross section per particle mass σ/m ≈ 0.1 cm2 g−1 is generally favored. Intriguingly, σ/m ≈ 0.1 cm2 g−1 can also be consistent with the strong lensing observations of early-type galaxies at redshift z ≈ 0.2, where adiabatic contraction is not observed overall.

Detectability of Supermassive Dark Stars with the Roman Space Telescope

Supermassive dark stars (SMDS) are luminous stellar objects formed in the early Universe at redshift z ∼ 10–20, made primarily of hydrogen and helium, yet powered by dark matter. We examine the capabilities of the Roman Space Telescope (RST), and find it able to identify ∼106 M ⊙ SMDSs at redshifts up to z ≃ 14. With a gravitational lensing factor of μ ∼ 100, RST could identify SMDS as small as ∼104 M ⊙ at z ∼ 12 with ∼106 s exposure. Differentiating SMDSs from early galaxies containing zero metallicity stars at similar redshifts requires spectral, photometric, and morphological comparisons. With only RST, the differentiation of SMDS, particularly those formed via adiabatic contraction with M ≳ 105 M ⊙ and lensed by μ ≳ 100, is possible due to their distinct photometric signatures from the first galaxies. Those formed via dark matter capture can be differentiated only by image morphology: i.e., point object (SMDSs) versus extended object (sufficiently magnified galaxies). By additionally employing James Webb Space Telescope (JWST) spectroscopy, we can identify the He ii λ1640 absorption line, a smoking gun for SMDS detection. Although RST does not cover the required wavelength band (for z emi ≳ 10), JWST does; hence, the two can be used in tandem to identify SMDS. The detection of SMDS would confirm a new type of star powered by dark matter and may shed light on the origins of the supermassive black holes powering bright quasars observed at z ≳ 6.

Modeling electron acceleration during the contraction of a magnetic island

Magnetic reconnection releases the magnetic energy through the contraction of multi-magnetic island leading to the electron acceleration as proposed by Drake et. al in 2006. However, how the released magnetic energy is converted into electron’s kinetic energy is still theoretically not well understood. We model in particular the kinetic process assuming the adiabatic contraction of magnetic island that induces electric field which is proportional to the vector potential of the magnetic island and approximate the magnetic island with an ellipse. Under this model, we show that the energy gain is achieved through the work of inductive electric field. We further show that the curvature drift which is along the inductive electric field dominates the energy gain. We compared our model with the magnetic island formed by tearing instability in a 2.5D particle-in-cell simulation of magnetic reconnection and found the results from the model consistent with that of the simulation.

MaNGA DynPop – VI. Matter density slopes from dynamical models of 6000 galaxies versus cosmological simulations: the interplay between baryonic and dark matter

ABSTRACT We try to understand the trends in the mass density slopes as a function of galaxy properties. We use the results from the best Jeans Anisotropic Modelling (JAM) of the integral-field stellar kinematics for near 6000 galaxies from the MaNGA DynPop project, with stellar masses $10^9\ {\rm {\rm M}_{\odot }}\lesssim M_*\lesssim 10^{12}\ {\rm {\rm M}_{\odot }}$, including both early-type and late-type galaxies. We use the mass-weighted density slopes for the stellar $\overline{\gamma }_*$, dark $\overline{\gamma }_{_{\rm DM}}$ and total $\overline{\gamma }_{_{\rm T}}$ mass from the MaNGA DynPop project. As previously reported, $\overline{\gamma }_{_{\rm T}}$ approaches a constant value of $\overline{\gamma }_{_{\rm T}}\approx 2.2$ for high σe galaxies, and flattens for $\lg (\sigma _{\rm e}/{\rm km\ s^{-1}})\lesssim 2.3$ galaxies, reaching $\overline{\gamma }_{_{\rm T}}\approx 1.5$ for $\lg (\sigma _{\rm e}/{\rm km\ s^{-1}})\approx 1.8$. We find that total and stellar slopes track each other tightly, with $\overline{\gamma }_{_{\rm T}}\approx \overline{\gamma }_*-0.174$ over the full σe range. This confirms the dominance of stellar matter within Re. We also show that there is no perfect conspiracy between baryonic and dark matter, as $\overline{\gamma }_*$ and $\overline{\gamma }_{_{\rm DM}}$ do not vary inversely within the σe range. We find that the central galaxies from TNG50 and TNG100 simulations do not reproduce the observed galaxy mass distribution, which we attribute to the overestimated dark matter fraction, possibly due to a constant IMF and excessive adiabatic contraction effects in the simulations. Finally, we present the stacked dark matter density profiles and show that they are slightly steeper than the pure dark matter simulation prediction of $\overline{\gamma }_{_{\rm DM}}\approx 1$, suggesting moderate adiabatic contraction in the central region of galaxies. Our work demonstrates the power of stellar dynamics modelling for probing the interaction between stellar and dark matter and testing galaxy formation theories.

The birth and early evolution of a low-mass protostar

Context. Understanding the collapse of dense molecular cloud cores to stellar densities and the subsequent evolution of the protostar is of importance to model the feedback effects such an object has on its surrounding environment, as well as describing the conditions with which it enters the stellar evolutionary track. This process is fundamentally multi-scale, both in density and in spatial extent, and requires the inclusion of complex physical processes such as self-gravity, turbulence, radiative transfer, and magnetic fields. As such, it necessitates the use of robust numerical simulations. Aims. We aim to model the birth and early evolution of a low-mass protostar. We also seek to describe the interior structure of the protostar and the radiative behavior of its accretion shock front. Methods. We carried out a high resolution numerical simulation of the collapse of a gravitationally unstable 1 M⊙ dense molecular cloud core to stellar densities using 3D radiation hydrodynamics under the gray flux-limited diffusion approximation. We followed the initial isothermal phase, the first adiabatic contraction, the second gravitational collapse triggered by the dissociation of H2 molecules, and ≈247 days of the subsequent main accretion phase. Results. We find that the subcritical radiative behavior of the protostar’s shock front causes it to swell as it accretes matter. We also find that the protostar is turbulent from the moment of its inception despite its radiative stability. This turbulence causes significant entropy mixing inside the protostar, which regulates the swelling. Furthermore, we find that the protostar is not fully ionized at birth, but the relative amount of ionized material within it increases as it accretes matter from its surroundings. Finally, we report in the appendix the results of the first 3D calculations involving a frequency-dependent treatment of radiative transfer, which has not produced any major differences with its gray counterpart.

The accretion history of the Milky Way: III. Hydrodynamical simulations of Galactic dwarf galaxies at first infall

ABSTRACT Most Milky Way dwarf galaxies are much less bound to their host than are relics of Gaia–Sausage–Enceladus and Sgr. These dwarfs are expected to have fallen into the Galactic halo less than 3 Gyr ago, and will therefore have undergone no more than one full orbit. Here, we have performed hydrodynamical simulations of this process, assuming that their progenitors are gas-rich, rotation-supported dwarfs. We follow their transformation through interactions with the hot corona and gravitational field of the galaxy. Our dedicated simulations reproduce the structural properties of three dwarf galaxies: Sculptor, Antlia II, and with somewhat a lower accuracy, Crater II. This includes reproducing their large velocity dispersions, which are caused by ram-pressure stripping and Galactic tidal shocks. Differences between dwarfs can be interpreted as due to different orbital paths, as well as to different initial conditions for their progenitor gas and stellar contents. However, we failed to suppress in a single orbit the rotational support of our Sculptor analogue if it is fully dark matter dominated. In addition, we have found that classical dwarf galaxies like Sculptor may have stellar cores sufficiently dense to survive the pericentre passage through adiabatic contraction. On the contrary, our Antlia II and Crater II analogues are tidally stripped, explaining their large sizes, extremely low surface brightnesses, and velocity dispersion. This modelling explains differences between dwarf galaxies by reproducing them as being at different stages of out-of-equilibrium stellar systems.

On the effect of angular momentum on the prompt cusp formation via the gravitational collapse

In this work, we extend the model proposed by White concerning the post-collapse evolution of density peaks while considering the role of angular momentum. On a timescale smaller than the peak collapse, $t_{0}$, the inner regions of the peak reach the equilibrium forming a cuspy profile, as in White's paper, but the power-law density profile is flatter, namely $\rho \propto r^{-1.52}$, using the specific angular momentum $J$ obtained in theoretical models of how it evolves in CDM universes, namely $J \propto M^{2/3}$. The previous result shows how angular momentum influences the slope of the density profile, and how a slightly flatter profile obtained in high-resolution numerical simulations, namely $\rho \propto r^{\alpha}$, $(\alpha \simeq -1.5)$ can be reobtained. Similarly to simulations, in our model adiabatic contraction was not taken into account. This means that more comprehensive simulations could give different values for the slope of the density profile, similar to an improvement of our model.

Slowly rotating Tolman VII solution

We present a model of a slowly rotating Tolman VII (T-VII) fluid sphere, at second order in the angular velocity. The structure of this configuration is obtained by integrating numerically the Hartle–Thorne equations for slowly rotating relativistic masses. We consider a sequence of models where we vary the parameter , where R is the radius of the configuration and is its Schwarzschild radius, representing an adiabatic and quasi-stationary contraction by progressively reducing the radius while keeping the angular momentum and gravitational mass constant. We determined the moment of inertia I, mass quadrupole moment Q, and the ellipticity ɛ, for various configurations. Similarly to previous results for Maclaurin and polytropic spheroids, in slow rotation, we found a change in the behavior of the ellipticity when reaches a certain critical value. Based on our analysis for the T-VII solution, we found variations of in the and relations, and variation in the I − Q relation, with respect to the universal fittings proposed for realistic neutron stars (NSs). Our results suggest that the T-VII solution can be considered a rather good approximation for the description of the interior of NSs.

A test of invariance of dark matter halo surface density using multiwavelength mock galaxy catalogues

ABSTRACT A large number of observations have shown that the dark matter halo surface density, given by the product of halo core radius and core density, is nearly constant for a diverse suite of galaxies. Although this invariance of the halo surface density is violated at galaxy cluster and group scales, it is still an open question on whether the aforementioned constancy on galactic scales can be explained within Lambda cold dark matter (ΛCDM). For this purpose, we probe the variation of halo surface density as a function of mass using multiwavelength mock galaxy catalogues from ΛCDM simulations, where the adiabatic contraction of dark matter haloes in the presence of baryons has been taken into account. We find that these baryonified ΛCDM haloes were best fitted with a generalized Navarro–Frenk–White profile, and the halo surface density from these haloes has a degeneracy with respect to both the halo mass and the virial concentration. We find that the correlation with mass when averaged over concentration is consistent with a constant halo surface density. However, a power-law dependence as a function of halo mass also cannot be ruled out.

A semi-analytic study of self-interacting dark-matter haloes with baryons

ABSTRACT We combine the isothermal Jeans model and the model of adiabatic halo contraction into a semi-analytic procedure for computing the density profile of self-interacting dark-matter (SIDM) haloes with the gravitational influence from the inhabitant galaxies. The model agrees well with cosmological SIDM simulations over the entire core-forming stage up to the onset of gravothermal core-collapse. Using this model, we show that the halo response to baryons is more diverse in SIDM than in CDM and depends sensitively on galaxy size, a desirable feature in the context of the structural diversity of bright dwarfs. The fast speed of the method facilitates analyses that would be challenging for numerical simulations – notably, we quantify the SIDM halo response as functions of the baryonic properties, on a fine mesh grid spanned by the baryon-to-total-mass ratio, Mb/Mvir, and galaxy compactness, r1/2/Rvir; we show with high statistical precision that for typical Milky-Way-like systems, the SIDM profiles are similar to their CDM counterparts; and we delineate the regime of core-collapse in the Mb/Mvir − r1/2/Rvir space, for a given cross section and concentration. Finally, we compare the isothermal Jeans model with the more sophisticated gravothermal fluid model, and show that the former yields faster core formation and agrees better with cosmological simulations. We attribute the difference to whether the target CDM halo is used as a boundary condition or as the initial condition for the gravothermal evolution, and thus comment on possible improvements of the fluid model. We have made our model publicly available at https://github.com/JiangFangzhou/SIDM.

Cosmic Reionization on Computers: Baryonic Effects on Halo Concentrations during the Epoch of Reionization

Baryons both increase halo concentration through adiabatic contraction and expel mass through feedback processes. However, it is not well understood how the radiation fields prevalent during the epoch of reionization affect the evolution of concentration in dark matter halos. We investigate how baryonic physics during the epoch of reionization modify the structure of dark matter halos in the Cosmic Reionization On Computers (CROC) simulations. We use two different measures of halo concentration to quantify the effects. We compare concentrations of halos matched between full-physics simulations and dark-matter-only simulations with identical initial conditions between 5 ≤ z ≤ 9. Baryons in full-physics simulations do pull matter toward the center, increasing the maximum circular velocity compared to dark-matter-only simulations. However, their overall effects are much less than if all the baryons were simply centrally concentrated indicating that heating processes efficiently counteract cooling effects. Finally, we show that the baryonic effects on halo concentrations at z ≈ 5 are relatively insensitive to environmental variations of reionization history. These results are pertinent to models of galaxy–halo connection during the epoch of reionization.

Splashback Radius in a Spherical Collapse Model

It was shown several years ago that dark matter halo outskirts are characterized by very steep density profiles in a very small radial range. This feature has been interpreted as a pile-up of different particle orbits at a similar location, namely, splashback material at half an orbit after collapse. Adhikari et al. (2014) obtained the location of the splashback radius through a very simple model by calculating a dark matter shell trajectory in the secondary infall model while it crosses a growing NFW profile-shaped dark matter halo. Because they imposed a halo profile instead of calculating it from the trajectories of the shells of dark matter, they were not able to find the dark matter profile around the splashback radius. In the present paper, we use an improved spherical infall model taking into account shell crossing as well as several physical effects such as ordered and random angular momentum, dynamical friction, adiabatic contraction, etc. This allows us to determine the density profile from the inner to the outer region and to study the behavior of the outer density profile. We compare the density profiles and their logarithmic slope of with the simulation results of Diemer and Kravtsov (2014), finding a good agreement between the prediction of the model and the simulations.

The impact of gas disc flaring on rotation curve decomposition and revisiting baryonic and dark matter relations for nearby galaxies

ABSTRACT Gas discs of late-type galaxies are flared, with scale heights increasing with the distance from the galaxy centres and often reaching kpc scales. We study the effects of gas disc flaring on the recovered dark matter halo parameters from rotation curve decomposition. For this, we carefully select a sample of 32 dwarf and spiral galaxies with high-quality neutral gas, molecular gas, and stellar mass profiles, robust H i rotation curves obtained via 3D kinematic modelling, and reliable bulge-disc decomposition. By assuming vertical hydrostatic equilibrium, we derive the scale heights of the atomic and molecular gas discs and fit dark matter haloes to the rotation curves self-consistently. We find that the effect of the gas flaring in the rotation curve decomposition can play an important role only for the smallest, gas-dominated dwarfs, while for most of the galaxies, the effect is minor and can be ignored. We revisit the stellar- and baryon-to-halo mass relations (M*–M200 and Mbar–M200). Both relations increase smoothly up to $M_{200} \approx 10^{12}~\rm { M_\odot }$, with galaxies at this end having high M*/M200 and Mbar/M200 ratios approaching the cosmological baryon fraction. At higher M200, the relations show a larger scatter. Most haloes of our galaxy sample closely follow the concentration–mass (c200–M200) relation resulting from N-body cosmological simulations. Interestingly, the galaxies deviating above and below the relation have the highest and lowest stellar and baryon factions, respectively, which suggests that the departures from the c200–M200 law are regulated by adiabatic contraction and an increasing importance of feedback.

Orbits and adiabatic contraction in scalar-field dark matter halos: revisiting the cusp-core problem in dwarf galaxies

ABSTRACT Bose–Einstein-condensed dark matter, also called scalar field dark matter (SFDM), has become a popular alternative to cold dark matter (CDM), because it predicts galactic cores, in contrast to the cusps of CDM halos (‘cusp-core problem’). We continue the study of SFDM with a strong, repulsive self-interaction; the Thomas–Fermi (TF) regime of SFDM (SFDM-TF). In this model, structure formation is suppressed below a scale related to the TF radius RTF, which is close to the radius of central cores in these halos. We investigate for the first time the impact of baryons onto realistic galactic SFDM-TF halo profiles by studying the process of adiabatic contraction (AC) in such halos. In doing so, we first analyse the underlying quantum Hamilton–Jacobi framework appropriate for SFDM and calculate dark matter orbits, in order to verify the validity of the assumptions usually required for AC. Then, we calculate the impact of AC onto SFDM-TF halos of mass $\sim 10^{11}\, {\rm M}_{\odot }$, with various baryon fractions and core radii, RTF ∼ (0.1–4) kpc, and compare our results with observational velocity data of dwarf galaxies. We find that AC-modified SFDM-TF halos with kpc-size core radii reproduce the data well, suggesting stellar feedback may not be required. On the other hand, halos with sub-kpc core radii face the same issue than CDM, in that they are not in accordance with galaxy data in the central halo parts.

Baryon-driven decontraction in Milky Way-mass haloes

ABSTRACT We select a sample of Milky Way (MW) mass haloes from a high-resolution version of the EAGLE simulation to study their inner dark matter (DM) content and how baryons alter it. As in previous studies, we find that all haloes are more massive at the centre compared to their dark matter-only (DMO) counterparts at the present day as a result of the dissipational collapse of baryons during the assembly of the galaxy. However, we identify two processes that can reduce the central halo mass during the evolution of the galaxy. First, gas blowouts induced by active galactic nuclei feedback can lead to a substantial decrease of the central DM mass. Secondly, the formation of a stellar bar and its interaction with the DM can induce a secular expansion of the halo; the rate at which DM is evacuated from the central region by this process is related to the average bar strength, and the time-scale on which it acts determines how much the halo has decontracted. Although the inner regions of the haloes we have investigated are still more massive than their DMO counterparts at z = 0, they are significantly less massive than in the past and less massive than expected from the classic adiabatic contraction model. Since the MW has both a central supermassive black hole and a bar, the extent to which its halo has contracted is uncertain. This may affect estimates of the mass of the MW halo and of the expected signals in direct and indirect DM detection experiments.

Constraining dark matter properties with the first generation of stars

Dark Matter (DM) can be trapped by the gravitational field of any star, since collisions with nuclei in dense environments can slow down the DM particle below the escape velocity ($v_{esc}$) at the surface of the star. If captured, the DM particles can self-annihilate, and, therefore, provide a new source of energy for the star. We investigate this phenomenon for capture of DM particles by the first generation of stars [Population III (Pop III) stars], by using the multiscatter capture formalism. Pop III stars are particularly good DM captors, since they form in DM-rich environments, at the center of$~\sim 10^6 M_\odot$ DM minihalos, at redshifts $z\sim 15$. Assuming a DM-proton scattering cross section ($\sigma)$ at the current deepest exclusion limits provided by the XENON1T experiment, we find that captured DM annihilations at the core of Pop III stars can lead, via the Eddington limit, to upper bounds in stellar masses that can be as low as a few $M_\odot$ if the ambient DM density ($\rho_X$) at the location of the Pop III star is sufficiently high. Conversely, when Pop III stars are identified, one can use their observed mass ($M_\star$) to place bounds on $\rho_X\sigma$. Using adiabatic contraction to estimate the ambient DM density in the environment surrounding Pop III stars, we place projected upper limits on $\sigma$, for $M_\star$ in the $100-1000~M_\odot$ range, and find bounds that are competitive with, or deeper than, those provided by the most sensitive current direct detection experiments for both spin independent and spin dependent interactions, for a wide range of DM masses. Most intriguingly, we find that Pop III stars with mass $M_\star \gtrsim 300 M_\odot$ could be used to probe the SD proton-DM cross section below the "neutrino floor," i.e. the region of parameter space where DM direct detection experiments will soon become overwhelmed by neutrino backgrounds.

Dark matter haloes of massive elliptical galaxies at z ∼ 0.2 are well described by the Navarro–Frenk–White profile

ABSTRACT We investigate the internal structure of elliptical galaxies at z ∼ 0.2 from a joint lensing–dynamics analysis. We model Hubble Space Telescope images of a sample of 23 galaxy–galaxy lenses selected from the Sloan Lens ACS (SLACS) survey. Whereas the original SLACS analysis estimated the logarithmic slopes by combining the kinematics with the imaging data, we estimate the logarithmic slopes only from the imaging data. We find that the distribution of the lensing-only logarithmic slopes has a median 2.08c ± 0.03 and intrinsic scatter 0.13 ± 0.02, consistent with the original SLACS analysis. We combine the lensing constraints with the stellar kinematics and weak lensing measurements, and constrain the amount of adiabatic contraction in the dark matter (DM) haloes. We find that the DM haloes are well described by a standard Navarro–Frenk–White halo with no contraction on average for both of a constant stellar mass-to-light ratio (M/L) model and a stellar M/L gradient model. For the M/L gradient model, we find that most galaxies are consistent with no M/L gradient. Comparison of our inferred stellar masses with those obtained from the stellar population synthesis method supports a heavy initial mass function (IMF) such as the Salpeter IMF. We discuss our results in the context of previous observations and simulations, and argue that our result is consistent with a scenario in which active galactic nucleus feedback counteracts the baryonic-cooling-driven contraction in the DM haloes.

Velocity-dependent self-interacting dark matter from groups and clusters of galaxies

We probe the self-interactions of dark matter with relaxed galaxy groups and clusters using observational data from strong and weak lensing and stellar kinematics. Our analysis uses the Jeans formalism and considers a wider range of systematic effects than in previous work, including adiabatic contraction and stellar anisotropy, to robustly constrain the self-interaction cross section. For both groups and clusters, our results show a mild preference for a nonzero cross section compared with cold collisionless dark matter. Our groups result, σ/m=0.5±0.2 cm2/g, places the first constraint on self-interacting dark matter (SIDM) at an intermediate scale M200∼ 1014 M⊙ between galaxies and massive clusters. For massive clusters with M200∼ 1015 M⊙, our result is σ/m=0.19±0.09 cm2/g, with an upper limit of σ / m < 0.35 cm2/g (95% CL) . Thus, our results disfavor a velocity-independent cross section of order 1 cm2/g or larger needed to impact small scale structure problems in galaxies, but are consistent with a velocity-dependent cross section that decreases with increasing scattering velocity. Comparing the cross sections with and without the effect of adiabatic contraction, we find that adiabatic contraction produces slightly larger values for our data sample, but they are consistent at the 1σ level. Finally, to validate our approach, we apply our Jeans analysis to a sample of mock data generated from SIDM-plus-baryons simulations with σ/m = 1 cm2/g. This is the first test of the Jeans model at the level of stellar and lensing observables directly measured from simulations. We find our analysis gives a robust determination of the cross section, as well as consistently inferring the true baryon and dark matter density profiles.

The orbital phase space of contracted dark matter haloes

ABSTRACT We study the orbital phase space of dark matter (DM) haloes in the auriga suite of cosmological hydrodynamics simulations of Milky Way (MW) analogues. We characterize haloes by their spherical action distribution, $F\left(J_{{r}},L\right)$, a function of the specific angular momentum, L, and the radial action, Jr, of the DM particles. By comparing DM-only and hydrodynamical simulations of the same haloes, we investigate the contraction of DM haloes caused by the accumulation of baryons at the centre. We find a small systematic suppression of the radial action in the DM haloes of the hydrodynamical simulations, suggesting that the commonly used adiabatic contraction approximation can result in an underestimate of the density by $\sim 8{{ \rm {per\ cent}}}$. We apply an iterative algorithm to contract the auriga DM haloes given a baryon density profile and halo mass, recovering the true contracted DM profiles with an accuracy of $\sim 15{{ \rm {per\ cent}}}$, that reflects halo-to-halo variation. Using this algorithm, we infer the total mass profile of the MW’s contracted DM halo. We derive updated values for the key astrophysical inputs to DM direct detection experiments: the DM density and velocity distribution in the Solar neighbourhood.