Source author record

Andrea V. Macciò

Andrea V. Macciò 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

40works
2topics
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

40 published item(s)

preprint2022arXiv

Cold-mode and hot-mode accretion in galaxy formation: an entropy approach

We have analysed two cosmological zoom simulations with $M_{\rm vir}\sim 10^{12}{\rm\,M}_\odot$ from the NIHAO series, both with and without feedback. We show that an entropy criterion based on the equation of state of the intergalactic medium can successfully separate cold- and hot-mode accretion. The shock-heated gas has non-negligible turbulent support and cools inefficiently. In the simulations without feedback, only a small fraction ($\sim 20$ per cent) of the stellar mass comes from baryons that have been in the hot circumgalactic medium, although quantitative conclusions should be taken with caution due to our small-number statistics. With feedback, the fraction is larger because of the reaccretion of gas heated by supernovae, which has lower entropies and shorter cooling times than the gas heated by accretion shocks. We have compared the results of NIHAO to predictions of the GalICS 2.1 semianalytic model of galaxy formation. The shock-stability criterion implemented in GalICS 2.1 successfully reproduces the transition from cold- to hot-mode accretion.

preprint2022arXiv

NIHAO XXVII: Crossing the green valley

The transition of high-mass galaxies from being blue and star forming to being red and dead is a crucial step in galaxy evolution, yet not fully understood. In this work, we use the NIHAO suite of galaxy simulations to investigate the relation between the transition time through the green valley and other galaxy properties. The typical green valley crossing time of our galaxies is approximately 400 Myr, somewhat shorter than observational estimates. The crossing of the green valley is triggered by the onset of AGN feedback and the subsequent shut down of star formation. Interestingly the time spent in the green valley is not related to any other galaxy properties, such as stellar age or metallicity, or the time at which the star formation quenching takes place. The crossing time is set by two main contributions: the ageing of the current stellar population and the residual star formation in the green valley. These effects are of comparable magnitude, while major and minor mergers have a negligible contribution. Most interestingly, we find the time that a galaxy spends to travel through the green valley is twice the $e$-folding time of the star formation quenching. This result is stable against galaxy properties and the exact numerical implementation of AGN feedback in the simulation. Assuming a typical crossing time of about one Gyr inferred from observations, our results imply that any mechanism or process aiming to quench star formation, must do it on a typical timescale of 500 Myr.

preprint2022arXiv

NIHAO XXVIII: Collateral effects of AGN on dark matter concentration and stellar kinematics

Although active galactic nuclei (AGN) feedback is required in simulations of galaxies to regulate star formation, further downstream effects on the dark matter distribution of the halo and stellar kinematics of the central galaxy can be expected. We combine simulations of galaxies with and without AGN physics from the Numerical Investigation of a Hundred Astrophysical Objects (NIHAO) to investigate the effect of AGN on the dark matter profile and central stellar rotation of the host galaxies. Specifically, we study how the concentration-halo mass ($c-M$) relation and the stellar spin parameter ($λ_R$) are affected by AGN feedback. We find that AGN physics is crucial to reduce the central density of simulated massive ($\gtrsim 10^{12}$ M$_\odot$) galaxies and bring their concentration to agreement with results from the Spitzer Photometry & Accurate Rotation Curves (SPARC) sample. Similarly, AGN feedback has a key role in reproducing the dichotomy between slow and fast rotators as observed by the ATLAS$^{3\text{D}}$ survey. Without star formation suppression due to AGN feedback, the number of fast rotators strongly exceeds the observational constraints. Our study shows that there are several collateral effects that support the importance of AGN feedback in galaxy formation, and these effects can be used to constrain its implementation in numerical simulations.

preprint2022arXiv

The diversity of spiral galaxies explained

An extensive catalog of spatially-resolved galaxy rotation curves and multi-band optical light profiles for 1752 observed spiral galaxies is assembled to explore the drivers of diversity in galaxy structural parameters, rotation curve shapes, and stellar mass profiles. Similar data were extracted from the NIHAO galaxy simulations to identify any differences between observations and simulations. Several parameters, including the inner slope "S" of a rotation curve (RC), were tested for diversity. Two distinct populations are found in observed and simulated galaxies; (i) blue, low mass spirals with stellar mass M* < 10^9.3 Msol and roughly constant "S", and (ii) redder, more massive and more diverse spirals with rapidly increasing "S". In all cases, the value of "S" seems equally contributed by the baryonic and non-baryonic (dark) matter. Diversity is shown to increase mildly with mass. Numerical simulations reproduce well most baryon-dominated galaxy parameter distributions, such as the inner stellar mass profile slope and baryonic scaling relations, but they struggle to match the full diversity of observed galaxy rotation curves (through "S") and most dark-matter-dominated parameters. To reproduce observations, the error broadening of the simulation's intrinsic spread of RC metrics would have to be tripled. The differences in various projections of observed and simulated scaling relations may reflect limitations of current sub-grid physics models to fully capture the complex nature of galaxies. For instance, AGNs are shown to have a significant effect on the shapes of simulated RCs. The inclusion of AGN feedback brings simulated and observed inner RC shapes into closer agreement.

preprint2022arXiv

Using Artificial Intelligence and real galaxy images to constrain parameters in galaxy formation simulations

Cosmological galaxy formation simulations are still limited by their spatial/mass resolution and cannot model from first principles some of the processes, like star formation, that are key in driving galaxy evolution. As a consequence they still rely on a set of 'effective parameters' that try to capture the scales and the physical processes that cannot be directly resolved in the simulation. In this study we show that it is possible to use Machine Learning techniques applied to real and simulated images of galaxies to discriminate between different values of these parameters by making use of the full information content of an astronomical image instead of collapsing it into a limited set of values like size, or stellar/ gas masses. In this work we apply our method to the NIHAO simulations and the THINGS and VLA-ANGST observations of HI maps in nearby galaxies to test the ability of different values of the star formation density threshold $n$ to reproduce observed HI maps. We show that observations indicate the need for a high value of $n \gtrsim 80$ ,cm$^{-3}$ (although the exact numerical value is model-dependent), which has important consequences for the dark matter distribution in galaxies. Our study shows that with innovative methods it is possible to take full advantage of the information content of galaxy images and compare simulations and observations in an interpretable, non-parametric and quantitative manner.

preprint2020arXiv

Abundance matching tested on small scales with galaxy dynamics

