Source author record

Michael Kuhlen

Michael Kuhlen appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

Catalog footprint

What is connected

20works
7topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

20 published item(s)

preprint2014arXiv

Measuring the power spectrum of dark matter substructure using strong gravitational lensing

In recent years, it has become possible to detect individual dark matter subhalos near images of strongly lensed extended background galaxies. Typically, only the most massive subhalos in the strong lensing region may be detected this way. In this work, we show that strong lenses may also be used to constrain the much more numerous population of lower mass subhalos that are too small to be detected individually. In particular, we show that the power spectrum of projected density fluctuations in galaxy halos can be measured using strong gravitational lensing. We develop the mathematical framework of power spectrum estimation, and test our method on mock observations. We use our results to determine the types of observations required to measure the substructure power spectrum with high significance. We predict that deep observations ($\sim10$ hours on a single target) with current facilities can measure this power spectrum at the $3σ$ level, with no apparent degeneracy with unknown clumpiness in the background source structure or fluctuations from detector noise. Upcoming ALMA measurements of strong lenses are capable of placing strong constraints on the abundance of dark matter subhalos and the underlying particle nature of dark matter.

preprint2014arXiv

The Distribution of Dark Matter in the Milky Way's Disk

We present an analysis of the effects of dissipational baryonic physics on the local dark matter (DM) distribution at the location of the Sun, with an emphasis on the consequences for direct detection experiments. Our work is based on a comparative analysis of two cosmological simulations of a Milky Way halo: Eris, a full hydrodynamic simulation, and ErisDark, its DM-only realization. We find that two distinct processes lead in Eris to a 30% enhancement of DM in the disk plane at the location of the Sun: the accretion and disruption of satellites resulting in a DM component with net angular momentum and the contraction of baryons pulling DM into the disk plane without forcing it to co-rotate. Owing to its particularly quiescent merger history for dark halos of Milky Way mass, the co-rotating dark disk in Eris is less massive than what has been suggested by previous work, contributing only 9% of the local DM density. The speed distribution in Eris is broadened and shifted to higher speeds compared to its DM-only twin simulation ErisDark. At high speeds $f(v)$ falls more steeply in Eris than in ErisDark or the SHM, easing the tension between recent results from the CDMS-II and XENON100 experiments. The non-Maxwellian aspects of $f(v)$ are still present, but much less pronounced in Eris than in DM-only runs. The weak dark disk increases the time-averaged scattering rate by only a few percent at low recoil energies. On the high velocity tail, however, the increase in typical speeds due to baryonic contraction results in strongly enhanced mean scattering rates compared to ErisDark, although they are still suppressed compared to the SHM. Similar trends are seen regarding the amplitude of the annual modulation, while the modulated fraction is increased compared to the SHM and decreased compared to ErisDark.

preprint2013arXiv

A Baryonic Solution to the Missing Satellites Problem

It has been demonstrated that the inclusion of baryonic physics can alter the dark matter densities in the centers of low-mass galaxies, making the central dark matter slope more shallow than predicted in pure cold dark matter simulations. This flattening of the dark matter profile can occur in the most luminous subhalos around Milky Way-mass galaxies. Zolotov et al. (2012) have suggested a correction to be applied to the central masses of dark matter-only satellites in order to mimic the affect of (1) the flattening of the dark matter cusp due to supernova feedback in luminous satellites, and (2) enhanced tidal stripping due to the presence of a baryonic disk. In this paper, we apply this correction to the z=0 subhalo masses from the high resolution, dark matter-only Via Lactea II (VL2) simulation, and find that the number of massive subhalos is dramatically reduced. After adopting a stellar mass to halo mass relationship for the VL2 halos, and identifying subhalos that are (1) likely to be destroyed by stripping and (2) likely to have star formation suppressed by photo-heating, we find that the number of massive, luminous satellites around a Milky Way-mass galaxy is in agreement with the number of observed satellites around the Milky Way or M31. We conclude that baryonic processes have the potential to solve the missing satellites problem.

preprint2013arXiv

An Off-center Density Peak in the Milky Way's Dark Matter Halo?

