Source author record

Andrey Kravtsov

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

ResearcherUnclaimed source record

Catalog footprint

What is connected

19works
5topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

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

Building this map preview

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

Published work

19 published item(s)

preprint2025arXiv

Picture an Astronomer: Best Practices for Retaining Talent in Astrophysics

Women are consistently underrepresented in astrophysics yet are simultaneously subject to disproportionate attrition at every career stage. This disparity between demonstrated efficacy in job performance and ultimate career outcome was the primary motivation for the Picture an Astronomer series, which included both targeted public outreach to increase representation of women in astrophysics and high-level, solution-oriented discussions among professional astronomers. In March 2025, more than 200 astronomers came together in a hybrid-format symposium focused on the state of the field for female scientists, combining scientific exchange with discussions of policies and practices to strengthen retention of talent in the field. This white paper is the result of those discussions, offering a wide range of recommendations developed in the context of gendered attrition in astrophysics but which ultimately support a healthier climate for all scientists alike.

preprint2022arXiv

Forward-modelling the Luminosity, Distance, and Size distributions of the Milky Way Satellites

We use \texttt{GRUMPY}, a simple regulator-type model for dwarf galaxy formation and evolution, to forward model the dwarf galaxy satellite population of the Milky Way (MW) using the Caterpillar zoom-in simulation suite. We show that luminosity and distance distributions of the model satellites are consistent with the distributions measured in the DES, PS1 and SDSS surveys, even without including a model for the orphan galaxies. We also show that our model for dwarf galaxy sizes can simultaneously reproduce the observed {\it distribution} of stellar half-mass radii, $r_{1/2}$, of the MW satellites and the overall $r_{1/2}-M_\star$ relation exhibited by observed dwarf galaxies. The model predicts that some of the observed faint stellar systems with $r_{1/2}<10$ pc are ultra-faint dwarf galaxies. Scaling of the stellar mass $M_\star$ and peak halo mass $M_{\rm peak}$ for the model satellites is not described by a power law, but has a clear flattening of $M_\star-M_{\rm peak}$ scaling at $M_{\rm peak}<10^8\,M_\odot$ imprinted by reionization. As a result, the fraction of low mass haloes ($M_{\rm peak} < 10^8\, M_\odot$) hosting galaxies with $M_V<0$ is predicted to be 50% at $M_{\rm peak} \sim 3.6 \times 10^7\,M_\odot$. We find that such high fraction at that halo mass helps to reproduce the number of dwarf galaxies discovered recently in the HSC-SSP survey. Using the model we forecast that there should be the total of $440^{+201}_{-147}$ (68\% confidence interval) MW satellites with $M_V < 0$ and $r_{1/2} > 10$ pc within 300 kpc and make specific predictions for the HSC-SSP, DELVE-WIDE and LSST surveys.

preprint2022arXiv

From dawn till disk: Milky Way's turbulent youth revealed by the APOGEE+Gaia data

We use accurate estimates of aluminium abundance provided as part of the APOGEE Data Release 17 and Gaia Early Data Release 3 astrometry to select a highly pure sample of stars with metallicity $-1.5\lesssim {\rm [Fe/H]}\lesssim 0.5$ born in-situ in the Milky Way proper. We show that the low-metallicity ([Fe/H]$\lesssim -1.3$) in-situ component that we dub Aurora is kinematically hot with an approximately isotropic velocity ellipsoid and a modest net rotation. Aurora stars exhibit large scatter in metallicity and in a number of element abundance ratios. The median tangential velocity of the in-situ stars increases sharply with increasing metallicity between [Fe/H]$=-1.3$ and $-0.9$, the transition that we call the spin-up. The observed and theoretically expected age-metallicity correlations imply that this increase reflects a rapid formation of the Milky Way disk over $\approx 1-2$ Gyrs. The transformation of the stellar kinematics as a function of [Fe/H] is accompanied by a qualitative change in chemical abundances: the scatter drops sharply once the Galaxy builds up a disk during later epochs corresponding to [Fe/H]$>-0.9$. Results of galaxy formation models presented in this and other recent studies strongly indicate that the trends observed in the Milky Way reflect generic processes during the early evolution of progenitors of MW-sized galaxies: a period of chaotic pre-disk evolution, when gas is accreted along cold narrow filaments and when stars are born in irregular configurations, and subsequent rapid disk formation. The latter signals formation of a stable hot gaseous halo around the MW progenitor, which changes the mode of gas accretion and allows development of coherently rotating disk.

