Source author record

Juerg Diemand

Juerg Diemand 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

19works
4topics
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

19 published item(s)

preprint2020arXiv

Globular Cluster Formation from Colliding Substructure

We investigate a scenario where the formation of Globular Clusters (GCs) is triggered by high-speed collisions between infalling atomic-cooling subhalos during the assembly of the main galaxy host, a special dynamical mode of star formation that operates at high gas pressures and is intimately tied to LCDM hierarchical galaxy assembly. The proposed mechanism would give origin to "naked" globulars, as colliding dark matter subhalos and their stars will simply pass through one another while the warm gas within them clashes at highly supersonic speed and decouples from the collisionless component, in a process reminiscent of the Bullet galaxy cluster. We find that the resulting shock-compressed layer cools on a timescale that is typically shorter than the crossing time, first by atomic line emission and then via fine-structure metal-line emission, and is subject to gravitational instability and fragmentation. Through a combination of kinetic theory approximation and high-resolution $N$-body simulations, we show that this model may produce: (a) a GC number-halo mass relation that is linear down to dwarf galaxy scales and agrees with the trend observed over five orders of magnitude in galaxy mass; (b) a population of old globulars with a median age of 12 Gyr and an age spread similar to that observed; (c) a spatial distribution that is biased relative to the overall mass profile of the host; and (d) a bimodal metallicity distribution with a spread similar to that observed in massive galaxies.

preprint2014arXiv

Milky Way Mass and Potential Recovery Using Tidal Streams in a Realistic Halo

We present a new method for determining the Galactic gravitational potential based on forward modeling of tidal stellar streams. We use this method to test the performance of smooth and static analytic potentials in representing realistic dark matter halos, which have substructure and are continually evolving by accretion. Our FAST-FORWARD method uses a Markov Chain Monte Carlo algorithm to compare, in 6D phase space, an "observed" stream to models created in trial analytic potentials. We analyze a large sample of streams evolved in the Via Lactea II (VL2) simulation, which represents a realistic Galactic halo potential. The recovered potential parameters are in agreement with the best fit to the global, present-day VL2 potential. However, merely assuming an analytic potential limits the dark matter halo mass measurement to an accuracy of 5 to 20%, depending on the choice of analytic parametrization. Collectively, mass estimates using streams from our sample reach this fundamental limit, but individually they can be highly biased. Individual streams can both under- and overestimate the mass, and the bias is progressively worse for those with smaller perigalacticons, motivating the search for tidal streams at galactocentric distances larger than 70 kpc. We estimate that the assumption of a static and smooth dark matter potential in modeling of the GD-1 and Pal5-like streams introduces an error of up to 50% in the Milky Way mass estimates.

preprint2014arXiv

Warm dark matter does not do better than cold dark matter in solving small-scale inconsistencies

Over the last decade, warm dark matter (WDM) has been repeatedly proposed as an alternative scenario to the standard cold dark matter (CDM) one, potentially resolving several disagreements between the CDM model and observations on small scales. Here, we reconsider the most important CDM small-scale discrepancies in the light of recent observational constraints on WDM. As a result, we find that a conventional thermal (or thermal-like) WDM cosmology with a particle mass in agreement with Lyman-$α$ is nearly indistinguishable from CDM on the relevant scales and therefore fails to alleviate any of the small-scale problems. The reason for this failure is that the power spectrum of conventional WDM falls off too rapidly. To maintain WDM as a significantly different alternative to CDM, more evolved production mechanisms leading to multiple dark matter components or a gradually decreasing small-scale power spectrum have to be considered.

preprint2013arXiv

Angular Momentum Acquisition in Galaxy Halos

We use high-resolution cosmological hydrodynamic simulations to study the angular momentum acquisition of gaseous halos around Milky Way sized galaxies. We find that cold mode accreted gas enters a galaxy halo with ~70% more specific angular momentum than dark matter averaged over cosmic time (though with a very large dispersion). In fact, we find that all matter has a higher spin parameter when measured at accretion than when averaged over the entire halo lifetime, and is well characterized by λ~0.1, at accretion. Combined with the fact that cold flow gas spends a relatively short time (1-2 dynamical times) in the halo before sinking to the center, this naturally explains why cold flow halo gas has a specific angular momentum much higher than that of the halo and often forms "cold flow disks". We demonstrate that the higher angular momentum of cold flow gas is related to the fact that it tends to be accreted along filaments.

preprint2013arXiv

Density Profiles of CDM Microhalos and their Implications for Annihilation Boost Factors

