Researcher profile

Andrey V. Kravtsov

Andrey V. Kravtsov contributes to research discovery and scholarly infrastructure.

ResearcherAffiliation not importedOpen to collaborate

astro-ph.CO astro-ph astro-ph.GA astro-ph.IM astro-ph.HE hep-ex hep-ph physics.comp-ph

Trust snapshot

Quick read

Trust 21 - EmergingVerification L1Unclaimed author

57works

0followers

8topics

4close collaborators

Actions

Decide how to stay connected

Follow researcher0

Claiming links this public author record to a researcher profile and unlocks direct collaboration workflows.

Direct collaboration

Open a focused conversation when the fit is right

Claim this author entity first to unlock direct invitations.

Inspect adjacent work, topics, institutions and collaborators without jumping out to a separate graph page.

Building this graph slice

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

preprint2022arXiv

Entropy-Conserving Scheme for Modeling Nonthermal Energies in Fluid Dynamics Simulations

We compare the performance of energy-based and entropy-conserving schemes for modeling nonthermal energy components, such as unresolved turbulence and cosmic rays, using idealized fluid dynamics tests and isolated galaxy simulations. While both methods are aimed to model advection and adiabatic compression or expansion of different energy components, the energy-based scheme numerically solves the nonconservative equation for the energy density evolution, while the entropy-conserving scheme uses a conservative equation for modified entropy. Using the standard shock tube and Zel'dovich pancake tests, we show that the energy-based scheme results in a spurious generation of nonthermal energy on shocks, while the entropy-conserving method evolves the energy adiabatically to machine precision. We also show that, in simulations of an isolated $L_\star$ galaxy, switching between the schemes results in $\approx 20-30\%$ changes of the total star formation rate and a significant difference in morphology, particularly near the galaxy center. We also outline and test a simple method that can be used in conjunction with the entropy-conserving scheme to model the injection of nonthermal energies on shocks. Finally, we discuss how the entropy-conserving scheme can be used to capture the kinetic energy dissipated by numerical viscosity into the subgrid turbulent energy implicitly, without explicit source terms that require calibration and can be rather uncertain. Our results indicate that the entropy-conserving scheme is the preferred choice for modeling nonthermal energy components, a conclusion that is equally relevant for Eulerian and moving-mesh fluid dynamics codes.

preprint2020arXiv

The Sheet of Giants: Unusual Properties of the Milky Way's Immediate Neighbourhood

We quantify the shape and overdensity of the galaxy distribution in the `Local Sheet' within a sphere of $R=8$ Mpc, and compare these properties with the expectations of the $Λ$CDM model. We measure ellipsoidal axis ratios of $c/a\approx0.16$ and $b/a\approx0.79$, indicating that the distribution of galaxies in the Local Volume can be approximated by a flattened oblate ellipsoid, consistent with the `sheet'-like configuration noted in previous studies. In contrast with previous estimates that the Local Sheet has a density close to average, we find that the number density of faint and bright galaxies in the Local Volume is $\approx1.7$ and $\approx5.2$ times denser, respectively, than the mean number density of galaxies of the same luminosity. Comparison with simulations shows that the number density contrasts of bright and faint galaxies within $8$ Mpc alone make the Local Volume a $\approx 2.5σ$ outlier in the $Λ$CDM cosmology. Our results indicate that the cosmic neighbourhood of the Milky Way may be unusual for galaxies of similar luminosity. The impact of the peculiar properties of our neighbourhood on the properties of the Milky Way and other nearby galaxies is not yet understood and warrants further study.

preprint2019arXiv

The Three Causes of Low-Mass Assembly Bias

We present a detailed analysis of the physical processes that cause halo assembly bias -- the dependence of halo clustering on proxies of halo formation time. We focus on the origin of assembly bias in the mass range corresponding to the hosts of typical galaxies and use halo concentration as our chief proxy of halo formation time. We also repeat our key analyses across a broad range of halo masses and for alternative formation time definitions. We show that splashback subhaloes are responsible for two thirds of the assembly bias signal, but do not account for the entire effect. After splashback subhaloes have been removed, we find that the remaining assembly bias signal is due to a relatively small fraction ($\lesssim 10\%$) of haloes in dense regions. We test a number of additional physical processes thought to contribute to assembly bias and demonstrate that the two key processes are the slowing of mass growth by large-scale tidal fields and by the high velocities of ambient matter in sheets and filaments. We also rule out several other proposed physical causes of halo assembly bias. Based on our results, we argue that there are three processes that contribute to assembly bias of low-mass halos: large-scale tidal fields, gravitational heating due to the collapse of large-scale structures, and splashback subhaloes located outside the virial radius.

preprint2015arXiv

A universal model for halo concentrations

We present a numerical study of dark matter halo concentrations in $Λ$CDM and self-similar cosmologies. We show that the relation between concentration, $c$, and peak height, $ν$, exhibits the smallest deviations from universality if halo masses are defined with respect to the critical density of the universe. These deviations can be explained by the residual dependence of concentration on the local slope of the matter power spectrum, $n$, which affects both the normalization and shape of the $c$-$ν$ relation. In particular, there is no well-defined floor in the concentration values. Instead, the minimum concentration depends on redshift: at fixed $ν$, halos at higher $z$ experience steeper slopes $n$, and thus have lower minimum concentrations. We show that the concentrations in our simulations can be accurately described by a universal seven-parameter function of only $ν$ and $n$. This model matches our $Λ$CDM results to $\lesssim 5\%$ accuracy up to $z = 6$, and matches scale-free $Ω_{\rm m} = 1$ models to $\lesssim 15\%$. The model also reproduces the low concentration values of Earth--mass halos at $z \approx 30$, and thus correctly extrapolates over $16$ orders of magnitude in halo mass. The predictions of our model differ significantly from all models previously proposed in the literature at high masses and redshifts. Our model is in excellent agreement with recent lensing measurements of cluster concentrations.

preprint2015arXiv

On the interplay between star formation and feedback in galaxy formation simulations

We investigate the star formation-feedback cycle in cosmological galaxy formation simulations, focusing on progenitors of Milky Way (MW)-sized galaxies. We find that in order to reproduce key properties of the MW progenitors, such as semi-empirically derived star formation histories and the shape of rotation curves, our implementation of star formation and stellar feedback requires 1) a combination of local early momentum feedback via radiation pressure and stellar winds and subsequent efficient supernovae feedback, and 2) efficacy of feedback that results in self-regulation of the global star formation rate on kiloparsec scales. We show that such feedback-driven self-regulation is achieved globally for a local star formation efficiency per free fall time of $ε_{\rm ff}\approx 10\%$. Although this value is larger that the $ε_{\rm ff}\sim 1\%$ value usually inferred from the Kennicutt-Schmidt (KS) relation, we show that it is consistent with direct observational estimates of $ε_{\rm ff}$ in molecular clouds. Moreover, we show that simulations with local efficiency of $ε_{\rm ff}\approx 10\%$ reproduce the global observed KS relation. Such simulations also reproduce the cosmic star formation history of the Milky Way sized galaxies and satisfy a number of other observational constraints. Conversely, we find that simulations that a priori assume an inefficient mode of star formation, instead of achieving it via stellar feedback regulation, fail to produce sufficiently vigorous outflows and do not reproduce observations. This illustrates the importance of understanding the complex interplay between star formation and feedback and the detailed processes that contribute to the feedback-regulated formation of galaxies.

preprint2014arXiv

Dependence of the outer density profiles of halos on their mass accretion rate

We present a systematic study of the density profiles of LCDM halos, focusing on the outer regions, 0.1 < r/Rvir < 9. We show that the median and mean profiles of halo samples of a given peak height exhibit significant deviations from the universal analytic profiles discussed previously in the literature, such as the Navarro-Frenk-White and Einasto profiles, at radii r > 0.5 R200m. In particular, at these radii the logarithmic slope of the median density profiles of massive or rapidly accreting halos steepens more sharply than predicted. The steepest slope of the profiles occurs at r ~ R200m, and its absolute value increases with increasing peak height or mass accretion rate, reaching slopes of -4 and steeper. Importantly, we find that the outermost density profiles at r > R200m are remarkably self-similar when radii are rescaled by R200m. This self-similarity indicates that radii defined with respect to the mean density are preferred for describing the structure and evolution of the outer profiles. However, the inner density profiles are most self-similar when radii are rescaled by R200c. We propose a new fitting formula that describes the median and mean profiles of halo samples selected by their peak height or mass accretion rate with accuracy < 10% at all radii, redshifts and masses we studied, r < 9 Rvir, 0 < z < 6 and Mvir > 1.7E10 Msun/h. We discuss observational signatures of the profile features described above, and show that the steepening of the outer profile should be detectable in future weak-lensing analyses of massive clusters. Such observations could be used to estimate the mass accretion rate of cluster halos.

preprint2014arXiv

On the mass of the Local Group

We use recent proper motion measurements of the tangential velocity of M31, along with its radial velocity and distance, to derive the likelihood of the sum of halo masses of the Milky Way and M31. This is done using a sample halo pairs in the Bolshoi cosmological simulation of $Λ$CDM cosmology selected to match properties and environment of the Local Group. The resulting likelihood gives estimate of the sum of masses of $M_{\rm MW,200}+M_{\rm M31,200}=$ $2.40_{-1.05}^{+1.95}\times10^{12}\,M_{\odot}$ ($90\%$ confidence interval). This estimate is consistent with individual mass estimates for the Milky Way and M31 and is consistent, albeit somewhat on the low side, with the mass estimated using the timing argument. We show that although the timing argument is unbiased on average for all pairs, for pairs constrained to have radial and tangential velocities similar to that of the Local Group the argument overestimates the sum of masses by a factor of $1.6$. Using similar technique we estimate the total dark matter mass enclosed within $1$ Mpc from the Local Group barycenter to be $M_{\rm LG}(r<1\, {\rm Mpc})=4.2_{-2.0}^{+3.4}\times10^{12}\,M_{\odot}$ ($90\%$ confidence interval).