preprint2022arXiv

GRUMPY: a simple framework for realistic forward-modelling of dwarf galaxies

We present a simple regulator-type framework designed specifically for modelling formation of dwarf galaxies. Despite its simplicity, when coupled with realistic mass accretion histories of haloes from simulations and reasonable choices for model parameter values, the framework can reproduce a remarkably broad range of observed properties of dwarf galaxies over seven orders of magnitude in stellar mass. In particular, we show that the model can simultaneously match observational constraints on the stellar mass--halo mass relation, as well as observed relations between stellar mass and gas phase and stellar metallicities, gas mass, size, and star formation rate, as well as general form and diversity of star formation histories (SFHs) of observed dwarf galaxies. The model can thus be used to predict photometric properties of dwarf galaxies hosted by dark matter haloes in $N$-body simulations, such as colors, surface brightnesses, and mass-to-light ratios and to forward model observations of dwarf galaxies. We present examples of such modelling and show that colors and surface brightness distributions of model galaxies are in good agreement with observed distributions for dwarfs in recent observational surveys. We also show that in contrast with the common assumption, the absolute magnitude-halo mass relation is generally predicted to have a non-power law form in the dwarf regime, and that the fraction of haloes that host detectable ultrafaint galaxies is sensitive to reionization redshift ($z_{\rm rei}$) and is predicted to be consistent with observations for $z_{\rm rei}\lesssim 9$.

preprint2022arXiv

Not so fuzzy: excluding FDM with sizes and stellar kinematics of ultra-faint dwarf galaxies

We use observations of ultra-faint dwarf (UFD) galaxies to constrain the particle mass of ultra-light dark matter. Potential fluctuations created by wave interference in virialized "fuzzy" dark matter (FDM) halos dynamically heat stellar orbits in UFDs, some of which exhibit velocity dispersions of $\lesssim$ 3 km/s and sizes $\lesssim$ 40 pc. Using simulations of FDM halos, and existing measurements of sizes and stellar radial velocities in Segue 1 and Segue 2 UFDs, we derive a lower limit on the dark matter particle mass of $m_{fdm} > 3\times 10^{-19}$ eV at 99% confidence, marginalized over host halo circular velocity. This constraint is conservative as it is derived under the assumption that soliton heating is negligible, and that no other sources of non-FDM dynamical heating of stars operate to increase velocity dispersion. It can potentially be strengthened by future spectroscopic observations of additional stars in ultra-faint galaxies and by tightening theoretical constraints on the soliton size-halo mass relation. However, even the current conservative lower limit on the FDM mass makes this model indistinguishable from Cold Dark Matter at the scales probed by existing astronomical observations.

preprint2019arXiv

Signatures of Self-Interacting dark matter on cluster density profile and subhalo distributions

Non-gravitational interactions between dark matter particles with strong scattering, but relatively small annihilation and dissipation, has been proposed to match various observables on cluster and group scales. In this paper, we present the results from large cosmological simulations which include the effects of different self-interaction scenarios. In particular we explore a model with the differential cross section that can depend on both the relative velocity of the interacting particles and the angle of scattering. We focus on how quantities, such as the stacked density profiles, subhalo counts and the splashback radius change as a function of different forms of self-interaction. We find that self-interactions not only affect the central region of the cluster, the effect well known from previous studies, but also significantly alter the distribution of subhalos and the density of particles out to the splashback radius. Our results suggest that current weak lensing data can already put constraints on the self-interaction cross-section that are only slightly weaker than the Bullet Cluster constraints ($σ/m \lesssim 2$ cm$^2/$g), and future lensing surveys should be able to tighten them even further making halo profiles on cluster scales a competitive probe for DM physics.