We show that the position of the central dark matter density peak may be expected to differ from the dynamical center of the Galaxy by several hundred parsec. In Eris, a high resolution cosmological hydrodynamics simulation of a realistic Milky-Way-analog disk galaxy, this offset is 300 - 400 pc (~3 gravitational softening lengths) after z=1. In its dissipationless dark-matter-only twin simulation ErisDark, as well as in the Via Lactea II and GHalo simulations, the offset remains below one softening length for most of its evolution. The growth of the DM offset coincides with a flattening of the central DM density profile in Eris inwards of ~1 kpc, and the direction from the dynamical center to the point of maximum DM density is correlated with the orientation of the stellar bar, suggesting a bar-halo interaction as a possible explanation. A dark matter density offset of several hundred parsec greatly affects expectations of the dark matter annihilation signals from the Galactic Center. It may also support a dark matter annihilation interpretation of recent reports by Weniger (2012) and Su & Finkbeiner (2012) of highly significant 130 GeV gamma-ray line emission from a region 1.5 degrees (~200 parsec projected) away from Sgr A* in the Galactic plane.

preprint2013arXiv

Dwarf Galaxy Formation with H2-Regulated Star Formation: II. Gas-Rich Dark Galaxies at Redshift 2.5

We present a cosmological hydrodynamic simulation of the formation of dwarf galaxies at redshifts z>~2.5 using a physically-motivated model for H2-regulated star formation. Our simulation, performed using the Enzo code and reaching a peak resolution of 109 proper parsecs at z=2.5, extends the results of Kuhlen et al. (2012) to significantly lower redshifts. We show that a star formation prescription regulated by the local H2 abundance leads to the suppression of star formation in dwarf galaxy halos with M_h <~ 10^10 Msun and to a large population of gas-rich "dark galaxies" at z=2.5 with low star formation efficiencies and gas depletion timescales >20 Gyr. The fraction of dark galaxies is 60% at M_h ~ 10^10 Msun and increases rapidly with decreasing halo mass. Dark galaxies form late and their gaseous disks never reach the surface densities, > ~5700 Msun / pc^2 (Z/10^-3 Zsun)^(-0.88), that are required to build a substantial molecular fraction. Despite this large population of dark galaxies, we show that our H2-regulated simulation is consistent with both the observed luminosity function of galaxies and the cosmological mass density of neutral gas at z>~2.5. Moreover, our results provide a theoretical explanation for the recent detection in fluorescent Ly-alpha emission of gaseous systems at high redshift with little or no associated star formation. We further propose that H2-regulation may offer a fresh solution to a number of outstanding "dwarf galaxy problems" in LambdaCDM. In particular, H2-regulation leads galaxy formation to become effectively stochastic on mass scales of M_h ~ 10^10 Msun, and thus these massive dwarfs are not "too big to fail".

preprint2013arXiv

Enzo: An Adaptive Mesh Refinement Code for Astrophysics

This paper describes the open-source code Enzo, which uses block-structured adaptive mesh refinement to provide high spatial and temporal resolution for modeling astrophysical fluid flows. The code is Cartesian, can be run in 1, 2, and 3 dimensions, and supports a wide variety of physics including hydrodynamics, ideal and non-ideal magnetohydrodynamics, N-body dynamics (and, more broadly, self-gravity of fluids and particles), primordial gas chemistry, optically-thin radiative cooling of primordial and metal-enriched plasmas (as well as some optically-thick cooling models), radiation transport, cosmological expansion, and models for star formation and feedback in a cosmological context. In addition to explaining the algorithms implemented, we present solutions for a wide range of test problems, demonstrate the code's parallel performance, and discuss the Enzo collaboration's code development methodology.

preprint2012arXiv

Dark Matter Substructure Detection Using Spatially Resolved Spectroscopy of Lensed Dusty Galaxies

We investigate how strong lensing of dusty, star-forming galaxies by foreground galaxies can be used as a probe of dark matter halo substructure. We find that spatially resolved spectroscopy of lensed sources allows dramatic improvements to measurements of lens parameters. In particular we find that modeling of the full, three-dimensional (angular position and radial velocity) data can significantly facilitate substructure detection, increasing the sensitivity of observables to lower mass subhalos. We carry out simulations of lensed dusty sources observed by early ALMA (Cycle 1) and use a Fisher matrix analysis to study the parameter degeneracies and mass detection limits of this method. We find that, even with conservative assumptions, it is possible to detect galactic dark matter subhalos of ~ 10^8 M_{\odot} with high significance in most lensed DSFGs. Specifically, we find that in typical DSFG lenses, there is a ~ 55 % probability of detecting a substructure with M>10^8 M_{\odot} with more than 5 sigma detection significance in each lens, if the abundance of substructure is consistent with previous lensing results. The full ALMA array, with its significantly enhanced sensitivity and resolution, should improve these estimates considerably. Given the sample of ~100 lenses provided by surveys like the South Pole Telescope, our understanding of dark matter substructure in typical galaxy halos is poised to improve dramatically over the next few years.

preprint2012arXiv

Direct Detection of Dark Matter Debris Flows