We present a comprehensive test of the relation between stellar and total mass in galaxies as predicted by popular models based on abundance matching (AM) techniques. We use the ``Spectroscopy and H-band Imaging of Virgo cluster galaxies'' (SHIVir) survey with photometric and dynamical profiles for 190 Virgo cluster galaxies to establish a relation between the stellar and dynamical masses measured within the isophotal radius $r_{23.5}$. Various dark matter and galaxy scaling relations are combined with results from the NIHAO suite of hydrodynamical simulations to recast AM predictions in terms of these observed quantities. Our results are quite insensitive to the exact choice of dark matter profile and halo response to baryon collapse. We find that theoretical models reproduce the slope and normalization of the observed stellar-to-halo mass relation (SHMR) over more than three orders of magnitude in stellar mass $(10^8 < M_*/M_{\odot} < 2 \times 10^{11})$. However, the scatter of the observed SHMR exceeds that of AM predictions by a factor $\sim$5. For systems with stellar masses exceeding $5 \times 10^{10}~M_{\odot}$, AM overpredicts the observed stellar masses for a given dynamical mass. The latter offset may support previous indications of a different stellar initial mass function in these massive galaxies. Overall our results support the validity of AM predictions on a wide dynamical range.

preprint2020arXiv

Addressing via N-body simulations the distribution of the satellite tidal debris in the Milky Way environment

We study the distribution of the Milky Way satellites stellar and dark matter debris. For the first time we address the question of the tidal disruption of satellites in simulations by utilising simultaneously a) a realistic set of orbits extracted from cosmological simulations, b) a three component host galaxy with live halo, disc and bulge components, and c) satellites from hydrodynamical simulations. We analyse the statistical properties of the satellite debris of all massive galaxies reaching the inner Milky Way on a timescale of 2 Gyr. Up to 80$\%$ of the dark matter is stripped from the satellites, while this happens for up to 30$\%$ of their stars. The stellar debris ends mostly in the inner Milky Way halo, whereas the dark matter debris shows a flat mass distribution over the full main halo. The dark matter debris follows a density profile with inner power law index $α_{\rm DM}=-0.66$ and outer index $β_{\rm DM}=2.94$, while for stars $α_{*}=-0.44$ and $β_{*}=6.17$. In the inner 25 kpc, the distribution of the stellar debris is flatter than that of the dark matter debris and the orientations of their short axes differ significantly. Changing the orientation of the stellar disc by 90$^{\rm{o}}$ has only a minor impact on the distribution of the satellite debris. Our results indicate that the dark matter is more easily stripped than stars from the Milky Way satellites. The structure of the debris is dominated by the satellite orbital properties. The radial profiles, the flattening and the orientation of the stellar and dark matter debris are significantly different, which prevents the prediction of the dark matter distribution from the observed stellar component.

preprint2020arXiv

NIHAO XXIV: Rotation or pressure supported systems? Simulated Ultra Diffuse Galaxies show a broad distribution in their stellar kinematics

In recent years a new window on galaxy evolution opened, thanks to the increasing discovery of galaxies with a low surface brightness, such as Ultra Diffuse Galaxies (UDGs). The formation mechanism of these systems is still a much debated question, and so are their kinematical properties. In this work, we address this topic by analyzing the stellar kinematics of isolated UDGs formed in the hydrodynamical simulation suite NIHAO. We construct projected line-of-sight velocity and velocity dispersion maps to compute the projected specific angular momentum, $λ_{\rm R}$, to characterize the kinematical support of the stars in these galaxies. We found that UDGs cover a broad distribution, ranging from dispersion to rotation supported galaxies, with similar abundances in both regimes. The degree of rotation support of simulated UDGs correlates with several properties such as galaxy morphology, higher HI fractions and larger effective radii with respect to the dispersion supported group, while the dark matter halo spin and mass accretion history are similar amongst the two populations. We demonstrate that the alignment of the infalling baryons into the protogalaxy at early $z$ is the principal driver of the $z$=0 stellar kinematic state: pressure supported isolated UDGs form via mis-aligned gas accretion while rotation supported ones build-up their baryons in an ordered manner. Accounting for random inclination effects, we predict that a comprehensive survey will find nearly half of field UDGs to have rotationally supported stellar disks.

preprint2020arXiv

NIHAO-UHD: The properties of MW-like stellar disks in high resolution cosmological simulations

Simulating thin and extended galactic disks has long been a challenge in computational astrophysics. We introduce the NIHAO-UHD suite of cosmological hydrodynamical simulations of Milky Way mass galaxies and study stellar disk properties such as stellar mass, size and rotation velocity which agree well with observations of the Milky Way and local galaxies. In particular, the simulations reproduce the age-velocity dispersion relation and a multi-component stellar disk as observed for the Milky Way. Half of our galaxies show a double exponential vertical profile, while the others are well described by a single exponential model which we link to the disk merger history. In all cases, mono-age populations follow a single exponential whose scale height varies monotonically with stellar age and radius. The scale length decreases with stellar age while the scale height increases. The general structure of the stellar disks is already set at time of birth as a result of the inside-out and upside-down formation. Subsequent evolution modifies this structure by increasing both the scale length and height of all mono-age populations. Thus, our results put tight constraints on how much dynamical memory stellar disks can retain over cosmological timescales. Our simulations demonstrate that it is possible to form thin galactic disks in cosmological simulations provided there are no significant stellar mergers at low redshifts. Most of the stellar mass is formed in-situ with only a few percent ($\lesssim5\%$) brought in by merging satellites at early times. Redshift zero snapshots and halo catalogues are publicly available.

preprint2020arXiv

NIHAO-XXIII: Dark Matter density shaped by Black Hole feedback

We present a systematic analysis of the reaction of dark matter distribution to galaxy formation across more than eight orders of magnitude in stellar mass. We extend the previous work presented in the NIHAO-IV paper (Tollet et al.) by adding 46 new high resolution simulations of massive galaxies performed with the inclusion of Black Hole feedback. We show that outflows generated by the AGN are able to partially counteract the dark matter contraction due to the large central stellar component in massive haloes. The net effect is to relax the central dark matter distribution that moves to a less cuspy density profiles at halo masses larger than $\approx 3 \times 10^{12}$ M$_{\odot}$. The scatter around the mean value of the density profile slope ($α$) is fairly constant ($Δα\approx 0.3$), with the exception of galaxies with halo masses around $10^{12}$ M$_{\odot}$, at the transition from stellar to AGN feedback dominated systems, where the scatter increases by almost a factor three. We provide useful fitting formulas for the slope of the dark matter density profiles at few percent of the virial radius for the whole stellar mass range: $10^5-10^{12}$ M$_{\odot}$ ($2 \times 10^9-5 \times 10^{13}$ M$_{\odot}$ in halo mass).

