Search

We present predictions of spectral energy distributions (SEDs), from the UV to the FIR, of simulated galaxies at z = 0. These were obtained by post-processing the results of an N-body+hydro simulation of a cosmological box of side 25 Mpc, which uses the Multi-Phase Particle Integrator (MUPPI) for star formation and stellar feedback, with the GRASIL-3D radiative transfer code that includes reprocessing of UV light by dust. Physical properties of our sample of ~500 galaxies resemble observed ones, though with some tension at small and large stellar masses. Comparing predicted SEDs of simulated galaxies with different samples of local galaxies, we find that these resemble observed ones, when normalized at 3.6 µm. A comparison with the Herschel Reference Survey shows that the average SEDs of galaxies, divided in bins of star formation rate (SFR), are reproduced in shape and absolute normalization to within a factor of ~2, while average SEDs of galaxies divided in bins of stellar mass show tensions that are an effect of the difference of simulated and observed galaxies in the stellar mass-SFR plane. We use our sample to investigate the correlation of IR luminosity in Spitzer and Herschel bands with several galaxy properties. SFR is the quantity that best correlates with IR light up to 160 µm, while at longer wavelengths better correlations are found with molecular mass and, at 500 μm, with dust mass. However, using the position of the FIR peak as a proxy for cold dust temperature, we assess that heating of cold dust is mostly determined by SFR, with stellar mass giving only a minor contribution. We finally show how our sample of simulated galaxies can be used as a guide to understand the physical properties and selection biases of observed samples. [ABSTRACT FROM AUTHOR]

We present a quantification of the properties of bars in two N-body+ smoothed particle hydrodynamics cosmological simulations of spiral galaxies, named GA and AqC. The initial conditions were obtained using the zoom-in technique and represent two dark matter haloes of 2-3 x 1012 M⊙, available at two different resolutions. The resulting galaxies are presented in the companion paper of Murante et al. We find that the GA galaxy has a bar of length 8.8 kpc, present at the two resolution levels even though with a slightly different strength. Classical bar signatures (e.g. pattern of streaming motions, high m = 2 Fourier mode with roughly constant phase) are consistently found at both resolutions. Though a close encounter with a merging satellite at z ~ 0.6 (mass ratio 1: 50) causes a strong, transient spiral pattern and some heating of the disc, we find that bar instability is due to secular process, caused by a low Toomre parameter Q ≲ 1 due to accumulation of mass in the disc. The AqC galaxy has a slightly different history: it suffers a similar tidal disturbance due to a merging satellite at z ~ 0.5 but with a mass ratio of 1: 32, that triggers a bar in the high-resolution simulation, while at low resolution the merging is found to take place at a later time, so that both secular evolution and merging are plausible triggers for bar instability. [ABSTRACT FROM AUTHOR]

We measure and quantify properties of galactic outflows and diffuse gas at z ≥ 1 in cosmological hydrodynamical simulations. Our novel subresolution model, Multi-Phase Particle Integrator (MUPPI), implements supernova feedback using fully local gas properties, where the wind velocity and mass loading are not given as input. We find the following trends at z = 2 by analysing central galaxies having a stellar mass higher than 109M⊙. The outflow velocity and mass outflow rate (Ṁout) exhibit positive correlations with galaxy mass and with the star formation rate (SFR). However, most of the relations present a large scatter. The outflow mass loading factor (η) is between 0.2 and 10. The comparison effective model generates a constant outflow velocity, and a negative correlation of η with halo mass. The number fraction of galaxies where outflow is detected decreases at lower redshifts, but remains more than 80 per cent over z = 1-5. The outflow velocity correlation with SFR becomes flatter at z = 1, and η displays a negative correlation with halo mass in massive galaxies. Our study demonstrates that both the MUPPI and effective models produce significant outflows at ~1/10 of the virial radius; at the same time shows that the properties of outflows generated can be different from the input speed and mass loading in the effective model. Our MUPPI model, using local properties of gas in the subresolution recipe, is able to develop galactic outflows whose properties correlate with global galaxy properties, and consistent with observations. [ABSTRACT FROM AUTHOR]

