Central Halo Research Articles

Dark Matter in Galaxies: Evidences and Challenges

The evidence of the phenomenon for which, in galaxies, the gravitating mass is distributed differently than the luminous mass, increases as new data become available. Furthermore, this discrepancy is well structured and it depends on the magnitude and the compactness of the galaxy and on the radius, in units of its luminous size $R_{opt}$, where the measurement is performed. For the disk systems with $-13\geq M_I\geq -24$ all this leads to an amazing scenario, revealed by the investigation of individual and coadded rotation curves, according to which, the circular velocity follows, from their centers out to their virial radii, an universal profile $V_{URC} (r/R_{opt}, M_I)$ function only of the properties of the luminous mass component. Moreover, from the Universal Rotation Curve, so as from many individual high-quality RCs, we discover that, in the innermost regions of galaxies, the DM halo density profiles are very shallow. Finally, the disk mass, the central halo density and its core radius, come out all related to each other and to two properties of the distribution of light in galaxies: the luminosity and the compactness. This phenomenology, being absent in the simplest $\Lambda CDM$ Cosmology scenario, poses serious challenges to the latter or, alternatively, it requires a substantial and tuned involvement of baryons in the formation of the galactic halos. On the other side, the URC helps to explain the two-accelerations relationship found by McGaugh et al 2016, in terms of only well known astrophysical processes, acting in a standard DM halos + luminous disks scenario

The impact of cored density profiles on the observable quantities of dwarf spheroidal galaxies

ABSTRACT We modify the chemo-dynamical code gear to simulate the impact of self-interacting dark matter (SIDM) on the observable quantities of 19 low-mass dwarf galaxies with a variety star-forming properties. We employ a relatively high, velocity independent cross-section of $\sigma /m=10\, \rm {cm^2\,g^{ -1}}$ and extract, in addition to integrated quantities, the total mass density profile, the luminosity profile, the line-of-sight velocities, the chemical abundance, and the star formation history. We find that despite the creation of large cores at the centre of the dark matter haloes, the impact of SIDM on the observable quantities of quenched galaxies is indiscernible, dominated mostly by the stochastic build up of the stellar matter. As such we conclude that it is impossible to make global statements on the density profile of dwarf galaxies from single or small samples. Although based mostly on quenched galaxies, this finding supports other recent work putting into question the reliability of inferred cored density profiles that are derived from observed line-of-sight velocities.

The different growth pathways of Brightest Cluster Galaxies and the Intra-Cluster Light

We study the growth pathways of Brightest Central Galaxies (BCGs) and Intra-Cluster Light (ICL) by means of a semi-analytic model. We assume that the ICL forms by stellar stripping of satellite galaxies and violent processes during mergers, and implement two independent models: (1) one considers both mergers and stellar stripping (named {\small STANDARD} model), and one considers only mergers (named {\small MERGERS} model). We find that BCGs and ICL form, grow and overall evolve at different times and with different timescales, but they show a clear co-evolution after redshift $z \sim 0.7-0.8$. Around 90\% of the ICL from stellar stripping is built-up in the innermost 150 Kpc from the halo centre and the dominant contribution comes from disk-like galaxies (B/T$<$0.4) through a large number of small/intermediate stripping events ($M_{strip}/M_{sat}<0.3$). The fractions of stellar mass in BCGs and in ICL over the total stellar mass within the virial radius of the halo evolve differently with time. At high redshift, the BCG accounts for the bulk of the mass, but its contribution gradually decreases with time and stays constant after $z\sim 0.4-0.5$. The ICL, instead, grows very fast and its contribution keeps increasing down to the present time. The {\small STANDARD} and the {\small MERGERS} models make very similar predictions in most of the cases, but predict different amounts of ICL associated to other galaxies within the virial radius of the group/cluster other than the BCG, at $z=0$. We then suggest that this quantity is a valid observable that can shed light on the relative importance of mergers and stellar stripping for the formation of the ICL.

Correlation between centre offsets and gas velocity dispersion of galaxy clusters in cosmological simulations