preprint2013arXiv

Modelling mass distribution in elliptical galaxies: mass profiles and their correlation with velocity dispersion profiles

We assemble a statistical set of global mass models for ~2,000 nearly spherical SDSS galaxies at a mean redshift of 0.12 based on their aperture velocity dispersions and newly derived luminosity profiles in conjunction with published velocity dispersion profiles and empirical properties and relations of galaxy and halo parameters. When two-component (i.e. stellar plus dark) mass models are fitted to the SDSS aperture velocity dispersions, the predicted velocity dispersion profile (VP) slopes within the effective radius R_eff match well the distribution in observed elliptical galaxies. In contrast, the single-component models cannot reproduce the VP slope distribution. From a number of input variations the models exhibit for the radial range 0.1 R_eff < r < R_eff a tight correlation <gamma_e>=(1.865+/-0.008)+(-4.93+/-0.15)<eta> where <gamma_e> is the mean slope absolute value of the total mass density and <eta> is the mean slope of the velocity dispersion profile, which leads to a super-isothermal <gamma_e> = 2.15+/-0.04 for <eta>=-0.058+/-0.008 in observed elliptical galaxies. Furthermore, the successful two-component models appear to imply a typical slope curvature pattern in the total mass profile because for the observed steep luminosity (stellar mass) profile and the weak lensing inferred halo profile at large radii a total mass profile with monotonically varying slope would require too high DM density in the optical region giving rise to too large aperture velocity dispersion and too shallow VP.

preprint2013arXiv

On the evolution of cluster scaling relations

Understanding the evolution of scaling relations between the observable properties of clusters and their total mass is key to realizing their potential as cosmological probes. In this study, we investigate whether the evolution of cluster scaling relations is affected by the spurious evolution of mass caused by the evolving reference density with respect to which halo masses are defined (pseudo-evolution). We use the relation between mass, M, and velocity dispersion, sigma, as a test case, and show that the deviation from the M-sigma relation of cluster-sized halos caused by pseudo-evolution is smaller than 10% for a wide range of mass definitions. The reason for this small impact is a tight relation between the velocity dispersion and mass profiles, sigma(<r) = const * (GM(<r) / r)^1/2, which holds across a wide range of radii. We show that such a relation is generically expected for a variety of density profiles, as long as halos are in approximate Jeans equilibrium. Thus, as the outer "virial" radius used to define the halo mass, R, increases due to pseudo-evolution, halos approximately preserve their M-sigma relation. This result highlights the fact that tight scaling relations are the result of tight equilibrium relations between radial profiles of physical quantities. We find exceptions at very small and very large radii, where the profiles deviate from the relations they exhibit at intermediate radii. We discuss the implications of these results for other cluster scaling relations, and argue that pseudo-evolution should have a small effect on most scaling relations, except for those that involve the stellar masses of galaxies. In particular, we show that the relation between stellar-mass fraction and total mass is affected by pseudo-evolution and is largely shaped by it for halo masses of less than 10^14 solar masses.

preprint2013arXiv

Satellites in MW-like hosts: Environment dependence and close pairs

Previous studies showed that an estimate of the likelihood distribution of the Milky Way halo mass can be derived using the properties of the satellites similar to the Large and Small Magellanic Clouds (LMC and SMC). However, it would be straightforward to interpret such an estimate only if the properties of the Magellanic Clouds (MCs) are fairly typical and are not biased by the environment. In this study we explore whether the environment of the Milky Way affects the properties of the SMC and LMC such as their velocities. To test for the effect of the environment, we compare velocity distributions for MC-sized subhalos around Milky Way hosts in a sample selected simply by mass and in the second sample of such halos selected with additional restrictions on the distance to the nearest cluster and the local galaxy density, designed to mimic the environment of the Local Group (LG). We find that satellites in halos in the LG-like environments do have somewhat larger velocities, as compared to the halos of similar mass in the sample without environmental constraints. We derive the host halo likelihood distribution for the samples in the LG-like envirionment and in the control sample and find that the environment does not significantly affect the derived likelihood. We use the updated properties of the SMC and LMC to derive the constraint on the MW halo mass $\log{({\rm M}_{200} /\msol)}=12.06^{+0.31}_{-0.19}$ (90% confidence interval). We also explore the incidence of close pairs with relative velocities and separations similar to those of the LMC and SMC and find that such pairs are quite rare among $Λ$CDM halos. Taking into account the close separation of the MCs in the Busha et al.\ 2011 method results in the shift of the MW halo mass estimate to smaller masses, with the peak shifting approximately by a factor of two.[Abridged]

preprint2013arXiv

Simulations of disk galaxies with cosmic ray driven galactic winds

We present results from high-resolution hydrodynamic simulations of isolated SMC- and Milky Way-sized galaxies that include a model for feedback from galactic cosmic rays (CRs). We find that CRs are naturally able to drive winds with mass loading factors of up to ~10 in dwarf systems. The scaling of the mass loading factor with circular velocity between the two simulated systems is consistent with \propto v_c^{1-2} required to reproduce the faint end of the galaxy luminosity function. In addition, simulations with CR feedback reproduce both the normalization and the slope of the observed trend of wind velocity with galaxy circular velocity. We find that winds in simulations with CR feedback exhibit qualitatively different properties compared to SN driven winds, where most of the acceleration happens violently in situ near star forming sites. In contrast, the CR-driven winds are accelerated gently by the large-scale pressure gradient established by CRs diffusing from the star-forming galaxy disk out into the halo. The CR-driven winds also exhibit much cooler temperatures and, in the SMC-sized system, warm (T~10^4 K) gas dominates the outflow. The prevalence of warm gas in such outflows may provide a clue as to the origin of ubiquitous warm gas in the gaseous halos of galaxies detected via absorption lines in quasar spectra.

preprint2012arXiv

A High Throughput Workflow Environment for Cosmological Simulations

The next generation of wide-area sky surveys offer the power to place extremely precise constraints on cosmological parameters and to test the source of cosmic acceleration. These observational programs will employ multiple techniques based on a variety of statistical signatures of galaxies and large-scale structure. These techniques have sources of systematic error that need to be understood at the percent-level in order to fully leverage the power of next-generation catalogs. Simulations of large-scale structure provide the means to characterize these uncertainties. We are using XSEDE resources to produce multiple synthetic sky surveys of galaxies and large-scale structure in support of science analysis for the Dark Energy Survey. In order to scale up our production to the level of fifty 10^10-particle simulations, we are working to embed production control within the Apache Airavata workflow environment. We explain our methods and report how the workflow has reduced production time by 40% compared to manual management.

preprint2012arXiv

Constraining Cluster Physics with the Shape of X-ray Clusters: Comparison of Local X-ray Clusters versus LCDM Clusters

Simulations of cluster formation have demonstrated that condensation of baryons into central galaxies during cluster formation can drive the shape of the gas distribution in galaxy clusters significantly rounder, even at radii as large as half of the virial radius. However, such simulations generally predict stellar fractions within cluster virial radii that are ~2 to 3 times larger than the stellar masses deduced from observations. In this work we compare ellipticity profiles of clusters simulated with and without baryonic cooling to the cluster ellipticity profiles derived from Chandra and ROSAT observations in an effort to constrain the fraction of gas that cools and condenses into the central galaxies within clusters. We find that the observed ellipticity profiles are fairly constant with radius, with an average ellipticity of 0.18 +/- 0.05. The observed ellipticity profiles are in good agreement with the predictions of non-radiative simulations. On the other hand, the ellipticity profiles of the clusters in simulations that include radiative cooling, star formation, and supernova feedback (but no AGN feedback) deviate significantly from the observed ellipticity profiles at all radii. The simulations with cooling overpredict (underpredict) ellipticity in the inner (outer) regions of galaxy clusters. By comparing the simulations with and without cooling, we show that the cooling of gas via cooling flows in the central regions of simulated clusters causes the gas distribution to be more oblate in the central regions, but makes the outer gas distribution more spherical. We find that late-time gas cooling and star formation are responsible for the significantly oblate gas distributions in cluster cores, but the gas shapes outside of cluster cores are set primarily by baryon dissipation at high redshift z > 2.

preprint2012arXiv

On the Origin of the High Column Density Turnover in the HI Column Density Distribution

We study the high column density regime of the HI column density distribution function and argue that there are two distinct features: a turnover at NHI ~ 10^21 cm^-2 which is present at both z=0 and z ~ 3, and a lack of systems above NHI ~ 10^22 cm^-2 at z=0. Using observations of the column density distribution, we argue that the HI-H2 transition does not cause the turnover at NHI ~ 10^21 cm^-2, but can plausibly explain the turnover at NHI > 10^22 cm^-2. We compute the HI column density distribution of individual galaxies in the THINGS sample and show that the turnover column density depends only weakly on metallicity. Furthermore, we show that the column density distribution of galaxies, corrected for inclination, is insensitive to the resolution of the HI map or to averaging in radial shells. Our results indicate that the similarity of HI column density distributions at z=3 and z=0 is due to the similarity of the maximum HI surface densities of high-z and low-z disks, set presumably by universal processes that shape properties of the gaseous disks of galaxies. Using fully cosmological simulations, we explore other candidate physical mechanisms that could produce a turnover in the column density distribution. We show that while turbulence within GMCs cannot affect the DLA column density distribution, stellar feedback can affect it significantly if the feedback is sufficiently effective in removing gas from the central 2-3 kpc of high-redshift galaxies. Finally, we argue that it is meaningful to compare column densities averaged over ~ kpc scales with those estimated from quasar spectra which probe sub-pc scales due to the steep power spectrum of HI column density fluctuations observed in nearby galaxies.

preprint2012arXiv

Semi-empirical catalog of early-type galaxy-halo systems: dark matter density profiles, halo contraction and dark matter annihilation strength