preprint2016arXiv

On Detecting Halo Assembly Bias with Galaxy Populations

The fact that the clustering of dark matter halos depends not only on their mass, but also the formation epoch, is a prominent, albeit subtle, feature of the cold dark matter structure formation theory, and is known as assembly bias. At low mass scales ($\sim 10^{12}\,h^{-1}M_\odot$), early-forming halos are predicted to be more strongly clustered than the late-forming ones. In this study we aim to robustly detect the signature of assembly bias observationally, making use of formation time indicators of central galaxies in low mass halos as a proxy for the halo formation history. Weak gravitational lensing is employed to ensure our early- and late-forming halo samples have similar masses, and are free of contamination of satellites from more massive halos. For the two formation time indicators used (resolved star formation history and current specific star formation rate), we do not find convincing evidence of assembly bias. For a pair of early- and late-forming galaxy samples with mean mass $M_{200c} \approx 9\times10^{11}\,h^{-1}M_\odot$, the relative bias is $1.00\pm 0.12$. We attribute the lack of detection to the possibilities that either the current measurements of these indicators are too noisy, or they do not correlate well with the halo formation history. Alternative proxies for the halo formation history that should perform better are suggested for future studies.

preprint2015arXiv

Preventing Star Formation in Early-Type Galaxies with Late-Time Stellar Heating

We revisit previous suggestions that the heating provided by the winds of dying low-mass stars plays an important role in preventing star formation in quiescent galaxies. At the end of their asymptotic giant branch phase, intermediate and low-mass stars eject their envelopes rapidly in a super-wind phase, usually giving rise to planetary nebulae. In spheroidal galaxies with high stellar velocity dispersions, the interaction of these ejected envelopes with the ambient diffuse gas can lead to significant, isotropic and steady-state heating that scales as $\dot{M}_\astσ_\ast^2$. We show that cooling of the central regions of the hot diffuse halo gas can be delayed for a Hubble time for halos more massive than $\sim10^{12.5}M_{\odot}$ at $0<z<2$, although stellar heating alone is unlikely to forestall cooling in the most massive clusters at $z=0$. This mechanism provides a natural explanation for the strong trend of galaxy quiescence with stellar surface density and velocity dispersion. In addition, since the ejected material will thermalize to $kT\simσ_\ast^2$, this mechanism provides an explanation for the observed similarity between the central temperature of the hot diffuse gas and $σ_\ast^2$, a result which is not trivial in light of the short inferred cooling times of the hot gas. The main uncertainty in this analysis is the ultimate fate of the stellar ejecta. Preventing accumulation of the ejecta in the central regions may require energy input from another source, such as Type Ia supernovae. Detailed simulations of the interaction of the stellar wind with the ambient gas are required to better quantify the net effect of AGB heating.

preprint2015arXiv

The splashback radius as a physical halo boundary and the growth of halo mass

The boundaries of cold dark matter halos are commonly defined to enclose a density contrast $Δ$ relative to a reference (mean or critical) density. We argue that a more physical boundary of halos is the radius at which accreted matter reaches its first orbital apocenter after turnaround. This splashback radius, $R_{sp}$, manifests itself as a sharp density drop in the halo outskirts, at a location that depends upon the mass accretion rate. We present calibrations of $R_{sp}$ and the enclosed mass, $M_{sp}$, as a function of the accretion rate and alternatively peak height. We find that $R_{sp}$ varies between $\approx0.8-1R_{200m}$ for rapidly accreting halos and $\approx1.5R_{200m}$ for slowly accreting halos. The extent of a halo and its associated environmental effects can thus extend well beyond the conventionally defined "virial" radius. We show that $M_{sp}$ and $R_{sp}$ evolve relatively strongly compared to other commonly used definitions. In particular, $M_{sp}$ evolves significantly even for the smallest dwarf-sized halos at $z=0$. We also contrast $M_{sp}$ with the mass enclosed within four scale radii of the halo density profile, $M_{<4rs}$, which characterizes the inner halo. During the early stages of halo assembly, $M_{sp}$ and $M_{<4rs}$ evolve similarly, but in the late stages $M_{<4rs}$ stops increasing while $M_{sp}$ continues to grow significantly. This illustrates that halos at low $z$ can have "quiet" interiors while continuing to accrete mass in their outskirts. We discuss potential observational estimates of the splashback radius and show that it may already have been detected in galaxy clusters.