The gas is the dominant component of baryonic matter in most galaxy groups and clusters. The spatial offsets of gas centre from the halo centre could be an indicator of the dynamical state of cluster. Knowledge of such offsets is important for estimate the uncertainties when using clusters as cosmological probes. In this paper, we study the centre offsets $r_{\rm off}$ between the gas and that of all the matter within halo systems in $\Lambda$CDM cosmological hydrodynamic simulations. We focus on two kinds of centre offsets: one is the three-dimensional PB offsets between the gravitational potential minimum of the entire halo and the barycentre of the ICM, and the other is the two-dimensional PX offsets between the potential minimum of the halo and the iterative centroid of the projected synthetic X-ray emission of the halo. Halos at higher redshifts tend to have larger values of rescaled offsets $r_{\rm off}/r_{200}$ and larger gas velocity dispersion $\sigma_{v}^{\rm gas}/\sigma_{200}$. For both types of offsets, we find that the correlation between the rescaled centre offsets $r_{\rm off}/r_{200}$ and the rescaled 3D gas velocity dispersion, $\sigma_v^{\rm gas}/\sigma_{200}$ can be approximately described by a quadratic function as $r_{{\rm off}}/r_{200} \propto (\sigma_{v}^{\rm gas}/\sigma_{200} - k_2)^{2}$. A Bayesian analysis with MCMC method is employed to estimate the model parameters. Dependence of the correlation relation on redshifts and the gas mass fraction are also investigated.

Painting galaxies into dark matter haloes using machine learning

We develop a machine learning (ML) framework to populate large dark matter-only simulations with baryonic galaxies. Our ML framework takes input halo properties including halo mass, environment, spin, and recent growth history, and outputs central galaxy and halo baryonic properties including stellar mass ($M_*$), star formation rate (SFR), metallicity ($Z$), neutral ($\rm HI$) and molecular ($\rm H_2$) hydrogen mass. We apply this to the MUFASA cosmological hydrodynamic simulation, and show that it recovers the mean trends of output quantities with halo mass highly accurately, including following the sharp drop in SFR and gas in quenched massive galaxies. However, the scatter around the mean relations is under-predicted. Examining galaxies individually, at $z=0$ the stellar mass and metallicity are accurately recovered ($\sigma\lesssim 0.2$~dex), but SFR and $\rm HI$ show larger scatter ($\sigma\gtrsim 0.3$~dex); these values improve somewhat at $z=1,2$. Remarkably, ML quantitatively recovers second parameter trends in galaxy properties, e.g. that galaxies with higher gas content and lower metallicity have higher SFR at a given $M_*$. Testing various ML algorithms, we find that none perform significantly better than the others, nor does ensembling improve performance, likely because none of the algorithms reproduce the large observed scatter around the mean properties. For the random forest algorithm, we find that halo mass and nearby ($\sim 200$~kpc) environment are the most important predictive variables followed by growth history, while halo spin and $\sim$Mpc scale environment are not important. Finally we study the impact of additionally inputting key baryonic properties $M_*$, SFR and $Z$, as would be available e.g. from an equilibrium model, and show that particularly providing the SFR enables $\rm HI$ to be recovered substantially more accurately.

How do stars affect ψDM haloes?

Wave dark matter ($\psi$DM) predicts a compact soliton core and a granular halo in every galaxy. This work presents the first simulation study of an elliptical galaxy by including both stars and $\psi$DM, focusing on the systematic changes of the central soliton and halo granules. With the addition of stars in the inner halo, we find the soliton core consistently becomes more prominent by absorbing mass from the host halo than that without stars, and the halo granules become "non-isothermal", "hotter" in the inner halo and "cooler" in the outer halo, as opposed to the isothermal halo in pure $\psi$DM cosmological simulations. Moreover, the composite (star+$\psi$DM) mass density is found to follow a $r^{-2}$ isothermal profile near the half-light radius in most cases. Most striking is the velocity dispersion of halo stars that increases rapidly toward the galactic center by a factor of at least 2 inside the half-light radius caused by the deepened soliton gravitational potential, a result that compares favorably with observations of elliptical galaxies and bulges in spiral galaxies. However in some rare situations we find a phase segregation turning a compact distribution of stars into two distinct populations with high and very low velocity dispersions; while the high-velocity component mostly resides in the halo, the very low-velocity component is bound to the interior of the soliton core, resembling stars in faint dwarf spheroidal galaxies.