With SDSS galaxy data and halo data from up-to-date N-body simulations we construct a semi-empirical catalog (SEC) of early-type systems by making a self-consistent bivariate statistical match of stellar mass (M_star) and velocity dispersion (sigma) with halo virial mass (M_vir). We then assign stellar mass profile and velocity dispersion profile parameters to each system in the SEC using their observed correlations with M_star and sigma. Simultaneously, we solve for dark matter density profile of each halo using the spherical Jeans equation. The resulting dark matter density profiles deviate in general from the dissipationless profile of NFW or Einasto and their mean inner density slope and concentration vary systematically with M_vir. Statistical tests of the distribution of profiles at fixed M_vir rule out the null hypothesis that it follows the distribution predicted by N-body simulations for M_vir ~< 10^{13.5-14.5} M_solar. These dark matter profiles imply that dark matter density is, on average, enhanced significantly in the inner region of halos with M_vir ~< 10^{13.5-14.5} M_solar supporting halo contraction. The main characteristics of halo contraction are: (1) the mean dark matter density within the effective radius has increased by a factor varying systematically up to ~ 3-4 at M_vir = 10^{12} M_solar, and (2) the inner density slope has a mean of <alpha> ~ 1.3 with rho(r) ~ r^{-alpha} and a halo-to-halo rms scatter of rms(alpha) ~ 0.4-0.5 for 10^{12} M_solar ~< M_vir ~< 10^{13-14} M_solar steeper than the NFW profile (alpha=1). Based on our results we predict that halos of nearby elliptical and lenticular galaxies can, in principle, be promising targets for gamma-ray emission from dark matter annihilation.

preprint2012arXiv

The X-factor in Galaxies: II. The molecular hydrogen -- star formation relation

There is ample observational evidence that the star formation rate (SFR) surface density, Sigma_SFR, is closely correlated with the surface density of molecular hydrogen, Sigma_H2. This empirical relation holds both for galaxy-wide averages and for individual >=kpc sized patches of the interstellar medium (ISM), but appears to degrade substantially at a sub-kpc scale. Identifying the physical mechanisms that determine the scale-dependent properties of the observed Sigma_H2-Sigma_SFR relation remains a challenge from a theoretical perspective. To address this question, we analyze the slope and scatter of the Sigma_H2-Sigma_SFR relation using a set of cosmological, galaxy formation simulations with a peak resolution of ~100 pc. These simulations include a chemical network for molecular hydrogen, a model for the CO emission, and a simple, stochastic prescription for star formation that operates on ~100 pc scales. Specifically, star formation is modeled as a Poisson process in which the average SFR is directly proportional to the present mass of H2. The predictions of our numerical model are in good agreement with the observed Kennicutt-Schmidt and Sigma_H2-Sigma_SFR relations. We show that observations based on CO emission are ill suited to reliably measure the slope of the latter relation at low (<20 M_sun pc^-2) H2 surface densities on sub-kpc scales. Our models also predict that the inferred Sigma_H2-Sigma_SFR relation steepens at high H2 surface densities as a result of the surface density dependence of the CO/H2 conversion factor. Finally, we show that on sub-kpc scales most of the scatter in the relation is a consequence of discreteness effects in the star formation process. In contrast, variations of the CO/H2 conversion factor are responsible for most of the scatter measured on super-kpc scales.

preprint2012arXiv

Towards a complete accounting of energy and momentum from stellar feedback in galaxy formation simulations

Stellar feedback plays a key role in galaxy formation by regulating star formation, driving interstellar turbulence and generating galactic scale outflows. Although modern simulations of galaxy formation can resolve scales of 10-100 pc, star formation and feedback operate on smaller, "subgrid" scales. Great care should therefore be taken in order to properly account for the effect of feedback on global galaxy evolution. We investigate the momentum and energy budget of feedback during different stages of stellar evolution, and study its impact on the interstellar medium using simulations of local star forming regions and galactic disks at the resolution affordable in modern cosmological zoom-in simulations. In particular, we present a novel subgrid model for the momentum injection due to radiation pressure and stellar winds from massive stars during early, pre-supernova evolutionary stages of young star clusters. Early injection of momentum acts to clear out dense gas in star forming regions, hence limiting star formation. The reduced gas density mitigates radiative losses of thermal feedback energy from subsequent supernova explosions, leading to an increased overall efficiency of stellar feedback. The detailed impact of stellar feedback depends sensitively on the implementation and choice of parameters. Somewhat encouragingly, we find that implementations in which feedback is efficient lead to approximate self-regulation of global star formation efficiency. We compare simulation results using our feedback implementation to other phenomenological feedback methods, where thermal feedback energy is allowed to dissipate over time scales longer than the formal gas cooling time. We find that simulations with maximal momentum injection suppress star formation to a similar degree as is found in simulations adopting adiabatic thermal feedback.

preprint2011arXiv

Fuel Efficient Galaxies: Sustaining Star Formation with Stellar Mass Loss

We examine the importance of secular stellar mass loss for fueling ongoing star formation in disk galaxies during the late stages of their evolution. For a galaxy of a given stellar mass, we calculate the total mass loss rate of its entire stellar population using star formation histories derived from the observed evolution of the M*-star formation rate relation, along with the predictions of standard stellar evolution models for stellar mass loss for a variety of initial stellar mass functions. Our model shows that recycled gas from stellar mass loss can provide most or all of the fuel required to sustain the current level of star formation in late type galaxies. Stellar mass loss can therefore remove the tension between the low gas infall rates that are derived from observations and the relatively rapid star formation occurring in disk galaxies. For galaxies where cold gas infall rates have been estimated, we demonstrate explicitly that stellar mass loss can account for most of the deficit between their star formation and infall rates.

preprint2011arXiv

Halo Contraction Effect in Hydrodynamic Simulations of Galaxy Formation

The condensation of gas and stars in the inner regions of dark matter halos leads to a more concentrated dark matter distribution. While this effect is based on simple gravitational physics, the question of its validity in hierarchical galaxy formation has led to an active debate in the literature. We use a collection of several state-of-the-art cosmological hydrodynamic simulations to study the halo contraction effect in systems ranging from dwarf galaxies to clusters of galaxies, at high and low redshift. The simulations are run by different groups with different codes and include hierarchical merging, gas cooling, star formation, and stellar feedback. We show that in all our cases the inner dark matter density increases relative to the matching simulation without baryon dissipation, at least by a factor of several. The strength of the contraction effect varies from system to system and cannot be reduced to a simple prescription. We present a revised analytical model that describes the contracted mass profile to an rms accuracy of about 10%. The model can be used to effectively bracket the response of the dark matter halo to baryon dissipation. The halo contraction effect is real and must be included in modeling of the mass distribution of galaxies and galaxy clusters.

preprint2011arXiv

Implementing the DC Mode in Cosmological Simulations with Supercomoving Variables

As emphasized by previous studies, proper treatment of the density fluctuation on the fundamental scale of a cosmological simulation volume - the "DC mode" - is critical for accurate modeling of spatial correlations on scales > 10% of simulation box size. We provide further illustration of the effects of the DC mode on the abundance of halos in small boxes and show that it is straightforward to incorporate this mode in cosmological codes that use the "supercomoving" variables. The equations governing evolution of dark matter and baryons recast with these variables are particularly simple and include the expansion factor, and hence the effect of the DC mode, explicitly only in the Poisson equation.

preprint2011arXiv

On determining the shape of matter distributions

A basic property of objects, like galaxies and halos that form in cosmological structure formation simulations, is their shape. Here, we critically investigate shape determination methods that are commonly used in the literature. It is found that using an enclosed integration volume and weight factors r^{-2} and r_{ell}^{-2} (elliptical radius) for the contribution of each particle or volume element in the shape tensor leads to biased axis ratios and smoothing of details when calculating the local shape as a function of distance from the center. To determine the local shape of matter distributions as a function of distance for well resolved objects (typically more than O(10^4) particles), we advocate a method that (1) uses an ellipsoidal shell (homoeoid) as an integration volume without any weight factors in the shape tensor and (2) removes subhalos.

preprint2011arXiv

On the Accuracy of Weak Lensing Cluster Mass Reconstructions

We study the bias and scatter in mass measurements of galaxy clusters resulting from fitting a spherically-symmetric Navarro, Frenk & White model to the reduced tangential shear profile measured in weak lensing observations. The reduced shear profiles are generated for ~10^4 cluster-sized halos formed in a LCDM cosmological N-body simulation of a 1 Gpc/h box. In agreement with previous studies, we find that the scatter in the weak lensing masses derived using this fitting method has irreducible contributions from the triaxial shapes of cluster-sized halos and uncorrelated large-scale matter projections along the line-of-sight. Additionally, we find that correlated large-scale structure within several virial radii of clusters contributes a smaller, but nevertheless significant, amount to the scatter. The intrinsic scatter due to these physical sources is ~20% for massive clusters, and can be as high as ~30% for group-sized systems. For current, ground-based observations, however, the total scatter should be dominated by shape noise from the background galaxies used to measure the shear. Importantly, we find that weak lensing mass measurements can have a small, ~5%-10%, but non-negligible amount of bias. Given that weak lensing measurements of cluster masses are a powerful way to calibrate cluster mass-observable relations for precision cosmological constraints, we strongly emphasize that a robust calibration of the bias requires detailed simulations which include more observational effects than we consider here. Such a calibration exercise needs to be carried out for each specific weak lensing mass estimation method, as the details of the method determine in part the expected scatter and bias. We present an iterative method for estimating mass M500c that can eliminate the bias for analyses of ground-based data.

preprint2011arXiv

Shapes of Gas, Gravitational Potential and Dark Matter in Lambda-CDM Clusters

