Search NASA⌕ Search

Engineering topics

Kuhlen, Michael

Publications and source records attributed to Kuhlen, Michael.

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 (∼10 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.

79 ASTRONOMY AND ASTROPHYSICS↗

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 with identical initial conditions of a Milky Way halo, one of which (Eris) is a full hydrodynamic simulation and the other (ErisDark) is a DM-only one. We find that in Eris two distinct processes lead 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 the 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. Yet, since the simulation results in a realistic Milky Way analog galaxy, its DM halo provides a plausible alternative to the Maxwellian standard halo model (SHM) commonly used in direct detection analyses. 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 the 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.

79 ASTRONOMY AND ASTROPHYSICS↗

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 one, two, and three 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.

79 ASTRONOMY AND ASTROPHYSICS↗

DWARF GALAXY FORMATION WITH H{sub 2}-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 H{sub 2}-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. to significantly lower redshifts. We show that a star formation prescription regulated by the local H{sub 2} abundance leads to the suppression of star formation in dwarf galaxy halos with M{sub h} ∼< 10{sup 10} M{sub ☉} 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{sub h} ≅ 10{sup 10} M{sub ☉} and increases rapidly with decreasing halo mass. Dark galaxies form late and their gaseous disks never reach the surface densities, ∼> 5700 M{sub ☉} pc{sup –2} (Z/10{sup –3} Z{sub ☉}){sup –0.88}, that are required to build a substantial molecular fraction. Despite this large population of dark galaxies, we show that our H{sub 2}-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α emission of gaseous systems at high redshift with little or no associated star formation. We further propose that H{sub 2}-regulation may offer a fresh solution to a number of outstanding 'dwarf galaxy problems' in ΛCDM. In particular, H{sub 2}-regulation leads galaxy formation to become effectively stochastic on mass scales of M{sub h} ∼ 10{sup 10} M{sub ☉}, and thus these massive dwarfs are not 'too big to fail'.

79 ASTRONOMY AND ASTROPHYSICS↗

DARK MATTER SUBSTRUCTURE DETECTION USING SPATIALLY RESOLVED SPECTROSCOPY OF LENSED DUSTY GALAXIES

We investigate how strong lensing of dusty, star-forming galaxies (DSFGs) 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 {approx}10{sup 8} M{sub Sun} with high significance in most lensed DSFGs. Specifically, we find that in typical DSFG lenses, there is a {approx}55% probability of detecting a substructure with M > 10{sup 8} M{sub Sun} 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 {approx}100 lenses provided by surveys such as the South Pole Telescope, our understanding of dark matter substructure in typical galaxy halos is poised to improve dramatically over the next few years.

79 ASTRONOMY AND ASTROPHYSICS↗

AN OFF-CENTER DENSITY PEAK IN THE MILKY WAY'S DARK MATTER HALO?

We show that the position of the central dark matter (DM) density peak may be expected to differ from the dynamical center of the Galaxy by several hundred parsecs. In Eris, a high-resolution cosmological hydrodynamics simulation of a realistic Milky-Way-analog disk galaxy, this offset is 300-400 pc ({approx}3 gravitational softening lengths) after z = 1. In its dissipationless DM-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 inward of {approx}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 DM density offset of several hundred parsecs greatly affects expectations of the DM annihilation signals from the Galactic center. It may also support a DM annihilation interpretation of recent reports by Weniger and Su and Finkbeiner of highly significant 130 GeV gamma-ray line emission from a region 1. Degree-Sign 5 ({approx}200 pc projected) away from Sgr A* in the Galactic plane.

79 ASTRONOMY AND ASTROPHYSICS↗

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. 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.

79 ASTRONOMY AND ASTROPHYSICS↗