Using velocity dispersion to estimate halo mass: Is the Local Group in tension with ΛCDM?

Satellite galaxies are commonly used as tracers to measure the line-of-sight velocity dispersion ($\sigma_{\rm LOS}$) of the dark matter halo associated with their central galaxy, and thereby to estimate the halo's mass. Recent observational dispersion estimates of the Local Group, including the Milky Way and M31, suggest $\sigma\sim$50 km/s, which is surprisingly low when compared to the theoretical expectation of $\sigma\sim$100s km/s for systems of their mass. Does this pose a problem for $\Lambda$CDM? We explore this tension using the {\small{SURFS}} suite of $N$-body simulations, containing over 10000 (sub)haloes with well tracked orbits. We test how well a central galaxy's host halo velocity dispersion can be recovered by sampling $\sigma_{\rm LOS}$ of subhaloes and surrounding haloes. Our results demonstrate that $\sigma_{\rm LOS}$ is biased mass proxy. We define an optimal window in $v_{\rm LOS}$ and projected distance ($D_p$) -- $0.5\lesssim D_p/R_{\rm vir}\lesssim1.0$ and $v_{\rm LOS} \lesssim0.5V_{\rm esc}$, where $R_{\rm vir}$ is the virial radius and $V_{\rm esc}$ is the escape velocity -- such that the scatter in LOS to halo dispersion is minimised - $\sigma_{\rm LOS}=(0.5\pm0.1)\sigma_{v,{\rm H}}$. We argue that this window should be used to measure line-of-sight dispersions as a proxy for mass, as it minimises scatter in the $\sigma_{\rm LOS}-M_{\rm vir}$ relation. This bias also naturally explains the results from \cite{mcconnachie2012a}, who used similar cuts when estimating $\sigma_{\rm LOS,LG}$, producing a bias of $\sigma_{\rm LG}=(0.44\pm0.14)\sigma_{v,{\rm H}}$. We conclude that the Local Group's velocity dispersion does not pose a problem for $\Lambda$CDM and has a mass of $\log M_{\rm LG, vir}/M_\odot=12.0^{+0.8}_{-2.0}$.

Why does Einasto profile index n∼ 6 occur so frequently?

We consider the behavior of spherically symmetric Einasto halos composed of gravitating particles in the Fokker-Planck approximation. This approach allows us to consider the undesirable influence of close encounters in the N-body simulations more adequately than the generally accepted criteria. The Einasto profile with index n ≈ 6 is a stationary solution of the Fokker-Planck equation in the halo center. There are some reasons to believe that the solution is an attractor. Then the Fokker-Planck diffusion tends to transform a density profile to the equilibrium one with the Einasto index n ≈ 6. We suggest this effect as a possible reason why the Einasto index n ≈ 6 occurs so frequently in the interpretation of N-body simulation results. The results obtained cast doubt on generally accepted criteria of N-body simulation convergence.

Constraints from microlensing experiments on clustered primordial black holes

Abstract It has recently been proposed that massive primordial black holes (PBH) could constitute all of the dark matter, providing a novel scenario of structure formation, with early reionization and a rapid growth of the massive black holes at the center of galaxies and dark matter halos. The scenario arises from broad peaks in the primordial power spectrum that give both a spatially clustered and an extended mass distribution of PBH. The constraints from the observed microlensing events on the extended mass function have already been addressed. Here we study the impact of spatial clustering on the microlensing constraints. We find that the bounds can be relaxed significantly for relatively broad mass distributions if the number of primordial black holes within each cluster is typically above one hundred. On the other hand, even if they arise from individual black holes within the cluster, the bounds from CMB anisotropies are less stringent due to the enhanced black hole velocity in such dense clusters. This way, the window between a few and ten solar masses has opened up for PBH to comprise the totality of the dark matter.

Incidental branch retinal artery occlusion on optical coherence tomography angiography presenting as segmental optic atrophy in a child: a case report