We present analysis of the three-dimensional shape of intracluster gas in clusters formed in cosmological simulations of the Lambda-CDM cosmology and compare it to the shape of dark matter distribution and the shape of the overall isopotential surfaces. We find that in simulations with radiative cooling, star formation and stellar feedback (CSF), intracluster gas outside the cluster core is more spherical compared to non-radiative (NR) simulations, while in the core the gas in the CSF runs is more triaxial and has a distinctly oblate shape. The latter reflects the ongoing cooling of gas, which settles into a thick oblate ellipsoid as it loses thermal energy. The shape of the gas in the inner regions of clusters can therefore be a useful diagnostic of gas cooling. We find that gas traces the shape of the underlying potential rather well outside the core, as expected in hydrostatic equilibrium. At smaller radii, however, the gas and potential shapes differ significantly. In the CSF runs, the difference reflects the fact that gas is partly rotationally supported. Interestingly, we find that in NR simulations the difference between gas and potential shape at small radii is due to random gas motions, which make the gas distribution more spherical than the equipotential surfaces. Finally, we use mock Chandra X-ray maps to show that the differences in shapes observed in three-dimensional distribution of gas are discernible in the ellipticity of X-ray isophotes. Contrasting the ellipticities measured in simulated clusters against observations can therefore constrain the amount of cooling of the intracluster medium and the presence of random gas motions in cluster cores.

preprint2011arXiv

The impact of baryon physics on the structure of high-redshift galaxies

We study the detailed structure of galaxies at redshifts z > 2 using cosmological simulations with improved modeling of the interstellar medium and star formation. The simulations follow the formation and dissociation of molecular hydrogen, and include star formation only in cold molecular gas. The molecular gas is more concentrated towards the center of galaxies than the atomic gas, and as a consequence, the resulting stellar distribution is very compact. For halos with total mass above 10^{11} Mo, the median half-mass radius of the stellar disks is 0.8 kpc at z = 3. The vertical structure of the molecular disk is much thinner than that of the atomic neutral gas. Relative to the non-radiative run, the inner regions of the dark matter halo change shape from prolate to mildly oblate and align with the stellar disk. However, we do not find evidence for a significant dark disk of dark matter around the stellar disk. The outer halo regions retain the orientation acquired during accretion and mergers, and are significantly misaligned with the inner regions. The radial profile of the dark matter halo contracts in response to baryon dissipation, establishing an approximately isothermal profile throughout most of the halo. This effect can be accurately described by a modified model of halo contraction. The angular momentum of a fixed amount of inner dark matter is approximately conserved over time, while in the dissipationless case most of it is transferred outward during mergers. The conservation of the dark matter angular momentum provides supporting evidence for the validity of the halo contraction model in a hierarchical galaxy formation process.

preprint2011arXiv

The X-factor in Galaxies: I. Dependence on Environment and Scale

Characterizing the conversion factor between CO emission and column density of molecular hydrogen, X_CO, is crucial in studying the gaseous content of galaxies, its evolution, and relation to star formation. In most cases the conversion factor is assumed to be close to that of giant molecular clouds (GMCs) in the Milky Way, except possibly for mergers and star-bursting galaxies. However, there are physical grounds to expect that it should also depend on the gas metallicity, surface density, and strength of the interstellar radiation field. The XCO factor may also depend on the scale on which CO emission is averaged due to effects of limited resolution. We study the dependence of X_CO on gas properties and averaging scale using a model that is based on a combination of results of sub-pc scale magneto-hydrodynamic simulations and on the gas distribution from self-consistent cosmological simulations of galaxy formation. Our model predicts a value of X_CO that is consistent with the Galactic value for interstellar medium conditions typical for the Milky Way. For such conditions the predicted X_CO varies by only a factor of two for gas surfaced densities in the range \sim 50 - 500 M_sun / pc^2. However, the model also predicts that more generally on the scale of GMCs, X_CO is a strong function of metallicity, and depends on the column density and the interstellar UV flux. We show explicitly that neglecting these dependencies in observational estimates can strongly bias the inferred distribution of H2 column densities of molecular clouds to have a narrower and offset range compared to the true distribution. We find that when averaged on \sim kpc scales the X-factor depends only weakly on radiation field and column density, but is still a strong function of metallicity. The predicted metallicity dependence can be approximated as X_CO \sim Z^{-γ} with γ ~ 0.5 - 0.8.

preprint2011arXiv

Ultra-faint dwarf galaxies as a test of early enrichment and metallicity-dependent star formation

The tight relation of star formation with molecular gas indicated by observations and assumed in recent models implies that the efficiency with which galaxies convert their gas into stars depends on gas metallicity. This is because the abundance of molecular hydrogen is sensitive to the abundance of dust, which catalyzes the formation of H_2 and helps to shield it from dissociating radiation. In this study we point out that in the absence of significant pre-enrichment by Population III stars forming out of zero metallicity gas, such H_2-based star formation is expected to leave an imprint in the form of bi-modality in the metallicity distribution among dwarf galaxies and in the metallicity distribution of stars within individual galaxies. The bi-modality arises because when gas metallicity (and dust abundance) is low, formation of molecular gas is inefficient, the gas consumption time scale is long, and star formation and metal enrichment proceed slowly. When metallicity reaches a critical threshold value star formation and enrichment accelerate, which leads to rapid increase in both stellar mass and metallicity of galaxies. We demonstrate this process both using a simple analytical model and full cosmological simulations. In contrast, observed metallicity distributions of dwarf galaxies or stars within them are not bi-modal. We argue that this discrepancy points to substantial early stochastic pre-enrichment by population III stars to levels Z ~ 0.01 Z_sun in dense, star forming regions of early galaxies.

preprint2010arXiv

Environmental Dependence of the Kennicutt-Schmidt Relation in Galaxies

We present a detailed description of a phenomenological H2 formation model and local star formation prescription based on the density of molecular (rather than total) gas. Such approach allows us to avoid the arbitrary density and temperature thresholds typically used in star formation recipes. We present results of the model based on realistic cosmological simulations of high-z galaxy formation for a grid of numerical models with varied dust-to-gas ratios and interstellar far UV (FUV) fluxes. Our results show that both the atomic-to-molecular transition on small, ~10 pc scales and the Kennicutt-Schmidt (KS) relation on ~kpc scales are sensititive to the dust-to-gas ratio and the FUV flux. The atomic-to-molecular transition as a function of gas density or column density has a large scatter but is rather sharp and shifts to higher densities with decreasing dust-to-gas ratio and/or increasing FUV flux. Consequently, star formation is concentrated to higher gas surface density regions, resulting in steeper slope and lower amplitude of the KS relation at a given gas surface density, in less dusty and/or higher FUV flux environments. These trends should have a particularly strong effect on the evolution of low-mass, low surface brightness galaxies which typically have low dust content and anemic star formation, but are also likely to be important for evolution of the Milky Way-sized systems. We parameterize the dependencies observed in our simulations in convenient fitting formulae, which can be used to model the dependence of the KS relation on the dust-to-gas ratio and FUV flux in semi-analytic models and in cosmological simulations that do not include radiative transfer and H2 formation.

preprint2010arXiv

The distribution of annihilation luminosities in dark matter substructure

We calculate the probability distribution function (PDF) of the expected annihilation luminosities of dark matter subhalos as a function of subhalo mass and distance from the Galactic center using a semi-analytical model of halo evolution. We find that the PDF of luminosities is relatively broad, exhibiting a spread of as much as an order of magnitude at fixed subhalo mass and halo-centric distance. The luminosity PDF allows for simple construction of mock samples of gamma-ray luminous subhalos and assessment of the variance in among predicted gamma-ray signals from dark matter annihilation. Other applications include quantifying the variance among the expected luminosities of dwarf spheroidal galaxies, assessing the level at which dark matter annihilation can be a contaminant in the expected gamma-ray signal from other astrophysical sources, as well as estimating the level at which nearby subhalos can contribute to the antimatter flux.

preprint2010arXiv

The Large Scale Bias of Dark Matter Halos: Numerical Calibration and Model Tests

We measure the clustering of dark matter halos in a large set of collisionless cosmological simulations of the flat LCDM cosmology. Halos are identified using the spherical overdensity algorithm, which finds the mass around isolated peaks in the density field such that the mean density is Delta times the background. We calibrate fitting functions for the large scale bias that are adaptable to any value of Delta we examine. We find a ~6% scatter about our best fit bias relation. Our fitting functions couple to the halo mass functions of Tinker et. al. (2008) such that bias of all dark matter is normalized to unity. We demonstrate that the bias of massive, rare halos is higher than that predicted in the modified ellipsoidal collapse model of Sheth, Mo, & Tormen (2001), and approaches the predictions of the spherical collapse model for the rarest halos. Halo bias results based on friends-of-friends halos identified with linking length 0.2 are systematically lower than for halos with the canonical Delta=200 overdensity by ~10%. In contrast to our previous results on the mass function, we find that the universal bias function evolves very weakly with redshift, if at all. We use our numerical results, both for the mass function and the bias relation, to test the peak-background split model for halo bias. We find that the peak-background split achieves a reasonable agreement with the numerical results, but ~20% residuals remain, both at high and low masses.

preprint2010arXiv

The Metal-Enriched Outer Disk of NGC 2915

We present optical emission-line spectra for outlying HII regions in the extended neutral gas disk surrounding the blue compact dwarf galaxy NGC 2915. Using a combination of strong-line R23 and direct oxygen abundance measurements, we report a flat, possibly increasing, metallicity gradient out to 1.2 times the Holmberg radius. We find the outer-disk of NGC 2915 to be enriched to a metallicity of 0.4 Z_solar. An analysis of the metal yields shows that the outer disk of NGC 2915 is overabundant for its gas fraction, while the central star-foming core is similarly under-abundant for its gas fraction. Star formation rates derived from very deep ~14 ks GALEX FUV exposures indicate that the low-level of star formation observed at large radii is not sufficient to have produced the measured oxygen abundances at these galactocentric distances. We consider 3 plausible mechanisms that may explain the metal-enriched outer gaseous disk of NGC 2915: radial redistribution of centrally generated metals, strong galactic winds with subsequent fallback, and galaxy accretion. Our results have implications for the physical origin of the mass-metallicity relation for gas-rich dwarf galaxies.