preprint2016arXiv

NIHAO IX: the role of gas inflows and outflows in driving the contraction and expansion of cold dark matter haloes

We use ~100 cosmological galaxy formation zoom-in simulations using the smoothed particle hydrodynamics code {\sc gasoline} to study the effect of baryonic processes on the mass profiles of cold dark matter haloes. The haloes in our study range from dwarf (M_{200}~10^{10}Msun) to Milky Way (M_{200}~10^{12}Msun) masses. Our simulations exhibit a wide range of halo responses, primarily varying with mass, from expansion to contraction, with up to factor ~10 changes in the enclosed dark matter mass at one per cent of the virial radius. Confirming previous studies, the halo response is correlated with the integrated efficiency of star formation: e_SF=(M_{star}/M_{200})/(Ω_b/Ω_m). In addition we report a new correlation with the compactness of the stellar system: e_R=r_{1/2}/R_{200}. We provide an analytic formula depending on e_SF and e_R for the response of cold dark matter haloes to baryonic processes. An observationally testable prediction is that, at fixed mass, larger galaxies experience more halo expansion, while the smaller galaxies more halo contraction. This diversity of dark halo response is captured by a toy model consisting of cycles of adiabatic inflow (causing contraction) and impulsive gas outflow (causing expansion). For net outflow, or equal inflow and outflow fractions, f, the overall effect is expansion, with more expansion with larger f. For net inflow, contraction occurs for small f (large radii), while expansion occurs for large f (small radii), recovering the phenomenology seen in our simulations. These regularities in the galaxy formation process provide a step towards a fully predictive model for the structure of cold dark matter haloes.

preprint2016arXiv

NIHAO project II: Halo shape, phase-space density and velocity distribution of dark matter in galaxy formation simulations

We use the NIHAO (Numerical Investigation of Hundred Astrophysical Objects) cosmological simulations to study the effects of galaxy formation on key properties of dark matter (DM) haloes. NIHAO consists of $\simeq 90$ high-resolution SPH simulations that include (metal-line) cooling, star formation, and feedback from massive stars and SuperNovae, and cover a wide stellar and halo mass range: $10^6 < M_* / M_{\odot} < 10^{11}$ ( $10^{9.5} < M_{\rm halo} / M_{\odot} < 10^{12.5}$). When compared to DM-only simulations, the NIHAO haloes have similar shapes at the virial radius, R_{\rm vir}, but are substantially rounder inside $\simeq 0.1R_{\rm vir}$. In NIHAO simulations $c/a$ increases with halo mass and integrated star formation efficiency, reaching $\sim 0.8$ at the Milky Way mass (compared to 0.5 in DM-only), providing a plausible solution to the long-standing conflict between observations and DM-only simulations. The radial profile of the phase-space $Q$ parameter ($ρ/σ^3$) is best fit with a single power law in DM-only simulations, but shows a flattening within $\simeq 0.1R_{\rm vir}$ for NIHAO for total masses $M>10^{11} M_{\odot}$. Finally, the global velocity distribution of DM is similar in both DM-only and NIHAO simulations, but in the solar neighborhood, NIHAO galaxies deviate substantially from Maxwellian. The distribution is more symmetric, roughly Gaussian, with a peak that shifts to higher velocities for Milky Way mass haloes. We provide the distribution parameters which can be used for predictions for direct DM detection experiments. Our results underline the ability of the galaxy formation processes to modify the properties of dark matter haloes.

preprint2016arXiv

NIHAO VIII: Circum-galactic medium and outflows - The puzzles of HI and OVI gas distributions

We study the hot and cold circum-galactic medium (CGM) of 86 galaxies of the cosmological, hydrodynamical simulation suite NIHAO. NIHAO allows a study of how the $z=0$ CGM varies across 5 orders of magnitude of stellar mass using OVI and HI as proxies for hot and cold gas. The cool HI covering fraction and column density profiles match observations well, particularly in the inner CGM. OVI shows increasing column densities with mass, a trend seemingly echoed in the observations. As in multiple previous simulations, the OVI column densities in simulations are lower than observed and optically thick HI does not extend as far out as in observations. We take a look at the collisional ionisation fraction of OVI as a function of halo mass. We make observable predictions of the bipolarity of outflows and their effect on the general shape of the CGM. Bipolar outflows can be seen out to around 40 kpc in intermediate and low mass halos ($M_{\mathrm{halo}}<10^{11}M_{\mathrm{sun}}$), but outside that radius, the CGM is too well mixed to detect an elongated shape. Larger halos have extended gas discs beyond the stellar disc that dominate the shape of the inner CGM. The simulated CGM is remarkably spherical even in low mass simulations. The chemical enrichment of both halo and disc gas follow expected increasing trends as a function of halo mass that are well fit with power laws. These relations can be used in non-hydrodynamic models, such as semi-analytic models.

preprint2016arXiv

NIHAO X: Reconciling the local galaxy velocity function with Cold Dark Matter via mock HI observations

We used 87 high resolution hydrodynamical cosmological simulations from the NIHAO suite to investigate the relation between the maximum circular velocity (Vmax) of a dark matter halo in a collisionless simulation and the velocity width of the HI gas in the same halo in the hydrodynamical simulation. These two quantities are normally used to compare theoretical and observational velocity functions and have led to a possible discrepancy between observations and predictions based on the Cold Dark Matter (CDM) model. We show that below 100 km/s, there is clear bias between HI based velocities and Vmax, that leads to an underestimation of the actual circular velocity of the halo. When this bias is taken into account the CDM model has no trouble in reproducing the observed velocity function and no lack of low velocity galaxies is actually present. Our simulations also reproduce the linewidth - stellar mass (Tully-Fisher) relation and HI sizes, indicating that the HI gas in our simulations is as extended as observed. The physical reason for the lower than expected linewidths is that, in contrast to high mass galaxies, low mass galaxies no longer have extended thin HI rotating disks, as is commonly assumed.

preprint2016arXiv

NIHAO XI: Formation of Ultra-Diffuse Galaxies by outflows

