Search

We investigate the impact of galactic outflow modelling on the formation and evolution of a disc galaxy, by performing a suite of cosmological simulations with zoomed-in initial conditions 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 initial conditions, 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 and Schaye (2012) and a scheme that is inspired by the Springel and 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$., Comment: Published in MNRAS. Figure 18 is a duplicate of Fig. 12 (left-hand panel) in the published version of the paper, due to an error occurred in the proof editing process. The real Fig. 18 is missing in the published version. This problem is not present here in the arxiv submission. Corrected a few typos to match the published version

We examine the properties of the galaxies and dark matter haloes residing in the cluster infall region surrounding the simulated $\Lambda$CDM galaxy cluster studied by Elahi et al. (2016) at z=0. The $1.1\times10^{15}h^{-1}\text{M}_{\odot}$ galaxy cluster has been simulated with eight different hydrodynamical codes containing a variety of hydrodynamic solvers and subgrid schemes. All models completed a dark-matter only, non-radiative and full-physics run from the same initial conditions. The simulations contain dark matter and gas with mass resolution $m_{\text{DM}}=9.01\times 10^8h^{-1}\text{M}_{\odot}$ and $m_{\text{gas}}=1.9\times 10^8h^{-1}\text{M}_{\odot}$ respectively. We find that the synthetic cluster is surrounded by clear filamentary structures that contain ~60% of haloes in the infall region with mass ~$10^{12.5} - 10^{14} h^{-1}\text{M}_{\odot}$, including 2-3 group-sized haloes ($> 10^{13}h^{-1}\text{M}_{\odot}$). However, we find that only ~10% of objects in the infall region are subhaloes residing in haloes, which may suggest that there is not much ongoing preprocessing occurring in the infall region at z=0. By examining the baryonic content contained within the haloes, we also show that the code-to-code scatter in stellar fraction across all halo masses is typically ~2 orders of magnitude between the two most extreme cases, and this is predominantly due to the differences in subgrid schemes and calibration procedures that each model uses. Models that do not include AGN feedback typically produce too high stellar fractions compared to observations by at least ~1 order of magnitude., Comment: 14 pages, 9 figures, 1 table; accepted for publication in MNRAS

Building on the initial results of the nIFTy simulated galaxy cluster comparison, we compare and contrast the impact of baryonic physics with a single massive galaxy cluster, run with 11 state-of-the-art codes, spanning adaptive mesh, moving mesh, classic and modern SPH approaches. For each code represented we have a dark matter only (DM) and non-radiative (NR) version of the cluster, as well as a full physics (FP) version for a subset of the codes. We compare both radial mass and kinematic profiles, as well as global measures of the cluster (e.g. concentration, spin, shape), in the NR and FP runs with that in the DM runs. Our analysis reveals good consistency (<= 20%) between global properties of the cluster predicted by different codes when integrated quantities are measured within the virial radius R200. However, we see larger differences for quantities within R2500, especially in the FP runs. The radial profiles reveal a diversity, especially in the cluster centre, between the NR runs, which can be understood straightforwardly from the division of codes into classic SPH and non-classic SPH (including the modern SPH, adaptive and moving mesh codes); and between the FP runs, which can also be understood broadly from the division of codes into those that include AGN feedback and those that do not. The variation with respect to the median is much larger in the FP runs with different baryonic physics prescriptions than in the NR runs with different hydrodynamics solvers., Comment: 24 pages, 13 figures, MNRAS submitted

We present a model for the seeding and evolution of magnetic fields in galaxies by supernovae (SN). SN explosions during galaxy assembly provide seed fields, which are subsequently amplified by compression, shear flows and random motions. Our model explains the origin of microG magnetic fields within galactic structures. We implement our model in the MHD version of the cosmological simulation code Gadget-3 and couple it with a multi-phase description of the interstellar medium. We perform simulations of Milky Way-like galactic halo formation and analyze the distribution and strength of the magnetic field. We investigate the intrinsic rotation measure (RM) evolution and find RM values exceeding 1000 rad/m*m at high redshifts and RM values around 10 rad/m*m at present-day. We compare our simulations to a limited set of observational data points and find encouraging similarities. In our model, galactic magnetic fields are a natural consequence of the very basic processes of star formation and galaxy assembly., Comment: 2 pages, proceedings of IAU symposium 315