preprint2009arXiv

Computational Eulerian Hydrodynamics and Galilean Invariance

Eulerian hydrodynamical simulations are a powerful and popular tool for modeling fluids in astrophysical systems. In this work, we critically examine recent claims that these methods violate Galilean invariance of the Euler equations. We demonstrate that Eulerian hydrodynamics methods do converge to a Galilean-invariant solution, provided a well-defined convergent solution exists. Specifically, we show that numerical diffusion, resulting from diffusion-like terms in the discretized hydrodynamical equations solved by Eulerian methods, accounts for the effects previously identified as evidence for the Galilean non-invariance of these methods. These velocity-dependent diffusive terms lead to different results for different bulk velocities when the spatial resolution of the simulation is kept fixed, but their effect becomes negligible as the resolution of the simulation is increased to obtain a converged solution. In particular, we find that Kelvin-Helmholtz instabilities develop properly in realistic Eulerian calculations regardless of the bulk velocity provided the problem is simulated with sufficient resolution (a factor of 2-4 increase compared to the case without bulk flows for realistic velocities). Our results reiterate that high-resolution Eulerian methods can perform well and obtain a convergent solution, even in the presence of highly supersonic bulk flows.

preprint2009arXiv

Cosmological Radiative Transfer Comparison Project II: The Radiation-Hydrodynamic Tests

The development of radiation hydrodynamical methods that are able to follow gas dynamics and radiative transfer self-consistently is key to the solution of many problems in numerical astrophysics. Such fluid flows are highly complex, rarely allowing even for approximate analytical solutions against which numerical codes can be tested. An alternative validation procedure is to compare different methods against each other on common problems, in order to assess the robustness of the results and establish a range of validity for the methods. Previously, we presented such a comparison for a set of pure radiative transfer tests (i.e. for fixed, non-evolving density fields). This is the second paper of the Cosmological Radiative Transfer (RT) Comparison Project, in which we compare 9 independent RT codes directly coupled to gasdynamics on 3 relatively simple astrophysical hydrodynamics problems: (5) the expansion of an H II region in a uniform medium; (6) an ionization front (I-front) in a 1/r^2 density profile with a flat core, and (7), the photoevaporation of a uniform dense clump. Results show a broad agreement between the different methods and no big failures, indicating that the participating codes have reached a certain level of maturity and reliability. However, many details still do differ, and virtually every code has showed some shortcomings and has disagreed, in one respect or another, with the majority of the results. This underscores the fact that no method is universal and all require careful testing of the particular features which are most relevant to the specific problem at hand.

preprint2009arXiv

Dark matter substructure and dwarf galactic satellites

A decade ago cosmological simulations of increasingly higher resolution were used to demonstrate that virialized regions of Cold Dark Matter (CDM) halos are filled with a multitude of dense, gravitationally-bound clumps. These dark matter subhalos are central regions of halos that survived strong gravitational tidal forces and dynamical friction during the hierarchical sequence of merging and accretion via which the CDM halos form. Comparisons with observations revealed that there is a glaring discrepancy between abundance of subhalos and luminous satellites of the Milky Way and Andromeda as a function of their circular velocity or bound mass within a fixed aperture. This large discrepancy, which became known as the ``substructure'' or the ``missing satellites'' problem, begs for an explanation. In this paper I review the progress made during the last several years both in quantifying the problem and in exploring possible scenarios in which it could be accommodated and explained in the context of galaxy formation in the framework of the CDM paradigm of structure formation. In particular, I show that the observed luminosity function, radial distribution, and the remarkable similarity of the inner density profiles of luminous satellites can be understood within hierarchical CDM framework using a simple model in which efficiency of star formation monotonically decreases with decreasing virial mass satellites had before their accretion without any actual sharp galaxy formation threshold.

preprint2009arXiv

Effects of Baryon Dissipation on the Dark Matter Virial Scaling Relation

We investigate effects of baryon dissipation on the dark matter virial scaling relation between total mass and velocity dispersion and the velocity bias of galaxies in groups and clusters using self-consistent cosmological simulations. We show that the baryon dissipation increases the velocity dispersion of dark matter within the virial radius by 5% - 10%. The effect is mainly driven by the change in density and gravitational potential in inner regions of cluster, and it is larger in lower mass systems where gas cooling and star formation are more efficient. We also show that the galaxy velocity bias depends on how galaxies are selected. Galaxies selected based on their stellar mass exhibit no velocity bias, while galaxies selected based on their total mass show positive bias of ~10%, consistent with previous results based on collisionless dark matter- only simulations. We further find that observational estimates of galaxy velocity dispersion are unbiased with respect to the velocity dispersion of dark matter, provided galaxies are selected using their stellar masses and and their velocity dispersions are computed with more than twenty most massive galaxies. Velocity dispersions estimated with fewer galaxies, on the other hand, can lead to significant underestimate of dynamical masses. Results presented in this paper should be useful in interpretating high-redshift groups and clusters as well as cosmological constraints derived from upcoming optical cluster surveys.

preprint2009arXiv

Evolution of Dark Matter Phase-Space Density Distributions in Equal-Mass Halo Mergers

We use dissipationless N-body simulations to investigate the evolution of the true coarse-grained phase-space density distribution f(x,v) in equal-mass mergers between dark matter (DM) halos. The halo models are constructed with various asymptotic power-law indices ranging from steep cusps to core-like profiles and we employ the phase-space density estimator ``Enbid'' developed by Sharma & Steinmetz to compute f(x,v). The adopted force resolution allows robust phase-space density profile estimates in the inner ~1% of the virial radii of the simulated systems. We confirm that mergers result in a decrease of the coarse-grained phase-space density in accordance with expectations from Mixing Theorems for collisionless systems. We demonstrate that binary mergers between identical DM halos produce remnants that retain excellent memories of the inner slopes and overall shapes of the phase-space density distribution of their progenitors. The robustness of the phase-space density profiles holds for a range of orbital energies, and a variety of encounter configurations including sequences of several consecutive merger events, designed to mimic hierarchical merging, and collisions occurring at different cosmological epochs. If the progenitor halos are constructed with appreciably different asymptotic power-law indices, we find that the inner slope and overall shape of the phase-space density distribution of the remnant are substantially closer to that of the initial system with the steepest central density cusp. These results explicitly demonstrate that mixing is incomplete in equal-mass mergers between DM halos, as it does not erase memory of the progenitor properties. Our results also confirm the recent analytical predictions of Dehnen (2005) regarding the preservation of merging self-gravitating central density cusps.

preprint2009arXiv

Modeling Molecular Hydrogen and Star Formation in Cosmological Simulations

We describe a phenomenological model for molecular hydrogen formation suited for applications in galaxy formation simulations, which includes on-equilibrium formation of molecular hydrogen on dust and approximate treatment of both its self-shielding and shielding by dust from the dissociating UV radiation. The model is applicable in simulations in which individual star forming regions - the giant molecular complexes - can be identified (resolution of tens of pc) and their mean internal density estimated reliably, even if internal structure is not resolved. In agreement with previous studies, calculations based on our model show that the transition from atomic to fully molecular phase depends primarily on the metallicity, which we assume is directly related to the dust abundance, and clumpiness of the interstellar medium. The clumpiness simply boosts the formation rate of molecular hydrogen, while dust serves both as a catalyst of molecular hydrogen formation and as an additional shielding from dissociating UV radiation. The upshot is that it is difficult to form fully-shielded giant molecular clouds while gas metallicity is low. However, once the gas is enriched to Z ~ 0.01-0.1 solar, the subsequent star formation and enrichment can proceed at a much faster rate. This may keep star formation efficiency in the low-mass, low-metallicity progenitors of galaxies very low for a certain period of time with the effect similar to a strong "feedback" mechanism. [abridged]

preprint2009arXiv

On the Kennicutt-Schmidt relation of low-metallicity high-redshift galaxies

We present results of self-consistent, high-resolution cosmological simulations of galaxy formation at z~3. The simulations employ recently developed recipe for star formation based on the local abundance of molecular hydrogen, which is tracked self-consistently during the course of simulation. The phenomenological H2 formation model accounts for the effects of dissociating UV radiation of stars in each galaxy, as well as self-shielding and shielding of H2 by dust, and therefore allows us to explore effects of lower metallicities and higher UV fluxes prevalent in high redshift galaxies on their star formation. We compare stellar masses, metallicities, and star formation rates of the simulated galaxies to available observations of the Lyman Break Galaxies (LBGs) and find a reasonable agreement. We find that the Kennicutt-Schmidt (KS) relation exhibited by our simulated galaxies at z~3 is substantially steeper and has a lower amplitude than the z=0 relation at Sigma_gas < 100 Msun/pc^2. The predicted relation, however, is consistent with existing observational constraints for the z~3 Damped Lyman $α$ (DLA) and LBGs. Our tests show that the main reason for the difference from the local KS relation is lower metallicity of the ISM in high redshift galaxies. We discuss several implications of the metallicity-dependence of the KS relation for galaxy evolution and interpretation of observations.

preprint2009arXiv

Residual Gas Motions in the Intracluster Medium and Bias in Hydrostatic Measurements of Mass Profiles of Clusters