We address the origin of Ultra-Diffuse Galaxies (UDGs), which have stellar masses typical of dwarf galaxies but effective radii of Milky Way-sized objects. Their formation mechanism, and whether they are failed $\rm L_{\star}$ galaxies or diffuse dwarfs, are challenging issues. Using zoom-in cosmological simulations from the NIHAO project, we show that UDG analogues form naturally in medium-mass haloes due to episodes of gas outflows associated with star formation. The simulated UDGs live in isolated haloes of masses $10^{10-11}\rm M_{\odot}$, have stellar masses of $10^{7-8.5}\rm M_{\odot}$, effective radii larger than 1 kpc and dark matter cores. They show a broad range of colors, an average Sérsic index of 0.83, a typical distribution of halo spin and concentration, and a non-negligible HI gas mass of $10^{7-9}\rm M_{\odot}$, which correlates with the extent of the galaxy. Gas availability is crucial to the internal processes that form UDGs: feedback driven gas outflows, and subsequent dark matter and stellar expansion, are the key to reproduce faint, yet unusually extended, galaxies. This scenario implies that UDGs represent a dwarf population of low surface brightness galaxies and should exist in the field. The largest isolated UDGs should contain more HI gas than less extended dwarfs of similar $\rm M_{\star}$.

preprint2016arXiv

Simulated ΛCDM analogues of the thin Plane of Satellites around the Andromeda galaxy are not kinematically coherent structures

A large fraction of the dwarf satellites orbiting the Andromeda galaxy are surprisingly aligned in a thin, extended and apparently kinematically coherent planar structure. Such a structure is not easily found in simulations based on the Cold Dark Matter model (ΛCDM). Using 21 high resolution cosmological simulations we analyse the kinematics of planes of satellites similar to the one around Andromeda. We find good agreement when co-rotation is characterized by the line-of-sight velocity. At the same time, when co-rotation is inferred by the angular momenta of the satellites, the planes are in agreement with the plane around our Galaxy. We find such planes to be common in our high concentration haloes. The number of co-rotating satellites obtained from the sign of the line-of-sight velocity shows large variations depending on the viewing angle and is consistent with that obtained from a sample with random velocities. We find that the clustering of angular momentum vectors of the satellites in the plane is a better measure of the kinematic coherence. Thus we conclude that the line-of- sight velocity is not well suited as a proxy for the kinematical coherence of the plane. Analysis of the kinematics of our planes shows a fraction of $\sim$30% chance aligned satellites. Tracking the satellites in the plane back in time reveals that these planes are a transient feature and not kinematically coherent as would appear at first sight. Thus we expect some of the satellites in the plane around Andromeda to have high velocities perpendicular to the plane.

preprint2015arXiv

A New Channel for the Formation of Kinematically Decoupled Cores in Early-type galaxies

We present the formation of a Kinematically Decoupled Core (KDC) in an elliptical galaxy, resulting from a major merger simulation of two disk galaxies. We show that although the two progenitor galaxies are initially following a prograde orbit, strong reactive forces during the merger can cause a short-lived change of their orbital spin; the two progenitors follow a retrograde orbit right before their final coalescence. This results in a central kinematic decoupling and the formation of a large-scale (~2 kpc radius) counter-rotating core (CRC) at the center of the final elliptical-like merger remnant (M*=1.3x10^11 Msun), while its outer parts keep the rotation direction of the initial orbital spin. The stellar velocity dispersion distribution of the merger remnant galaxy exhibits two symmetrical off-centered peaks, comparable to the observed "2-sigma galaxies". The KDC/CRC consists mainly of old, pre-merger population stars (older than 5 Gyr), remaining prominent in the center of the galaxy for more than 2 Gyr after the coalescence of its progenitors. Its properties are consistent with KDCs observed in massive elliptical galaxies. This new channel for the formation of KDCs from prograde mergers is in addition to previously known formation scenarios from retrograde mergers and can help towards explaining the substantial fraction of KDCs observed in early-type galaxies.

preprint2015arXiv

Effects of Coupled Dark Energy on the Milky Way and its Satellites

We present the first numerical simulations in coupled dark energy cosmologies with high enough resolution to investigate the effects of the coupling on galactic and sub-galactic scales. We choose two constant couplings and a time-varying coupling function and we run simulations of three Milky-Way-size halos ($\sim$10$^{12}$M$_{\odot}$), a lower mass halo (6$\times$10$^{11}$M$_{\odot}$) and a dwarf galaxy halo (5$\times$10$^{9}$M$_{\odot}$). We resolve each halo with several millions dark matter particles. On all scales the coupling causes lower halo concentrations and a reduced number of substructures with respect to LCDM. We show that the reduced concentrations are not due to different formation times, but they are related to the extra terms that appear in the equations describing the gravitational dynamics. On the scale of the Milky Way satellites, we show that the lower concentrations can help in reconciling observed and simulated rotation curves, but the coupling values necessary to have a significant difference from LCDM are outside the current observational constraints. On the other hand, if other modifications to the standard model allowing a higher coupling (e.g. massive neutrinos) are considered, coupled dark energy can become an interesting scenario to alleviate the small-scale issues of the LCDM model.

preprint2015arXiv

Evidence for Early Filamentary Accretion from the Andromeda Galaxy's Thin Plane of Satellites

Recently it has been shown that a large fraction of the dwarf satellite galaxies orbiting the Andromeda galaxy are surprisingly aligned in a thin, extended and kinematically coherent planar structure. The presence of such a structure seems to challenge the current Cold Dark Matter paradigm of structure formation, which predicts a more uniform distribution of satellites around central objects. We show that it is possible to obtain a thin, extended, rotating plane of satellites resembling the one in Andromeda in cosmological collisionless simulations based on the Cold Dark Matter model. Our new high resolution simulations show a correlation between the formation time of the dark matter halo and the thickness of the plane of satellites. Our simulations have a high incidence of satellite planes as thin, extended, and as rich as the one in Andromeda and with a very coherent kinematic structure when we select high concentration/early forming halos. By tracking the formation of the satellites in the plane we show that they have been mainly accreted onto the main object along thin dark matter filaments at high redshift. Our results show that the presence of a thin, extended, rotating plane of satellites is not a challenge for the Cold Dark Matter paradigm, but actually supports one of the predictions of this paradigm related to the presence of filaments of dark matter around galaxies at high redshift.

preprint2015arXiv

NIHAO project I: Reproducing the inefficiency of galaxy formation across cosmic time with a large sample of cosmological hydrodynamical simulations