We examine subhaloes and galaxies residing in a simulated LCDM galaxy cluster ($M^{\rm crit}_{200}=1.1\times10^{15}M_\odot/h$) produced by hydrodynamical codes ranging from classic Smooth Particle Hydrodynamics (SPH), newer SPH codes, adaptive and moving mesh codes. These codes use subgrid models to capture galaxy formation physics. We compare how well these codes reproduce the same subhaloes/galaxies in gravity only, non-radiative hydrodynamics and full feedback physics runs by looking at the overall subhalo/galaxy distribution and on an individual objects basis. We find the subhalo population is reproduced to within $\lesssim10\%$ for both dark matter only and non-radiative runs, with individual objects showing code-to-code scatter of $\lesssim0.1$ dex, although the gas in non-radiative simulations shows significant scatter. Including feedback physics significantly increases the diversity. Subhalo mass and $V_{max}$ distributions vary by $\approx20\%$. The galaxy populations also show striking code-to-code variations. Although the Tully-Fisher relation is similar in almost all codes, the number of galaxies with $10^{9}M_\odot/h\lesssim M_*\lesssim 10^{12}M_\odot/h$ can differ by a factor of 4. Individual galaxies show code-to-code scatter of $\sim0.5$ dex in stellar mass. Moreover, strong systematic differences exist, with some codes producing galaxies $70\%$ smaller than others. The diversity partially arises from the inclusion/absence of AGN feedback. Our results combined with our companion papers demonstrate that subgrid physics is not just subject to fine-tuning, but the complexity of building galaxies in all environments remains a challenge. We argue even basic galaxy properties, such as the stellar mass to halo mass, should be treated with errors bars of $\sim0.2-0.4$ dex., Comment: 17 pages (+4 page appendix), 16 figures, 2 tables; accepted for publication in MNRAS

We present an on-the-fly geometrical approach for shock detection and Mach number calculation in simulations employing smoothed particle hydrodynamics (SPH). We utilize pressure gradients to select shock candidates and define up- and downstream positions. We obtain hydrodynamical states in the up- and downstream regimes with a series of normal and inverted kernel weightings parallel and perpendicular to the shock normals. Our on-the-fly geometrical Mach detector incorporates well within the SPH formalism and has low computational cost. We implement our Mach detector into the simulation code GADGET and alongside many SPH improvements. We test our shock finder in a sequence of shock-tube tests with successively increasing Mach numbers exceeding by far the typical values inside galaxy clusters. For all shocks, we resolve the shocks well and the correct Mach numbers are assigned. An application to a strong magnetized shock-tube gives stable results in full magnetohydrodynamic set-ups. We simulate a merger of two idealized galaxy clusters and study the shock front. Shock structures within the merging clusters as well as the cluster shock are well-captured by our algorithm and assigned correct Mach numbers., Comment: 9 pages, 9 figures, accepted to MNRAS

We have simulated the formation of a massive galaxy cluster (M$_{200}^{\rm crit}$ = 1.1$\times$10$^{15}h^{-1}M_{\odot}$) in a $\Lambda$CDM universe using 10 different codes (RAMSES, 2 incarnations of AREPO and 7 of GADGET), modeling hydrodynamics with full radiative subgrid physics. These codes include Smoothed-Particle Hydrodynamics (SPH), spanning traditional and advanced SPH schemes, adaptive mesh and moving mesh codes. Our goal is to study the consistency between simulated clusters modeled with different radiative physical implementations - such as cooling, star formation and AGN feedback. We compare images of the cluster at $z=0$, global properties such as mass, and radial profiles of various dynamical and thermodynamical quantities. We find that, with respect to non-radiative simulations, dark matter is more centrally concentrated, the extent not simply depending on the presence/absence of AGN feedback. The scatter in global quantities is substantially higher than for non-radiative runs. Intriguingly, adding radiative physics seems to have washed away the marked code-based differences present in the entropy profile seen for non-radiative simulations in Sembolini et al. (2015): radiative physics + classic SPH can produce entropy cores. Furthermore, the inclusion/absence of AGN feedback is not the dividing line -as in the case of describing the stellar content- for whether a code produces an unrealistic temperature inversion and a falling central entropy profile. However, AGN feedback does strongly affect the overall stellar distribution, limiting the effect of overcooling and reducing sensibly the stellar fraction., Comment: 20 pages, 13 figures, submitted to MNRAS