We present results of cosmological simulations of disc galaxies carried out with the gadget-3 TreePM+SPH code, where star formation and stellar feedback are described using our MUlti Phase Particle Integrator model. This description is based on a simple multiphase model of the interstellar medium at unresolved scales, where mass and energy flows among the components are explicitly followed by solving a system of ordinary differential equations. Thermal energy from supernovae is injected into the local hot phase, so as to avoid that it is promptly radiated away. A kinetic feedback prescription generates the massive outflows needed to avoid the overproduction of stars. We use two sets of zoomed-in initial conditions of isolated cosmological haloes with masses (2-3) × 1012M⊙, both available at several resolution levels. In all cases we obtain spiral galaxies with small bulge-over-total stellar mass ratios (B/T ~ 0.2), extended stellar and gas discs, flat rotation curves and realistic values of stellar masses. Gas profiles are relatively flat, molecular gas is found to dominate at the centre of galaxies, with star formation rates following the observed Schmidt-Kennicutt relation. Stars kinematically belonging to the bulge form early, while disc stars show a clear inside-out formation pattern and mostly form after redshift z = 2. However, the baryon conversion efficiencies in our simulations differ from the relation given by Moster et al. at a 3σ level, thus indicating that our stellar discs are still too massive for the dark matter halo in which they reside. Results are found to be remarkably stable against resolution. This further demonstrates the feasibility of carrying out simulations producing a realistic population of galaxies within representative cosmological volumes, at a relatively modest resolution. [ABSTRACT FROM AUTHOR]

We compare predictions of cooled masses and cooling rates from three stripped-down semi-analytic models (SAMs) of galaxy formation with the results of N-body+Smoothed Particle Hydrodynamics (SPH) simulations with gas particle mass of 3.9 × 106 h−1 M⊙, where radiative cooling of a gas of primordial composition is implemented. We also run a simulation where cooling is switched on at redshift ∼2, in order to test cooling models in a regime in which their approximations are expected to be valid. We confirm that cooling models implemented in SAMs are able to predict the amount of cooled mass at z = 0 to within ∼20 per cent. However, some relevant discrepancies are found. (i) When the contribution from poorly resolved haloes is subtracted out, SAMs tend to underpredict by ∼30 per cent the mass that cools in the infall-dominated regime. (ii) At large halo masses, SAMs tend to overpredict cooling rates, though the numerical result may be affected by the use of a standard version of SPH. (iii) As found in our previous work, cooling rates are found to be significantly affected by model details: simulations disfavour models with large cores and with quenching of cooling at major mergers. (iv) When cooling is switched on at z ∼ 2, cold gas accumulates very quickly in the simulated haloes. This accumulation is reproduced by SAMs with varying degrees of accuracy. [ABSTRACT FROM AUTHOR]

We present a new multiphase subresolution model for star formation and feedback in smoothed particle hydrodynamics (SPH) numerical simulations of galaxy formation. Our model, called MUPPI (MUlti-Phase Particle Integrator), describes each gas particle as a multiphase system, with cold and hot gas phases, coexisting in pressure equilibrium and a stellar component. Cooling of the hot tenuous gas phase feeds the cold gas phase. We compute the cold gas molecular fraction using the phenomenological relation of Blitz & Rosolowsky between this fraction and the external disc pressure, which we identify with the SPH pressure. Stars are formed out of molecular gas with a given efficiency, which scales with the dynamical time of the cold phase. Our prescription for star formation is not based on imposing the Schmidt–Kennicutt relation, which is instead naturally produced by MUPPI. Energy from supernova explosions is deposited partly into the hot phase of the gas particles, and partly to that of neighbouring particles. Mass and energy flows among the different phases of each particle are described by a set of ordinary differential equations which we explicitly integrate for each gas particle, instead of relying on equilibrium solutions. This system of equations also includes the response of the multiphase structure to energy changes associated to the thermodynamics of the gas. Our model has an intrinsically runaway behaviour: energy from supernovae increases gas pressure which increases in turn the star formation rate through the molecular fraction. This runaway is stabilized in simulations by the hydrodynamic response of the gas: when it receives enough energy, it expands thereby decreasing its pressure. We apply our model to two isolated disc galaxy simulations and two spherical cooling flows. MUPPI is able to reproduce the Schmidt–Kennicutt relation for disc galaxies. It also reproduces the basic properties of the interstellar medium in disc galaxies, the surface densities of cold and molecular gas, of stars and of star formation rate, the vertical velocity dispersion of cold clouds and the flows connected to the galactic fountains. Quite remarkably, MUPPI also provides efficient stellar feedback without the need to include a scheme of kinetic energy feedback. [ABSTRACT FROM AUTHOR]