We introduce project NIHAO (Numerical Investigation of a Hundred Astrophysical Objects), a set of 100 cosmological zoom-in hydrodynamical simulations performed using the gasoline code, with an improved implementation of the SPH algorithm. The haloes in our study range from dwarf to Milky Way masses, and represent an unbiased sampling of merger histories, concentrations and spin parameters. The particle masses and force softenings are chosen to resolve the mass profile to below 1% of the virial radius at all masses, ensuring that galaxy half-light radii are well resolved. Using the same treatment of star formation and stellar feedback for every object, the simulated galaxies reproduce the observed inefficiency of galaxy formation across cosmic time as expressed through the stellar mass vs halo mass relation, and the star formation rate vs stellar mass relation. We thus conclude that stellar feedback is the chief piece of physics required to limit the efficiency of star formation in galaxies less massive than the Milky Way.

preprint2015arXiv

NIHAO V: Too big doesn't fail -- reconciling the conflict between LCDM predictions and the circular velocities of nearby field galaxies

We compare the half-light circular velocities, V_{1/2}, of dwarf galaxies in the Local Group to the predicted circular velocity curves of galaxies in the NIHAO suite of LCDM simulations. We use a subset of 34 simulations in which the central galaxy has a stellar luminosity in the range 0.5 x 10^5 < L_V < 2 x 10^8 L_{sun}. The NIHAO galaxy simulations reproduce the relation between stellar mass and halo mass from abundance matching, as well as the observed half-light size vs luminosity relation. The corresponding dissipationless simulations over-predict the V_{1/2}, recovering the problem known as too big to fail (TBTF). By contrast, the NIHAO simulations have expanded dark matter haloes, and provide an excellent match to the distribution of V_{1/2} for galaxies with L_V > 2 x 10^6 L_{sun}. For lower luminosities our simulations predict very little halo response, and tend to over predict the observed circular velocities. In the context of LCDM, this could signal the increased stochasticity of star formation in haloes below M_{halo} \sim 10^{10} M_{sun}, or the role of environmental effects. Thus, haloes that are "too big to fail", do not fail LCDM, but haloes that are "too small to pass" (the galaxy formation threshold) provide a future test of LCDM.

preprint2015arXiv

Strongly Coupled Dark Energy Cosmologies: preserving LCDM success and easing low scale problems I - Linear theory revisited

In this first paper we discuss the linear theory and the background evolution of a new class of models we dub SCDEW: Strongly Coupled DE, plus WDM. In these models, WDM dominates today's matter density; like baryons, WDM is uncoupled. Dark Energy is a scalar field $Φ$; its coupling to ancillary CDM, whose today's density is $\ll 1\, \%$, is an essential model feature. Such coupling, in fact, allows the formation of cosmic structures, in spite of very low WDM particle masses ($\sim 100$ eV). SCDEW models yields Cosmic Microwave Background and linear Large Scale features substantially undistinguishable from $Λ$CDM, but thanks to the very low WDM masses they strongly alleviate $Λ$CDM issues on small scales, as confirmed via numerical simulations in the II associated paper. Moreover SCDEW cosmologies significantly ease the coincidence and fine tuning problems of $Λ$CDM and, by using a field theory approach, we also outline possible links with inflationary models. We also discuss a possible fading of the coupling at low redshifts which prevents non linearities on the CDM component to cause computational problems. The (possible) low-$z$ coupling suppression, its mechanism, and its consequences are however still open questions -not necessarily problems- for SCDEW models. The coupling intensity and the WDM particle mass, although being extra parameters in respect to $Λ$CDM, are found to be substantially constrained a priori so that, if SCDEW is the underlying cosmology, we expect most data to fit also $Λ$CDM predictions.

preprint2015arXiv

Strongly Coupled Dark Energy Cosmologies: preserving LCDM success and easing low scale problems II - Cosmological simulations

In this second paper we present the first Nbody cosmological simulations of strongly coupled Dark Energy models (SCDEW), a class of models that alleviates theoretical issues related to the nature of dark energy. SCDEW models assume a strong coupling between Dark Energy (DE) and an ancillary Cold Dark Matter (CDM) component together with the presence of an uncoupled Warm Dark Matter component. The strong coupling between CDM and DE allows us to preserve small scale fluctuations even if the warm particle is quite light ($\approx 100$ eV). Our large scale simulations show that, for $10^{11}<M/M_\odot<10^{14}$, SCDEW haloes exhibit a number density and distribution similar to a standard Lambda Cold Dark Matter (LCDM) model, even though they have lower concentration parameters. High resolution simulation of a galactic halo ($M\sim 10^{12} M_{\odot} $) shows $\sim 60\%$ less substructures than its LCDM counterpart, but the same cuspy density profile. On the scale of galactic satellites ($M\sim 10^{9} M_{\odot}$) SCDEW haloes dramatically differ from LCDM. Due to the high thermal velocities of the WDM component they are almost devoid of any substructures and present strongly cored dark matter density profiles. These density cores extend for several hundreds of parsecs, in very good agreement with Milky Way satellites observations. Strongly coupled models, thanks to their ability to match observations on both large and small scales might represent a valid alternative to a simple LCDM model.

preprint2015arXiv

The response of dark matter haloes to elliptical galaxy formation: a new test for quenching scenarios

We use cosmological hydrodynamical zoom-in simulations with the SPH code gasoline of four haloes of mass M_{200} \sim 10^{13}\Msun to study the response of the dark matter to elliptical galaxy formation. Our simulations include metallicity dependent gas cooling, star formation, and feedback from massive stars and supernovae, but not active galactic nuclei (AGN). At z=2 the progenitor galaxies have stellar to halo mass ratios consistent with halo abundance matching, assuming a Salpeter initial mass function. However by z=0 the standard runs suffer from the well known overcooling problem, overpredicting the stellar masses by a factor of > 4. To mimic a suppressive halo quenching scenario, in our forced quenching (FQ) simulations, cooling and star formation are switched off at z=2. The resulting z=0 galaxies have stellar masses, sizes and circular velocities close to what is observed. Relative to the control simulations, the dark matter haloes in the FQ simulations have contracted, with central dark matter density slopes d\logρ/d\log r \sim -1.5, showing that dry merging alone is unable to fully reverse the contraction that occurs at z>2. Simulations in the literature with AGN feedback however, have found expansion or no net change in the dark matter halo. Thus the response of the dark matter halo to galaxy formation may provide a new test to distinguish between ejective and suppressive quenching mechanisms.