BackgroundRetinal artery occlusion is extremely rare in the pediatric population and most patients have risk factors. We report a case of a healthy child with segmental optic atrophy, complicated by incidental branch retinal artery occlusion (BRAO).Case presentationA 10-year-old boy who had a history of his mother’s gestational diabetes presented with an inferonasal visual field defect in the left eye. His best-corrected visual acuities were 20/20 in both eyes (OU). Fundoscopic examination revealed segmental pallor of the left optic disc, thinning of the superotemporal rim, a relative superior entrance of the central retinal artery and superior peripapillary scleral halo. Fluorescein angiography showed patchy filling delays in the corresponding disc area without retinal vascular abnormalities. Spectral domain optical coherence tomography (SD OCT) via automated segmentation analysis demonstrated sectoral absence of the ganglion cell layer and retinal nerve fiber layer with thinning of the inner plexiform layer, inner nuclear layer and outer plexiform layer in the corresponding retina. OCT angiography (OCTA) showed focal attenuation of superficial and intermediate/deep capillary plexuses in the corresponding areas. Systemic evaluation was unremarkable. The patient was diagnosed with segmental optic atrophy caused by incidental BRAO.ConclusionsRetinal vascular occlusions are rare in childhood, and may present as segmental optic atrophy mimicking congenital anomalies. OCTA allows the detection of previous microvascular abnormalities in the chronic phase. To the best of our knowledge, this is the first report of a child with segmental optic atrophy presumably caused by BRAO, which was documented by SD OCT and OCTA in detail.

Gas flows in the circumgalactic medium around simulated high-redshift galaxies

We analyse the properties of circumgalactic gas around simulated galaxies in the redshift range z >= 3, utilising a new sample of cosmological zoom simulations. These simulations are intended to be representative of the observed samples of Lyman-alpha emitters recently obtained with the MUSE instrument (halo masses ~10^10-10^11 solar masses). We show that supernova feedback has a significant impact on both the inflowing and outflowing circumgalactic medium by driving outflows, reducing diffuse inflow rates, and by increasing the neutral fraction of inflowing gas. By temporally stacking simulation outputs we find that significant net mass exchange occurs between inflowing and outflowing phases: none of the phases are mass-conserving. In particular, we find that the mass in neutral outflowing hydrogen declines exponentially with radius as gas flows outwards from the halo centre. This is likely caused by a combination of both fountain-like cycling processes and gradual photo/collisional ionization of outflowing gas. Our simulations do not predict the presence of fast-moving neutral outflows in the CGM. Neutral outflows instead move with modest radial velocities (~ 50 kms^-1), and the majority of the kinetic energy is associated with tangential rather than radial motion.

Tau neutrinos from ultracompact dark matter minihalos and constraints on the primordial curvature perturbations

The observations and research on the neutrinos provide a kind of indirect way of revealing the properties of dark matter particles. For the detection of muon neutrinos, the main issue is the large atmospheric background, which is caused by the interactions between the cosmic rays and atoms within the atmosphere. Compared with muon neutrinos, tau neutrinos have a smaller atmospheric background especially for the downward-going direction. Except for the classical neutrino sources, dark matter particles can also annihilate into the neutrinos and are the potential high energy astrophysical sources. The annihilation rate of dark matter particles is proportional to the square of number density; therefore, the annihilation rate is large near the center of dark matter halos especially for the new kind of dark matter structures named ultracompact dark matter minihalos (UCMHs). In previous works, we have investigated the potential muon neutrino flux from UCMHs due to dark matter annihilation. Moreover, since the formation of UCMHs is related to the primordial density perturbations of small scales, we get the constraints on the amplitude of the primordial curvature perturbations of small scales, $1 \lesssim k \lesssim 10^{7} ~\rm Mpc^{-1}$. In this work, we focus on the downward-going tau neutrinos from UCMHs due to dark matter annihilation. Compared with the background of tau neutrino flux we get the constraints on the mass fraction of UCMHs. Then using the limits on the mass fraction of UCMHs we got the constraints on the amplitude of the primordial curvature perturbations which are extended to the scale $k \sim 10^{8} ~ \rm Mpc^{-1}$ compared with previous results.

The structure and assembly history of cluster-sized haloes in self-interacting dark matter