preprint2014arXiv

On the intermediate-redshift central stellar mass-halo mass relation, and implications for the evolution of the most massive galaxies since z~1

The stellar mass-halo mass relation is a key constraint in all semi-analytic, numerical, and semi-empirical models of galaxy formation and evolution. However, its exact shape and redshift dependence remain debated. Several recent works support a relation in the local Universe steeper than previously thought. Based on the comparisons with a variety of data on massive central galaxies, we show that this steepening holds up to z~1, for stellar masses Mstar>2e11 Msun. Specifically, we find significant evidence for a high-mass end slope of β>0.35-0.70, instead of the usual β~0.20-0.30 reported by a number of previous results. When including the independent constraints from the recent BOSS clustering measurements, the data, independent of any systematic errors in stellar masses, tend to favor a model with a very small scatter (< 0.15 dex) in stellar mass at fixed halo mass, in the redshift range z < 0.8 and for Mstar>3e11 Msun, suggesting a close connection between massive galaxies and host halos even at relatively recent epochs. We discuss the implications of our results with respect to the evolution of the most massive galaxies since z~1.

preprint2013arXiv

The pseudo-evolution of halo mass

A dark matter halo is commonly defined as a spherical overdensity of matter with respect to a reference density, such as the critical density or the mean matter density of the Universe. Such definitions can lead to a spurious pseudo-evolution of halo mass simply due to redshift evolution of the reference density, even if its physical density profile remains constant over time. We estimate the amount of such pseudo-evolution of mass between z=1 to 0 for halos identified in a large N-body simulation, and show that it accounts for almost the entire mass evolution of the majority of halos with M200 of about 1E12 solar masses and can be a significant fraction of the apparent mass growth even for cluster-sized halos. We estimate the magnitude of the pseudo-evolution assuming that halo density profiles remain static in physical coordinates, and show that this simple model predicts the pseudo-evolution of halos identified in numerical simulations to good accuracy, albeit with significant scatter. We discuss the impact of pseudo-evolution on the evolution of the halo mass function and show that the non-evolution of the low-mass end of the halo mass function is the result of a fortuitous cancellation between pseudo-evolution and the absorption of small halos into larger hosts. We also show that the evolution of the low mass end of the concentration-mass relation observed in simulations is almost entirely due to the pseudo-evolution of mass. Finally, we discuss the implications of our results for the interpretation of the evolution of various scaling relations between the observable properties of galaxies and galaxy clusters and their halo masses.

preprint2012arXiv

Formation of Galaxy Clusters

In this review, we describe our current understanding of cluster formation: from the general picture of collapse from initial density fluctuations in an expanding Universe to detailed simulations of cluster formation including the effects of galaxy formation. We outline both the areas in which highly accurate predictions of theoretical models can be obtained and areas where predictions are uncertain due to uncertain physics of galaxy formation and feedback. The former includes the description of the structural properties of the dark matter halos hosting cluster, their mass function and clustering properties. Their study provides a foundation for cosmological applications of clusters and for testing the fundamental assumptions of the standard model of structure formation. The latter includes the description of the total gas and stellar fractions, the thermodynamical and non-thermal processes in the intracluster plasma. Their study serves as a testing ground for galaxy formation models and plasma physics. In this context, we identify a suitable radial range where the observed thermal properties of the intra-cluster plasma exhibit the most regular behavior and thus can be used to define robust observational proxies for the total cluster mass. We put particular emphasis on examining assumptions and limitations of the widely used self-similar model of clusters. Finally, we discuss the formation of clusters in non-standard cosmological models, such as non-Gaussian models for the initial density field and models with modified gravity, along with prospects for testing these alternative scenarios with large cluster surveys in the near future.