preprint2015arXiv

The Star Formation and AGN luminosity relation: Predictions from a semi-analytical model

In a Universe where AGN feedback regulates star formation in massive galaxies, a strong correlation between these two quantities is expected. If the gas causing star formation is also responsible for feeding the central black hole, then a positive correlation is expected. If powerful AGNs are responsible for the star formation quenching, then a negative correlation is expected. Observations so far have mainly found a mild correlation or no correlation at all (i.e. a flat relation between star formation rate (SFR) and AGN luminosity), raising questions about the whole paradigm of "AGN feedback". In this paper, we report the predictions of the GALFORM semi-analytical model, which has a very strong coupling between AGN activity and quenching of star formation. The predicted SFR-AGN luminosity correlation appears negative in the low AGN luminosity regime, where AGN feedback acts, but becomes strongly positive in the regime of the brightest AGN. Our predictions reproduce reasonably well recent observations by Rosario et al., yet there is some discrepancy in the normalisation of the correlation at low luminosities and high redshifts. Though this regime could be strongly influenced by observational biases, we argue that the disagreement could be ascribed to the fact that GALFORM neglects AGN variability effects. Interestingly, the galaxies that dominate the regime where the observations imply a weak correlation are massive early-type galaxies that are subject to AGN feedback. Nevertheless, these galaxies retain high enough molecular hydrogen contents to maintain relatively high star formation rates and strong infrared emission.

preprint2014arXiv

A mass-dependent density profile for dark matter haloes including the influence of galaxy formation

We introduce a mass dependent density profile to describe the distribution of dark matter within galaxies, which takes into account the stellar-to-halo mass dependence of the response of dark matter to baryonic processes. The study is based on the analysis of hydrodynamically simulated galaxies from dwarf to Milky Way mass, drawn from the MaGICC project, which have been shown to match a wide range of disk scaling relationships. We find that the best fit parameters of a generic double power-law density profile vary in a systematic manner that depends on the stellar-to-halo mass ratio of each galaxy. Thus, the quantity Mstar/Mhalo constrains the inner ($γ$) and outer ($β$) slopes of dark matter density, and the sharpness of transition between the slopes($α$), reducing the number of free parameters of the model to two. Due to the tight relation between stellar mass and halo mass, either of these quantities is sufficient to describe the dark matter halo profile including the effects of baryons. The concentration of the haloes in the hydrodynamical simulations is consistent with N-body expectations up to Milky Way mass galaxies, at which mass the haloes become twice as concentrated as compared with pure dark matter runs. This mass dependent density profile can be directly applied to rotation curve data of observed galaxies and to semi analytic galaxy formation models as a significant improvement over the commonly used NFW profile.

preprint2014arXiv

Cold dark matter haloes in the Planck era: evolution of structural parameters for Einasto and NFW profiles

We present the evolution of the structure of relaxed cold dark matter haloes in the cosmology from the Planck satellite. Our simulations cover 5 decades in halo mass, from dwarf galaxies to galaxy clusters. Due to the increased matter density and power spectrum normalization the concentration mass relation in the Planck cosmology has a 20 percent higher normalization at redshift z=0 compared to WMAP cosmology. We confirm that CDM haloes are better described by the Einasto profile; for example, at scales near galaxy half-light radii CDM haloes have significantly steeper density profiles than implied by NFW fits. There is a scatter of 0.2 dex in the Einasto shape parameter at fixed halo mass, adding further to the diversity of CDM halo profiles. The evolution of the concentration mass relation in our simulations is not reproduced by any of the analytic models in the literature. We thus provide a simple fitting formula that accurately describes the evolution between redshifts z=5 to z=0 for both NFW and Einasto fits. Finally, the observed concentrations and halo masses of spiral galaxies, groups and clusters of galaxies at low redshifts are in good agreement with our simulations, suggesting only mild halo response to galaxy formation on these scales.

preprint2014arXiv

Cold stream stability during minor mergers

We use high-resolution Eulerian simulations to study the stability of cold gas flows in a galaxy size dark matter halo (10^12 Msun) at redshift z=2. Our simulations show that a cold stream penetrating a hot gaseous halo is stable against thermal convection and Kelvin-Helmholtz instability. We then investigate the effect of a satellite orbiting the main halo in the plane of the stream. The satellite is able to perturb the stream and to inhibit cold gas accretion towards the center of the halo for 0.5 Gyr. However, if the supply of cold gas at large distances is kept constant, the cold stream is able to re-establish itself after 0.3 Gyr. We conclude that cold streams are very stable against a large variety of internal and external perturbations.

preprint2014arXiv

Dark MaGICC: the effect of Dark Energy on galaxy formation. Cosmology does matter

We present the Dark MaGICC project, which aims to investigate the effect of Dark Energy (DE) modeling on galaxy formation via hydrodynamical cosmological simulations. Dark MaGICC includes four dynamical Dark Energy scenarios with time varying equations of state, one with a self-interacting Ratra-Peebles model. In each scenario we simulate three galaxies with high resolution using smoothed particle hydrodynamics (SPH). The baryonic physics model is the same used in the Making Galaxies in a Cosmological Context (MaGICC) project, and we varied only the background cosmology. We find that the Dark Energy parameterization has a surprisingly important impact on galaxy evolution and on structural properties of galaxies at z=0, in striking contrast with predictions from pure Nbody simulations. The different background evolutions can (depending on the behavior of the DE equation of state) either enhance or quench star formation with respect to a LCDM model, at a level similar to the variation of the stellar feedback parameterization, with strong effects on the final galaxy rotation curves. While overall stellar feedback is still the driving force in shaping galaxies, we show that the effect of the Dark Energy parameterization plays a larger role than previously thought, especially at lower redshifts. For this reason, the influence of Dark Energy parametrization on galaxy formation must be taken into account, especially in the era of precision cosmology.

preprint2014arXiv

The Distribution of Satellites Around Central Galaxies in a Cosmological Hydrodynamical Simulation