In a standard cold dark matter (CDM) cosmology, microhalos at the CDM cutoff scale are the first and smallest objects expected to form in the universe. Here we present results of high resolution simulations of three representative roughly Earth-mass microhalos in order to determine their inner density profile. We find that CDM microhalos in simulations without a cutoff in the power spectrum roughly follow the NFW density profile, just like the much larger CDM halos on galaxy and galaxy cluster scales. But having a cutoff in the initial power spectrum at a typical neutralino free streaming scale of $10^{-7} M_{\odot}$ makes their inner density profiles considerably steeper, i.e. $ρ\propto r^{-(1.3-1.4)}$, in good agreement with the results from Ishiyama et al. (2010). An extrapolation of the halo and subhalo mass functions down to the cutoff scale indicates that microhalos are extremely abundant throughout the present day dark matter distribution and might contribute significantly to indirect dark matter detection signals. Assuming a transition from a NFW to a steeper inner profile ($ρ\propto r^{-1.4}$) two orders of magnitude above the cutoff scale, the total boost factor for a Milky Way sized dark matter halo increases from about 3.5 to 4. We further find that CDM microhalo concentrations are consistent with the Bullock et al. (2001) model and clearly rule out simplistic power law models for the mass dependence of concentrations and subhalo annihilation, which would erroneously lead to very large boost factors (a few hundred for galaxy halos and over 1000 for clusters).

preprint2013arXiv

Structure Finding in Cosmological Simulations: The State of Affairs

The ever increasing size and complexity of data coming from simulations of cosmic structure formation demands equally sophisticated tools for their analysis. During the past decade, the art of object finding in these simulations has hence developed into an important discipline itself. A multitude of codes based upon a huge variety of methods and techniques have been spawned yet the question remained as to whether or not they will provide the same (physical) information about the structures of interest. Here we summarize and extent previous work of the "halo finder comparison project": we investigate in detail the (possible) origin of any deviations across finders. To this extent we decipher and discuss differences in halo finding methods, clearly separating them from the disparity in definitions of halo properties. We observe that different codes not only find different numbers of objects leading to a scatter of up to 20 per cent in the halo mass and Vmax function, but also that the particulars of those objects that are identified by all finders differ. The strength of the variation, however, depends on the property studied, e.g. the scatter in position, bulk velocity, mass, and the peak value of the rotation curve is practically below a few per cent, whereas derived quantities such as spin and shape show larger deviations. Our study indicates that the prime contribution to differences in halo properties across codes stems from the distinct particle collection methods and -- to a minor extent -- the particular aspects of how the procedure for removing unbound particles is implemented. We close with a discussion of the relevance and implications of the scatter across different codes for other fields such as semi-analytical galaxy formation models, gravitational lensing, and observables in general.

preprint2013arXiv

The Galactic Halo in Mixed Dark Matter Cosmologies

A possible solution to the small scale problems of the cold dark matter (CDM) scenario is that the dark matter consists of two components, a cold and a warm one. We perform a set of high resolution simulations of the Milky Way halo varying the mass of the WDM particle ($m_{\rm WDM}$) and the cosmic dark matter mass fraction in the WDM component ($\bar{f}_{\rm W}$). The scaling ansatz introduced in combined analysis of LHC and astroparticle searches postulates that the relative contribution of each dark matter component is the same locally as on average in the Universe (e.g. $f_{\rm W,\odot} = \bar{f}_{\rm W}$). Here we find however, that the normalised local WDM fraction ($f_{\rm W,\odot}$ / $\bar{f}_{\rm W}$) depends strongly on $m_{\rm WDM}$ for $m_{\rm WDM} <$ 1 keV. Using the scaling ansatz can therefore introduce significant errors into the interpretation of dark matter searches. To correct this issue a simple formula that fits the local dark matter densities of each component is provided.

preprint2012arXiv

Gamma-Ray Constraints on the First Stars from Annihilation of Light WIMPs

We calculate the limits on the fraction of viable dark matter minihalos in the early universe to host Population III.1 stars, surviving today as dark matter spikes in our Milky Way halo. Motivated by potential hints of light dark matter from the DAMA and CoGeNT direct dark matter searches, we consider thermal relic WIMP dark matter with masses of 5, 10, and 20 GeV, and annihilation to mu^+ mu^-, tau^+ tau^-, and q bar{q}. From this brief study we conclude that, if dark matter is light, either the typical black hole size is \lesssim 100 M_\odot (i.e. there is no significant Dark Star phase), and/or dark matter annihilates primarily to mu^+ mu^- or other final states that result in low gamma-ray luminosity, and/or that an extremely small fraction of minihalos in the early universe that seem suitable to host the formation of the first stars actually did.

