We present a simple model for the response of a dissipationless spherical system to an instantaneous mass change at its center, describing the formation of flat cores in dark matter haloes and ultra-diffuse galaxies (UDGs) from feedback-driven outflow episodes in a specific mass range. This model generalizes an earlier simplified analysis of an isolated shell into a system with continuous density, velocity and potential profiles. The response is divided into an instantaneous change of potential at constant velocities due to a given mass loss or gain, followed by energy-conserving relaxation to a new Jeans equilibrium. The halo profile is modeled by a two-parameter function with a variable inner slope and an analytic potential profile (Dekel et al. 2017), which enables determining the associated kinetic energy at equilibrium. The model is tested against NIHAO cosmological zoom-in simulations, where it successfully predicts the evolution of the inner dark-matter profile between successive snapshots in about 75% of the cases, failing mainly in merger situations. This model provides a simple understanding of the formation of dark-matter halo cores and UDGs by supernova-driven outflows, and a useful analytic tool for studying such processes.

We study ultra-diffuse galaxies (UDGs) in zoom in cosmological simulations, seeking the origin of UDGs in the field versus galaxy groups. We find that while field UDGs arise from dwarfs in a characteristic mass range by multiple episodes of supernova feedback (Di Cintio et al. 2017), group UDGs may also form by tidal puffing up and they become quiescent by ram-pressure stripping. The field and group UDGs share similar properties, independent of distance from the group centre. Their dark-matter haloes have ordinary spin parameters and centrally dominant dark-matter cores. Their stellar components tend to have a prolate shape with a Sersic index n~1 but no significant rotation. Ram pressure removes the gas from the group UDGs when they are at pericentre, quenching star formation in them and making them redder. This generates a colour/star-formation-rate gradient with distance from the centre, as observed in clusters. We find that ~20 per cent of the field UDGs that fall into a massive halo survive as satellite UDGs. In addition, normal field dwarfs on highly eccentric orbits can become UDGs near pericentre due to tidal puffing up, contributing about half of the group-UDG population. We interpret our findings using simple toy models, showing that gas stripping is mostly due to ram pressure rather than tides. We estimate that the energy deposited by tides in the bound component of a satellite over one orbit can cause significant puffing up provided that the orbit is sufficiently eccentric.

We analyze 88 independent high-resolution cosmological zoom-in simulations of disk galaxies in the NIHAO simulations suite to explore the connection between the atomic gas fraction and angular momentum of baryons throughout cosmic time. The study is motivated by the analytic model of \citet{obreschkow16}, which predicts a relation between the atomic gas fraction $f_{\rm atm}$ and the global atomic stability parameter $q \equiv j\sigma / (GM)$, where $M$ and $j$ are the mass and specific angular momentum of the galaxy (stars+cold gas) and $\sigma$ is the velocity dispersion of the atomic gas. We show that the simulated galaxies follow this relation from their formation ($z\simeq4$) to present within $\sim 0.5$ dex. To explain this behavior, we explore the evolution of the local Toomre stability and find that $90\%$--$100\%$ of the atomic gas in all simulated galaxies is stable at any time. In other words, throughout the entire epoch of peak star formation until today, the timescale for accretion is longer than the timescale to reach equilibrium, thus resulting in a quasi-static equilibrium of atomic gas at any time. Hence, the evolution of $f_{\rm atm}$ depends on the complex hierarchical growth history primarily via the evolution of $q$. An exception are galaxies subject to strong environmental effects.