We present cosmological zoom-in hydrodynamical simulations for the formation of disc galaxies, implementing dust evolution and dust promoted cooling of hot gas. We couple an improved version of our previous treatment of dust evolution, which adopts the two-size approximation to estimate the grain-size distribution, with the MUPPI star formation and feedback subresolution model. Our dust evolution model follows carbon and silicate dust separately. To distinguish differences induced by the chaotic behaviour of simulations from those genuinely due to different simulation set-up, we run each model six times, after introducing tiny perturbations in the initial conditions. With this method, we discuss the role of various dust-related physical processes and the effect of a few possible approximations adopted in the literature. Metal depletion and dust cooling affect the evolution of the system, causing substantial variations in its stellar, gas, and dust content. We discuss possible effects on the Spectral Energy Distribution of the significant variations of the size distribution and chemical composition of grains, as predicted by our simulations during the evolution of the galaxy. We compare dust surface density, dust-to-gas ratio, and small-to-large grain mass ratio as a function of galaxy radius and gas metallicity predicted by our fiducial run with recent observational estimates for three disc galaxies of different masses. The general agreement is good, in particular taking into account that we have not adjusted our model for this purpose. [ABSTRACT FROM AUTHOR]

We investigate the accuracy of weak lensing simulations by comparing the results of five independently developed lensing simulation codes run on the same input N -body simulation. Our comparison focuses on the lensing convergence maps produced by the codes, and in particular on the corresponding PDFs, power spectra, and peak counts. We find that the convergence power spectra of the lensing codes agree to |$\lesssim 2{{\ \rm per\ cent}}$| out to scales ℓ ≈ 4000. For lensing peak counts, the agreement is better than |$5{{\ \rm per\ cent}}$| for peaks with signal-to-noise ≲ 6. We also discuss the systematic errors due to the Born approximation, line-of-sight discretization, particle noise, and smoothing. The lensing codes tested deal in markedly different ways with these effects, but they none-the-less display a satisfactory level of agreement. Our results thus suggest that systematic errors due to the operation of existing lensing codes should be small. Moreover their impact on the convergence power spectra for a lensing simulation can be predicted given its numerical details, which may then serve as a validation test. [ABSTRACT FROM AUTHOR]

We present simulations of galaxy formation, based on the gadget -3 code, in which a sub-resolution model for star formation and stellar feedback is interfaced with a new model for active galactic nucleus (AGN) feedback. Our sub-resolution model describes a multiphase interstellar medium (ISM), accounting for hot and cold gas within the same resolution element: we exploit this feature to investigate the impact of coupling AGN feedback energy to the different phases of the ISM over cosmic time. Our fiducial model considers that AGN feedback energy coupling is driven by the covering factors of the hot and cold phases. We perform a suite of cosmological hydrodynamical simulations of disc galaxies (⁠|$M_{\rm halo,DM} \simeq 2 \times 10^{12}$|  M⊙, at z  = 0), to investigate: (i) the effect of different ways of coupling AGN feedback energy to the multiphase ISM; (ii) the impact of different prescriptions for gas accretion (i.e. only cold gas, both cold and hot gas, with the additional possibility of limiting gas accretion from cold gas with high angular momentum); (iii) how different models of gas accretion and coupling of AGN feedback energy affect the coevolution of supermassive black holes (BHs) and their host galaxy. We find that at least a share of the AGN feedback energy has to couple with the diffuse gas, in order to avoid an excessive growth of the BH mass. When the BH only accretes cold gas, it experiences a growth that is faster than in the case in which both cold and hot gas are accreted. If the accretion of cold gas with high angular momentum is reduced, the BH mass growth is delayed, the BH mass at z  = 0 is reduced by up to an order of magnitude, and the BH is prevented from accreting below z ≲ 2, when the galaxy disc forms. [ABSTRACT FROM AUTHOR]

We study the accuracy of several approximate methods for gravitational dynamics in terms of halo power spectrum multipoles and their estimated covariance matrix. We propagate the differences in covariances into parameter constraints related to growth rate of structure, Alcock–Paczynski distortions, and biasing. We consider seven methods in three broad categories: algorithms that solve for halo density evolution deterministically using Lagrangian trajectories (ICE–COLA, pinocchio, and peakpatch), methods that rely on halo assignment schemes on to dark matter overdensities calibrated with a target N-body run (halogen, patchy), and two standard assumptions about the full density probability distribution function (Gaussian and lognormal). We benchmark their performance against a set of three hundred N-body simulations, running similar sets of approximate simulations with matched initial conditions, for each method. We find that most methods reproduce the monopole to within |$5{{\ \rm per\ cent}}$|⁠, while residuals for the quadrupole are sometimes larger and scale dependent. The variance of the multipoles is typically reproduced within |$10{{\ \rm per\ cent}}$|⁠. Overall, we find that covariances built from approximate simulations yield errors on model parameters within |$10{{\ \rm per\ cent}}$| of those from the N-body-based covariance. [ABSTRACT FROM AUTHOR]