We perform dark-matter-only simulations of 28 relaxed massive cluster-sized haloes for Cold Dark Matter (CDM) and Self-Interacting Dark Matter (SIDM) models, to study structural differences between the models at large radii, where the impact of baryonic physics is expected to be very limited. We find that the distributions for the radial profiles of the density, ellipsoidal axis ratios, and velocity anisotropies ($\beta$) of the haloes differ considerably between the models (at the $\sim1\sigma$ level), even at $\gtrsim10\%$ of the virial radius, if the self-scattering cross section is $\sigma/m_\chi=1$ cm$^2$ gr$^{-1}$. Direct comparison with observationally inferred density profiles disfavours SIDM for $\sigma/m_\chi=1$ cm$^2$ gr$^{-1}$, but in an intermediate radial range ($\sim3\%$ of the virial radius), where the impact of baryonic physics is uncertain. At this level of the cross section, we find a narrower $\beta$ distribution in SIDM, clearly skewed towards isotropic orbits, with no SIDM (90\% of CDM) haloes having $\beta>0.12$ at $7\%$ of the virial radius. We estimate that with an observational sample of $\sim30$ ($\sim10^{15}$ M$_\odot$) relaxed clusters, $\beta$ can potentially be used to put competitive constraints on SIDM, once observational uncertainties improve by a factor of a few. We study the suppression of the memory of halo assembly history in SIDM clusters. For $\sigma/m_\chi=1$ cm$^2$ gr$^{-1}$, we find that this happens only in the central halo regions ($\sim1/4$ of the scale radius of the halo), and only for haloes that assembled their mass within this region earlier than a formation redshift $z_f\sim2$. Otherwise, the memory of assembly remains and is reflected in ways similar to CDM, albeit with weaker trends.

Brightest galaxies as halo centre tracers in SDSS DR7

Abstract Determining the positions of halo centres in large-scale structure surveys is crucial for many cosmological studies. A common assumption is that halo centres correspond to the location of their brightest member galaxies. In this paper, we study the dynamics of brightest galaxies with respect to other halo members in the Sloan Digital Sky Survey DR7. Specifically, we look at the line-of-sight velocity and spatial offsets between brightest galaxies and their neighbours. We compare those to detailed mock catalogues, constructed from high-resolution, dark-matter-only N-body simulations, in which it is assumed that satellite galaxies trace dark matter subhaloes. This allows us to place constraints on the fraction fBNC of haloes in which the brightest galaxy is not the central. Compared to previous studies, we explicitly take into account the unrelaxed state of the host haloes, velocity offsets of halo cores and correlations between fBNC and the satellite occupation. We find that fBNC strongly decreases with the luminosity of the brightest galaxy and increases with the mass of the host halo. Overall, in the halo mass range 1013–1014.5 h− 1M⊙ we find fBNC ∼ 30 per cent, in good agreement with a previous study by Skibba et al. We discuss the implications of these findings for studies inferring the galaxy–halo connection from satellite kinematics, models of the conditional luminosity function and galaxy formation in general.

Observations of the galaxy cluster CIZA J2242.8+5301 with the Sardinia Radio Telescope

We observed the galaxy cluster CIZA J2242.8+5301 with the Sardinia Radio Telescope to provide new constraints on its spectral properties at high frequency. We conducted observations in three frequency bands centred at 1.4 GHz, 6.6 GHz and 19 GHz, resulting in beam resolutions of 14$^{\prime}$, 2.9$^{\prime}$ and 1$^{\prime}$ respectively. These single-dish data were also combined with archival interferometric observations at 1.4 and 1.7 GHz. From the combined images, we measured a flux density of ${\rm S_{1.4GHz}=(158.3\pm9.6)\,mJy}$ for the central radio halo and ${\rm S_{1.4GHz}=(126\pm8)\,mJy}$ and ${\rm S_{1.4GHz}=(11.7\pm0.7)\,mJy}$ for the northern and the southern relic respectively. After the spectral modelling of the discrete sources, we measured at 6.6 GHz ${\rm S_{6.6GHz}=(17.1\pm1.2)\,mJy}$ and ${\rm S_{6.6GHz}=(0.6\pm0.3)\,mJy}$ for the northern and southern relic respectively. Assuming simple diffusive shock acceleration, we interpret measurements of the northern relic with a continuous injection model represented by a broken power-law. This yields an injection spectral index ${\rm \alpha_{inj}=0.7\pm0.1}$ and a Mach number ${\rm M=3.3\pm0.9}$, consistent with recent X-ray estimates. Unlike other studies of the same object, no significant steepening of the relic radio emission is seen in data up to 8.35 GHz. By fitting the southern relic spectrum with a simple power-law (${\rm S_{\nu}\propto\nu^{-\alpha}}$) we obtained a spectral index ${\rm \alpha\approx1.9}$ corresponding to a Mach number (${\rm M\approx1.8}$) in agreement with X-ray estimates. Finally, we evaluated the rotation measure of the northern relic at 6.6 GHz. These results provide new insights on the magnetic structure of the relic, but further observations are needed to clarify the nature of the observed Faraday rotation.