Tidal stripping of dark matter from subhalos falling into the Milky Way produces narrow, cold tidal streams as well as more spatially extended "debris flows" in the form of shells, sheets, and plumes. Here we focus on the debris flow in the Via Lactea II simulation, and show that this incompletely phase-mixed material exhibits distinctive high velocity behavior. Unlike tidal streams, which may not necessarily intersect the Earth's location, debris flow is spatially uniform at 8 kpc and thus guaranteed to be present in the dark matter flux incident on direct detection experiments. At Earth-frame speeds greater than 450 km/s, debris flow comprises more than half of the dark matter at the Sun's location, and up to 80% at even higher speeds. Therefore, debris flow is most important for experiments that are particularly sensitive to the high speed tail of the dark matter distribution, such as searches for light or inelastic dark matter or experiments with directional sensitivity. We show that debris flow yields a distinctive recoil energy spectrum and a broadening of the distribution of incidence direction.

preprint2012arXiv

Identifying Local Group Field Galaxies which have interacted with the Milky Way

We distinguish between Local Group field galaxies which may have passed through the virial volume of the Milky Way, and those which have not, via a statistical compari- son against populations of dark matter haloes in the Via Lactea II (VLII) simulation with known orbital histories. Analysis of VLII provides expectations for this escaped population: they contribute 13 per cent of the galactic population between 300 and 1500 kpc from the Milky Way, and hence we anticipate that about 7 of the 54 known Local Group galaxies in that distance range are likely to be Milky Way escapees. These objects can be of any mass below that of the Milky Way, and they are expected to have positive radial velocities with respect to the Milky Way. Comparison of the radius-velocity distributions of VLII populations and measurements of Local Group galaxies presents a strong likelihood that Tucana, Cetus, NGC3109, SextansA, SextansB, Antlia, NGC6822, Phoenix, LeoT, and NGC185 have passed through the Milky Way. Most of these dwarfs have a lower HI mass fraction than the majority of dwarfs lying at similar distances to either the Milky Way or M31. Indeed, several of these galaxies - especially those with lower masses - contain signatures in their morphology, star formation history and/or gas content indicative of evolution seen in simulations of satellite/parent galactic interactions. Our results offer strong support for scenarios in which dwarfs of different types form a sequence in morphology and gas content, with evolution along the sequence being driven by interaction history.

preprint2012arXiv

Numerical Simulations of the Dark Universe: State of the Art and the Next Decade

We present a review of the current state of the art of cosmological dark matter simulations, with particular emphasis on the implications for dark matter detection efforts and studies of dark energy. This review is intended both for particle physicists, who may find the cosmological simulation literature opaque or confusing, and for astro-physicists, who may not be familiar with the role of simulations for observational and experimental probes of dark matter and dark energy. Our work is complementary to the contribution by M. Baldi in this issue, which focuses on the treatment of dark energy and cosmic acceleration in dedicated N-body simulations. Truly massive dark matter-only simulations are being conducted on national supercomputing centers, employing from several billion to over half a trillion particles to simulate the formation and evolution of cosmologically representative volumes (cosmic scale) or to zoom in on individual halos (cluster and galactic scale). These simulations cost millions of core-hours, require tens to hundreds of terabytes of memory, and use up to petabytes of disk storage. The field is quite internationally diverse, with top simulations having been run in China, France, Germany, Korea, Spain, and the USA. Predictions from such simulations touch on almost every aspect of dark matter and dark energy studies, and we give a comprehensive overview of this connection. We also discuss the limitations of the cold and collisionless DM-only approach, and describe in some detail efforts to include different particle physics as well as baryonic physics in cosmological galaxy formation simulations, including a discussion of recent results highlighting how the distribution of dark matter in halos may be altered. We end with an outlook for the next decade, presenting our view of how the field can be expected to progress. (abridged)

preprint2011arXiv

On the assembly of the Milky Way dwarf satellites and their common mass scale