We present analysis of bulk and random gas motions in the intracluster medium using high-resolution Eulerian cosmological simulations of sixteen simulated clusters, including both very relaxed and unrelaxed systems and spanning a virial mass range of 5*10^13 - 2*10^15 Msun/h. We investigate effects of the residual subsonic gas motions on the hydrostatic estimates of mass profiles and concentrations of galaxy clusters. In agreement with previous studies we find that the gas motions contribute up to ~ 5%-15% of the total pressure support in relaxed clusters with contribution increasing with cluster-centric radius. The fractional pressure support is higher in unrelaxed systems. This contribution would not be accounted for in hydrostatic estimates of the total mass profile and would lead to systematic underestimate of mass. We demonstrate that total mass can be recovered accurately if pressure due to gas motions measured in simulations is explicitly taken into account in the equation of hydrostatic equilibrium. Given that the underestimate of mass is increasing at larger radii, where gas is less relaxed and contribution of gas motions to pressure is larger, the total density profile derived from hydrostatic analysis is more concentrated than the true profile. This may at least partially explain some high values of concentrations of clusters estimated from hydrostatic analysis of X-ray data.

preprint2008arXiv

Cold Dark Matter Substructure and Galactic Disks I: Morphological Signatures of Hierarchical Satellite Accretion

(Abridged) We conduct a series of high-resolution, dissipationless N-body simulations to investigate the cumulative effect of substructure mergers onto thin disk galaxies in the context of the LCDM paradigm of structure formation. Our simulation campaign is based on a hybrid approach. Substructure properties are culled directly from cosmological simulations of galaxy-sized cold dark matter (CDM) halos. In contrast to what can be inferred from statistics of the present-day substructure populations, accretions of massive subhalos onto the central regions of host halos, where the galactic disk resides, since z~1 should be common occurrences. One host halo merger history is subsequently used to seed controlled numerical experiments of repeated satellite impacts on an initially-thin Milky Way-type disk galaxy. We show that these accretion events produce several distinctive observational signatures in the stellar disk including: a ring-like feature in the outskirts; a significant flare; a central bar; and faint filamentary structures that (spuriously) resemble tidal streams. The final distribution of disk stars exhibits a complex vertical structure that is well-described by a standard ``thin-thick'' disk decomposition. We conclude that satellite-disk encounters of the kind expected in LCDM models can induce morphological features in galactic disks that are similar to those being discovered in the Milky Way, M31, and in other disk galaxies. These results highlight the significant role of CDM substructure in setting the structure of disk galaxies and driving galaxy evolution. Upcoming galactic structure surveys and astrometric satellites may be able to distinguish between competing cosmological models by testing whether the detailed structure of galactic disks is as excited as predicted by the CDM paradigm.

preprint2007arXiv

Effects of Galaxy Formation on Thermodynamics of the Intracluster Medium

We present detailed comparisons of the intracluster medium (ICM) in cosmological Eulerian cluster simulations with deep Chandra observations of nearby relaxed clusters. To assess the impact of galaxy formation, we compare two sets of simulations, one performed in the non-radiative regime and another with radiative cooling and several physical processes critical to various aspects of galaxy formation: star formation, metal enrichment and stellar feedback. We show that the observed ICM properties outside cluster cores are well-reproduced in the simulations that include cooling and star formation, while the non-radiative simulations predict an overall shape of the ICM profiles inconsistent with observations. In particular, we find that the ICM entropy in our runs with cooling is enhanced to the observed levels at radii as large as half of the virial radius. We also find that outside cluster cores entropy scaling with the mean ICM temperature in both simulations and Chandra observations is consistent with being self-similar within current error bars. We find that the pressure profiles of simulated clusters are also close to self-similar and exhibit little cluster-to-cluster scatter. The X-ray observable-total mass relations for our simulated sample agree with the Chandra measurements to \~10%-20% in normalization. We show that this systematic difference could be caused by the subsonic gas motions, unaccounted for in X-ray hydrostatic mass estimates. The much improved agreement of simulations and observations in the ICM profiles and scaling relations is encouraging and the existence of tight relations of X-ray observables, such as Yx, and total cluster mass and the simple redshift evolution of these relations hold promise for the use of clusters as cosmological probes.

preprint2007arXiv

Escape of Ionizing Radiation from High Redshift Galaxies

We model the escape of ionizing radiation from high-redshift galaxies using high-resolution Adaptive Mesh Refinement N-body + hydrodynamics simulations. Our simulations include time-dependent and spatially-resolved transfer of ionizing radiation in three dimensions, including effects of dust absorption. For galaxies of total mass M > 10^11 Msun and star formation rates SFR ~ 1-5 Msun/yr, we find angular averaged escape fractions of 0.01-0.03 over the entire redshift interval studied (3<z<9). In addition, we find that the escape fraction varies by more than an order of magnitude along different lines-of-sight within individual galaxies, from the largest values near galactic poles to the smallest along the galactic disk. The escape fraction declines steeply at lower masses and SFR. We show that the low values of escape fractions are due to a small fraction of young stars located just outside the edge of HI disk. We compare our predicted escape fraction of ionizing photons with previous results, and find a general agreement with both other simulation results and available direct detection measurements at z ~ 3. We also compare our simulations with a novel method to estimate the escape fraction in galaxies from the observed distribution of neutral hydrogen column densities along the lines of sights to long duration gamma-ray bursts. Using this method we find escape fractions of the GRB host galaxies of 2-3%, consistent with our theoretical predictions. [abridged]

preprint2007arXiv

Scaling Relations of Dwarf Galaxies without Supernova-Driven Winds

Nearby dwarf galaxies exhibit tight correlations between their global stellar and dynamical properties, such as circular velocity, mass-to-light ratio, stellar mass, surface brightness, and metallicity. Such correlations have often been attributed to gas or metal-rich outflows driven by supernova energy feedback to the interstellar medium. We use high-resolution cosmological simulations of high-redshift galaxies with and without energy feedback, as well as analytic modeling, to investigate whether the observed correlations can arise without supernova-driven outflows. We find that the simulated dwarf galaxies exhibit correlations similar to those observed as early as z~10, regardless of whether supernova feedback is included. We also show that the correlations can be well reproduced by our analytic model that accounts for realistic gas inflow but assumes no outflows, and star formation rate obeying the Kennicutt-Schmidt law with a critical density threshold. We argue that correlations in simulated galaxies arise due to the increasingly inefficient conversion of gas into stars in low-mass dwarf galaxies rather than supernova-driven outflows. We also show that the decrease of the observed effective yield in low-mass objects, often used as an indicator of gas and metal outflows, can be reasonably reproduced in our simulations without outflows. We show that this trend can arise if a significant fraction of metals in small galaxies is spread to the outer regions of the halo outside the stellar extent via mixing. In this case the effective yield can be significantly underestimated if only metals within the stellar radius are taken into account. Measurements of gas metallicity in the outskirts of gaseous disks of dwarfs would thus provide a key test of such explanation.

preprint2006arXiv

A Large Dark Matter Core in the Fornax Dwarf Spheroidal Galaxy?

We use measurements of the stellar velocity dispersion profile of the Fornax dwarf spheroidal galaxy to derive constraints on its dark matter distribution. Though the data are unable to distinguish between models with small cores and those with cusps, we show that a large > 1 kpc dark matter core in Fornax is highly implausible. Irrespective of the origin of the core, reasonable dynamical limits on the mass of the Fornax halo constrain its core radius to be no larger than 700 pc. We derive an upper limit core radius of 300 pc by demanding that the central phase space density of Fornax not exceed that directly inferred from the rotation curves of low-mass spiral galaxies. Further, if the halo is composed of warm dark matter then phase-space constraints force the core to be quite small in order to avoid conservative limits from the Ly alpha forest power spectrum, implying a core radius < 85 pc. We discuss our results in the context of the idea that the extended globular cluster distribution in Fornax can be explained by the presence of a large 1.5 kpc core. A self-consistent core of this size would be drastically inconsistent with the expectations of standard warm or cold dark matter models, and would also require an unreasonably massive dark matter halo, with a maximum circular velocity of 200 km/s.

preprint2006arXiv

A New Robust Low-Scatter X-ray Mass Indicator for Clusters of Galaxies

We present comparison of X-ray proxies for the total cluster mass, including the spectral temperature (Tx), gas mass measured within r500 (Mg), and the new proxy, Yx, which is a simple product of Tx and Mg and is related to the total thermal energy of the ICM. We use mock Chandra images constructed for a sample of clusters simulated with the eulerian N-body+gasdynamics adaptive mesh refinement ART code in the concordance LCDM cosmology. The simulations achieve high spatial and mass resolution and include radiative cooling, star formation, and other processes accompanying galaxy formation. Our analysis shows that simulated clusters exhibit a high degree of regularity and tight correlations between the considered observables and total mass. The normalizations of the M-Tx, Mg-Tx, and M-Yx relations agree to better than 10-15% with the current observational measurements of these relations. Our results show that Yx is the best mass proxy with a remarkably low scatter of only ~5-7% in M500 for a fixed Yx, at both low and high redshifts and regardless of whether clusters are relaxed or not. In addition, we show that redshift evolution of the Yx-M500 relation is close to the self-similar prediction, which makes Yx a very attractive mass indicator for measurements of the cluster mass function from X-ray selected samples.

preprint2006arXiv

Modeling Chandra X-ray observations of Galaxy Clusters using Cosmological Simulations

X-ray observations of galaxy clusters potentially provide powerful cosmological probes if systematics due to our incomplete knowledge of the intracluster medium (ICM) physics are understood and controlled. In this paper, we study the effects of galaxy formation on the properties of the ICM and X-ray observable-mass relations using high-resolution self-consistent cosmological simulations of galaxy clusters and comparing their results with recent Chandra X-ray observations. We show that despite complexities of their formation and uncertainties in their modeling, clusters of galaxies both in observations and numerical simulations are remarkably regular outside of their cores, which holds great promise for their use as cosmological probes.

preprint2006arXiv

On relaxation processes in collisionless mergers