preprint2011arXiv

Baryon Content of Massive Galaxy Clusters at z=0-0.6

We study the relationship between two major baryonic components in galaxy clusters, namely the stars in galaxies, and the ionized gas in the intracluster medium (ICM), using 94 clusters that span the redshift range 0-0.6. Accurately measured total and ICM masses from Chandra observations, and stellar masses derived from the Wide-field Infrared Survey Explorer and the Two-Micron All-Sky Survey allow us to trace the evolution of cluster baryon content in a self-consistent fashion. We find that, within r_{500}, the evolution of the ICM mass--total mass relation is consistent with the expectation of self-similar model, while there is no evidence for redshift evolution in the stellar mass--total mass relation. This suggests that the stellar mass and ICM mass in the inner parts of clusters evolve differently.

preprint2011arXiv

The overdensity and masses of the friends-of-friends halos and universality of the halo mass function

The friends-of-friends algorithm (hereafter, FOF) is a percolation algorithm which is routinely used to identify dark matter halos from N-body simulations. We use results from percolation theory to show that the boundary of FOF halos does not correspond to a single density threshold but to a range of densities close to a critical value that depends upon the linking length parameter, b. We show that for the commonly used choice of b = 0.2, this critical density is equal to 81.62 times the mean matter density. Consequently, halos identified by the FOF algorithm enclose an average overdensity which depends on their density profile (concentration) and therefore changes with halo mass contrary to the popular belief that the average overdensity is ~180. We derive an analytical expression for the overdensity as a function of the linking length parameter b and the concentration of the halo. Results of tests carried out using simulated and actual FOF halos identified in cosmological simulations show excellent agreement with our analytical prediction. We also find that the mass of the halo that the FOF algorithm selects crucially depends upon mass resolution. We find a percolation theory motivated formula that is able to accurately correct for the dependence on number of particles for the mock realizations of spherical and triaxial Navarro-Frenk-White halos. However, we show that this correction breaks down when applied to the real cosmological FOF halos due to presence of substructures. Given that abundance of substructure depends on redshift and cosmology, we expect that the resolution effects due to substructure on the FOF mass and halo mass function will also depend on redshift and cosmology and will be difficult to correct for in general. Finally, we discuss the implications of our results for the universality of the mass function.

preprint2009arXiv

Collapse Barriers and Halo Abundance: Testing the Excursion Set Ansatz

Our heuristic understanding of the abundance of dark matter halos centers around the concept of a density threshold, or "barrier", for gravitational collapse. If one adopts the ansatz that regions of the linearly evolved density field smoothed on mass scale M with an overdensity that exceeds the barrier will undergo gravitational collapse into halos of mass M, the corresponding abundance of such halos can be estimated simply as a fraction of the mass density satisfying the collapse criterion divided by the mass M. The key ingredient of this ansatz is therefore the functional form of the collapse barrier as a function of mass M or, equivalently, of the variance sigma^2(M). Several such barriers based on the spherical, Zel'dovich, and ellipsoidal collapse models have been extensively discussed. Using large scale cosmological simulations, we show that the relation between the linear overdensity and the mass variance for regions that collapse to form halos by the present epoch resembles expectations from dynamical models of ellipsoidal collapse. However, we also show that using such a collapse barrier with the excursion set ansatz predicts a halo mass function inconsistent with that measured directly in cosmological simulations. This inconsistency demonstrates a failure of the excursion set ansatz as a physical model for halo collapse. We discuss implications of our results for understanding the collapse epoch for halos as a function of mass, and avenues for improving consistency between analytical models for the collapse epoch and the results of cosmological simulations.