We use a particle tagging technique to dynamically populate the N-body Via Lactea II high-resolution simulation with stars. The method is calibrated using the observed luminosity function of Milky Way satellites and the concentration of their stellar populations, and self-consistently follows the accretion and disruption of progenitor dwarfs and the build-up of the stellar halo in a cosmological "live host". Simple prescriptions for assigning stellar populations to collisionless particles are able to reproduce many properties of the observed Milky Way halo and its surviving dwarf satellites, like velocity dispersions, sizes, brightness profiles, metallicities, and spatial distribution. Our model predicts the existence of approximately 1,850 subhalos harboring "extremely faint" satellites (with mass-to-light ratios >5,000) lying beyond the Sloan Digital Sky Survey detection threshold. Of these, about 20 are "first galaxies", i.e. satellites that formed a stellar mass above 10 Msun before redshift 9. The ten most luminous satellites (L> 1e6 Lsun) in the simulation are hosted by subhalos with peak circular velocities today in the range V_max=10-40 km/s that have shed between 80% and 99% of their dark mass after being accreted at redshifts 1.7< z <4.6. The satellite maximum circular velocity and stellar line-of-sight velocity dispersion today follow the relation V_max=2.2 sigma_los. We apply a standard mass estimation algorithm based on Jeans modelling of the line-of-sight velocity dispersion profiles to the simulated dwarf spheroidals, and test the accuracy of this technique. The inner (within 300 pc) mass-luminosity relation for currently detectable satellites is nearly flat in our model, in qualitative agreement with the "common mass scale" found in Milky Way dwarfs. We do, however, predict a weak, but significant positive correlation for these objects: M_300 ~L^{0.088 \pm 0.024}.

preprint2010arXiv

Fermi-LAT Sensitivity to Dark Matter Annihilation in Via Lactea II Substructure

We present a study of the ability of the Fermi Gamma-ray Space Telescope to detect dark-matter annihilation signals from the Galactic subhalos predicted by the Via Lactea II N-body simulation. We implement an improved formalism for estimating the boost factor needed to account for the effect of dark-matter clumping on scales below the resolution of the simulation, and we incorporate a detailed Monte Carlo simulation of the response of the Fermi-LAT telescope, including a simulation of its all-sky observing mode integrated over a ten year mission. We find that for WIMP masses up to about 150 GeV in standard supersymmetric models with velocity-averaged cross section 3*10^-26 cm^3 s^-1, a few subhalos could be detectable with >5 standard deviations significance and would likely deviate significantly from the appearance of a point source.

preprint2010arXiv

Galactic Substructure and Dark Matter Annihilation in the Milky Way Halo

We study the effects of substructure on the rate of dark-matter annihilation in the Galactic halo. We use an analytic model for substructure that can extend numerical simulation results to scales too small to be resolved by the simulations. We first calibrate the analytic model to numerical simulations, and then determine the annihilation boost factor, for standard WIMP models as well as those with Sommerfeld (or other) enhancements, as a function of Galactocentric radius in the Milky Way. We provide an estimate of the dependence of the gamma-ray intensity of WIMP annihilation as a function of angular distance from the Galactic center. This methodology, coupled with future numerical simulation results can be a powerful tool that can be used to constrain WIMP properties using Fermi all-sky data.

preprint2010arXiv

The Origin of Dark Matter Halo Profiles

A longstanding puzzle of fundamental importance in modern cosmology has been the origin of the nearly universal density profiles of dark matter halos found in N-body simulations -- the so-called NFW profile. We show how this behavior may be understood, simply, by applying adiabatic contraction to peaks of Gaussian random fields. We argue that dynamical friction acts to reduce enormously the effect of random scatter in the properties of initial peaks, providing a key simplification. We compare our model predictions with results of the ultra-high resolution Via Lactea-II N-body simulation, and find superb agreement. We show how our model may be used to predict the distribution of halo properties like concentration. Our results suggest that many of the basic properties of halo structure may be understood using extremely simple physics.

preprint2009arXiv

The Dark Matter Annihilation Signal from Dwarf Galaxies and Subhalos

Dark Matter annihilation holds great potential for directly probing the clumpiness of the Galactic halo that is one of the key predictions of the Cold Dark Matter paradigm of hierarchical structure formation. Here we review the gamma-ray signal arising from dark matter annihilation in the centers of Galactic subhalos. We consider both known Galactic dwarf satellite galaxies and dark clumps without a stellar component as potential sources. Utilizing the Via Lactea II numerical simulation, we estimate fluxes for 18 Galactic dwarf spheroidals with published central densities. The most promising source is Segue 1, followed by Ursa Major II, Ursa Minor, Draco, and Carina. We show that if any of the known Galactic satellites can be detected, then at least ten times more subhalos should be visible, with a significant fraction of them being dark clumps.

preprint2008arXiv

Formation and evolution of galaxy dark matter halos and their substructure