We analyze N-body simulations of halo mergers to investigate the mechanisms responsible for driving mixing in phase-space and the evolution to dynamical equilibrium. We focus on mixing in energy and angular momentum and show that mixing occurs in step-like fashion following pericenter passages of the halos. This makes mixing during a merger unlike other well known mixing processes such as phase mixing and chaotic mixing whose rates scale with local dynamical time. We conclude that the mixing process that drives the system to equilibrium is primarily a response to energy and angular momentum redistribution that occurs due to impulsive tidal shocking and dynamical friction rather than a result of chaotic mixing in a continuously changing potential. We also analyze the merger remnants to determine the degree of mixing at various radii by monitoring changes in radius, energy and angular momentum of particles. We confirm previous findings that show that the majority of particles retain strong memory of their original kinetic energies and angular momenta but do experience changes in their potential energies owing to the tidal shocks they experience during pericenter passages. Finally, we show that a significant fraction of mass (~ 40%) in the merger remnant lies outside its formal virial radius and that this matter is ejected roughly uniformly from all radii outside the inner regions. This highlights the fact that mass, in its standard virial definition, is not additive in mergers. We discuss the implications of these results for our understanding of relaxation in collisionless dynamical systems.

preprint2005arXiv

Impact of Dark Matter Substructure on the Matter and Weak Lensing Power Spectra

We explore the effect of substructure in dark matter halos on the power spectrum and bispectrum of matter fluctuations and weak lensing shear. By experimenting with substructure in a cosmological N = 512^3 simulation, we find that when a larger fraction of the host halo mass is in subhalos, the resulting power spectrum has less power at 1 < k < 100 h Mpc^{-1} and more power at k > 100 h Mpc^{-1}. We explain this effect using an analytic halo model including subhalos, which shows that the 1 < k < 100 h Mpc^{-1} regime depends sensitively on the radial distribution of subhalo centers while the interior structure of subhalos is important at k > 100 h Mpc^{-1}. The corresponding effect due to substructures on the weak lensing power spectrum is up to about 11% at angular scale l < 10^4. Predicting the nonlinear power spectrum to a few percent accuracy for future surveys would therefore require large cosmological simulations that also have exquisite numerical resolution to model accurately the survivals of dark matter subhalos in the tidal fields of their hosts.

preprint2005arXiv

Supersonic Motions of Galaxies in Clusters

We study motions of galaxies in galaxy clusters formed in the concordance LCDM cosmology. We use high-resolution cosmological simulations that follow dynamics of dark matter and gas and include various physical processes critical for galaxy formation: gas cooling, heating and star formation. Analysing motions of galaxies and the properties of intracluster gas in the sample of eight simulated clusters at z=0, we study velocity dispersion profiles of the dark matter, gas, and galaxies. We measure the mean velocity of galaxy motions and gas sound speed as a function of radius and calculate the average Mach number of galaxy motions. The simulations show that galaxies, on average, move supersonically with the average Mach number of ~1.4, approximately independent of the cluster-centric radius. The supersonic motions of galaxies may potentially provide an important source of heating for the intracluster gas by driving weak shocks and via dynamical friction, although these heating processes appear to be inefficient in our simulations. We also find that galaxies move faster than the dark matter particles in clusters. The magnitude of the velocity bias, b_v~1.1, is, however, smaller than the bias estimated for subhalos in dissipationless simulations. Interestingly, we find velocity bias in the tangential component of the velocity dispersion, but not in the radial component. Finally, we find significant random bulk motions of gas. The typical gas velocities are of order ~20-30% of the gas sound speed. These random motions provide about 10% of the total pressure support in our simulated clusters. The non-thermal pressure support, if neglected, will bias measurements of the total mass in the hydrostatic analyses of the X-ray cluster observations.

preprint2004arXiv

Modeling Galaxy-Mass Correlations in Dissipationless Simulations

We use high-resolution dissipationless simulations of the concordance flat LCDM model to make predictions for the galaxy--mass (GM) correlations and compare them to the recent SDSS weak lensing measurements.We use a simple observationally motivated scheme to assign luminosities and colors to the halos.This allows us to closely match the selection criteria used to define observational samples.The simulations reproduce the observed GM correlation function and its observed dependencies on luminosity and color.The luminosity dependence of the correlation function is primarily determined by the changing relative contribution of central and satellite galaxies at different luminosities. The color dependence of the GM correlations reflects the difference in the typical environments of blue and red galaxies. We also find agreement between the predicted and observed cross-bias, b_x=b/r,at all probed scales.The GM correlation coefficient, r, is close to unity on scales >1/h Mpc.The cross bias is thus expected to measure the actual bias of galaxy clustering on these scales.The aperture mass-to-light ratio is independent of galaxy color.The aperture mass scales approximately linearly with luminosity at L_r>10^{10}h^{-2} Lsun, while at lower luminosities the scaling is shallower: L_r^{0.5}. We show that most of the luminous galaxies (M_r<-21) are near the centers of their halos and their GM correlation function at r<100/h kpc can therefore be interpreted as the average dark matter density profile of these galaxies. We find that for galaxies in a given narrow luminosity range, there is a broad and possibly non-gaussian distribution of halo virial masses. Therefore, the average relation between mass and luminosity derived from the weak lensing analyses should be interpreted with caution.

preprint2004arXiv

Response of dark matter halos to condensation of baryons: cosmological simulations and improved adiabatic contraction model

The cooling of gas in the centers of dark matter halos is expected to lead to a more concentrated dark matter distribution. The response of dark matter to the condensation of baryons is usually calculated using the model of adiabatic contraction, which assumes spherical symmetry and circular orbits. In contrast, halos in the hierarchical structure formation scenarios grow via multiple violent mergers and accretion along filaments, and particle orbits in the halos are highly eccentric. We study the effects of the cooling of gas in the inner regions of halos using high-resolution cosmological simulations which include gas dynamics, radiative cooling, and star formation. We find that the dissipation of gas indeed increases the density of dark matter and steepens its radial profile in the inner regions of halos compared to the case without cooling. For the first time, we test the adiabatic contraction model in cosmological simulations and find that the standard model systematically overpredicts the increase of dark matter density in the inner 5% of the virial radius. We show that the model can be improved by a simple modification of the assumed invariant from M(r)r to M(r_av)r, where r and r_av are the current and orbit-averaged particle positions. This modification approximately accounts for orbital eccentricities of particles and reproduces simulation profiles to within 10-20%. We present analytical fitting functions that accurately describe the transformation of the dark matter profile in the modified model and can be used for interpretation of observations.

preprint2004arXiv

The tumultuous lives of galactic dwarfs and the missing satellites problem

Hierarchical Cold Dark Matter (CDM) models predict that Milky Way (MW) sized halos contain hundreds of dense low-mass dark satellites, an order of magnitude more than the number of observed satellites in the Local Group (LG). If the CDM paradigm is correct, we need to understand why most of these halos failed to form stars and become galaxies. We analyze the dynamical evolution of the satellite halos in a high-resolution cosmological simulation of MW sized halos in the LCDM cosmology. We find that about 10% of the substructure halos with the present masses <10^8-10^9 Msun (circular velocities Vmax<30 km/s) had considerably larger masses and circular velocities when they formed at redshifts z>2. After the initial period of mass accretion in isolation, these objects experience dramatic mass loss due to tidal stripping, in some cases even before they are accreted by their host halo. This can explain how the smallest dwarf spheroidal galaxies of the LG were able to build up a sizable stellar mass in their seemingly shallow potential wells. We propose a new model in which all of the luminous dwarf spheroidals in the Local Group are descendants of the relatively massive (>10^9 Msun) high-redshift systems, in which the gas could cool efficiently by atomic line emission and which were not significantly affected by the extragalactic ultraviolet radiation. We present a simple galaxy formation model based on the trajectories extracted from the simulation, which accounts for the bursts of star formation after strong tidal shocks and the inefficiency of gas cooling in halos with virial temperatures Tvir<~10^4 K. Our model reproduces the abundance, spatial distribution, and morphological segregation of the observed Galactic satellites. The results are insensitive to the redshift of reionization.

preprint2003arXiv

On the Origin of the Global Schmidt Law of Star Formation

One of the most puzzling properties of observed galaxies is the universality of the empirical correlation between the star formation rate and average gas surface density on kiloparsec scales (the Schmidt law). In this study I present results of self-consistent cosmological simulations of high-redshift galaxy formation that reproduce the Schmidt law naturally, without assuming it, and provide some clues to this puzzle. The simulations incorporate the main physical processes critical to various aspects of galaxy formation and have a dynamic range high enough to identify individual star forming regions. The results indicate that the global Schmidt law is a manifestation of the overall density distribution of the interstellar medium (ISM). In particular, the density probability distribution function (PDF) in the simulated disks is similar to that observed in recent state-of-the-art modeling of the turbulent ISM and has a well-defined generic shape. The shape of the PDF in a given region of the disk depends on the local average surface density Sigma_g. The dependence is such that the fraction of gas mass in the high-density tail of the distribution scales as Sigma_g^{n-1} with n~1.4, which gives rise to the Schmidt-like correlation. The high-density tail of the PDF is remarkably insensitive to the inclusion of feedback and details of the cooling and heating processes. This indicates that the global star formation rate is determined by the supersonic turbulence driven by gravitational instabilities on large scales, rather than stellar feedback or thermal instabilities on small scales.

preprint2002arXiv

Effect of Internal Flows on Sunyaev-Zeldovich Measurements of Cluster Peculiar Velocities

Galaxy clusters are potentially powerful probes of the large-scale velocity field in the Universe because their peculiar velocity can be estimated directly via the kinematic Sunyaev-Zeldovich effect (kSZ). Using high-resolution cosmological simulations of an evolving cluster of galaxies, we evaluate how well the average velocity obtained via a kSZ measurement reflects the actual cluster peculiar velocity. We find that the internal velocities in the intracluster gas are comparable to the overall cluster peculiar velocity, 20 to 30% of the sound speed even when a cluster is relatively relaxed. Nevertheless, the velocity averaged over the kSZ map inside a circular aperture matched to the cluster virial region provides an unbiased estimate of a cluster's radial peculiar velocity with a dispersion of 50 to 100 km/s, depending on the line of sight and dynamical state of the cluster. This dispersion puts a lower limit on the accuracy with which cluster peculiar velocity can be measured. Although the dispersion of the average is relatively small, the velocity distribution is broad; regions of low signal must be treated with care to avoid bias. We discuss the extent to which systematic errors might be modelled, and the resulting limitations on using galaxy clusters as cosmological velocity tracers.