preprint2009arXiv

Cosmological simulations of galaxy clusters

We review recent progress in the description of the formation and evolution of galaxy clusters in a cosmological context by using numerical simulations. We focus our presentation on the comparison between simulated and observed X-ray properties, while we will also discuss numerical predictions on properties of the galaxy population in clusters. Many of the salient observed properties of clusters, such as X-ray scaling relations, radial profiles of entropy and density of the intracluster gas, and radial distribution of galaxies are reproduced quite well. In particular, the outer regions of cluster at radii beyond about 10 per cent of the virial radius are quite regular and exhibit scaling with mass remarkably close to that expected in the simplest case in which only the action of gravity determines the evolution of the intra-cluster gas. However, simulations generally fail at reproducing the observed cool-core structure of clusters: simulated clusters generally exhibit a significant excess of gas cooling in their central regions, which causes an overestimate of the star formation and incorrect temperature and entropy profiles. The total baryon fraction in clusters is below the mean universal value, by an amount which depends on the cluster-centric distance and the physics included in the simulations, with interesting tensions between observed stellar and gas fractions in clusters and predictions of simulations. Besides their important implications for the cosmological application of clusters, these puzzles also point towards the important role played by additional physical processes, beyond those already included in the simulations. We review the role played by these processes, along with the difficulty for their implementation, and discuss the outlook for the future progress in numerical modeling of clusters.

preprint2009arXiv

Evolution of the Dark Matter Phase-Space Density Distributions of LCDM Halos

We study the evolution of phase-space density during the hierarchical structure formation of LCDM halos. We compute both a spherically-averaged surrogate for phase-space density (Q) and the coarse-grained distribution function f(x,v) for dark matter particles that lie within~2 virial radii of four Milky-Way-sized dark matter halos. The estimated f(x,v) spans over four decades at any radius. Dark matter particles that end up within two virial radii of a Milky-Way-sized DM halo at $z=0$ have an approximately Gaussian distribution in log(f) at early redshifts, but the distribution becomes increasingly skewed at lower redshifts. The value corresponding to the peak of the Gaussian decreases as the evolution progresses and is well described by a power-law in (1+z). The highest values of f are found at the centers of dark matter halos and subhalos, where f can be an order of magnitude higher than in the center of the main halo. The power-law Q(r) profile likely reflects the distribution of entropy (K = sigma^2/rho^{2/3} \propto r^{1.2}), which dark matter acquires as it is accreted onto a growing halo. The estimated f(x, v), on the other hand, exhibits a more complicated behavior. Although the median coarse-grained phase-space density profile F(r) can be approximated by a power-law in the inner regions of halos and at larger radii the profile flattens significantly. This is because phase-space density averaged on small scales is sensitive to the high-f material associated with surviving subhalos, as well as relatively unmixed material (probably in streams) resulting from disrupted subhalos, which contribute a sizable fraction of matter at large radii. (ABRIDGED)

preprint2006arXiv

The Impact of Baryonic Cooling on Giant Arc Abundances

Using ray tracing for simple analytic profiles, we demonstrate that the lensing cross section for producing giant arcs has distinct contributions due to arcs formed through image distortion only, and arcs form from the merging of two or three images. We investigate the dependence of each of these contributions on halo ellipticity and on the slope of the density profile, and demonstrate that at fixed Einstein radius, the lensing cross section increases as the halo profile becomes steeper. We then compare simulations with and without baryonic cooling of the same cluster for a sample of six clusters, and demonstrate that cooling can increase the overall abundance of giant arcs by factors of a few. The net boost to the lensing probability for individual clusters is mass dependent, and can lower the effective low mass limit of lensing clusters. This last effect can potentially increase the number of lensing clusters by an extra 50%. While the magnitude of these effects may be overestimated due to the well known overcooling problem in simulations, it is evident that baryonic cooling has a non-negligible impact on the expected abundance of giant arcs, and hence cosmological constraints from giant arc abundances may be subject to large systematic errors.