Hydrodynamic cosmological simulations at present usually employ either the Lagrangian smoothed particle hydrodynamics (SPH) technique or Eulerian hydrodynamics on a Cartesian mesh with (optional) ...adaptive mesh refinement (AMR). Both of these methods have disadvantages that negatively impact their accuracy in certain situations, for example the suppression of fluid instabilities in the case of SPH, and the lack of Galilean invariance and the presence of overmixing in the case of AMR. We here propose a novel scheme which largely eliminates these weaknesses. It is based on a moving unstructured mesh defined by the Voronoi tessellation of a set of discrete points. The mesh is used to solve the hyperbolic conservation laws of ideal hydrodynamics with a finite-volume approach, based on a second-order unsplit Godunov scheme with an exact Riemann solver. The mesh-generating points can in principle be moved arbitrarily. If they are chosen to be stationary, the scheme is equivalent to an ordinary Eulerian method with second-order accuracy. If they instead move with the velocity of the local flow, one obtains a Lagrangian formulation of continuum hydrodynamics that does not suffer from the mesh distortion limitations inherent in other mesh-based Lagrangian schemes. In this mode, our new method is fully Galilean invariant, unlike ordinary Eulerian codes, a property that is of significant importance for cosmological simulations where highly supersonic bulk flows are common. In addition, the new scheme can adjust its spatial resolution automatically and continuously, and hence inherits the principal advantage of SPH for simulations of cosmological structure growth. The high accuracy of Eulerian methods in the treatment of shocks is also retained, while the treatment of contact discontinuities improves. We discuss how this approach is implemented in our new code arepo, both in 2D and in 3D, and is parallelized for distributed memory computers. We also discuss techniques for adaptive refinement or de-refinement of the unstructured mesh. We introduce an individual time-step approach for finite-volume hydrodynamics, and present a high-accuracy treatment of self-gravity for the gas that allows the new method to be seamlessly combined with a high-resolution treatment of collisionless dark matter. We use a suite of test problems to examine the performance of the new code and argue that the hydrodynamic moving-mesh scheme proposed here provides an attractive and competitive alternative to current SPH and Eulerian techniques.
This review discusses smoothed particle hydrodynamics (SPH) in the astrophysical context, with a focus on inviscid gas dynamics. The particle-based SPH technique allows an intuitive and simple ...formulation of hydrodynamics that has excellent conservation properties and can be coupled to self-gravity with high accuracy. The Lagrangian character of SPH allows it to automatically adjust its resolution to the clumping of matter, a property that makes the scheme ideal for many application areas in astrophysics, where often a large dynamic range in density is encountered. We discuss the derivation of the basic SPH equations in their modern formulation, and give an overview about extensions of SPH developed to treat physics such as radiative transfer, thermal conduction, relativistic dynamics, or magnetic fields. We also briefly describe some of the most important applications areas of SPH in astrophysical research. Finally, we provide a critical discussion of the accuracy of SPH for different hydrodynamical problems, including measurements of its convergence rate for important classes of problems. PUBLICATION ABSTRACT
ABSTRACT
An ultralight bosonic particle of mass around $10^{-22}\, \mathrm{eV}/c^2$ is of special interest as a dark matter candidate, as it both has particle physics motivations, and may give rise ...to notable differences in the structures on highly non-linear scales due to the manifestation of quantum-physical wave effects on macroscopic scales, which could address a number of contentious small-scale tensions in the standard cosmological model, ΛCDM. Using a spectral technique, we here discuss simulations of such fuzzy dark matter (FDM), including the full non-linear wave dynamics, with a comparatively large dynamic range and for larger box sizes than considered previously. While the impact of suppressed small-scale power in the initial conditions associated with FDM has been studied before, the characteristic FDM dynamics are often neglected; in our simulations, we instead show the impact of the full non-linear dynamics on physical observables. We focus on the evolution of the matter power spectrum, give first results for the FDM halo mass function directly based on full FDM simulations, and discuss the computational challenges associated with the FDM equations. FDM shows a pronounced suppression of power on small scales relative to cold dark matter (CDM), which can be understood as a damping effect due to ‘quantum pressure’. In certain regimes, however, the FDM power can exceed that of CDM, which may be interpreted as a reflection of order-unity density fluctuations occurring in FDM. In the halo mass function, FDM shows a significant abundance reduction below a characteristic mass scale only. This could in principle alleviate the need to invoke very strong feedback processes in small galaxies to reconcile ΛCDM with the observed galaxy luminosity function, but detailed studies that also include baryons will be needed to ultimately judge the viability of FDM.
ABSTRACT
It is well established that the properties of supermassive black holes (BHs) and their host galaxies are correlated through scaling relations. While hydrodynamical cosmological simulations ...have begun to account for the coevolution of BHs and galaxies, they typically have neglected the BH spin, even though it may play an important role in modulating the growth and feedback of BHs. Here we introduce a new sub-grid model for the BH spin evolution in the moving-mesh code arepo in order to improve the physical faithfulness of the BH modelling in galaxy formation simulations. We account for several different channels of spin evolution, in particular gas accretion through a Shakura–Sunyaev α-disc, chaotic accretion, and BH mergers. For BH feedback, we extend the IllustrisTNG model, which considers two different BH feedback modes, a thermal quasar mode for high accretion states and a kinetic mode for low Eddington ratios, with a self-consistent accounting of spin-dependent radiative efficiencies and thus feedback strength. We find that BHs with a mass $M_{\mbox{{bh}}}\lesssim 10^{8}\, {\rm M}_{\odot }$ reach high spin values as they typically evolve in the coherent gas accretion regime, in which consecutive accretion episodes are aligned. On the other hand, BHs with a mass $M_{\mbox{{bh}}}\gtrsim 10^{8}\, {\rm M}_{\odot }$ have lower spins as BH mergers become more frequent, and their accretion discs fragment due to self-gravity, inducing chaotic accretion. We also explore the hypothesis that the transition between the quasar and kinetic feedback modes is mediated by the accretion mode of the BH disc itself, i.e. the kinetic feedback mode is activated when the disc enters the self-gravity regime instead of by an ad hoc switch tied to the BH mass. We find excellent agreement between the galaxy and BH populations for this approach and the fiducial TNG model with no spin evolution. Furthermore, our new approach alleviates a tension in the galaxy morphology–colour relation of the original TNG model.
We describe a new iterative approach for the realization of equilibrium N-body systems for given density distributions. Our method uses elements of Schwarzschild's technique and of the ...made-to-measure method, but is based on a different principle. Starting with some initial assignment of particle velocities, the difference of the time-averaged density response produced by the particle orbits with respect to the initial density configuration is characterized through a merit function, and a stationary solution of the collisionless Boltzmann equation is found by minimizing this merit function directly by iteratively adjusting the initial velocities. Because the distribution function is in general not unique for a given density structure, we augment the merit function with additional constraints that single out a desired target solution. The velocity adjustment is carried out with a stochastic process in which new velocities are randomly drawn from an approximate solution of the distribution function, but are kept only when they improve the fit. Our method converges rapidly and is flexible enough to allow the construction of solutions with third integrals of motion, including disc galaxies in which radial and vertical dispersions are different. A parallel code for the calculation of compound galaxy models with this new method is made publicly available.
We present a new massively parallel code for N-body and cosmological hydrodynamical simulations of modified gravity models. The code employs a multigrid-accelerated Newton-Gauss-Seidel relaxation ...solver on an adaptive mesh to efficiently solve for perturbations in the scalar degree of freedom of the modified gravity model. As this new algorithm is implemented as a module for the p-gadget3 code, it can at the same time follow the baryonic physics included in p-gadget3, such as hydrodynamics, radiative cooling and star formation. We demonstrate that the code works reliably by applying it to simple test problems that can be solved analytically, as well as by comparing cosmological simulations to results from the literature. Using the new code, we perform the first non-radiative and radiative cosmological hydrodynamical simulations of an f (R)-gravity model. We also discuss the impact of active galactic nucleus feedback on the matter power spectrum, as well as degeneracies between the influence of baryonic processes and modifications of gravity.
ABSTRACT
We study the evolution of magnetic fields generated by charge segregation ahead of ionization fronts during the Epoch of Reionization, and their effects on galaxy formation. We compare this ...magnetic seeding process with the Biermann battery, injection from supernovae, and an imposed seed field at redshift z ≳ 127. Using a suite of self-consistent cosmological and zoom-in simulations based on the Auriga galaxy-formation model, we determine that all mechanisms produce galactic magnetic fields that equally affect galaxy formation, and are nearly indistinguishable at z ≲ 1.5. The former is compatible with observed values, while the latter is correlated with the gas metallicity below a seed-dependent redshift. Low-density gas and haloes below a seed-dependent mass threshold retain memory of the initial magnetic field. We produce synthetic Faraday rotation measure maps, showing that they have the potential to constrain the seeding process, although current observations are not yet sensitive enough. Our results imply that the ad-hoc assumption of a primordial seed field – widely used in galaxy formation simulations but of uncertain physical origin – can be replaced by physically motivated mechanisms for magnetogenesis with negligible impact on galactic properties. Additionally, magnetic fields generated ahead of ionization fronts appear very similar but weaker than those produced by the Biermann battery. Hence, in a realistic scenario where both mechanisms are active, the former will be negligible compared to the latter. Finally, our results highlight that the high-redshift Universe is a fruitful testing ground for our understanding of magnetic fields generation.
ABSTRACT
We introduce the Stars and MUltiphase Gas in GaLaxiEs – SMUGGLE model, an explicit and comprehensive stellar feedback model for the moving-mesh code arepo. This novel sub-resolution model ...resolves the multiphase gas structure of the interstellar medium and self-consistently generates gaseous outflows. The model implements crucial aspects of stellar feedback including photoionization, radiation pressure, energy, and momentum injection from stellar winds and from supernovae. We explore this model in high-resolution isolated simulations of Milky Way like disc galaxies. Stellar feedback regulates star formation to the observed level and naturally captures the establishment of a Kennicutt–Schmidt relation. This result is achieved independent of the numerical mass and spatial resolution of the simulations. Gaseous outflows are generated with average mass loading factors of the order of unity. Strong outflow activity is correlated with peaks in the star formation history of the galaxy with evidence that most of the ejected gas eventually rains down on to the disc in a galactic fountain flow that sustains late-time star formation. Finally, the interstellar gas in the galaxy shows a distinct multiphase distribution with a coexistence of cold, warm, and hot phases.
Abstract
The IllustrisTNG project is a new suite of cosmological magnetohydrodynamical simulations of galaxy formation performed with the arepo code and updated models for feedback physics. Here, we ...introduce the first two simulations of the series, TNG100 and TNG300, and quantify the stellar mass content of about 4000 massive galaxy groups and clusters (1013 ≤ M200c/M⊙ ≤ 1015) at recent times (z ≤ 1). The richest clusters have half of their total stellar mass bound to satellite galaxies, with the other half being associated with the central galaxy and the diffuse intracluster light. Haloes more massive than about 5 × 1014 M⊙ have more diffuse stellar mass outside 100 kpc than within 100 kpc, with power-law slopes of the radial mass density distribution as shallow as the dark matter's ( − 3.5 ≲ α3D ≲ −3). Total halo mass is a very good predictor of stellar mass, and vice versa: at z = 0, the 3D stellar mass measured within 30 kpc scales as ∝(M500c)0.49 with a ∼0.12 dex scatter. This is possibly too steep in comparison to the available observational constraints, even though the abundance of The Next Generation less-massive galaxies ( ≲ 1011 M⊙ in stars) is in good agreement with the measured galaxy stellar mass functions at recent epochs. The 3D sizes of massive galaxies fall too on a tight (∼0.16 dex scatter) power-law relation with halo mass, with $r^{\rm stars}_{\rm 0.5} \propto (M_{\rm 200c})^{0.53}$. Even more fundamentally, halo mass alone is a good predictor for the whole stellar mass profiles beyond the inner few kiloparsecs, and we show how on average these can be precisely recovered given a single-mass measurement of the galaxy or its halo.