Observations have shown that the spatial distribution of satellite galaxies is not random, but rather is aligned with the major axes of central galaxies (CGs). The strength of the alignment is dependent on the properties of both the satellites and centrals. Theoretical studies using dissipationless N-body simulations are limited by their inability to directly predict the shape of CGs. Using hydrodynamical simulations including gas cooling, star formation, and feedback, we carry out a study of galaxy alignment and its dependence on the galaxy properties predicted directly from the simulations.We found that the observed alignment signal is well produced, as is the color dependence: red satellites and red centrals both show stronger alignments than their blue counterparts. The reason for the stronger alignment of red satellites is that most of them stay in the inner region of the dark matter halo where the shape of the CG better traces the dark matter distribution. The dependence of alignment on the color of CGs arises from the halo mass dependence, since the alignment between the shape of the central stellar component and the inner halo increases with halo mass. We also find that the alignment of satellites is most strongly dependent on their metallicity, suggesting that the metallicity of satellites, rather than color, is a better tracer of galaxy alignment on small scales. This could be tested in future observational studies.

preprint2013arXiv

MaGICC-WDM: the effects of warm dark matter in hydrodynamical simulations of disc galaxy formation

We study the effect of warm dark matter (WDM) on hydrodynamic simulations of galaxy formation as part of the Making Galaxies in a Cosmological Context (MaGICC) project. We simulate three different galaxies using three WDM candidates of 1, 2 and 5 keV and compare results with pure cold dark matter simulations. WDM slightly reduces star formation and produces less centrally concentrated stellar profiles. These effects are most evident for the 1 keV candidate but almost disappear for $m_{\mathrm{WDM}}>2$ keV. All simulations form similar stellar discs independent of WDM particle mass. In particular, the disc scale length does not change when WDM is considered. The reduced amount of star formation in the case of 1 keV particles is due to the effects of WDM on merging satellites which are on average less concentrated and less gas rich. The altered satellites cause a reduced starburst during mergers because they trigger weaker disc instabilities in the main galaxy. Nevertheless we show that disc galaxy evolution is much more sensitive to stellar feedback than it is to WDM candidate mass. Overall we find that WDM, especially when restricted to current observational constraints ($m_{\mathrm{WDM}}>2$ keV), has a minor impact on disc galaxy formation.

preprint2013arXiv

The effect of Warm Dark Matter on galaxy properties: constraints from the stellar mass function and the Tully-Fisher relation

In this paper we combine high resolution N-body simulations with a semi analytical model of galaxy formation to study the effects of a possible Warm Dark Matter (WDM) component on the observable properties of galaxies. We compare three WDM models with a dark matter mass of 0.5, 0.75 and 2.0 keV, with the standard Cold Dark Matter case. For a fixed set of parameters describing the baryonic physics the WDM models predict less galaxies at low (stellar) masses, as expected due to the suppression of power on small scales, while no substantial difference is found at the high mass end. However these differences in the stellar mass function, vanish when different set of parameters are used to describe the (largely unknown) galaxy formation processes. We show that is possible to break this degeneracy between DM properties and the parameterization of baryonic physics by combining observations on the stellar mass function with the Tully-Fisher relation (the relation between stellar mass and the rotation velocity at large galactic radii as probed by resolved HI rotation curves). WDM models with a too warm candidate (m<0.75 keV) cannot simultaneously reproduce the stellar mass function and the Tully-Fisher relation. We conclude that accurate measurements of the galaxy stellar mass function and the link between galaxies and dark matter haloes down to the very low-mass end can give very tight constraints on the nature of DM candidates.

preprint2013arXiv

The MaGICC volume: reproducing statistical properties of high redshift galaxies

We present a cosmological hydrodynamical simulation of a representative volume of the Universe, as part of the Making Galaxies in a Cosmological Context (MaGICC) project. MaGICC uses a thermal implementation for supernova and early stellar feedback. This work tests the feedback model at lower resolution across a range of galaxy masses, morphologies and merger histories. The simulated sample compares well with observations of high redshift galaxies ($z \ge 2$) including the stellar mass - halo mass ($M_\star - M_h$ ) relation, the Galaxy Stellar Mass Function (GSMF) at low masses ($M_\star \lt 5 \times 10^{10} M_\odot$ ) and the number density evolution of low mass galaxies. The poor match of $M_\star - M_h$ and the GSMF at high masses ($M_\star \ge 5 \times 10^{10} M_\odot$ ) indicates supernova feedback is insufficient to limit star formation in these haloes. At $z = 0$, our model produces too many stars in massive galaxies and slightly underpredicts the stellar mass around $L_\star$ mass galaxy. Altogether our results suggest that early stellar feedback, in conjunction with supernovae feedback, plays a major role in regulating the properties of low mass galaxies at high redshift.

preprint2012arXiv

An attractor for the dynamical state of the intracluster medium

Galaxy clusters provide us with important information about the cosmology of our universe. Observations of the X-ray radiation or of the SZ effect allow us to measure the density and temperature of the hot intergalactic medium between the galaxies in a cluster, which then allow us to calculate the total mass of the galaxy cluster. However, no simple connection between the density and the temperature profiles has been identified. Here we use controlled high-resolution numerical simulations to identify a relation between the density and temperature of the gas in equilibrated galaxy clusters. We demonstrate that the temperature-density relation is a real attractor, by showing that a wide range of equilibrated structures all move towards the attractor when perturbed and subsequently allowed to relax. For structures which have undergone sufficient perturbations for this connection to hold, one can therefore extract the mass profile directly from the X-ray intensity profile.

preprint2012arXiv

MAGICC haloes: confronting simulations with observations of the circumgalactic medium at z=0

We explore the circumgalactic medium (CGM) of two simulated star-forming galaxies with luminosities L ~ 0.1 and 1 L* generated using the smooth particle hydrodynamic code GASOLINE. These simulations are part of the Making Galaxies In a Cosmological Context (MAGICC) program in which the stellar feedback is tuned to match the stellar mass-halo mass relationship. For comparison, each galaxy was also simulated using a 'lower feedback' (LF) model which has strength comparable to other implementations in the literature. The 'MAGICC feedback' (MF) model has a higher incidence of massive stars and an approximately two times higher energy input per supernova. Apart from the low-mass halo using LF, each galaxy exhibits a metal-enriched CGM that extends to approximately the virial radius. A significant fraction of this gas has been heated in supernova explosions in the disc and subsequently ejected into the CGM where it is predicted to give rise to substantial O VI absorption. The simulations do not yet address the question of what happens to the O VI when the galaxies stop forming stars. Our models also predict a reservoir of cool H I clouds that show strong Lyαabsorption to several hundred kpc. Comparing these models to recent surveys with the Hubble Space Telescope, we find that only the MF models have sufficient O VI and H I gas in the CGM to reproduce the observed distributions. In separate analyses, these same MF models also show better agreement with other galaxy observables (e.g. rotation curves, surface brightness profiles and H I gas distribution). We infer that the CGM is the dominant reservoir of baryons for galaxy haloes.