preprint2011arXiv

Black Holes in our Galactic Halo: Compatibility with FGST and PAMELA Data and Constraints on the First Stars

10 to 10^5 solar mass black holes with dark matter spikes that formed in early minihalos and still exist in our Milky Way Galaxy today are examined in light of recent data from the Fermi Gamma-Ray Space Telescope (FGST). The dark matter spikes surrounding black holes in our Galaxy are sites of significant dark matter annihilation. We examine the signatures of annihilations into gamma-rays, electrons and positrons, and neutrinos. We find that some significant fraction of the point sources detected by FGST might be due to dark matter annihilation near black holes in our Galaxy. We obtain limits on the properties of dark matter annihilations in the spikes using the information in the FGST First Source Catalog as well as the diffuse gamma-ray flux measured by FGST. We determine the maximum fraction of high redshift minihalos that could have hosted the formation of the first generation of stars and, subsequently, their black hole remnants. The strength of the limits depends on the choice of annihilation channel and black hole mass; limits are strongest for the heaviest black holes and annhilation to $b \bar{b}$ and $W^+W^-$ final states. The larger black holes considered in this paper may arise as the remnants of Dark Stars after the dark matter fuel is exhausted and thermonuclear burning runs its course; thus FGST observations may be used to constrain the properties of Dark Stars. Additionally, we comment on the excess positron flux found by PAMELA and its possible interpretation in terms of dark matter annihilation around these black hole spikes.

preprint2011arXiv

Haloes gone MAD: The Halo-Finder Comparison Project

[abridged] We present a detailed comparison of fundamental dark matter halo properties retrieved by a substantial number of different halo finders. These codes span a wide range of techniques including friends-of-friends (FOF), spherical-overdensity (SO) and phase-space based algorithms. We further introduce a robust (and publicly available) suite of test scenarios that allows halo finder developers to compare the performance of their codes against those presented here. This set includes mock haloes containing various levels and distributions of substructure at a range of resolutions as well as a cosmological simulation of the large-scale structure of the universe. All the halo finding codes tested could successfully recover the spatial location of our mock haloes. They further returned lists of particles (potentially) belonging to the object that led to coinciding values for the maximum of the circular velocity profile and the radius where it is reached. All the finders based in configuration space struggled to recover substructure that was located close to the centre of the host halo and the radial dependence of the mass recovered varies from finder to finder. Those finders based in phase space could resolve central substructure although they found difficulties in accurately recovering its properties. Via a resolution study we found that most of the finders could not reliably recover substructure containing fewer than 30-40 particles. However, also here the phase space finders excelled by resolving substructure down to 10-20 particles. By comparing the halo finders using a high resolution cosmological volume we found that they agree remarkably well on fundamental properties of astrophysical significance (e.g. mass, position, velocity, and peak of the rotation curve).

preprint2011arXiv

The Total Mass of Dark Matter Haloes

The simple, conventional dark matter halo mass definitions commonly used in cosmological simulations ("virial" mass, FoF mass, $M_{50,100,200,...}$) only capture part of the collapsed material and are therefore inconsistent with the halo mass concept used in analytical treatments of structure formation. Simulations have demonstrated that typical dark matter particle orbits extend out to about 90 per cent of their turnaround radius, which results in apocenter passages outside of the current "virial" radius on the first and also on the second orbit. Here we describe how the formation history of haloes can be used to identify those particles which took part in the halo collapse, but are missed by conventional group-finders because of their remote present location. These particles are added to the part of the halo already identified by FoF. The corrected masses of dark haloes are significantly higher (the median mass increase is 25 per cent) and there is a considerable shift of the halo mass function towards the Press & Schechter form. We conclude that meaningful quantitative comparisons between (semi-)analytic predictions of halo properties (e.g. mass functions, mass accretion rates, merger rates, spatial clustering, etc.) and simulation results will require using the same halo definition in both approaches.

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

Signatures of Dark Star Remnants in the Galactic Halo

The very first stars likely formed from metal-free, molecular hydrogen-cooled gas at the centers of dark matter minihalos. Prior to nuclear fusion, these stars may have been supported by dark matter heating from annihilations in the star, in which case they could have grown to be quite massive before collapsing to black holes. Many remnant black holes and their surrounding dark matter density spikes may be part of our Milky Way halo today. Here we explore the gamma-ray signatures of dark matter annihilations in the dark matter spikes surrounding these black holes for a range of star formation scenarios, black hole masses, and dark matter annihilation modes. Data from the Fermi Gamma-Ray Space Telescope are used to constrain models of dark matter annihilation and the formation of the first stars.

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.

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.