Fire in the field: simulating the threshold of galaxy formation

We present a suite of 15 cosmological zoom-in simulations of isolated dark matter halos, all with masses of $M_{\rm halo} \approx 10^{10}\,{\rm M}_\odot$ at $z=0$, in order to understand the relationship between halo assembly, galaxy formation, and feedback's effects on the central density structure in dwarf galaxies. These simulations are part of the Feedback in Realistic Environments (FIRE) project and are performed at extremely high resolution. The resultant galaxies have stellar masses that are consistent with rough abundance matching estimates, coinciding with the faintest galaxies that can be seen beyond the virial radius of the Milky Way ($M_\star/{\rm M}_\odot\approx 10^5-10^7$). This non-negligible spread in stellar mass at $z=0$ in halos within a narrow range of virial masses is strongly correlated with central halo density or maximum circular velocity $V_{\rm max}$. Much of this dependence of $M_\star$ on a second parameter (beyond $M_{\rm halo}$) is a direct consequence of the $M_{\rm halo}\sim10^{10}\,{\rm M}_\odot$ mass scale coinciding with the threshold for strong reionization suppression: the densest, earliest-forming halos remain above the UV-suppression scale throughout their histories while late-forming systems fall below the UV-suppression scale over longer periods and form fewer stars as a result. In fact, the latest-forming, lowest-concentration halo in our suite fails to form any stars. Halos that form galaxies with $M_\star\gtrsim2\times10^{6}\,{\rm M}_\odot$ have reduced central densities relative to dark-matter-only simulations, and the radial extent of the density modifications is well-approximated by the galaxy half-mass radius $r_{1/2}$. This apparent stellar mass threshold of $M_\star \approx 2\times 10^{6} \approx 2\times 10^{-4} \,M_{\rm halo}$ is broadly consistent with previous work and provides a testable prediction of FIRE feedback models in LCDM.

The origin of the mass discrepancy–acceleration relation in ΛCDM

We examine the origin of the mass discrepancy--radial acceleration relation (MDAR) of disk galaxies. This is a tight empirical correlation between the disk centripetal acceleration and that expected from the baryonic component. The MDAR holds for most radii probed by disk kinematic tracers, regardless of galaxy mass or surface brightness. The relation has two characteristic accelerations; $a_0$, above which all galaxies are baryon-dominated; and $a_{\rm min}$, an effective minimum aceleration probed by kinematic tracers in isolated galaxies. We use a simple model to show that these trends arise naturally in $\Lambda$CDM. This is because: (i) disk galaxies in $\Lambda$CDM form at the centre of dark matter haloes spanning a relatively narrow range of virial mass; (ii) cold dark matter halo acceleration profiles are self-similar and have a broad maximum at the centre, reaching values bracketed precisely by $a_{\rm min}$ and $a_0$ in that mass range; and (iii) halo mass and galaxy size scale relatively tightly with the baryonic mass of a galaxy in any successful $\Lambda$CDM galaxy formation model. Explaining the MDAR in $\Lambda$CDM does not require modifications to the cuspy inner mass profiles of dark haloes, although these may help to understand the detailed rotation curves of some dwarf galaxies and the origin of extreme outliers from the main relation. The MDAR is just a reflection of the self-similar nature of cold dark matter haloes and of the physical scales introduced by the galaxy formation process.

Not so lumpy after all: modelling the depletion of dark matter subhaloes by Milky Way-like galaxies