We use the ``Via Lactea'' simulation to study the co-evolution of a Milky Way-size LambdaCDM halo and its subhalo population. While most of the host halo mass is accreted over the first 6 Gyr in a series of major mergers, the physical mass distribution [not M_vir(z)] remains practically constant since z=1. The same is true in a large sample of LambdaCDM galaxy halos. Subhalo mass loss peaks between the turnaround and virialization epochs of a given mass shell, and declines afterwards. 97% of the z=1 subhalos have a surviving bound remnant at the present epoch. The retained mass fraction is larger for initially lighter subhalos: satellites with maximum circular velocities Vmax=10 km/s at z=1 have today about 40% of their mass back then. At the first pericenter passage a larger average mass fraction is lost than during each following orbit. Tides remove mass in substructure from the outside in, leading to higher concentrations compared to field halos of the same mass. This effect, combined with the earlier formation epoch of the inner satellites, results in strongly increasing subhalo concentrations towards the Galactic center. We present individual evolutionary tracks and present-day properties of the likely hosts of the dwarf satellites around the Milky Way. The formation histories of ``field halos'' that lie today beyond the Via Lactea host are found to strongly depend on the density of their environment. This is caused by tidal mass loss that affects many field halos on eccentric orbits.

preprint2008arXiv

Quantifying the heart of darkness with GHALO - a multi-billion particle simulation of our galactic halo

We perform a series of simulations of a Galactic mass dark matter halo at different resolutions, our largest uses over three billion particles and has a mass resolution of 1000 M_sun. We quantify the structural properties of the inner dark matter distribution and study how they depend on numerical resolution. We can measure the density profile to a distance of 120 pc (0.05% of R_vir) where the logarithmic slope is -0.8 and -1.4 at (0.5% of R_vir). We propose a new two parameter fitting function that has a linearly varying logarithmic density gradient which fits the GHALO and VL2 density profiles extremely well. Convergence in the density profile and the halo shape scales as N^(-1/3), but the shape converges at a radius three times larger at which point the halo becomes more spherical due to numerical resolution. The six dimensional phase-space profile is dominated by the presence of the substructures and does not follow a power law, except in the smooth under-resolved inner few kpc.

preprint2007arXiv

Redefining the Missing Satellites Problem

Numerical simulations of Milky-Way size Cold Dark Matter (CDM) halos predict a steeply rising mass function of small dark matter subhalos and a substructure count that greatly outnumbers the observed satellites of the Milky Way. Several proposed explanations exist, but detailed comparison between theory and observation in terms of the maximum circular velocity (Vmax) of the subhalos is hampered by the fact that Vmax for satellite halos is poorly constrained. We present comprehensive mass models for the well-known Milky Way dwarf satellites, and derive likelihood functions to show that their masses within 0.6 kpc (M_0.6) are strongly constrained by the present data. We show that the M_0.6 mass function of luminous satellite halos is flat between ~ 10^7 and 10^8 M_\odot. We use the ``Via Lactea'' N-body simulation to show that the M_0.6 mass function of CDM subhalos is steeply rising over this range. We rule out the hypothesis that the 11 well-known satellites of the Milky Way are hosted by the 11 most massive subhalos. We show that models where the brightest satellites correspond to the earliest forming subhalos or the most massive accreted objects both reproduce the observed mass function. A similar analysis with the newly-discovered dwarf satellites will further test these scenarios and provide powerful constraints on the CDM small-scale power spectrum and warm dark matter models.

preprint2006arXiv

Dark matter substructure and gamma-ray annihilation in the Milky Way halo

We present initial results from ``Via Lactea'', the highest resolution simulation to date of Galactic CDM substructure. It follows the formation of a Milky Way-size halo with Mvir=1.8x10^12 Msun in a WMAP 3-year cosmology, using 234 million particles. Over 10,000 subhalos can be identified at z=0: Their cumulative mass function is well-fit by N(>Msub)= 0.0064 (Msub/Mvir)^(-1) down to Msun=4x10^6 Msun. The total mass fraction in subhalos is 5.3%, while the fraction of surface mass density in substructure within a projected distance of 10 kpc from the halo center is 0.3%. Because of the significant contribution from the smallest resolved subhalos, these fractions have not converged yet. Sub-substructure is apparent in all the larger satellites, and a few dark matter lumps are resolved even in the solar vicinity. The number of dark satellites with peak circular velocities above 10 km/s (5 km/s) is 124 (812): of these, 5 (26) are found within 0.1 Rvir, a region that appeared practically smooth in previous simulations. The neutralino self-annihilation gamma-ray emission from dark matter clumps is approximately constant per subhalo mass decade. Therefore, while in our run the contribution of substructure to the gamma-ray luminosity of the Galactic halo amounts to only 40% of the total spherically-averaged smooth signal, we expect this fraction to grow significantly as resolution is increased further. An all-sky map of the expected annihilation gamma-ray flux reaching a fiducial observer at 8 kpc from the Galactic center shows that at the current resolution a small number of subhalos start to be bright enough to be visible against the background from the smooth density field surrounding the observer.