By means of zoom-in hydrodynamic simulations we quantify the amount of neutral hydrogen (HI) hosted by groups and clusters of galaxies. Our simulations, which are based on an improved formulation of smoothed particle hydrodynamics (SPH), include radiative cooling, star formation, metal enrichment and supernova feedback, and can be split in two different groups, depending on whether feedback from active galactic nuclei (AGN) is turned on or off. Simulations are analyzed to account for HI self-shielding and the presence of molecular hydrogen. We find that the mass in neutral hydrogen of dark matter halos monotonically increases with the halo mass and can be well described by a power-law of the form $M_{\rm HI}(M,z)\propto M^{3/4}$. Our results point out that AGN feedback reduces both the total halo mass and its HI mass, although it is more efficient in removing HI. We conclude that AGN feedback reduces the neutral hydrogen mass of a given halo by $\sim50\%$, with a weak dependence on halo mass and redshift. The spatial distribution of neutral hydrogen within halos is also affected by AGN feedback, whose effect is to decrease the fraction of HI that resides in the halo inner regions. By extrapolating our results to halos not resolved in our simulations we derive astrophysical implications from the measurements of $\Omega_{\rm HI}(z)$: halos with circular velocities larger than $\sim25~{\rm km/s}$ are needed to host HI in order to reproduce observations. We find that only the model with AGN feedback is capable of reproducing the value of $\Omega_{\rm HI}b_{\rm HI}$ derived from available 21cm intensity mapping observations., Comment: 19 pages, 12 figures. Accepted for publication in MNRAS

We have simulated the formation of a galaxy cluster in a $\Lambda$CDM universe using twelve different codes modeling only gravity and non-radiative hydrodynamics (\art, \arepo, \hydra\ and 9 incarnations of GADGET). This range of codes includes particle based, moving and fixed mesh codes as well as both Eulerian and Lagrangian fluid schemes. The various GADGET implementations span traditional and advanced smoothed-particle hydrodynamics (SPH) schemes. The goal of this comparison is to assess the reliability of cosmological hydrodynamical simulations of clusters in the simplest astrophysically relevant case, that in which the gas is assumed to be non-radiative. We compare images of the cluster at $z=0$, global properties such as mass, and radial profiles of various dynamical and thermodynamical quantities. The underlying gravitational framework can be aligned very accurately for all the codes allowing a detailed investigation of the differences that develop due to the various gas physics implementations employed. As expected, the mesh-based codes ART and AREPO form extended entropy cores in the gas with rising central gas temperatures. Those codes employing traditional SPH schemes show falling entropy profiles all the way into the very centre with correspondingly rising density profiles and central temperature inversions. We show that methods with modern SPH schemes that allow entropy mixing span the range between these two extremes and the latest SPH variants produce gas entropy profiles that are essentially indistinguishable from those obtained with grid based methods., Comment: 21 pages, 13 figures, 4 tables - submitted to MNRAS

We present an implementation of thermal conduction including the anisotropic effects of magnetic fields for SPH. The anisotropic thermal conduction is mainly proceeding parallel to magnetic fields and suppressed perpendicular to the fields. We derive the SPH formalism for the anisotropic heat transport and solve the corresponding equation with an implicit conjugate gradient scheme. We discuss several issues of unphysical heat transport in the cases of extreme ansiotropies or unmagnetized regions and present possible numerical workarounds. We implement our algorithm into the GADGET code and study its behaviour in several test cases. In general, we reproduce the analytical solutions of our idealised test problems, and obtain good results in cosmological simulations of galaxy cluster formations. Within galaxy clusters, the anisotropic conduction produces a net heat transport similar to an isotropic Spitzer conduction model with an efficiency of one per cent. In contrast to isotropic conduction our new formalism allows small-scale structure in the temperature distribution to remain stable, because of their decoupling caused by magnetic field lines. Compared to observations, isotropic conduction with more than 10 per cent of the Spitzer value leads to an oversmoothed temperature distribution within clusters, while the results obtained with anisotropic thermal conduction reproduce the observed temperature fluctuations well. A proper treatment of heat transport is crucial especially in the outskirts of clusters and also in high density regions. It's connection to the local dynamical state of the cluster also might contribute to the observed bimodal distribution of cool core and non cool core clusters. Our new scheme significantly advances the modelling of thermal conduction in numerical simulations and overall gives better results compared to observations., Comment: 25 pages, 18 figures