Galaxy formation entails the hierarchical assembly of mass, along with the condensation of baryons and the ensuing, self-regulating star formation. The stars form a collisionless system whose orbit distribution retains dynamical memory that can constrain a galaxy's formation history. The ordered-rotation dominated orbits with near maximum circularity $\lambda_z \simeq1$ and the random-motion dominated orbits with low circularity $\lambda_z \simeq0$ are called kinematically cold and kinematically hot, respectively. The fraction of stars on `cold' orbits, compared to the fraction of stars on `hot' orbits, speaks directly to the quiescence or violence of the galaxies' formation histories. Here we present such orbit distributions, derived from stellar kinematic maps via orbit-based modelling for a well defined, large sample of 300 nearby galaxies. The sample, drawn from the CALIFA survey, includes the main morphological galaxy types and spans the total stellar mass range from $10^{8.7}$ to $10^{11.9}$ solar masses. Our analysis derives the orbit-circularity distribution as a function of galaxy mass, $p(\lambda_z~|~M_\star)$, and its volume-averaged total distribution, $p(\lambda_z)$. We find that across most of the considered mass range and across morphological types, there are more stars on `warm' orbits defined as $0.25\le \lambda_z \le 0.8$ than on either `cold' or `hot' orbits. This orbit-based "Hubble diagram" provides a benchmark for galaxy formation simulations in a cosmological context.

We present useful functions for the profiles of dark-matter (DM) haloes with a free inner slope, from cusps to cores, where the profiles of density, mass-velocity and potential are simple analytic expressions. Analytic velocity is obtained by expressing the mean density as a simple functional form, and deriving the local density by differentiation. The function involves four shape parameters, with only two or three free: a concentration parameter $c$, inner and outer asymptotic slopes $\alpha$ and $\bar{\gamma}$, and a middle shape parameter $\beta$. Analytic expressions for the potential and velocity dispersion exist for $\bar{\gamma}=3$ and for $\beta$ a natural number. We match the models to the DM haloes in cosmological simulations, with and without baryons, ranging from steep cusps to flat cores. Excellent fits are obtained with three free parameters ($c$, $\alpha$, $\bar{\gamma}$) and $\beta=2$. For an analytic potential, similar fits are obtained for $\bar{\gamma}=3$ and $\beta=2$ with only two free parameters ($c$, $\alpha$); this is our favorite model. A linear combination of two such profiles, with an additional free concentration parameter, provides excellent fits also for $\beta=1$, where the expressions are simpler. The fit quality is comparable to non-analytic popular models. An analytic potential is useful for modeling the inner-halo evolution due to gas inflows and outflows, studying environmental effects on the outer halo, and generating halo potentials or initial conditions for simulations. The analytic velocity can quantify simulated and observed rotation curves without numerical integrations.

We present a useful function for describing the profiles of dark-matter haloes with a varying asymptotic inner slope -\alpha, ranging from a cusp to a core, where the profiles of density, mass-velocity and potential are simple analytic expressions for any \alpha. The idea is to express the mean-density profile as a simple functional form, and obtain the local density by derivative. The model involves a concentration parameter c analogous to the NFW profile. More flexibility is provided by a sum of two such functions, with concentrations c_1 and c_2. An optional additional parameter is the asymptotic outer slope -\gamma, which allows more flexibility in the outskirts of haloes. For \gamma different than 3, there are analytic expressions for the profiles of density and mass-velocity but not for the potential. We match the proposed function, in different variants, to the dark-matter profiles of haloes in cosmological simulations, with and without baryons, spanning inner profiles that range from steep cusps to flat cores, and find excellent fits. The analytic potential profile is useful in modeling the evolution of the inner halo due to episodes of gas inflow and outflow. The analytic profile, especially with a free \gamma, is useful for a study of environmental effects in the outer halo. In general, these analytic profiles can serve for straightforwardly quantifying the shapes of simulated and observed rotation curves, without the need for numerical integrations.

We use semi-analytic models and cosmological merger trees to provide the initial conditions for multi-merger numerical hydrodynamic simulations, and exploit these simulations to explore the effect of galaxy interaction and merging on star formation (SF). We compute numerical realisations of twelve merger trees from z=1.5 to z=0. We include the effects of the large hot gaseous halo around all galaxies, following recent obervations and predictions of galaxy formation models. We find that including the hot gaseous halo has a number of important effects. Firstly, as expected, the star formation rate on long timescales is increased due to cooling of the hot halo and refuelling of the cold gas reservoir. Secondly, we find that interactions do not always increase the SF in the long term. This is partially due to the orbiting galaxies transferring gravitational energy to the hot gaseous haloes and raising their temperature. Finally we find that the relative size of the starburst, when including the hot halo, is much smaller than previous studies showed. Our simulations also show that the order and timing of interactions are important for the evolution of a galaxy. When multiple galaxies interact at the same time, the SF enhancement is less than when galaxies interact in series. All these effects show the importance of including hot gas and cosmologically motivated merger trees in galaxy evolution models.

We study how X-rays from stellar binary systems and the hot intracluster medium (ICM) affect the radiative cooling rates of gas in galaxies. Our study uses a novel implementation of gas cooling in the moving-mesh hydrodynamics code \textsc{arepo}. X-rays from stellar binaries do not affect cooling at all as their emission spectrum is too hard to effectively couple with galactic gas. In contrast, X-rays from the ICM couple well with gas in the temperature range $10^4 - 10^6$ K. Idealised simulations show that the hot halo radiation field has minimal impact on the dynamics of cooling flows in clusters because of the high virial temperature ($> 10^7$K), making the interaction between the gas and incident photons very ineffective. Satellite galaxies in cluster environments, on the other hand, experience a high radiation flux due to the emission from the host halo. Low mass satellites ($< 10^{12}\rm{M_\odot}$) in particular have virial temperatures that are exactly in the regime where the effect of the radiation field is maximal. Idealised simulations of satellite galaxies including only the effect of host halo radiation (no ram pressure stripping or tidal effects) fields show a drastic reduction in the amount of cool gas formed ($\sim 40\%$) on a short timescale of about $0.5$ Gyrs. A galaxy merger simulation including all the other environmental quenching mechanisms, shows about $20\%$ reduction in the stellar mass of the satellite and about $\sim 30\%$ reduction in star formation rate after $1$ Gyr due to the host hot halo radiation field. These results indicate that the hot halo radiation fields potentially play an important role in quenching galaxies in cluster environments.

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

We use 31 simulated galaxies from the MaGICC project to investigate the effects of baryonic feedback on the density profiles of dark matter (DM) haloes. The sample covers a wide mass range: 9.4e9<Mhalo/Msun<7.8e11, hosting galaxies with stellar masses: 5.0e5<M*/Msun<8.3e10, i.e. from dwarf to L*. The galaxies are simulated with several baryonic prescriptions, including a range of stellar feedbacks. The main result is a clear dependence of the inner slope of the DM density profile, \alpha\ in \rho r^\alpha, on the ratio between stellar-to-halo mass (M*/Mhalo). This relation is independent of the stellar feedback scheme, allowing a prediction for cusp vs core formation. When M*/Mhalo is low, ~0.01%, energy from stellar feedback is insufficient to significantly alter the inner DM density and the galaxy retains a cuspy profile. At higher M*/Mhalo, feedback drives the expansion of the DM and generates cored profiles. The flattest profiles form where M*/Mhalo~0.5%. Above this ratio, stars formed in the central regions deepen the gravitational potential enough to oppose this supernova-driven expansion process, resulting in smaller cores and cuspier profiles. Combining the dependence of \alpha\ on M*/Mhalo with the abundance matching relation between M* and Mhalo provides a prediction for how \alpha\ varies with M*. Further, using the Tully-Fisher relation allows a prediction for the dependence of the DM inner slope on the observed rotation velocity of galaxies. The most cored galaxies are expected to have Vrot~50km/s, with \alpha\ decreasing for more massive disc galaxies: spirals with Vrot~150km/s have central slopes \alpha<-0.8, approaching the NFW profile. This novel prediction for the dependence of \alpha\ on disc galaxy mass can be tested using current observational data sets, and can be applied to theoretical modeling of mass profiles and populations of disc galaxies.

Quasars represent the brightest Active Galactic Nuclei (AGN) in the Universe and are thought to indicate the location of prodigiously growing Black Holes (BHs), with luminosities as high as 10^48 erg/sec. It is often expected though that such an extremely energetic process will take place in the most massive bound structures in the dark matter (DM) distribution. We show that in contrast to this expectation, in a galaxy formation model which includes AGN feedback, quasars are predicted to live in average DM halo environments with typical masses of a few times 10^12 Msun. This fundamental prediction arises from the fact that quasar activity (i.e., BH accretion with luminosity greater than 10^46 erg/sec) is inhibited in DM haloes where AGN feedback operates. The galaxy hosts of quasars in our simulations are identified with over massive (in gas and stars) spheroidal galaxies, in which BH accretion is triggered via a galaxy merger or secular processes. We further show that the z=0 descendants of high redshift (z~6) QSOs span a wide range of morphologies, galaxy and halo masses. The z~6 BHs typically grow only by a modest factor by the present day. Remarkably, high redshift QSOs never inhabit the largest DM haloes at that time and their descendants are very seldom found in the most massive haloes at z=0. We also show that observationally it is very likely to find an enhancement in the abundance of galaxies around quasars at z~5. However, these enhancements are considerably weaker compared to the overdensities expected at the extreme peaks of the DM distribution. Thus, it is very unlikely that a quasar detected in the $z\gtrsim5$ Universe pinpoints the location of the progenitors of superclusters in the local Universe.

The nature of dark matter is still unknown and one of the most fundamental scientific mysteries. Although successfully describing large scales, the standard cold dark matter model (CDM) exhibits possible shortcomings on galactic and sub-galactic scales. It is exactly at these highly non-linear scales where strong astrophysical constraints can be set on the nature of the dark matter particle. While observations of the Lyman-$\alpha$ forest probe the matter power spectrum in the mildly non-linear regime, satellite galaxies of the Milky Way provide an excellent laboratory as a test of the underlying cosmology on much smaller scales. Here we present results from a set of high resolution simulations of a Milky Way sized dark matter halo in eight distinct cosmologies: CDM, warm dark matter (WDM) with a particle mass of 2 keV and six different cold plus warm dark matter (C+WDM) models, varying the fraction, $f_{\rm wdm}$, and the mass, $m_{\rm wdm}$, of the warm component. We used three different observational tests based on Milky Way satellite observations: the total satellite abundance, their radial distribution and their mass profile. We show that the requirement of simultaneously satisfying all three constraints sets very strong limits on the nature of dark matter. This shows the power of a multi-dimensional small scale approach in ruling out models which would be still allowed by large scale observations.

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

Recent studies have found a dramatic difference between the observed number density evolution of low mass galaxies and that predicted by semi-analytic models. While models accurately reproduce the z=0 number density, they require that the evolution occurs rapidly at early times, which is incompatible with the strong late evolution found observationally. We report here the same discrepancy in two state-of-the-art cosmological hydrodynamical simulations, which is evidence that the problem is fundamental. We search for the underlying cause of this problem using two complementary methods. Firstly, we look for evidence of a different history of today's low mass galaxies in models and observations and we find that the models yield too few young, strongly star-forming galaxies. Secondly, we construct a toy model to link the observed evolution of specific star formation rates (sSFR) with the evolution of the galaxy stellar mass function. We infer from this model that a key problem in both semi-analytic and hydrodynamical models is the presence of a positive instead of a negative correlation between sSFR and stellar mass. A similar positive correlation is found between the specific dark matter halo accretion rate and the halo mass, indicating that model galaxies are growing in a way that follows the growth of their host haloes too closely. It therefore appears necessary to find a mechanism that decouples the growth of low mass galaxies, which occurs primarily at late times, from the growth of their host haloes, which occurs primarily at early times. We argue that the current form of star-formation driven feedback implemented in most galaxy formation models is unlikely to achieve this goal, owing to its fundamental dependence on host halo mass and time. [Abridged]

We examine the chemical properties of 5 cosmological hydrodynamical simulations of an M33-like disc galaxy which have been shown to be consistent with the morphological characteristics and bulk scaling relations expected of late-type spirals. These simulations are part of the Making Galaxies In a Cosmological Context (MaGICC) Project, in which stellar feedback is tuned to match the stellar mass -- halo mass relationship. Each realisation employed identical initial conditions and assembly histories, but differed from one another in their underlying baryonic physics prescriptions, including (a) the efficiency with which each supernova energy couples to the ISM, (b) the impact of feedback associated with massive star radiation pressure, (c) the role of the minimum shut-off time for radiative cooling of Type II SNe remnants, (d) the treatment of metal diffusion, and (e) varying the IMF. Our analysis focusses on the resulting stellar metallicity distribution functions (MDFs) in each simulated (analogous) `solar neighbourhood' and central `bulge' region. We compare the simulated MDFs' skewness, kurtosis, and dispersion (inter-quartile, inter-decile, inter-centile, and inter-tenth-percentile regions) with that of the empirical solar neighbourhood MDF and Local Group dwarfs. We find that the MDFs of the simulated discs are more negatively skewed, with higher kurtosis, than those observed locally. We can trace this difference to the simulations' tight and correlated age-metallicity relations (compared with that of the Milky Way), suggesting that these relations within `dwarf' discs might be steeper than in L* discs and/or the degree of stellar orbital re-distribution and migration inferred locally has not been captured in their entirety, at the resolution of our simulations. The important role of metal diffusion in ameliorating the over-production of extremely metal-poor stars is highlighted.

The dark energy dominated warm dark matter (WDM) model is a promising alternative cosmological scenario. We explore large-scale structure formation in this paradigm. We do this in two different ways: with the halo model approach and with the help of an ensemble of high resolution N-body simulations. Combining these quasi-independent approaches, leads to a physical understanding of the important processes which shape the formation of structures. We take a detailed look at the halo mass function, the concentrations and the linear halo bias of WDM. In all cases we find interesting deviations with respect to CDM. In particular, the concentration-mass relation displays a turnover for group scale dark matter haloes, for the case of WDM particles with masses of the order ~0.25 keV. This may be interpreted as a hint for top-down structure formation on small scales. We implement our results into the halo model and find much better agreement with simulations. On small scales the WDM halo model now performs as well as its CDM counterpart.

A candidate diffuse stellar substructure was previously reported in the halo of the nearby dwarf starburst galaxy NGC 4449 by Karachentsev et al. We map and analyze this feature using a unique combination of deep integrated-light images from the Black Bird 0.5-meter telescope, and high-resolution wide-field images from the 8-meter Subaru telescope, which resolve the nebulosity into a stream of red giant branch stars, and confirm its physical association with NGC 4449. The properties of the stream imply a massive dwarf spheroidal progenitor, which after complete disruption will deposit an amount of stellar mass that is comparable to the existing stellar halo of the main galaxy. The ratio between luminosity or stellar-mass between the two galaxies is ~1:50, while the indirectly measured dynamical mass-ratio, when including dark matter, may be ~1:10-1:5. This system may thus represent a "stealth" merger, where an infalling satellite galaxy is nearly undetectable by conventional means, yet has a substantial dynamical influence on its host galaxy. This singular discovery also suggests that satellite accretion can play a significant role in building up the stellar halos of low-mass galaxies, and possibly in triggering their starbursts.

We show that the black hole-bulge mass scaling relations observed from the local to the high-z Universe can be largely or even entirely explained by a non-causal origin, i.e. they do not imply the need for any physically coupled growth of black hole and bulge mass, for example through feedback by active galactic nuclei (AGN). Provided some physics for the absolute normalisation, the creation of the scaling relations can be fully explained by the hierarchical assembly of black hole and stellar mass through galaxy merging, from an initially uncorrelated distribution of BH and stellar masses in the early Universe. We show this with a suite of dark matter halo merger trees for which we make assumptions about (uncorrelated) black hole and stellar mass values at early cosmic times. We then follow the halos in the presence of global star formation and black hole accretion recipes that (i) work without any coupling of the two properties per individual galaxy and (ii) correctly reproduce the observed star formation and black hole accretion rate density in the Universe. With disk-to-bulge conversion in mergers included, our simulations even create the observed slope of ~1.1 for the M_BH-M_bulge-relations at z=0. This also implies that AGN feedback is not a required (though still a possible) ingredient in galaxy evolution. In light of this, other mechanisms that can be invoked to truncate star formation in massive galaxies are equally justified.

We investigate the effects of a coupled Dark Energy (cDE) scalar field on the alignment between satellites and matter distributions in galaxy clusters. Using high-resolution N-body simulations for LCDM and cDE cosmological models, we compute the probability density distribution for the alignment angle between the satellite galaxies and underlying matter distributions, finding a difference between the two scenarios. With respect to LCDM, in cDE cosmologies the satellite galaxies are less preferentially located along the major axis of the matter distribution, possibly reducing the tension with obersevational data. A physical explanation is that the coupling between dark matter and dark energy acts as an additional tidal force on the satellite galaxies diminishing the alignments between their distribution and the matter one. Through a Wald test based on the generalized chi-square statistics, the null hypothesis that the two probability distributions come from the same parent population is rejected at the 99 % confidence level. It is concluded that the galaxy-matter alignment in clusters may provide a unique probe of dark sector interactions as well as the nature of dark energy.

We present non-parametric radial mass profiles for ten QSO strong lensing galaxies. Five of the galaxies have profiles close to $\rho(r)\propto r^{-2}$, while the rest are closer to r^{-1}, consistent with an NFW profile. The former are all relatively isolated early-types and dominated by their stellar light. The latter --though the modeling code did not know this-- are either in clusters, or have very high mass-to-light, suggesting dark-matter dominant lenses (one is a actually pair of merging galaxies). The same models give $H_0^{-1} = 15.2_{-1.7}^{+2.5}\Gyr$ ($H_0 = 64_{-9}^{+8} \legacy$), consistent with a previous determination. When tested on simulated lenses taken from a cosmological hydrodynamical simulation, our modeling pipeline recovers both H_0 and $\rho(r)$ within estimated uncertainties. Our result is contrary to some recent claims that lensing time delays imply either a low H_0 or galaxy profiles much steeper than r^{-2}. We diagnose these claims as resulting from an invalid modeling approximation: that small deviations from a power-law profile have a small effect on lensing time-delays. In fact, as we show using using both perturbation theory and numerical computation from a galaxy-formation simulation, a first-order perturbation of an isothermal lens can produce a zeroth-order change in the time delays.

We use the conditional luminosity function (CLF) and data from the 2dFGRS to constrain the average relation between light and mass in a LCDM cosmology with Omega_m=0.23 and sigma_8=0.74 (hereafter WMAP3 cosmology). Reproducing the observed luminosity dependence of the galaxy two-point correlation function results in average mass-to-light ratios that are about 35 percent lower than in a LCDM cosmology with Omega_m=0.3 and sigma_8=0.9 (hereafter WMAP1 cosmology). This removes an important problem with previous halo occupation models which had a tendency to predict cluster mass-to-light ratios that were too high. For the WMAP3 cosmology our model yields average mass-to-light ratios, central galaxy luminosities, halo occupation numbers, satellite fractions, and luminosity-gap statistics, that are all in excellent agreement with those obtained from a 2dFGRS group catalogue and from other independent studies. We also use our CLF model to compute the probability distribution P(M|L_cen), that a central galaxy of luminosity L_cen resides in a halo of mass M. We find this distribution to be much broader than what is typically assumed in HOD models, which has important implications for the interpretation of galaxy-galaxy lensing data. Finally, reproducing the luminosity dependence of the pairwise velocity dispersions in the 2dFGRS requires relatively low mass-to-light ratios for clusters and a satellite fraction that decreases strongly with increasing luminosity. This is only marginally consistent with our CLF constraints. We argue that a cosmology with parameters between those of the WMAP1 and WMAP3 cosmologies is likely to yield results with a higher level of consistency.