preprint2006arXiv

Empirical models for Dark Matter Halos. I. Nonparametric Construction of Density Profiles and Comparison with Parametric Models

We use techniques from nonparametric function estimation theory to extract the density profiles, and their derivatives, from a set of N-body dark matter halos. We consider halos generated from LCDM simulations of gravitational clustering, as well as isolated, spherical collapses. The logarithmic density slopes gamma = d(log rho)/d(log r) of the LCDM halos are found to vary as power-laws in radius, reaching values of gamma ~ -1 at the innermost resolved radii (~0.01 r_virial). This behavior is significantly different from that of broken power-law models like the NFW profile, but similar to that of models like de Vaucouleurs'. Accordingly, we compare the N-body density profiles with various parametric models to find which provide the best fit. We consider an NFW-like model with arbitrary inner slope; Dehnen & McLaughlin's anisotropic model; Einasto's model (identical in functional form to Sersic's model but fit to the space density); and the density model of Prugniel & Simien that was designed to match the deprojected form of Sersic's R^{1/n} law. Overall, the best-fitting model to the LCDM halos is Einasto's, although the Prugniel-Simien and Dehnen-McLaughlin models also perform well. With regard to the spherical collapse halos, both the Prugniel-Simien and Einasto models describe the density profiles well, with an rms scatter some four times smaller than that obtained with either the NFW-like model or the 3-parameter Dehnen-McLaughlin model. Finally, we confirm recent claims of a systematic variation in profile shape with halo mass.

preprint2006arXiv

Empirical Models for Dark Matter Halos. II. Inner profile slopes, dynamical profiles, and rho/sigma^3

We have recently shown that both the Prugniel-Simien model and Sersic's function (hereafter referred to as the Einasto model when applied to internal density profiles) describe simulated dark matter halos better than an NFW-like model with equal number of parameters. Here we provide analytical expressions for the logarithmic slopes of these models, and compare them with data from real galaxies. Depending on the Einasto parameters of the dark matter halo, one can expect an extrapolated, inner (0.01--1 kpc), logarithmic profile slope ranging from -0.2 to -1.5, with a typical value at 0.1 kpc around -0.7. Application of this (better fitting) model therefore alleviates some of the past disagreement with observations on this issue. We additionally provide useful expressions for the concentration and assorted scale radii: r_s, r_{-2}, r_e, R_e, r_virial, and r_max -- the radius where the circular velocity profile has its maximum value. We also present the circular velocity profiles and the radial behavior of rho(r)/sigma(r)^3 for both the Einasto and Prugniel-Simien models. We find this representation of the phase-space density profile to be well approximated by a power-law with slope slightly shallower than -2 near r=r_{-2}.

preprint2006arXiv

Empirical Models for Dark Matter Halos. III. The Kormendy relation and the log(rho_e)-log(R_e) relation

We have recently shown that the 3-parameter density-profile model from Prugniel & Simien provides a better fit to simulated, galaxy- and cluster-sized, dark matter halos than an NFW-like model with arbitrary inner profile slope gamma (Paper I). By construction, the parameters of the Prugniel-Simien model equate to those of the Sersic R^{1/n} function fitted to the projected distribution. Using the Prugniel-Simien model, we are therefore able to show that the location of simulated (10^{12} M_sun) galaxy-sized dark matter halos in the <mu>_e-log(R_e) diagram coincides with that of brightest cluster galaxies, i.e., the dark matter halos appear consistent with the Kormendy relation defined by luminous elliptical galaxies. These objects are also seen to define the new, and equally strong, relation log(rho_e) = 0.5 - 2.5log(R_e), in which rho_e is the internal density at r=R_e. Simulated (10^{14.5} M_sun) cluster-sized dark matter halos and the gas component of real galaxy clusters follow the relation log(rho_e) = 2.5[1 - log(R_e)]. Given the shapes of the various density profiles, we are able to conclude that while dwarf elliptical galaxies and galaxy clusters can have dark matter halos with effective radii of comparable size to the effective radii of their baryonic component, luminous elliptical galaxies can not. For increasingly large elliptical galaxies, with increasingly large profile shapes `n', to be dark matter dominated at large radii requires dark matter halos with increasingly large effective radii compared to the effective radii of their stellar component.