We extend the description of the isotropic and anisotropic random component of the small-scale magnetic field within the existing magnetic field model of the Milky Way from Jansson&Farrar, by including random realizations of the small-scale component. Using a magnetic-field power spectrum with Gaussian random fields, the NE2001 model for the thermal electrons and the Galactic cosmic-ray electron distribution from the current GALPROP model we derive full-sky maps for the total and polarized synchrotron intensity as well as the Faraday rotation-measure distribution. While previous work assumed that small-scale fluctuations average out along the line-of-sight or which only computed ensemble averages of random fields, we show that these fluctuations need to be carefully taken into account. Comparing with observational data we obtain not only good agreement with 408 MHz total and WMAP7 22 GHz polarized intensity emission maps, but also an improved agreement with Galactic foreground rotation-measure maps and power spectra, whose amplitude and shape strongly depend on the parameters of the random field. We demonstrate that a correlation length of ~220pc (50pc being a 5\sigma lower limit) is needed to match the slope of the observed power spectrum of Galactic foreground rotation-measure maps. Using multiple realizations allows us also to infer errors on individual observables. We find that previously-used amplitudes for random and anisotropic random magnetic field components need to be rescaled by factors of ~0.3 and 0.6 to account for the new small-scale contributions. Our model predicts a rotation measure of -2.8\pm7.1 rad/m^2 and 4.4\pm11.0 rad/m^2 for the north and south Galactic poles respectively, in good agreement with observations. Applying our model to deflections of ultra-high-energy cosmic rays we infer a mean deflection of ~3.5\pm1.1 degree for 60 EeV protons arriving from CenA., Comment: 28 pages, 12 figures, 2 tables

We present a model for the seeding and evolution of magnetic fields in protogalaxies. Supernova (SN) explosions during the assembly of a protogalaxy provide magnetic seed fields, which are subsequently amplified by compression, shear flows and random motions. We implement the model into the MHD version of the cosmological N-body / SPH simulation code GADGET and we couple the magnetic seeding directly to the underlying multi-phase description of star formation. We perform simulations of Milky Way-like galactic halo formation using a standard LCDM cosmology and analyse the strength and distribution of the subsequent evolving magnetic field. A dipole-shape divergence-free magnetic field is injected at a rate of 10^{-9}G / Gyr within starforming regions, given typical dimensions and magnetic field strengths in canonical SN remnants. Subsequently, the magnetic field strength increases exponentially on timescales of a few ten million years. At redshift z=0, the entire galactic halo is magnetized and the field amplitude is of the order of a few $\mu$G in the center of the halo, and 10^{-9} G at the virial radius. Additionally, we analyse the intrinsic rotation measure (RM) of the forming galactic halo over redshift. The mean halo intrinsic RM peaks between redshifts z=4 and z=2 and reaches absolute values around 1000 rad m^{-2}. While the halo virializes towards redshift z=0, the intrinsic RM values decline to a mean value below 10 rad m^{-2}. At high redshifts, the distribution of individual starforming, and thus magnetized regions is widespread. In our model for the evolution of galactic magnetic fields, the seed magnetic field amplitude and distribution is no longer a free parameter, but determined self-consistently by the star formation process occuring during the formation of cosmic structures., Comment: 13 pages, 14 figures, accepted to MNRAS after moderate revision

We study the possible magnetization of cosmic voids by void galaxies. Recently, observations revealed isolated starforming galaxies within the voids. Furthermore, a major fraction of a voids volume is expected to be filled with magnetic fields of a minimum strength of about $10^{-15}$ G on Mpc scales. We estimate the transport of magnetic energy by cosmic rays (CR) from the void galaxies into the voids. We assume that CRs and winds are able to leave small isolated void galaxies shortly after they assembled, and then propagate within the voids. For a typical void, we estimate the magnetic field strength and volume filling factor depending on its void galaxy population and possible contributions of strong active galactic nuclei (AGN) which border the voids. We argue that the lower limit on the void magnetic field can be recovered, if a small fraction of the magnetic energy contained in the void galaxies or void bordering AGNs is distributed within the voids., Comment: 5 pages, 2 figures, accepted to MNRAS Letters after moderate revision