preprint2012arXiv

Numerical hydrodynamic simulations based on semi-analytic galaxy merger trees: method and Milky-Way like galaxies

We present a new approach to study galaxy evolution in a cosmological context. We combine cosmological merger trees and semi-analytic models of galaxy formation to provide the initial conditions for multi-merger hydrodynamic simulations. In this way we exploit the advantages of merger simulations (high resolution and inclusion of the gas physics) and semi-analytic models (cosmological background and low computational cost), and integrate them to create a novel tool. This approach allows us to study the evolution of various galaxy properties, including the treatment of the hot gaseous halo from which gas cools and accretes onto the central disc, which has been neglected in many previous studies. This method shows several advantages over other methods. As only the particles in the regions of interest are included, the run time is much shorter than in traditional cosmological simulations, leading to greater computational efficiency. Using cosmological simulations, we show that multiple mergers are expected to be more common than sequences of isolated mergers, and therefore studies of galaxy mergers should take this into account. In this pilot study, we present our method and illustrate the results of simulating ten Milky Way-like galaxies since z=1. We find good agreement with observations for the total stellar masses, star formation rates, cold gas fractions and disc scale length parameters. We expect that this novel numerical approach will be very useful for pursuing a number of questions pertaining to the transformation of galaxy internal structure through cosmic time.

preprint2011arXiv

Interaction between dark matter sub-halos & galactic gaseous disk

We investigate the idea that the interaction of Dark Matter (DM) sub-halos with the gaseous disk of galaxies can be the origin for the observed holes and shells found in their neutral hydrogen (HI) distributions. We use high resolution hydrodynamic simulations to show that pure DM sub-halos impacting a galactic disk are not able to produce holes; on the contrary, they result in high density regions in the disk. However, sub-halos containing a small amount of gas (a few percent of the total DM mass of the sub-halo) are able to displace the gas in the disk and form holes and shells. The size and lifetime of these holes depend on the sub-halo gas mass, density and impact velocity. A DM sub-halo of mass $10^8$ $M_{\odot}$ and a gas mass fraction of $\sim 3$%, is able to create a kpc scale hole, with a lifetime similar to those observed in nearby galaxies. We also register an increase in the star formation rate at the rim of the hole, again in agreement with observations. Even though the properties off these simulated structures resemble those found in observations we find that the number of predicted holes (based on mass and orbital distributions of DM halos derived from cosmological N-body simulations) falls short compared to the observations. Only a handful of holes are produced per Giga year. This leads us to conclude that DM halo impact is not the major channel through which these holes are formed.

preprint2011arXiv

Non-linear weak lensing forecasts

We investigate the impact of non-linear corrections on dark energy parameter estimation from weak lensing probes. We find that using halofit expressions, suited to LCDM models, implies substantial discrepancies with respect to results directly obtained from N-body simulations, when w(z)\neq-1. Discrepancies appear strong when using models with w'(z=0)>0, as fiducial models; they are however significant even in the neighborhood of LCDM, where neglecting the degrees of freedom associated with the DE state equation can lead to a misestimate of the matter density parameter Ω_m.

preprint2011arXiv

What Sets the Sizes of the Faintest Galaxies?

We provide a comprehensive description and offer an explanation for the sizes of the faintest known galaxies in the universe, the dwarf spheroidal (dSph) satellites of the Milky Way and Andromeda. After compiling a consistent data set of half-light radii (r_{1/2}) and luminosities, we describe the size-luminosity relation of dSphs by a log-normal distribution in r_{1/2} with a mean size that varies as a function of luminosity. Accounting for modest number statistics, measurement uncertainties and surface brightness limitations, we find that the size-luminosity relations of the Milky Way and Andromeda dSph populations are statistically indistinguishable, and also very similar: their mean sizes at a given stellar luminosity differ by no more than 30%. In addition, we find that the mean size, slope and scatter of this log-normal size description of Local Group dSphs matches onto the relation of more massive low-concentration galaxies. This suggests that the stellar sizes of dSphs are ultimately related to their overall initial baryonic angular momentum. To test this hypothesis we perform a series of high resolution N-body simulations that we couple with a semi-analytic model of galaxy formation. These predict the same mean size and slope as observed in dSph satellites. At the same time, these models predict that the size-luminosity distributions for satellite galaxies around similar host-halos must be similar providing a natural explanation as to why the size distributions of Milky Way and Andromeda satellites are similar. Although strong rotation is currently not observed in dSphs, this may well be consistent with our angular-momentum-based explanation for their sizes if the disks of these galaxies have become sufficiently stirred through tidal interaction.

preprint2009arXiv

A quantitative explanation of the observed population of Milky Way satellite galaxies

We revisit the well known discrepancy between the observed number of Milky Way (MW) dwarf satellite companions and the predicted population of cold dark matter (CDM) sub-halos, in light of the dozen new low luminosity satellites found in SDSS imaging data and our recent calibration of the SDSS satellite detection efficiency, which implies a total population far larger than these dozen discoveries. We combine a dynamical model for the CDM sub-halo population with simple, physically motivated prescriptions for assigning stellar content to each sub-halo, then apply observational selection effects and compare to the current observational census. As expected, models in which the stellar mass is a constant fraction F(Omega_b/Omega_m) of the sub-halo mass M_sat at the time it becomes a satellite fail for any choice of F. However, previously advocated models that invoke suppression of gas accretion after reionization in halos with circular velocity v_c <~ 35 km/s can reproduce the observed satellite counts for -15 < M_V < 0, with F ~ 10^{-3}. Successful models also require strong suppression of star formation BEFORE reionization in halos with v_c <~ 10 km/s; models without pre-reionization suppression predict far too many satellites with -5 < M_V < 0. Our models also reproduce the observed stellar velocity dispersions ~ 5-10 km/s of the SDSS dwarfs given the observed sizes of their stellar distributions, and model satellites have M(<300 pc) ~ 10^7 M_sun as observed even though their present day total halo masses span more than two orders of magnitude. Our modeling shows that natural physical mechanisms acting within the CDM framework can quantitatively explain the properties of the MW satellite population as it is presently known, thus providing a convincing solution to the `missing satellite' problem.