Among the most important goals in cosmology is detecting and quantifying small ($M_{\rm halo}\simeq10^{6-9}~\mathrm{M}_\odot$) dark matter (DM) subhalos. Current probes around the Milky Way (MW) are most sensitive to such substructure within $\sim20$ kpc of the halo center, where the galaxy contributes significantly to the potential. We explore the effects of baryons on subhalo populations in $\Lambda$CDM using cosmological zoom-in baryonic simulations of MW-mass halos from the Latte simulation suite, part of the Feedback In Realistic Environments (FIRE) project. Specifically, we compare simulations of the same two halos run using (1) DM-only (DMO), (2) full baryonic physics, and (3) DM with an embedded disk potential grown to match the FIRE simulation. Relative to baryonic simulations, DMO simulations contain $\sim2\times$ as many subhalos within 100 kpc of the halo center; this excess is $\gtrsim5\times$ within 25 kpc. At $z=0$, the baryonic simulations are completely devoid of subhalos down to $3\times10^6~\mathrm{M}_\odot$ within $15$ kpc of the MW-mass galaxy, and fewer than 20 surviving subhalos have orbital pericenters <20 kpc. Despite the complexities of baryonic physics, the simple addition of an embedded central disk potential to DMO simulations reproduces this subhalo depletion, including trends with radius, remarkably well. Thus, the additional tidal field from the central galaxy is the primary cause of subhalo depletion. Subhalos on radial orbits that pass close to the central galaxy are preferentially destroyed, causing the surviving subhalo population to have tangentially biased orbits compared to DMO predictions. Our method of embedding a disk potential in DMO simulations provides a fast and accurate alternative to full baryonic simulations, thus enabling suites of cosmological simulations that can provide accurate and statistical predictions of substructure populations.

Exponentially growing bubbles around early supermassive black holes

We addressed the so far unexplored issue of outflows induced by exponentially growing power sources, focusing on early supermassive black holes (BHs). We assumed that these objects grow to $10^9\;M_{\odot}$ by z=6 by Eddington-limited accretion and convert 5% of their bolometric output into a wind. We first considered the case of energy-driven and momentum-driven outflows expanding in a region where the gas and total mass densities are uniform and equal to the average values in the Universe at $z>6$. We derived analytic solutions for the evolution of the outflow, finding that, for an exponentially growing power with e-folding time $t_{Sal}$, the late time expansion of the outflow radius is also exponential, with e-folding time of $5t_{Sal}$ and $4t_{Sal}$ in the energy-driven and momentum-driven limit, respectively. We then considered energy-driven outflows produced by QSOs at the center of early dark matter halos of different masses and powered by BHs growing from different seeds. We followed the evolution of the source power and of the gas and dark matter density profiles in the halos from the beginning of the accretion until $z=6$. The final bubble radius and velocity do not depend on the seed BH mass but are instead smaller for larger halo masses. At z=6, bubble radii in the range 50-180 kpc and velocities in the range 400-1000 km s$^{-1}$ are expected for QSOs hosted by halos in the mass range $3\times10^{11}-10^{13}\;M_{\odot}$. By the time the QSO is observed, we found that the total thermal energy injected within the bubble in the case of an energy-driven outflow is $E_{th}\sim5 \times 10^{60}$ erg. This is in excellent agreement with the value of $E_{th}=(6.2\pm 1.7)\times 10^{60}$ erg measured through the detection of the thermal Sunyaev-Zeldovich effect around a large population of luminous QSOs at lower redshift. [abridged]

Properties and Origin of Galaxy Velocity Bias in the Illustris Simulation

Abstract We use hydrodynamical galaxy formation simulations from the Illustris suite to study the origin and properties of galaxy velocity bias, i.e., the difference between the velocity distributions of galaxies and dark matter inside halos. We find that galaxy velocity bias decreases with increasing ratio of galaxy stellar mass to host halo mass. In general, central galaxies are not at rest with respect to dark matter halos or the core of halos, with a velocity dispersion above 0.04 times that of the dark matter. The central galaxy velocity bias is found to be mostly caused by close interactions between the central and satellite galaxies. For satellite galaxies, the velocity bias is related to their dynamical and tidal evolution history after being accreted onto the host halos. It depends on the time after the accretion and their distances from the halo centers, with massive satellites generally moving more slowly than the dark matter. The results are in broad agreement with those inferred from modeling small-scale redshift-space galaxy clustering data, and the study can help improve models of redshift-space galaxy clustering.