We present simulations of the compact galaxy group Stephan's Quintet (SQ) including magnetic fields, performed with the N-body/smoothed particle hydrodynamics (SPH) code \textsc{Gadget}. The simulations include radiative cooling, star formation and supernova feedback. Magnetohydrodynamics (MHD) is implemented using the standard smoothed particle magnetohydrodynamics (SPMHD) method. We adapt two different initial models for SQ based on Renaud et al. and Hwang et al., both including four galaxies (NGC 7319, NGC 7320c, NGC 7318a and NGC 7318b). Additionally, the galaxies are embedded in a magnetized, low density intergalactic medium (IGM). The ambient IGM has an initial magnetic field of $10^{-9}$ G and the four progenitor discs have initial magnetic fields of $10^{-9} - 10^{-7}$ G. We investigate the morphology, regions of star formation, temperature, X-ray emission, magnetic field structure and radio emission within the two different SQ models. In general, the enhancement and propagation of the studied gaseous properties (temperature, X-ray emission, magnetic field strength and synchrotron intensity) is more efficient for the SQ model based on Renaud et al., whose galaxies are more massive, whereas the less massive SQ model based on Hwang et al. shows generally similar effects but with smaller efficiency. We show that the large shock found in observations of SQ is most likely the result of a collision of the galaxy NGC 7318b with the IGM. This large group-wide shock is clearly visible in the X-ray emission and synchrotron intensity within the simulations of both SQ models. The order of magnitude of the observed synchrotron emission within the shock front is slightly better reproduced by the SQ model based on Renaud et al., whereas the distribution and structure of the synchrotron emission is better reproduced by the SQ model based on Hwang et al.., Comment: 20 pages, 15 figures, accepted to MNRAS

An analytical model predicting the growth rates, the absolute growth times and the saturation values of the magnetic field strength within galactic haloes is presented. The analytical results are compared to cosmological MHD simulations of Milky-Way like galactic halo formation performed with the N-body / \textsc{Spmhd} code \textsc{Gadget}. The halo has a mass of $\approx{}3\cdot{}10^{12}$ $M_{\odot}$ and a virial radius of $\approx{}$270 kpc. The simulations in a $\Lambda$CDM cosmology also include radiative cooling, star formation, supernova feedback and the description of non-ideal MHD. A primordial magnetic seed field ranging from $10^{-10}$ to $10^{-34}$ G in strength agglomerates together with the gas within filaments and protohaloes. There, it is amplified within a couple of hundred million years up to equipartition with the corresponding turbulent energy. The magnetic field strength increases by turbulent small-scale dynamo action. The turbulence is generated by the gravitational collapse and by supernova feedback. Subsequently, a series of halo mergers leads to shock waves and amplification processes magnetizing the surrounding gas within a few billion years. At first, the magnetic energy grows on small scales and then self-organizes to larger scales. Magnetic field strengths of $\approx{}10^{-6}$ G are reached in the center of the halo and drop to $\approx{}10^{-9}$ G in the IGM. Analyzing the saturation levels and growth rates, the model is able to describe the process of magnetic amplification notably well and confirms the results of the simulations., Comment: 13 pages, 12 figures, accepted to MNRAS after minor revision

We investigate the magnetic field evolution in a series of galaxy minor mergers using the N-body/smoothed particle hydrodynamics (SPH) code \textsc{Gadget}. The simulations include the effects of radiative cooling, star formation and supernova feedback. Magnetohydrodynamics (MHD) is implemented using the SPH method. We present 32 simulations of binary mergers of disc galaxies with mass ratios of 2:1 up to 100:1, whereby we have additionally varied the initial magnetic field strengths, disc orientations and resolutions. We investigate the amplification of a given initial magnetic field within the galaxies and an ambient intergalactic medium (IGM) during the interaction. We find that the magnetic field strengths of merger remnants with mass ratios up to 10:1 saturate at a common value of several $\mu$G. For higher mass ratios, the field strength saturates at lower values. The saturation values correspond to the equipartition value of magnetic and turbulent energy density. The initial magnetization, disc orientation and numerical resolution show only minor effects on the saturation value of the magnetic field. We demonstrate that a higher impact energy of the progenitor galaxies leads to a more efficient magnetic field amplification. The magnetic and turbulent energy densities are higher for larger companion galaxies, consistent with the higher impact energy supplied to the system. We present a detailed study of the evolution of the temperature and the bolometric X-ray luminosity within the merging systems. Thereby we find that magnetic fields cause a more efficient increase of the IGM temperature and the corresponding IGM X-ray luminosity after the first encounter. However, the presence of magnetic fields does not enhance the total X-ray luminosity. Generally, the final value of the X-ray luminosity is even clearly lower for higher initial magnetic fields., Comment: 20 pages, 12 figures. Submitted to MNRAS