preprint2002arXiv

Sample Variance Considerations for Cluster Surveys

We present a general statistical framework for describing the effect of sample variance in the number counts of virialized objects and examine its effect on cosmological parameter estimation. Specifically, we consider effects of sample variance on the power spectrum normalization and properties of dark energy extracted from current and future local and high-redshift samples of clusters. We show that for future surveys that probe ever lower cluster masses and temperatures, sample variance is generally comparable to or greater than shot noise and thus cannot be neglected in deriving precision cosmological constraints. For example, sample variance is usually more important than shot variance in constraints on the equation of state of the dark energy from z < 1 clusters. Although we found that effects of sample variance on the sigma_8-Omega_m constraints from the current flux and temperature limited X-ray surveys are not significant, they may be important for future studies utilizing the shape of the temperature function to break the sigma_8-Omega_m degeneracy. We also present numerical tests clarifying the definition of cluster mass employed in cosmological modelling and accurate fitting formula for the conversion between different definitions of halo mass (e.g., virial vs. fixed overdensity).

preprint2000arXiv

A universal angular momentum profile for galactic halos

[Abridged] We study the angular-momentum profiles of a statistical sample of halos drawn from a high-resolution N-body simulation of the LCDM cosmology. We find that the cumulative mass distribution of specific angular momentum, j, in a halo of mass Mv is well fit by a universal function, M(<j) = Mv μj/(j_0+j). This profile is defined by one shape parameter (μor j_0) in addition to the global spin parameter λ. It follows a power-law over most of the mass, and flattens at large j, with the flattening more pronounced for small values of μ. Compared to a uniform sphere in solid-body rotation, most halos have a higher fraction of their mass in the low- and high-j tails of the distribution. The spatial distribution of angular momentum in halos tends to be cylindrical and is well-aligned within each halo for ~80% of the halos. We investigate two ideas for the origin of this profile. The first is based on a revised version of linear tidal-torque theory combined with extended Press-Schechter mass accretion, and the second focuses on j transport in minor mergers. Finally, we briefly explore implications of the M(<j) profile on the formation of galactic disks assuming that j is conserved during an adiabatic baryonic infall. The implied gas density profile deviates from an exponential disk, with a higher density at small radii and a tail extending to large radii. The steep central density profiles may imply disk scale lengths that are smaller than observed. This is reminiscent of the "angular-momentum problem" seen in hydrodynamic simulations, even though we have assumed perfect j conservation. A possible solution is to associate the central excesses with bulge components and the outer regions with extended gaseous disks.

preprint2000arXiv

On the supernovae heating of intergalactic medium

We present estimates of the energy input from supernovae (SNe) into the intergalactic medium using (i) recent measurements of Si and Fe abundances in the intracluster medium (ICM) and (ii) self-consistent gasdynamical galaxy formation simulations that include processes of cooling, star formation, SNe feedback, and a multi-phase model of the interstellar medium. We estimate the energy input from observed abundances using two different assumptions: (i) spatial uniformity of metal abundances in the ICM and (ii) radial abundance gradients. We show that these two cases lead to energy input estimates which are different by an order of magnitude, highlighting a need for observational data on large-scale abundance gradients in clusters. Analysis of galaxy formation results and estimates from observed Fe and Si abundances indicates that the SNe energy input can be important for heating of the entire ICM (providing energy of ~1 keV per particle) only if the ICM abundances are uniform and the efficiency of gas heating by SN explosions is close to 100% (implying that all of the initial kinetic energy of the explosion goes into heating of the ICM). We conclude that unless these most favorable conditions are met, SNe alone are unlikely to provide sufficient energy input to heat all of the cluster ICM and may need to be supplemented or even substituted by some other heating process(es). (Abridged)

preprint2000arXiv

Profiles of dark haloes: evolution, scatter, and environment

We study dark-matter halo density profiles in a high-resolution N-body simulation of an LCDM cosmology. Our statistical sample contains ~5000 haloes in the range 10^{11}-10^{14} M_sun. The profiles are parameterized by an NFW form with two parameters, an inner radius r_s and a virial radius r_v, and we define the halo concentration c_v = r_v/r_s. First, we find that, for a given halo mass, the redshift dependence of the median concentration is c_v ~ 1/(1+z), corresponding to a roughly constant r_s with redshift. We present an improved analytic treatment of halo formation that fits the measured relations between halo parameters and their redshift dependence. The implications are that high-redshift galaxies are predicted to be more extended and dimmer than expected before. Second, we find that the scatter in log(c_v) is ~0.18, corresponding to a scatter in maximum rotation velocities of dV/V ~ 0.12. We discuss implications for modelling the Tully-Fisher relation, which has a smaller reported intrinsic scatter. Third, haloes in dense environments tend to be more concentrated than isolated haloes. These results suggest that c_v is an essential parameter for the theory of galaxy modelling, and we briefly discuss implications for the universality of the Tully-Fisher relation, the formation of low surface brightness galaxies, and the origin of the Hubble sequence.

Andrey V. Kravtsov

Quick read

Decide how to stay connected

How to connect with this researcher

Open a focused conversation when the fit is right

See the researcher in context

Building this graph slice

57 published item(s)

Entropy-Conserving Scheme for Modeling Nonthermal Energies in Fluid Dynamics Simulations

The Sheet of Giants: Unusual Properties of the Milky Way&#39;s Immediate Neighbourhood

The Three Causes of Low-Mass Assembly Bias

A universal model for halo concentrations

On the interplay between star formation and feedback in galaxy formation simulations

Dependence of the outer density profiles of halos on their mass accretion rate

On the mass of the Local Group

Modelling mass distribution in elliptical galaxies: mass profiles and their correlation with velocity dispersion profiles

On the evolution of cluster scaling relations

Satellites in MW-like hosts: Environment dependence and close pairs

Simulations of disk galaxies with cosmic ray driven galactic winds

A High Throughput Workflow Environment for Cosmological Simulations

Constraining Cluster Physics with the Shape of X-ray Clusters: Comparison of Local X-ray Clusters versus LCDM Clusters

On the Origin of the High Column Density Turnover in the HI Column Density Distribution

Semi-empirical catalog of early-type galaxy-halo systems: dark matter density profiles, halo contraction and dark matter annihilation strength

The X-factor in Galaxies: II. The molecular hydrogen -- star formation relation

Towards a complete accounting of energy and momentum from stellar feedback in galaxy formation simulations

Fuel Efficient Galaxies: Sustaining Star Formation with Stellar Mass Loss

Halo Contraction Effect in Hydrodynamic Simulations of Galaxy Formation

Implementing the DC Mode in Cosmological Simulations with Supercomoving Variables

On determining the shape of matter distributions

On the Accuracy of Weak Lensing Cluster Mass Reconstructions

Shapes of Gas, Gravitational Potential and Dark Matter in Lambda-CDM Clusters

The impact of baryon physics on the structure of high-redshift galaxies

The X-factor in Galaxies: I. Dependence on Environment and Scale

Ultra-faint dwarf galaxies as a test of early enrichment and metallicity-dependent star formation

Environmental Dependence of the Kennicutt-Schmidt Relation in Galaxies

The distribution of annihilation luminosities in dark matter substructure

The Large Scale Bias of Dark Matter Halos: Numerical Calibration and Model Tests

The Metal-Enriched Outer Disk of NGC 2915

Computational Eulerian Hydrodynamics and Galilean Invariance

Cosmological Radiative Transfer Comparison Project II: The Radiation-Hydrodynamic Tests

Dark matter substructure and dwarf galactic satellites

Effects of Baryon Dissipation on the Dark Matter Virial Scaling Relation

Evolution of Dark Matter Phase-Space Density Distributions in Equal-Mass Halo Mergers

Modeling Molecular Hydrogen and Star Formation in Cosmological Simulations

On the Kennicutt-Schmidt relation of low-metallicity high-redshift galaxies

Residual Gas Motions in the Intracluster Medium and Bias in Hydrostatic Measurements of Mass Profiles of Clusters

Cold Dark Matter Substructure and Galactic Disks I: Morphological Signatures of Hierarchical Satellite Accretion

Effects of Galaxy Formation on Thermodynamics of the Intracluster Medium

Escape of Ionizing Radiation from High Redshift Galaxies

Scaling Relations of Dwarf Galaxies without Supernova-Driven Winds

A Large Dark Matter Core in the Fornax Dwarf Spheroidal Galaxy?

A New Robust Low-Scatter X-ray Mass Indicator for Clusters of Galaxies

Modeling Chandra X-ray observations of Galaxy Clusters using Cosmological Simulations

On relaxation processes in collisionless mergers

Impact of Dark Matter Substructure on the Matter and Weak Lensing Power Spectra

Supersonic Motions of Galaxies in Clusters

Modeling Galaxy-Mass Correlations in Dissipationless Simulations

Response of dark matter halos to condensation of baryons: cosmological simulations and improved adiabatic contraction model

The tumultuous lives of galactic dwarfs and the missing satellites problem

On the Origin of the Global Schmidt Law of Star Formation

Effect of Internal Flows on Sunyaev-Zeldovich Measurements of Cluster Peculiar Velocities

Sample Variance Considerations for Cluster Surveys

A universal angular momentum profile for galactic halos

On the supernovae heating of intergalactic medium

Profiles of dark haloes: evolution, scatter, and environment

The Sheet of Giants: Unusual Properties of the Milky Way's Immediate Neighbourhood