Search

Cosmological galaxy formation simulations are powerful tools to understand the complex processes that govern the formation and evolution of galaxies. However, evaluating the realism of these simulations remains a challenge. The two common approaches for evaluating galaxy simulations is either through scaling relations based on a few key physical galaxy properties, or through a set of pre-defined morphological parameters based on galaxy images. This paper proposes a novel image-based method for evaluating the quality of galaxy simulations using unsupervised deep learning anomaly detection techniques. By comparing full galaxy images, our approach can identify and quantify discrepancies between simulated and observed galaxies. As a demonstration, we apply this method to SDSS imaging and NIHAO simulations with different physics models, parameters, and resolution. We further compare the metric of our method to scaling relations as well as morphological parameters. We show that anomaly detection is able to capture similarities and differences between real and simulated objects that scaling relations and morphological parameters are unable to cover, thus indeed providing a new point of view to validate and calibrate cosmological simulations against observed data., Comment: 13 pages, 16 figures, published in MNRAS

We append two additional black hole (BH) accretion models, namely viscous disc and gravitational torque-driven accretion, into the Numerical Investigation of a Hundred Astrophysical Objects (NIHAO) project of galaxy simulations. We show that these accretion models, characterized by a weaker dependence on the BH mass compared to the commonly used Bondi-Hoyle accretion, naturally create a common evolutionary track (co-existence) between the mass of the BH and the stellar mass of the galaxy, even without any direct coupling via feedback (FB). While FB is indeed required to control the final BH and stellar mass of the galaxies, our results suggest that FB might not be the leading driver of the cosmic co-evolution between these two quantities; in these models, co-evolution is simply determined by the shared central gas supply. Conversely, simulations using Bondi-Hoyle accretion show a two-step evolution, with an early growth of stellar mass followed by exponential growth of the central supermassive black hole (SMBH). Our results show that the modelling of BH accretion (sometimes overlooked) is an extremely important part of BH evolution and can improve our understanding of how scaling relations emerge and evolve, and whether SMBH and stellar mass co-exist or co-evolve through cosmic time., Comment: 13 pages, 10 figures

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 ($\lambda_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., Comment: Accepted for publication in MNRAS. 15 pages, 10 figures

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., Comment: Accepted by MNRAS

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., Comment: 22 pages, 14 figures, submitted to MNRAS, comments welcome

A central question regarding Ultra Diffuse Galaxies (UDGs) is whether they are a separate category to Low Surface Brightness (LSB) galaxies, or just their natural continuation towards low stellar masses. In this letter, we show that the rotation curve of the gas rich UDG AGC 242019 is well fit by a dark matter halo with inner slope that asymptotes to -0.54, and that such fit provides a concentration parameter that matches theoretical expectations. This finding, together with previously works in which shallow inner profiles are derived for UDGs, shows that the structural properties of these galaxies are like other observed LSBs. UDGs show slowly rising rotation curves and this favours formation scenarios in which internal processes, such as SNae driven gas outflows, are acting to modify UDGs profiles., Comment: Accepted for publication in ApJL

We use cosmological hydrodynamical galaxy formation simulations from the NIHAO project to investigate the response of cold dark matter (CDM) haloes to baryonic processes. Previous work has shown that the halo response is primarily a function of the ratio between galaxy stellar mass and total virial mass, and the density threshold above which gas is eligible to form stars, $n [{\rm cm}^{-3}]$. At low $n$ all simulations in the literature agree that dwarf galaxy haloes are cuspy, but at high $n\ge 100$ there is no consensus. We trace halo contraction in dwarf galaxies with $n\ge 100$ reported in some previous simulations to insufficient spatial resolution. Provided the adopted star formation threshold is appropriate for the resolution of the simulation, we show that the halo response is remarkably stable for $n\ge 5$, up to the highest star formation threshold that we test, $n=500$. This free parameter can be calibrated using the observed clustering of young stars. Simulations with low thresholds $n\le 1$ predict clustering that is too weak, while simulations with high star formation thresholds $n\ge 5$, are consistent with the observed clustering. Finally, we test the CDM predictions against the circular velocities of nearby dwarf galaxies. Low thresholds predict velocities that are too high, while simulations with $n\sim 10$ provide a good match to the observations. We thus conclude that the CDM model provides a good description of the structure of galaxies on kpc scales provided the effects of baryons are properly captured., Comment: 15 pages, 10 figures, published in MNRAS

We investigate how the NIHAO galaxies match the observed star formation main sequence (SFMS) and what the origin of its scatter is. The NIHAO galaxies reproduce the SFMS and generally agree with observations, but the slope is about unity and thus significantly larger than observed values. This is because observed galaxies at large stellar masses, although still being part of the SFMS, are already influenced by quenching. This partial suppression of star formation by AGN feedback leads to lower star formation rates and therefore to lower observed slopes. We confirm that including the effects of AGN in our galaxies leads to slopes in agreement with observations. We find the deviation of a galaxy from the SFMS is correlated with its $z=0$ dark matter halo concentration and thus with its halo formation time. This means galaxies with a higher-than-average star formation rate (SFR) form later and vice versa. We explain this apparent correlation with the SFR by re-interpreting galaxies that lie above the SFMS (higher-than-average SFR) as lying to the left of the SFMS (lower-than-average stellar mass) and vice versa. Thus later forming haloes have a lower-than-average stellar mass, this is simply because they have had less-than-average time to form stars, and vice versa. It is thus the nature, i.e. how and when these galaxies form, that sets the path of a galaxy in the SFR versus stellar mass plane., Comment: 8 pages, 11 figures, accepted by MNRAS

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 ($\alpha$) is fairly constant ($\Delta\alpha \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)., Comment: 5 pages, 4 figures, accepted for publication in MNRAS Letters

We implement an optically thin approximation for the effects of the local radiation field from stars and hot gas on the gas heating and cooling in the N-body SPH code GASOLINE2. We resimulate three galaxies from the NIHAO project: one dwarf, one Milky Way-like and one massive spiral, and study what are the local radiation field effects on various galaxy properties. We also study the effects of varying the Ultra Violet Background (UVB) model, by running the same galaxies with two different UVBs. Galaxy properties at $z=0$ like stellar mass, stellar effective mass radius, HI mass, and radial extent of the HI disc, show significant changes between the models with and without the local radiation field, and smaller differences between the two UVB models. The intrinsic effect of the local radiation field through cosmic time is to increase the equilibrium temperature at the interface between the galaxies and their circumgalactic media (CGM), moving this boundary inwards, while leaving relatively unchanged the gas inflow rate. Consequently, the temperature of the inflow increases when considering the local radiation sources. This temperature increase is a function of total galaxy mass, with a median CGM temperature difference of one order of magnitude for the massive spiral. The local radiation field suppresses the stellar mass growth by $\sim$20 per cent by $z=0$ for all three galaxies, while the HI mass is roughly halfed. The differences in the gas phase diagrams, significantly impact the HI column densities, shifting their peaks in the distributions towards lower $N_{\rm HI}$., Comment: resubmitted to MNRAS after minor revision

We introduce algorithms for black hole physics, i.e., black hole formation, accretion and feedback, into the NIHAO (Numerical Investigation of a Hundred Astrophysical Objects) project of galaxy simulations. This enables us to study high mass, elliptical galaxies, where feedback from the central black hole is generally thought to have a significant effect on their evolution. We furthermore extend the NIHAO suite by 45 simulations that encompass $z=0$ halo masses from $1 \times 10^{12}$ to $4 \times 10^{13}\,\mathrm{M}_{\odot}$, and resimulate five galaxies from the original NIHAO sample with black hole physics, which have $z=0$ halo masses from $8 \times 10^{11}$ to $3 \times 10^{12}\,\mathrm{M}_{\odot}$. Now NIHAO contains 144 different galaxies and thus has the largest sample of zoom-in simulations of galaxies, spanning $z=0$ halo masses from $9 \times 10^{8}$ to $4 \times 10^{13}\,\mathrm{M}_{\odot}$. In this paper we focus on testing the algorithms and calibrating their free parameters against the stellar mass versus halo mass relation and the black hole mass versus stellar mass relation. We also investigate the scatter of these relations, which we find is a decreasing function with time and thus in agreement with observations. For our fiducial choice of parameters we successfully quench star formation in objects above a $z=0$ halo mass of $10^{12}\,\mathrm{M}_{\odot}$, thus transforming them into red and dead galaxies., Comment: 14 pages, 11 figures, accepted by MNRAS

We use thirty-eight high-resolution simulations of galaxy formation between redshift 10 and 5 to study the impact of a 3 keV warm dark matter (WDM) candidate on the high-redshift Universe. We focus our attention on the stellar mass function and the global star formation rate and consider the consequences for reionization, namely the neutral hydrogen fraction evolution and the electron scattering optical depth. We find that three different effects contribute to differentiate warm and cold dark matter (CDM) predictions: WDM suppresses the number of haloes with mass less than few $10^9$ M$_{\odot}$; at a fixed halo mass, WDM produces fewer stars than CDM; and finally at halo masses below $10^9$ M$_{\odot}$, WDM has a larger fraction of dark haloes than CDM post-reionization. These three effects combine to produce a lower stellar mass function in WDM for galaxies with stellar masses at and below $\sim 10^7$ M$_{\odot}$. For $z > 7$, the global star formation density is lower by a factor of two in the WDM scenario, and for a fixed escape fraction, the fraction of neutral hydrogen is higher by 0.3 at $z \sim 6$. This latter quantity can be partially reconciled with CDM and observations only by increasing the escape fraction from 23 per cent to 34 per cent. Overall, our study shows that galaxy formation simulations at high redshift are a key tool to differentiate between dark matter candidates given a model for baryonic physics., Comment: 11 pages, 8 figures, submitted to MNRAS

In this third paper of the series, we investigate the effects of warm dark matter with a particle mass of $m_\mathrm{WDM}=3\,\mathrm{keV}$ on the smallest galaxies in our Universe. We present a sample of 21 hydrodynamical cosmological simulations of dwarf galaxies and 20 simulations of satellite-host galaxy interaction that we performed both in a Cold Dark Matter (CDM) and Warm Dark Matter (WDM) scenario. In the WDM simulations, we observe a higher critical mass for the onset of star formation. Structure growth is delayed in WDM, as a result WDM haloes have a stellar population on average two Gyrs younger than their CDM counterparts. Nevertheless, despite this delayed star formation, CDM and WDM galaxies are both able to reproduce the observed scaling relations for velocity dispersion, stellar mass, size, and metallicity at $z=0$. WDM satellite haloes in a Milky Way mass host are more susceptible to tidal stripping due to their lower concentrations, but their galaxies can even survive longer than the CDM counterparts if they live in a dark matter halo with a steeper central slope. In agreement with our previous CDM satellite study we observe a steepening of the WDM satellites' central dark matter density slope due to stripping. The difference in the average stellar age for satellite galaxies, between CDM and WDM, could be used in the future for disentangling these two models., Comment: 10 pages, 11 figures, accepted for publication on MNRAS

We use cosmological hydrodynamical galaxy formation simulations from the NIHAO project to investigate the impact of the threshold for star formation on the response of the dark matter (DM) halo to baryonic processes. The fiducial NIHAO threshold, $n=10\, {\rm cm}^{-3}$, results in strong expansion of the DM halo in galaxies with stellar masses in the range $10^{7.5} < M_{star} < 10^{9.5} M_{\odot}$. We find that lower thresholds such as $n=0.1$ (as employed by the EAGLE/APOSTLE and Illustris/AURIGA projects) do not result in significant halo expansion at any mass scale. Halo expansion driven by supernova feedback requires significant fluctuations in the local gas fraction on sub-dynamical times (i.e., < 50 Myr at galaxy half-light radii), which are themselves caused by variability in the star formation rate. At one per cent of the virial radius, simulations with $n=10$ have gas fractions of $\simeq 0.2$ and variations of $\simeq 0.1$, while $n=0.1$ simulations have order of magnitude lower gas fractions and hence do not expand the halo. The observed DM circular velocities of nearby dwarf galaxies are inconsistent with CDM simulations with $n=0.1$ and $n=1$, but in reasonable agreement with $n=10$. Star formation rates are more variable for higher $n$, lower galaxy masses, and when star formation is measured on shorter time scales. For example, simulations with $n=10$ have up to 0.4 dex higher scatter in specific star formation rates than simulations with $n=0.1$. Thus observationally constraining the sub-grid model for star formation, and hence the nature of DM, should be possible in the near future., Comment: 18 pages, 13 figures, accepted to MNRAS

In the first paper we presented 27 hydrodynamical cosmological simulations of galaxies with total masses between $5 \times 10^8$ and $10^{10}\,\mathrm{M}_\odot$. In this second paper we use a subset of these cosmological simulations as initial conditions (ICs) for more than forty hydrodynamical simulations of satellite and host galaxy interaction. Our cosmological ICs seem to suggest that galaxies on these mass scales have very little rotational support and are velocity dispersion ($\sigma$) dominated. Accretion and environmental effects increase the scatter in the galaxy scaling relations (e.g. size - velocity dispersion) in very good agreement with observations. Star formation is substantially quenched after accretion. Mass removal due to tidal forces has several effects: it creates a very flat stellar velocity dispersion profiles, and it reduces the dark matter content at all scales (even in the centre), which in turn lowers the stellar velocity on scales around 0.5 kpc even when the galaxy does not lose stellar mass. Satellites that start with a cored dark matter profile are more prone to either be destroyed or to end up in a very dark matter poor galaxy. Finally, we found that tidal effects always increase the "cuspyness" of the dark matter profile, even for haloes that infall with a core., Comment: 13 pages, 15 figures, accepted to be published in MNRAS (with minor changes to the published version). This is the second paper of a series of two. The first paper can be found here arXiv:1707.01106

The satellites of the Milky Way and Andromeda represent the smallest galaxies we can observe in our Universe. In this series of papers we aim to shed light on their formation and evolution using cosmological hydrodynamical simulations. In this first paper we focus on the galaxy properties before accretion, by simulating twenty seven haloes with masses between $5\times 10^8$ and $10^{10} M_\odot$. Out of this set nineteen haloes successfully form stars, while eight remain dark. The simulated galaxies match quite well present day observed scaling relations between stellar mass, size and metallicity, showing that such relations are in place before accretion. Our galaxies show a large variety of star formation histories, from extended star formation periods to single bursts. As in more massive galaxies, large star formation bursts are connected with major mergers events, which greatly contribute to the overall stellar mass build up. The intrinsic stochasticity of mergers induces a large scatter in the stellar mass halo mass relation, up to two orders of magnitude. Despite the bursty star formation history, on these mass scales baryons are very ineffective in modifying the dark matter profiles, and galaxies with a stellar mass below $\approx 10^6 M_\odot$ retain their cuspy central dark matter distribution, very similar to results from pure N-body simulations., Comment: 12 pages, 15 figures. Two new figures added, discussion improved, several new references, conclusions unchanged. Accepted to be published in MNRAS. the second paper can found at arXiv:1707.01102

There is strong observational evidence indicating a time lag of order of some 100 Myr between the onset of starburst and AGN activity in galaxies. Dynamical time lags have been invoked to explain this. We extend this approach by introducing a viscous time lag the gas additionally needs to flow through the AGN's accretion disc before it reaches the central black hole. Our calculations reproduce the observed time lags and are in accordance with the observed correlation between black hole mass and stellar velocity dispersion., Comment: 11 pages, 8 figures, accepted by MNRAS

The circumnuclear disc (CND) orbiting the Galaxy's central black hole is a reservoir of material that can ultimately provide energy through accretion, or form stars in the presence of the black hole, as evidenced by the stellar cluster that is presently located at the CND's centre. In this paper, we report the results of a computational study of the dynamics of the CND. The results lead us to question two paradigms that are prevalent in previous research on the Galactic Centre. The first is that the disc's inner cavity is maintained by the interaction of the central stellar cluster's strong winds with the disc's inner rim, and second, that the presence of unstable clumps in the disc implies that the CND is a transient feature. Our simulations show that, in the absence of a magnetic field, the interaction of the wind with the inner disc rim actually leads to a filling of the inner cavity within a few orbital time-scales, contrary to previous expectations. However, including the effects of magnetic fields stabilizes the inner disc rim against rapid inward migration. Furthermore, this interaction causes instabilities that continuously create clumps that are individually unstable against tidal shearing. Thus the occurrence of such unstable clumps does not necessarily mean that the disc is itself a transient phenomenon. The next steps in this investigation are to explore the effect of the magnetorotational instability on the disc evolution and to test whether the results presented here persist for longer time-scales than those considered here., Comment: 18 pages, 11 figures, accepted by MNRAS

We derive analytical expressions for the mass flow rates in and from accretion discs, taking into account the Eddington limit. This allows us to connect the basic properties of outflows with those of the accretion flow. As an example, we derive the radial surface density distribution in the accretion disc of Mrk 231., Comment: 6 pages, 2 figures, submitted to MNRAS