We perform a suite of cosmological hydrodynamical simulations of disc galaxies, with zoomed-in initial conditions leading to the formation of a halo of mass |$M_{\rm halo, \, DM} \simeq 2 \times 10^{12}$|  M⊙ at redshift z  = 0. These simulations aim at investigating the chemical evolution and the distribution of metals in a disc galaxy, and at quantifying the effect of (i) the assumed IMF, (ii) the adopted stellar yields, and (iii) the impact of binary systems originating SNe Ia on the process of chemical enrichment. We consider either a Kroupa, Tout & Gilmore (1993) or a more top-heavy Kroupa (2001) IMF, two sets of stellar yields and different values for the fraction of binary systems suitable to give rise to SNe Ia. We investigate stellar ages, SN rates, stellar and gas metallicity gradients, and stellar α-enhancement in simulations, and compare predictions with observations. We find that a Kroupa et al. (1993) IMF has to be preferred when modelling late-type galaxies in the local Universe. On the other hand, the comparison of stellar metallicity profiles and α-enhancement trends with observations of Milky Way stars shows a better agreement when a Kroupa (2001) IMF is assumed. Comparing the predicted SN rates and stellar α-enhancement with observations supports a value for the fraction of binary systems producing SNe Ia of 0.03, at least for late-type galaxies and for the considered IMFs. Adopted stellar yields are crucial in regulating cooling and star formation, and in determining patterns of chemical enrichment for stars, especially for those located in the galaxy bulge. [ABSTRACT FROM AUTHOR]

We investigate the impact of galactic outflowmodelling on the formation and evolution of a disc galaxy, by performing a suite of cosmological simulations with zoomed-in initial conditions (ICs) of a Milky Way-sized halo. We verify how sensitive the general properties of the simulated galaxy are to the way in which stellar feedback triggered outflows are implemented, keeping ICs, simulation code and star formation (SF) model all fixed. We present simulations that are based on a version of the GADGET3 code where our sub-resolution model is coupled with an advanced implementation of smoothed particle hydrodynamics that ensures a more accurate fluid sampling and an improved description of gas mixing and hydrodynamical instabilities. We quantify the strong interplay between the adopted hydrodynamic scheme and the sub-resolution model describing SF and feedback. We consider four different galactic outflow models, including the one introduced by Dalla Vecchia & Schaye (2012) and a scheme that is inspired by the Springel & Hernquist (2003) model. We find that the sub-resolution prescriptions adopted to generate galactic outflows are the main shaping factor of the stellar disc component at low redshift. The key requirement that a feedback model must have to be successful in producing a disc-dominated galaxy is the ability to regulate the high-redshift SF (responsible for the formation of the bulge component), the cosmological infall of gas from the large-scale environment, and gas fall-back within the galactic radius at low redshift, in order to avoid a too high SF rate at z = 0. [ABSTRACT FROM AUTHOR]

We study the luminosity function (LF) and the radial distribution of satellite galaxies within Milky Way (MW) sized haloes as predicted in cold dark matter based models of galaxy formation, making use of numerical N-body techniques as well as three different semi-analytic models (SAMs) galaxy formation codes. We extract merger trees from very high-resolution dissipationless simulations of four Galaxy-sized DM haloes, and use these as common input for the SAMs. We present a detailed comparison of our predictions with the observational data recently obtained on the MW satellite LF. We find that SAMs with rather standard astrophysical ingredients are able to reproduce the observed LF over six orders of magnitude in luminosity, down to magnitudes as faint as . We also perform a comparison with the actual observed number of satellites as a function of luminosity, by applying the selection criteria of the SDSS survey to our simulations instead of correcting the observations for incompleteness. Using this approach, we again find good agreement for both the luminosity and radial distributions of MW satellites. We investigate which physical processes in our models are responsible for shaping the predicted satellite LF, and find that tidal destruction, suppression of gas infall by a photoionizing background, and supernova feedback all make important contributions. We conclude that the number and luminosity of MW satellites can be naturally accounted for within the (Λ)cold dark matter paradigm, and this should no longer be considered a problem. [ABSTRACT FROM AUTHOR]