-
Prevention is better than cure? Feedback from high specific energy winds in cosmological simulations with Arkenstone
Authors:
Jake S. Bennett,
Matthew C. Smith,
Drummond B. Fielding,
Greg L. Bryan,
Chang-Goo Kim,
Volker Springel,
Lars Hernquist
Abstract:
We deploy the new Arkenstone galactic wind model in cosmological simulations for the first time, allowing us to robustly resolve the evolution and impact of high specific energy winds. In a (25$\,h^{-1}\,$Mpc)$^3$ box we perform a set of numerical experiments that systematically vary the mass and energy loadings of such winds, finding that their energy content is the key parameter controlling the…
▽ More
We deploy the new Arkenstone galactic wind model in cosmological simulations for the first time, allowing us to robustly resolve the evolution and impact of high specific energy winds. In a (25$\,h^{-1}\,$Mpc)$^3$ box we perform a set of numerical experiments that systematically vary the mass and energy loadings of such winds, finding that their energy content is the key parameter controlling the stellar to dark matter mass ratio. Increasing the mass loading, at fixed energy, actually results in mildly enhanced star formation, counter to prevailing wisdom but in agreement with recent analytic models. Of the simple parametrisations that we test, we find that an energy loading that scales inversely with halo mass best matches a wide range of observations, and can do so with mass loadings drastically lower than those in most previous cosmological simulations. In this scenario, much less material is ejected from the interstellar medium. Instead, winds both heat gas in the circumgalactic medium, slowing infall onto the galaxy, and also drive shocks beyond the virial radius, preventing accretion onto the halo in the first place. We have not yet tied the mass and energy loadings to high-resolution simulations (a key goal of the Learning the Universe collaboration); however, we can already report that a much lower fraction of the available supernova energy is needed in preventative galaxy regulation than required by ejective wind feedback models such as IllustrisTNG.
△ Less
Submitted 16 October, 2024;
originally announced October 2024.
-
Anisotropic thermal conduction on a moving mesh for cosmological simulations
Authors:
Rosie Y. Talbot,
Rüdiger Pakmor,
Christoph Pfrommer,
Volker Springel,
Maria Werhahn,
Rebekka Bieri,
Freeke van de Voort
Abstract:
In weakly collisional, strongly magnetised plasmas such as the intracluster medium (ICM), hot accretion flows and the solar corona, the transport of heat and momentum occurs primarily along magnetic field lines. In this paper we present a new scheme for modelling anisotropic thermal conduction which we have implemented in the moving mesh code AREPO. Our implementation uses a semi-implicit time int…
▽ More
In weakly collisional, strongly magnetised plasmas such as the intracluster medium (ICM), hot accretion flows and the solar corona, the transport of heat and momentum occurs primarily along magnetic field lines. In this paper we present a new scheme for modelling anisotropic thermal conduction which we have implemented in the moving mesh code AREPO. Our implementation uses a semi-implicit time integration scheme which works accurately and efficiently with individual timestepping, making the scheme highly suitable for use in cosmological simulations. We apply the scheme to a number of test-problems including the diffusion of a hot patch of gas in a circular magnetic field, the progression of a point explosion in the presence of thermal conduction, and the evolution and saturation of buoyancy instabilities in anisotropically conducting plasmas. We use these idealised tests to demonstrate the accuracy and stability of the solver and highlight the ways in which anisotropic conduction can fundamentally change the behaviour of the system. Finally, we demonstrate the solver's capability when applied to highly non-linear problems with deep timestep hierarchies by performing high-resolution cosmological zoom-in simulations of a galaxy cluster with conduction. We show that anisotropic thermal conduction can have a significant impact on the temperature distribution of the ICM and that whistler suppression may be relevant on cluster scales. The new scheme is, therefore, well suited for future work which will explore the role of anisotropic thermal conduction in a range of astrophysical contexts including the ICM of clusters and the circumgalactic medium of galaxies.
△ Less
Submitted 9 October, 2024;
originally announced October 2024.
-
Moving-mesh non-ideal magnetohydrodynamical simulations of the collapse of cloud cores to protostars
Authors:
Alexander C. Mayer,
Oliver Zier,
Thorsten Naab,
Rüdiger Pakmor,
Paola Caselli,
Alexei V. Ivlev,
Volker Springel,
Stefanie Walch
Abstract:
Magnetic fields have been shown both observationally and through theoretical work to be an important factor in the formation of protostars and their accretion disks. Accurate modelling of the evolution of the magnetic field in low-ionization molecular cloud cores requires the inclusion of non-ideal magnetohydrodynamics (MHD) processes, specifically Ohmic and ambipolar diffusion and the Hall effect…
▽ More
Magnetic fields have been shown both observationally and through theoretical work to be an important factor in the formation of protostars and their accretion disks. Accurate modelling of the evolution of the magnetic field in low-ionization molecular cloud cores requires the inclusion of non-ideal magnetohydrodynamics (MHD) processes, specifically Ohmic and ambipolar diffusion and the Hall effect. These have a profound influence on the efficiency of magnetic removal of angular momentum from protostellar disks and simulations that include them can avoid the `magnetic-braking catastrophe' in which disks are not able to form. However, the impact of the Hall effect, in particular, is complex and remains poorly studied. In this work, we perform a large suite of simulations of the collapse of cloud cores to protostars with several non-ideal MHD chemistry models and initial core geometries using the moving-mesh code {\small AREPO}. We find that the efficiency of angular momentum removal is significantly reduced with respect to ideal MHD, in line with previous results. The Hall effect has a varied influence on the evolution of the disk which depends on the initial orientation of the magnetic field. This extends to the outflows seen in a subset of the models, where this effect can act to enhance or suppress them and open up new outflow channels. We conclude, in agreement with a subset of the previous literature, that the Hall effect is the dominant non-ideal MHD process in some collapse scenarios and thus should be included in simulations of protostellar disk formation.
△ Less
Submitted 9 October, 2024;
originally announced October 2024.
-
Environment matters: stronger magnetic fields in satellite galaxies
Authors:
Maria Werhahn,
Rüdiger Pakmor,
Rebekka Bieri,
Freeke van de Voort,
Rosie Y. Talbot,
Volker Springel
Abstract:
Magnetic fields are ubiquitous in the universe and an important component of the interstellar medium. It is crucial to accurately model and understand their properties in different environments and across all mass ranges to interpret observables related to magnetic fields correctly. However, the assessment of the role of magnetic fields in galaxy evolution is often hampered by limited numerical re…
▽ More
Magnetic fields are ubiquitous in the universe and an important component of the interstellar medium. It is crucial to accurately model and understand their properties in different environments and across all mass ranges to interpret observables related to magnetic fields correctly. However, the assessment of the role of magnetic fields in galaxy evolution is often hampered by limited numerical resolution in cosmological simulations, in particular for satellite galaxies. To this end, we study the magnetic fields in high-resolution cosmological zoom simulations of disk galaxies (with $M_{200}\approx10^{10}$ to $10^{13}\,\mathrm{M_\odot}$) and their satellites within the Auriga galaxy formation model including cosmic rays. We find significantly higher magnetic field strengths in satellite galaxies compared to isolated dwarfs with a similar mass or star-formation rate, in particular after they had their first close encounter with their host galaxy. While this result is ubiquitous and independent of resolution in the satellites that are past their first infall, there seems to be a wide range of amplification mechanisms acting together. Our result highlights the importance of considering the environment of dwarf galaxies when interpreting their magnetic field properties as well as related observables such as their gamma-ray and radio emission, the latter being particularly relevant for future observations such as with the SKA observatory.
△ Less
Submitted 25 September, 2024;
originally announced September 2024.
-
The impact of baryons on the internal structure of dark matter haloes from dwarf galaxies to superclusters in the redshift range 0<z<7
Authors:
Daniele Sorini,
Sownak Bose,
Rüdiger Pakmor,
Lars Hernquist,
Volker Springel,
Boryana Hadzhiyska,
César Hernández-Aguayo,
Rahul Kannan
Abstract:
We investigate the redshift evolution of the concentration-mass relationship of dark matter haloes in state-of-the-art cosmological hydrodynamic simulations and their dark-matter-only counterparts. By combining the IllustrisTNG suite and the novel MillenniumTNG simulation, our analysis encompasses a wide range of box size ($50 - 740 \: \rm cMpc$) and mass resolution (…
▽ More
We investigate the redshift evolution of the concentration-mass relationship of dark matter haloes in state-of-the-art cosmological hydrodynamic simulations and their dark-matter-only counterparts. By combining the IllustrisTNG suite and the novel MillenniumTNG simulation, our analysis encompasses a wide range of box size ($50 - 740 \: \rm cMpc$) and mass resolution ($8.5 \times 10^4 - 3.1 \times 10^7 \: \rm M_{\odot}$ per baryonic mass element). This enables us to study the impact of baryons on the concentration-mass relationship in the redshift interval $0<z<7$ over an unprecedented halo mass range, extending from dwarf galaxies to superclusters ($\sim 10^{9.5}-10^{15.5} \, \rm M_{\odot}$). We find that the presence of baryons increases the steepness of the concentration-mass relationship at higher redshift, and demonstrate that this is driven by adiabatic contraction of the profile, due to gas accretion at early times, which promotes star formation in the inner regions of haloes. At lower redshift, when the effects of feedback start to become important, baryons decrease the concentration of haloes below the mass scale $\sim 10^{11.5} \, \rm M_{\odot}$. Through a rigorous information criterion test, we show that broken power-law models accurately represent the redshift evolution of the concentration-mass relationship, and of the relative difference in the total mass of haloes induced by the presence of baryons. We provide the best-fit parameters of our empirical formulae, enabling their application to models that mimic baryonic effects in dark-matter-only simulations over six decades in halo mass in the redshift range $0<z<7$.
△ Less
Submitted 3 September, 2024;
originally announced September 2024.
-
Probing galaxy evolution from $z=0$ to $z\simeq10$ through galaxy scaling relations in three L-Galaxies flavours
Authors:
Akash Vani,
Mohammadreza Ayromlou,
Guinevere Kauffmann,
Volker Springel
Abstract:
We present a comprehensive examination of the three most recent versions of the L-Galaxies semi-analytic galaxy formation model, focusing on the evolution of galaxy properties across a broad stellar mass range ($10^7\:{\rm M}_{\odot}\lesssim{M_\star}\lesssim10^{12}\:{\rm M}_{\odot}$) from $z=0$ to $z\simeq10$. We compare the predictions with the latest multiband data from key astronomical surveys,…
▽ More
We present a comprehensive examination of the three most recent versions of the L-Galaxies semi-analytic galaxy formation model, focusing on the evolution of galaxy properties across a broad stellar mass range ($10^7\:{\rm M}_{\odot}\lesssim{M_\star}\lesssim10^{12}\:{\rm M}_{\odot}$) from $z=0$ to $z\simeq10$. We compare the predictions with the latest multiband data from key astronomical surveys, including SDSS, CANDELS, and COSMOS along with HST, JWST, and ALMA. We assess the models' ability to reproduce various time-dependent galaxy scaling relations for star-forming and quenched galaxies. Key focus areas include global galaxy properties such as stellar mass functions, cosmic star formation rate density, and the evolution of the main sequence of star-forming galaxies. Additionally, we examine resolved morphological properties such as the galaxy mass-size relation, alongside core $(R<1\,{\rm{kpc}})$ and effective $(R<R_{\rm{e}})$ stellar-mass surface densities as a function of stellar mass. This analysis reveals that the \textsc{L-Galaxies} models are in qualitatively good agreement with observed global galaxy scaling relations up to $z\simeq10$. However, significant discrepancies exist at both low and high redshifts in accurately reproducing the number density, size, and surface density evolution of quenched galaxies. These issues are most pronounced for massive galaxies, where the simulations underpredict the abundance of quenched galaxies at $z\geq1.5$, reaching a discrepancy of a factor of 60 by $z\approx3$, with sizes several times larger than observed. Therefore, we propose that the physical prescriptions governing the cessation of star formation in galaxies, such as AGN feedback and processes related to merging, require improvement to be better supported by observational data.
△ Less
Submitted 1 August, 2024;
originally announced August 2024.
-
The MillenniumTNG Project: Impact of massive neutrinos on the cosmic large-scale structure and the distribution of galaxies
Authors:
César Hernández-Aguayo,
Volker Springel,
Sownak Bose,
Carlos Frenk,
Adrian Jenkins,
Monica Barrera,
Fulvio Ferlito,
Rüdiger Pakmor,
Simon D. M. White,
Lars Hernquist,
Ana Maria Delgado,
Rahul Kannan,
Boryana Hadzhiyska
Abstract:
We discuss the cold dark matter plus massive neutrinos simulations of the MillenniumTNG (MTNG) project, which aim to improve understanding of how well ongoing and future large-scale galaxy surveys will measure neutrino masses. Our largest simulations, $3000\,{\rm Mpc}$ on a side, use $10240^3$ particles of mass $m_{p} = 6.66\times 10^{8}\,h^{-1}{\rm M}_\odot$ to represent cold dark matter, and…
▽ More
We discuss the cold dark matter plus massive neutrinos simulations of the MillenniumTNG (MTNG) project, which aim to improve understanding of how well ongoing and future large-scale galaxy surveys will measure neutrino masses. Our largest simulations, $3000\,{\rm Mpc}$ on a side, use $10240^3$ particles of mass $m_{p} = 6.66\times 10^{8}\,h^{-1}{\rm M}_\odot$ to represent cold dark matter, and $2560^3$ to represent a population of neutrinos with summed mass $M_ν= 100\,{\rm meV}$. Smaller volume runs with $\sim 630\,{\rm Mpc}$ also include cases with $M_ν= 0\,\textrm{and}\, 300\,{\rm meV}$. All simulations are carried out twice using the paired-and-fixed technique for cosmic variance reduction. We evolve the neutrino component using the particle-based $δf$ importance sampling method, which greatly reduces shot noise in the neutrino density field. In addition, we modify the GADGET-4 code to account both for the influence of relativistic and mildly relativistic components on the expansion rate and for non-Newtonian effects on the largest represented simulation scales. This allows us to quantify accurately the impact of neutrinos on basic statistical measures of nonlinear structure formation, such as the matter power spectrum and the halo mass function. We use semi-analytic models of galaxy formation to predict the galaxy population and its clustering properties as a function of summed neutrino mass, finding significant ($\sim 10\%$) impacts on the cosmic star formation rate history, the galaxy mass function, and the clustering strength. This offers the prospect of identifying combinations of summary statistics that are optimally sensitive to the neutrino mass.
△ Less
Submitted 30 July, 2024;
originally announced July 2024.
-
Percolation Statistics in the MillenniumTNG Simulations
Authors:
Eniko Regos,
Volker Springel,
Sownak Bose,
Boryana Hadzhiyska,
Cesar Hernandez-Aguayo
Abstract:
The statistical analysis of cosmic large-scale structure is most often based on simple two-point summary statistics, like the power spectrum or the two-point correlation function of a sample of galaxies or other types of tracers. In contrast, topological measures of clustering are also sensitive to higher-order correlations, and thus offer the prospect to access additional information that may har…
▽ More
The statistical analysis of cosmic large-scale structure is most often based on simple two-point summary statistics, like the power spectrum or the two-point correlation function of a sample of galaxies or other types of tracers. In contrast, topological measures of clustering are also sensitive to higher-order correlations, and thus offer the prospect to access additional information that may harbor important constraining power. We here revisit one such geometric measure of the cosmic web in the form of the so-called percolation analysis, using the recent MillenniumTNG simulation suite of the LCDM paradigm. We analyze continuum percolation statistics both for high resolution dark matter particle distributions, as well as for galaxy mock catalogues from a semi-analytic galaxy formation model within a periodic simulation volume of 3000 Mpc on a side. For comparison, we also investigate the percolation statistics of random particle sets and neutrino distributions with two different summed particle masses. We find that the percolation statistics of the dark matter distribution evolves strongly with redshift and thus clustering strength, yielding progressively lower percolation threshold towards later times. However, there is a sizable residual dependence on numerical resolution which we interpret as a residual influence of different levels of shot noise. This is corroborated by our analysis of galaxy mock catalogues whose results depend on sampling density more strongly than on galaxy selection criteria. While this limits the discriminative power of percolation statistics, our results suggest that it still remains useful as a complementary cosmological test when controlled for sampling density.
△ Less
Submitted 21 August, 2024; v1 submitted 2 July, 2024;
originally announced July 2024.
-
Improving the Accuracy of Halo Mass Based Statistics For Fast Approximate N-body Simulations
Authors:
Yiheng Wu,
Hong Guo,
Volker Springel
Abstract:
Approximate N-body methods, such as FastPM and COLA, have been successful in modelling halo and galaxy clustering statistics, but their low resolution on small scales is a limitation for applications that require high precision. Full N-body simulations can provide better accuracy but are too computationally expensive for a quick exploration of cosmological parameters. This paper presents a method…
▽ More
Approximate N-body methods, such as FastPM and COLA, have been successful in modelling halo and galaxy clustering statistics, but their low resolution on small scales is a limitation for applications that require high precision. Full N-body simulations can provide better accuracy but are too computationally expensive for a quick exploration of cosmological parameters. This paper presents a method for correcting distinct haloes identified in fast N-body simulations, so that various halo statistics improve to a percent level accuracy. The scheme seeks to find empirical corrections to halo properties such that the virial mass is the same as that of a corresponding halo in a full N-body simulation. The modified outer density contour of the corrected halo is determined on the basis of the FastPM settings and the number of particles inside the halo. This method only changes some parameters of the halo finder, and does not require any extra CPU-cost. We demonstrate that the adjusted halo catalogues of FastPM simulations significantly improve the precision of halo mass-based statistics from redshifts $z=0.0$ to $1.0$, and that our calibration can be applied to different cosmologies without needing to be recalibrated.
△ Less
Submitted 14 June, 2024;
originally announced June 2024.
-
Bursty Star Formation in Dwarfs is Sensitive to Numerical Choices in Supernova Feedback Models
Authors:
Eric Zhang,
Laura V Sales,
Federico Marinacci,
Paul Torrey,
Mark Vogelsberger,
Volker Springel,
Hui Li,
Rüdiger Pakmor,
Thales A Gutcke
Abstract:
Simulations of galaxy formation are mostly unable to resolve the energy-conserving phase of individual supernova events, having to resort to subgrid models to distribute the energy and momentum resulting from stellar feedback. However, the properties of these simulated galaxies, including the morphology, stellar mass formed and the burstiness of the star formation history, are highly sensitive to…
▽ More
Simulations of galaxy formation are mostly unable to resolve the energy-conserving phase of individual supernova events, having to resort to subgrid models to distribute the energy and momentum resulting from stellar feedback. However, the properties of these simulated galaxies, including the morphology, stellar mass formed and the burstiness of the star formation history, are highly sensitive to numerical choices adopted in these subgrid models. Using the {\small SMUGGLE} stellar feedback model, we compute idealized simulations of a $M_{\rm vir} \sim 10^{10} \, \msun$ dwarf galaxy, a regime where most simulation codes predict significant burstiness in star formation, resulting in strong gas flows that lead to the formation of dark matter cores. We find that by varying only the directional distribution of momentum imparted from supernovae to the surrounding gas, while holding the total momentum per supernova constant, bursty star formation may be amplified or completely suppressed, and the total stellar mass formed can vary by as much as a factor of $\sim 3$. In particular, when momentum is primarily directed perpendicular to the gas disk, less bursty and lower overall star formation rates result, yielding less gas turbulence, more disky morphologies and a retention of cuspy dark matter density profiles. An improved understanding of the non-linear coupling of stellar feedback into inhomogeneous gaseous media is thus needed to make robust predictions for stellar morphologies and dark matter core formation in dwarfs independent of uncertain numerical choices in the baryonic treatment.
△ Less
Submitted 14 June, 2024;
originally announced June 2024.
-
Bar formation and evolution in the cosmological context: Inputs from the Auriga simulations
Authors:
Francesca Fragkoudi,
Robert Grand,
Rüdiger Pakmor,
Facundo Gómez,
Federico Marinacci,
Volker Springel
Abstract:
Galactic bars drive the internal evolution of spiral galaxies, while their formation is tightly coupled to the properties of their host galaxy and dark matter halo. To explore what drives bar formation in the cosmological context and how these structures evolve throughout cosmic history, we use the Auriga suite of magneto-hydrodynamical cosmological zoom-in simulations. We find that bars are robus…
▽ More
Galactic bars drive the internal evolution of spiral galaxies, while their formation is tightly coupled to the properties of their host galaxy and dark matter halo. To explore what drives bar formation in the cosmological context and how these structures evolve throughout cosmic history, we use the Auriga suite of magneto-hydrodynamical cosmological zoom-in simulations. We find that bars are robust and long-lived structures, and we recover a decreasing bar fraction with increasing redshift which plateaus around $\sim20\%$ at $z\sim3$. We find that bars which form at low and intermediate redshifts grow longer with time, while bars that form at high redshifts are born `saturated' in length, likely due to their merger-induced formation pathway. This leads to a larger bar-to-disc size ratio at high redshifts as compared to the local Universe. We subsequently examine the multi-dimensional parameter space thought to drive bar formation. We find that barred galaxies tend to have lower Toomre $Q$ values at the time of their formation, while we do not find a difference in the gas fraction of barred and unbarred populations when controlling for stellar mass. Barred galaxies tend to be more baryon-dominated at all redshifts, assembling their stellar mass earlier, while galaxies that are baryon-dominated but that do not host a bar, have a higher ex-situ bulge fraction. We explore the implications of the baryon-dominance of barred galaxies on the Tully-Fisher relation, finding an offset from the unbarred relation; confirming this in observations would serve as additional evidence for dark matter, as this behaviour is not readily explained in modified gravity scenarios.
△ Less
Submitted 19 June, 2024; v1 submitted 12 June, 2024;
originally announced June 2024.
-
Ray-tracing vs. Born approximation in full-sky weak lensing simulations of the MillenniumTNG project
Authors:
Fulvio Ferlito,
Christopher T. Davies,
Volker Springel,
Martin Reinecke,
Alessandro Greco,
Ana Maria Delgado,
Simon D. M. White,
César Hernández-Aguayo,
Sownak Bose,
Lars Hernquist
Abstract:
Weak gravitational lensing is a powerful tool for precision tests of cosmology. As the expected deflection angles are small, predictions based on non-linear N-body simulations are commonly computed with the Born approximation. Here we examine this assumption using ${\small DORIAN}$, a newly developed full-sky ray-tracing scheme applied to high-resolution mass-shell outputs of the two largest simul…
▽ More
Weak gravitational lensing is a powerful tool for precision tests of cosmology. As the expected deflection angles are small, predictions based on non-linear N-body simulations are commonly computed with the Born approximation. Here we examine this assumption using ${\small DORIAN}$, a newly developed full-sky ray-tracing scheme applied to high-resolution mass-shell outputs of the two largest simulations in the MillenniumTNG suite, each with a 3000 Mpc box containing almost 1.1 trillion cold dark matter particles in addition to 16.7 billion particles representing massive neutrinos. We examine simple two-point statistics like the angular power spectrum of the convergence field, as well as statistics sensitive to higher order correlations such as peak and minimum statistics, void statistics, and Minkowski functionals of the convergence maps. Overall, we find only small differences between the Born approximation and a full ray-tracing treatment. While these are negligibly small at power-spectrum level, some higher order statistics show more sizable effects; ray-tracing is necessary to achieve percent level precision. At the resolution reached here, full-sky maps with 0.8 billion pixels and an angular resolution of 0.43 arcmin, we find that interpolation accuracy can introduce appreciable errors in ray-tracing results. We therefore implemented an interpolation method based on nonuniform fast Fourier transforms (NUFFT) along with more traditional methods. Bilinear interpolation introduces significant smoothing, while nearest grid point sampling agrees well with NUFFT, at least for our fiducial source redshift, $z_s=1.0$, and for the 1 arcmin smoothing we use for higher-order statistics.
△ Less
Submitted 12 June, 2024;
originally announced June 2024.
-
A calibrated model for N-body dynamical friction acting on supermassive black holes
Authors:
Anna Genina,
Volker Springel,
Antti Rantala
Abstract:
Black holes are believed to be crucial in regulating star formation in massive galaxies, which makes it essential to faithfully represent the physics of these objects in cosmological hydrodynamics simulations. Limited spatial and mass resolution and the associated discreteness noise make following the dynamics of black holes especially challenging. In particular, dynamical friction, which is respo…
▽ More
Black holes are believed to be crucial in regulating star formation in massive galaxies, which makes it essential to faithfully represent the physics of these objects in cosmological hydrodynamics simulations. Limited spatial and mass resolution and the associated discreteness noise make following the dynamics of black holes especially challenging. In particular, dynamical friction, which is responsible for driving massive black holes towards the centres of galaxies, cannot be accurately modelled with softened $N$-body interactions. A number of subgrid models have been proposed to mimic dynamical friction or directly include its full effects in simulations. Each of these methods has its individual benefits and shortcomings, while all suffer from a common issue of being unable to represent black holes with masses below a few times the simulated dark matter particle mass. In this paper, we propose a correction for unresolved dynamical friction, which has been calibrated on simulations run with the code KETJU, in which gravitational interactions of black holes are not softened. We demonstrate that our correction is able to sink black holes with masses greater than the dark matter particle mass at the correct rate. We show that the impact of stochasticity is significant for low-mass black holes ($M_{\rm BH} \leq 5 M_{\rm DM}$) and propose a correction for stochastic heating. Combined, this approach is applicable to next generation cosmological hydrodynamics simulations that jointly track galaxy and black hole growth with realistic black hole orbits.
△ Less
Submitted 12 September, 2024; v1 submitted 14 May, 2024;
originally announced May 2024.
-
The effect of local photoionization on the galaxy properties and the circumgalactic medium in simulations of Milky Way-sized galaxies
Authors:
Bocheng Zhu,
Volker Springel
Abstract:
In this study, we investigate the impact of local stellar radiation in cosmological zoom simulations of the formation of Milky Way-sized galaxies. We include the radiation field as an additional feedback component that is computed alongside gravity with a tree code in an optically thin approximation. We resimulate the initial conditions of five Milk Way-like systems taken from the Auriga project w…
▽ More
In this study, we investigate the impact of local stellar radiation in cosmological zoom simulations of the formation of Milky Way-sized galaxies. We include the radiation field as an additional feedback component that is computed alongside gravity with a tree code in an optically thin approximation. We resimulate the initial conditions of five Milk Way-like systems taken from the Auriga project with and without stellar radiation, and study the effects of local stellar radiation on several properties of the galaxies and the circumgalactic medium (CGM). Similar to previous findings, we observe with our current model that local stellar radiation can modify gas cooling in the CGM and thus suppress star formation and the surface densities of young stars and HI gas, while having little impact on the total gas content. In particular, it also suppresses the peak of the rotation curve and reduces the mass of the stellar bulge. In the CGM region, the young stellar radiation exceeds the external UVB and dominates the radiation field within the virial halo at all redshifts. Nevertheless, we find that the local stellar radiation, as implemented in the current study, has overall little impact on the radial density and temperature profile of the CGM gas. However, for the ion species HI and MgII the column densities within $\sim 0.3\,R_{\rm vir}$ are reduced, while the OVI column density is hardly impacted by the radiation field due to a lack of soft X-ray components in our current model. Additional effects can be expected from the radiation of the central AGN during phases of quasar activity and from soft X-ray sources, which have not yet been included in the simulations of the present study.
△ Less
Submitted 3 September, 2024; v1 submitted 21 April, 2024;
originally announced April 2024.
-
Cosmological gas accretion history onto the stellar discs of Milky Way-like galaxies in the Auriga simulations -- (II) The inside-out growth of discs
Authors:
Federico G. Iza,
Sebastián E. Nuza,
Cecilia Scannapieco,
Robert J. J. Grand,
Facundo A. Gómez,
Volker Springel,
Rüdiger Pakmor,
Federico Marinacci,
Francesca Fragkoudi
Abstract:
We investigate the growth of stellar discs in Milky Way-mass galaxies using the magnetohydrodynamical simulations of the Auriga Project in a full cosmological context. We focus on the gas accretion process along the discs, calculating the net, infall and outflow rates as a function of galactocentric distance, and investigate the relation between them and the star formation activity. The stellar di…
▽ More
We investigate the growth of stellar discs in Milky Way-mass galaxies using the magnetohydrodynamical simulations of the Auriga Project in a full cosmological context. We focus on the gas accretion process along the discs, calculating the net, infall and outflow rates as a function of galactocentric distance, and investigate the relation between them and the star formation activity. The stellar distributions of around 70% of the simulated galaxies exhibit an ``inside-out'' pattern, with older (younger) stellar populations preferentially located in the inner (outer) disc regions. In all cases, we find a very tight correlation between the infall, outflow and net accretion rates, as well as between these three quantities and the star formation rate. This is because the amount of gas which is ultimately available for star formation in each radial ring depends not only on the infall rates, but also on the amount of gas leaving the disc in outflows, which directly relates to the local star formation level. Therefore, any of these rates can be used to identify galaxies with inside-out growth. For these galaxies, the correlation between the dominant times of accretion/star formation and disc radius is well fitted by a linear function. We also find that, when averaged over galaxies with formation histories similar to the Milky Way, the simulated accretion rates show a similar evolution (both temporally- and radially-integrated) to the usual accretion prescriptions used in chemical evolution models, although some major differences arise at early times and in the inner disc regions.
△ Less
Submitted 10 April, 2024;
originally announced April 2024.
-
The influence of baryons on low-mass haloes
Authors:
Haonan Zheng,
Sownak Bose,
Carlos S. Frenk,
Liang Gao,
Adrian Jenkins,
Shihong Liao,
Volker Springel,
Jie Wang,
Simon D. M. White
Abstract:
The Voids-within-Voids-within-Voids (VVV) project used dark-matter-only simulations to study the abundance and structure of dark matter haloes over the full mass range populated in the standard $Λ\mathrm{CDM}$ cosmology. Here we explore how baryonic effects modify these results for $z=0$ halo masses in the range $10^4$ to $10^7~\mathrm{M_\odot}$, below the threshold for galaxy formation. Our main…
▽ More
The Voids-within-Voids-within-Voids (VVV) project used dark-matter-only simulations to study the abundance and structure of dark matter haloes over the full mass range populated in the standard $Λ\mathrm{CDM}$ cosmology. Here we explore how baryonic effects modify these results for $z=0$ halo masses in the range $10^4$ to $10^7~\mathrm{M_\odot}$, below the threshold for galaxy formation. Our main study focuses on three simulations from identical initial conditions at $z=127$, one following dark matter only, one including non-radiative gas, and one additionally including the baryonic physics relevant in this halo mass range (cooling and photoheating). In the non-radiative simulation, above $10^{5.5}~\mathrm{M_\odot}$, halo abundance and internal structure are very similar to the dark-matter-only simulation, and the baryon to dark matter ratio is everywhere close to the cosmic value. At lower mass, this ratio drops and haloes are less concentrated and less massive in the non-radiative case. Test simulations at higher resolution show this to be mainly a resolution effect; the expected drop in baryon content due to residual pressure effects only becomes substantial for $z=0$ haloes below $\sim 10^{2.7}~\mathrm{M_\odot}$. However, gas is heated by reionization at $z=6$ in our "full physics" run, and this results in almost complete expulsion of gas from all haloes in our simulated mass range. This suppresses the halo mass function by $\sim 30 \%$, lowers halo concentration, and consequently weakens the dark matter annihilation signal by $\sim 40-60 \%$.
△ Less
Submitted 9 August, 2024; v1 submitted 25 March, 2024;
originally announced March 2024.
-
A deep learning model for the density profiles of subhaloes in IllustrisTNG
Authors:
Luisa Lucie-Smith,
Giulia Despali,
Volker Springel
Abstract:
We present a machine-learning-based model for the total density profiles of subhaloes with masses $M \gtrsim 7\times 10^8\,h^{-1}{\rm M}_\odot$ in the IllustrisTNG100 simulation. The model is based on an interpretable variational encoder (IVE) which returns the independent factors of variation in the density profiles within a low-dimensional representation, as well as the predictions for the densi…
▽ More
We present a machine-learning-based model for the total density profiles of subhaloes with masses $M \gtrsim 7\times 10^8\,h^{-1}{\rm M}_\odot$ in the IllustrisTNG100 simulation. The model is based on an interpretable variational encoder (IVE) which returns the independent factors of variation in the density profiles within a low-dimensional representation, as well as the predictions for the density profiles themselves. The IVE returns accurate and unbiased predictions on all radial ranges, including the outer region profile where the subhaloes experience tidal stripping; here its fit accuracy exceeds that of the commonly used Einasto profile. The IVE discovers three independent degrees of freedom in the profiles, which can be interpreted in terms of the formation history of the subhaloes. In addition to the two parameters controlling the normalization and inner shape of the profile, the IVE discovers a third parameter that accounts for the impact of tidal stripping onto the subhalo outer profile; this parameter is sensitive to the mass loss experienced by the subhalo after its infall onto its parent halo. Baryonic physics in the IllustrisTNG galaxy formation model does not impact the number of degrees of freedom identified in the profile compared to the pure dark matter expectations, nor their physical interpretation. Our newly proposed profile fit can be used in strong lensing analyses or other observational studies which aim to constrain cosmology from small-scale structures.
△ Less
Submitted 18 June, 2024; v1 submitted 18 March, 2024;
originally announced March 2024.
-
The origin of lopsided satellite galaxy distribution around isolated systems in MillenniumTNG
Authors:
Yikai Liu,
Peng Wang,
Hong Guo,
Volker Springel,
Sownak Bose,
Rüdiger Pakmor,
Lars Hernquist
Abstract:
Dwarf satellites in galaxy groups are distributed in an anisotropic and asymmetric manner, which is called the ``lopsided satellite distribution''. This lopsided signal has been observed not only in galaxy pairs but also in isolated systems. However, the physical origin of the lopsided signal in isolated systems is still unknown. In this work, we investigate this in the state-of-the-art hydrodynam…
▽ More
Dwarf satellites in galaxy groups are distributed in an anisotropic and asymmetric manner, which is called the ``lopsided satellite distribution''. This lopsided signal has been observed not only in galaxy pairs but also in isolated systems. However, the physical origin of the lopsided signal in isolated systems is still unknown. In this work, we investigate this in the state-of-the-art hydrodynamical simulation of the MillenniumTNG Project by tracing each system back to high redshift. We find that the lopsided signal is dominated by satellites located in the outer regions of the halo and is also dominated by recently accreted satellites. The lopsided signal originates from the anisotropic accretion of galaxies from the surrounding large-scale structure and that, after accretion, the nonlinear evolution of satellites inside the dark-matter halo weakens the lopsidedness. The signal decreases as cosmic time passes because of a competition between anisotropic accretion and internal evolution within dark matter halos. Our findings provide a useful perspective for the study of galaxy evolution, especially for the origin of the spatial satellite galaxy distributions.
△ Less
Submitted 2 March, 2024;
originally announced March 2024.
-
The boundary of cosmic filaments
Authors:
Wei Wang,
Peng Wang,
Hong Guo,
Xi Kang,
Noam I. Libeskind,
Daniela Galarraga-Espinosa,
Volker Springel,
Rahul Kannan,
Lars Hernquist,
Rudiger Pakmor,
Haoran Yu,
Sownak Bose,
Quan Guo,
Luo Yu,
Cesar Hernandez-Aguayo
Abstract:
For decades, the boundary of cosmic filaments have been a subject of debate. In this work, we determine the physically-motivated radii of filaments by constructing stacked galaxy number density profiles around the filament spines. We find that the slope of the profile changes with distance to the filament spine, reaching its minimum at approximately 1 Mpc at z = 0 in both state-of-the-art hydrodyn…
▽ More
For decades, the boundary of cosmic filaments have been a subject of debate. In this work, we determine the physically-motivated radii of filaments by constructing stacked galaxy number density profiles around the filament spines. We find that the slope of the profile changes with distance to the filament spine, reaching its minimum at approximately 1 Mpc at z = 0 in both state-of-the-art hydrodynamical simulations and observational data. This can be taken as the average value of the filament radius. Furthermore, we note that the average filament radius rapidly decreases from z = 4 to z = 1, and then slightly increases. Moreover, we find that the filament radius depends on the filament length, the distance from connected clusters, and the masses of the clusters. These results suggest a two-phase formation scenario of cosmic filaments. The filaments experience rapid contraction before z = 1, but their density distribution has remained roughly stable since then. The subsequent mass transport along the filaments to the connected clusters is likely to have contributed to the formation of the clusters themselves.
△ Less
Submitted 3 August, 2024; v1 submitted 18 February, 2024;
originally announced February 2024.
-
Brightest Cluster Galaxy Offsets in Cold Dark Matter
Authors:
Cian Roche,
Michael McDonald,
Josh Borrow,
Mark Vogelsberger,
Xuejian Shen,
Volker Springel,
Lars Hernquist,
Ruediger Pakmor,
Sownak Bose,
Rahul Kannan
Abstract:
The distribution of offsets between the brightest cluster galaxies of galaxy clusters and the centroid of their dark matter distributions is a promising probe of the underlying dark matter physics. In particular, since this distribution is sensitive to the shape of the potential in galaxy cluster cores, it constitutes a test of dark matter self-interaction on the largest mass scales in the univers…
▽ More
The distribution of offsets between the brightest cluster galaxies of galaxy clusters and the centroid of their dark matter distributions is a promising probe of the underlying dark matter physics. In particular, since this distribution is sensitive to the shape of the potential in galaxy cluster cores, it constitutes a test of dark matter self-interaction on the largest mass scales in the universe. We examine these offsets in three suites of modern cosmological simulations; IllustrisTNG, MillenniumTNG and BAHAMAS. For clusters above $10^{14}\rm{M_\odot}$, we examine the dependence of the offset distribution on gravitational softening length, the method used to identify centroids, redshift, mass, baryonic physics, and establish the stability of our results with respect to various nuisance parameter choices. We find that offsets are overwhelmingly measured to be smaller than the minimum converged length scale in each simulation, with a median offset of $\sim1\rm{kpc}$ in the highest resolution simulation considered, TNG300-1, which uses a gravitational softening length of $1.48\rm{kpc}$. We also find that centroids identified via source extraction on smoothed dark matter and stellar particle data are consistent with the potential minimum, but that observationally relevant methods sensitive to cluster strong gravitational lensing scales, or those using gas as a tracer for the potential can overestimate offsets by factors of $\sim10$ and $\sim30$, respectively. This has the potential to reduce tensions with existing offset measurements which have served as evidence for a nonzero dark matter self-interaction cross section.
△ Less
Submitted 5 August, 2024; v1 submitted 1 February, 2024;
originally announced February 2024.
-
The SRG/eROSITA All-Sky Survey: Constraints on AGN Feedback in Galaxy Groups
Authors:
Y. E. Bahar,
E. Bulbul,
V. Ghirardini,
J. S. Sanders,
X. Zhang,
A. Liu,
N. Clerc,
E. Artis,
F. Balzer,
V. Biffi,
S. Bose,
J. Comparat,
K. Dolag,
C. Garrel,
B. Hadzhiyska,
C. Hernández-Aguayo,
L. Hernquist,
M. Kluge,
S. Krippendorf,
A. Merloni,
K. Nandra,
R. Pakmor,
P. Popesso,
M. Ramos-Ceja,
R. Seppi
, et al. (3 additional authors not shown)
Abstract:
We investigate the impact of AGN feedback, on the entropy and characteristic temperature measurements of galaxy groups detected in the SRG/eROSITA's first All-Sky Survey (eRASS1) to shed light on the characteristics of the feedback mechanisms. We analyze deeper eROSITA observations of 1178 galaxy groups detected in eRASS1. We divide the sample into 271 subsamples and extract average thermodynamic…
▽ More
We investigate the impact of AGN feedback, on the entropy and characteristic temperature measurements of galaxy groups detected in the SRG/eROSITA's first All-Sky Survey (eRASS1) to shed light on the characteristics of the feedback mechanisms. We analyze deeper eROSITA observations of 1178 galaxy groups detected in eRASS1. We divide the sample into 271 subsamples and extract average thermodynamic properties, including electron density, temperature, and entropy at three characteristic radii along with the integrated temperature by jointly analyzing X-ray images and spectra following a Bayesian approach. We present the tightest constraints on the impact of AGN feedback through our average entropy and characteristic temperature measurements of the largest group sample used in X-ray studies, incorporating major systematics in our analysis. We find that entropy shows an increasing trend with temperature in the form of a power-law-like relation at the higher intra-group medium temperatures, while for the low mass groups, a slight flattening is observed on the average entropy. Overall, the observed entropy measurements agree well with the earlier measurements in the literature. The comparisons with the state-of-the-art cosmological hydrodynamic simulations (MillenniumTNG, Magneticum, OWL simulations) after the applications of the selection function calibrated for our galaxy groups reveal that observed entropy profiles in the cores are below the predictions of simulations. At the mid-region, the entropy measurements agree well with the Magneticum simulations, whereas the predictions of MillenniumTNG and OWL simulations fall below observations. At the outskirts, the overall agreement between the observations and simulations improves, with Magneticum simulations reproducing the observations the best. Our measurements will pave the way for more realistic AGN feedback implementations in simulations.
△ Less
Submitted 30 January, 2024;
originally announced January 2024.
-
Overview and public data release of the augmented Auriga Project: cosmological simulations of dwarf and Milky Way-mass galaxies
Authors:
Robert J. J. Grand,
Francesca Fragkoudi,
Facundo A. Gómez,
Adrian Jenkins,
Federico Marinacci,
Rüdiger Pakmor,
Volker Springel
Abstract:
We present an extended suite of the Auriga cosmological gravo-magnetohydrodynamical "zoom-in" simulations of 40 Milky Way-mass halos and 26 dwarf galaxy-mass halos run with the moving-mesh code Arepo. Auriga adopts the $Λ$ Cold Dark Matter ($Λ$CDM) cosmogony and includes a comprehensive galaxy formation physics model following the coupled cosmic evolution of dark matter, gas, stars, and supermassi…
▽ More
We present an extended suite of the Auriga cosmological gravo-magnetohydrodynamical "zoom-in" simulations of 40 Milky Way-mass halos and 26 dwarf galaxy-mass halos run with the moving-mesh code Arepo. Auriga adopts the $Λ$ Cold Dark Matter ($Λ$CDM) cosmogony and includes a comprehensive galaxy formation physics model following the coupled cosmic evolution of dark matter, gas, stars, and supermassive black holes which has been shown to produce numerically well-converged galaxy properties for Milky Way-mass systems. We describe the first public data release of this augmented suite of Auriga simulations, which includes raw snapshots, group catalogues, merger trees, initial conditions, and supplementary data, as well as public analysis tools with worked examples of how to use the data. To demonstrate the value and robustness of the simulation predictions, we analyse a series of low-redshift global properties that compare well with many observed scaling relations, such as the Tully-Fisher relation, the star-forming main sequence, and HI gas fraction/disc thickness. Finally, we show that star-forming gas discs appear to build rotation and velocity dispersion rapidly for $z\gtrsim 3$ before they "settle" into ever-increasing rotation-dispersion ratios ($V/σ$). This evolution appears to be in rough agreement with some kinematic measurements from H$α$ observations, and demonstrates an application of how to utilise the released data.
△ Less
Submitted 3 July, 2024; v1 submitted 16 January, 2024;
originally announced January 2024.
-
Supersonic turbulence simulations with GPU-based high-order Discontinuous Galerkin hydrodynamics
Authors:
Miha Cernetic,
Volker Springel,
Thomas Guillet,
Rüdiger Pakmor
Abstract:
We investigate the numerical performance of a Discontinuous Galerkin (DG) hydrodynamics implementation when applied to the problem of driven, isothermal supersonic turbulence. While the high-order element-based spectral approach of DG is known to efficiently produce accurate results for smooth problems (exponential convergence with expansion order), physical discontinuities in solutions, like shoc…
▽ More
We investigate the numerical performance of a Discontinuous Galerkin (DG) hydrodynamics implementation when applied to the problem of driven, isothermal supersonic turbulence. While the high-order element-based spectral approach of DG is known to efficiently produce accurate results for smooth problems (exponential convergence with expansion order), physical discontinuities in solutions, like shocks, prove challenging and may significantly diminish DG's applicability to practical astrophysical applications. We consider whether DG is able to retain its accuracy and stability for highly supersonic turbulence, characterized by a network of shocks. We find that our new implementation, which regularizes shocks at sub-cell resolution with artificial viscosity, still performs well compared to standard second-order schemes for moderately high Mach number turbulence, provided we also employ an additional projection of the primitive variables onto the polynomial basis to regularize the extrapolated values at cell interfaces. However, the accuracy advantage of DG diminishes significantly in the highly supersonic regime. Nevertheless, in turbulence simulations with a wide dynamic range that start with supersonic Mach numbers and can resolve the sonic point, the low numerical dissipation of DG schemes still proves advantageous in the subsonic regime. Our results thus support the practical applicability of DG schemes for demanding astrophysical problems that involve strong shocks and turbulence, such as star formation in the interstellar medium. We also discuss the substantial computational cost of DG when going to high order, which needs to be weighted against the resulting accuracy gain. For problems containing shocks, this favours the use of comparatively low DG order.
△ Less
Submitted 1 October, 2024; v1 submitted 12 January, 2024;
originally announced January 2024.
-
The evolutionary path of void galaxies in TNG300 simulation
Authors:
Agustín Rodríguez-Medrano,
Volker Springel,
Federico Stasyszyn,
Dante Paz
Abstract:
The properties of galaxies in low-density regions of the universe suggest an interplay between galaxy formation and environment. However, the specific reason why this particular large-scale environment influences the evolution of galaxies remains unclear. This paper examines the properties and evolutionary paths of galaxies within cosmic voids using the Illustris TNG300 simulation. The population…
▽ More
The properties of galaxies in low-density regions of the universe suggest an interplay between galaxy formation and environment. However, the specific reason why this particular large-scale environment influences the evolution of galaxies remains unclear. This paper examines the properties and evolutionary paths of galaxies within cosmic voids using the Illustris TNG300 simulation. The population of void galaxies at z = 0 has a higher star formation rate, a smaller stellar-to-halo-mass ratio, higher gas metallicity, and lower stellar metallicity in comparison with non-void galaxies at fixed stellar mass. Our analysis shows that these differences are mainly due to the characteristics of galaxies classified as satellites, for which the largest differences between void and non-void samples are found. Although the mean number of mergers is similar between void and non-void samples at a fixed stellar mass, void galaxies tend to experience mergers at later times, resulting in a more recent accumulation of accreted stellar mass. While the mean net accreted mass is comparable for high mass galaxies, low mass void galaxies tend to exhibit higher fractions of accreted stars than non-void galaxies. This finding challenges the common notion that void galaxies predominantly experience growth with infrequent mergers or interactions.
△ Less
Submitted 15 January, 2024; v1 submitted 14 December, 2023;
originally announced December 2023.
-
The THESAN project: connecting ionized bubble sizes to their local environments during the Epoch of Reionization
Authors:
Meredith Neyer,
Aaron Smith,
Rahul Kannan,
Mark Vogelsberger,
Enrico Garaldi,
Daniela Galárraga-Espinosa,
Josh Borrow,
Lars Hernquist,
Rüdiger Pakmor,
Volker Springel
Abstract:
An important characteristic of cosmic hydrogen reionization is the growth of ionized gas bubbles surrounding early luminous objects. Ionized bubble sizes are beginning to be probed using Lyman-$α$ emission from high-redshift galaxies, and will also be probed by upcoming 21-cm maps. We present results from a study of bubble sizes using the state-of-the-art THESAN radiation-hydrodynamics simulation…
▽ More
An important characteristic of cosmic hydrogen reionization is the growth of ionized gas bubbles surrounding early luminous objects. Ionized bubble sizes are beginning to be probed using Lyman-$α$ emission from high-redshift galaxies, and will also be probed by upcoming 21-cm maps. We present results from a study of bubble sizes using the state-of-the-art THESAN radiation-hydrodynamics simulation suite, which self-consistently models radiation transport and realistic galaxy formation. We employ the mean-free path method, and track the evolution of the effective ionized bubble size at each point ($R_{\rm eff}$) throughout the Epoch of Reionization. We show there is a slow growth period for regions ionized early, but a rapid "flash ionization" process for regions ionized later as they immediately enter a large, pre-existing bubble. We also find that bright sources are preferentially in larger bubbles, and find consistency with recent observational constraints at $z \gtrsim 9$, but tension with idealized Lyman-$α$ damping-wing models at $z \approx 7$. We find that high overdensity regions have larger characteristic bubble sizes, but the correlation decreases as reionization progresses, likely due to runaway formation of large percolated bubbles. Finally, we compare the redshift at which a region transitions from neutral to ionized ($z_{\rm reion}$) with the time it takes to reach a given bubble size and conclude that $z_{\rm reion}$ is a reasonable local probe of small-scale bubble size statistics ($R_\text{eff} \lesssim 1\,\rm{cMpc}$). However, for larger bubbles, the correspondence between $z_{\rm reion}$ and size statistics weakens due to the time delay between the onset of reionization and the expansion of large bubbles, particularly at high redshifts.
△ Less
Submitted 5 June, 2024; v1 submitted 5 October, 2023;
originally announced October 2023.
-
Statistics of thermal gas pressure as a probe of cosmology and galaxy formation
Authors:
Ziyang Chen,
Drew Jamieson,
Eiichiro Komatsu,
Sownak Bose,
Klaus Dolag,
Boryana Hadzhiyska,
César Hernández-Aguayo,
Lars Hernquist,
Rahul Kannan,
Rüediger Pakmor,
Volker Springel
Abstract:
The statistics of thermal gas pressure are a new and promising probe of cosmology and astrophysics. The large-scale cross-correlation between galaxies and the thermal Sunyaev-Zeldovich effect gives the bias-weighted mean electron pressure, $\langle b_\mathrm{h}P_e\rangle$. In this paper, we show that $\langle b_\mathrm{h}P_e\rangle$ is sensitive to the amplitude of fluctuations in matter density,…
▽ More
The statistics of thermal gas pressure are a new and promising probe of cosmology and astrophysics. The large-scale cross-correlation between galaxies and the thermal Sunyaev-Zeldovich effect gives the bias-weighted mean electron pressure, $\langle b_\mathrm{h}P_e\rangle$. In this paper, we show that $\langle b_\mathrm{h}P_e\rangle$ is sensitive to the amplitude of fluctuations in matter density, for example $\langle b_\mathrm{h}P_e\rangle\propto \left(σ_8Ω_\mathrm{m}^{0.81}h^{0.67}\right)^{3.14}$ at redshift $z=0$. We find that at $z<0.5$ the observed $\langle b_\mathrm{h}P_e\rangle$ is smaller than that predicted by the state-of-the-art hydrodynamical simulations of galaxy formation, MillenniumTNG, by a factor of $0.93$. This can be explained by a lower value of $σ_8$ and $Ω_\mathrm{m}$, similar to the so-called "$S_8$ tension'' seen in the gravitational lensing effect, although the influence of astrophysics cannot be completely excluded. The difference between Magneticum and MillenniumTNG at $z<2$ is small, indicating that the difference in the galaxy formation models used by these simulations has little impact on $\langle b_\mathrm{h}P_e\rangle$ at this redshift range. At higher $z$, we find that both simulations are in a modest tension with the existing upper bounds on $\langle b_\mathrm{h}P_e\rangle$. We also find a significant difference between these simulations there, which we attribute to a larger sensitivity to the galaxy formation models in the high redshift regime. Therefore, more precise measurements of $\langle b_\mathrm{h}P_e\rangle$ at all redshifts will provide a new test of our understanding of cosmology and galaxy formation.
△ Less
Submitted 28 September, 2023;
originally announced September 2023.
-
Non-ideal magnetohydrodynamics on a moving mesh II: Hall effect
Authors:
Oliver Zier,
Alexander C. Mayer,
Volker Springel
Abstract:
In this work we extend the non-ideal magnetohydrodynamics (MHD) solver in the moving mesh code AREPO to include the Hall effect. The core of our algorithm is based on an estimation of the magnetic field gradients by a least-square reconstruction on the unstructured mesh, which we also used in the companion paper for Ohmic and ambipolar diffusion. In an extensive study of simulations of a magnetic…
▽ More
In this work we extend the non-ideal magnetohydrodynamics (MHD) solver in the moving mesh code AREPO to include the Hall effect. The core of our algorithm is based on an estimation of the magnetic field gradients by a least-square reconstruction on the unstructured mesh, which we also used in the companion paper for Ohmic and ambipolar diffusion. In an extensive study of simulations of a magnetic shock, we show that without additional magnetic diffusion our algorithm for the Hall effect becomes unstable at high resolution. We can however stabilise it by artificially increasing the Ohmic resistivity, $η_{\rm OR}$, so that it satisfies the condition $η_{\rm OR} \geq η_{\rm H} /5$, where $η_{\rm H}$ is the Hall diffusion coefficient. Adopting this solution we find second order convergence for the C-shock and are also able to accurately reproduce the dispersion relation of the whistler waves. As a first application of the new scheme, we simulate the collapse of a magnetised cloud with constant Hall parameter $η_{\rm H}$ and show that, depending on the sign of $η_{\rm H}$, the magnetic braking can either be weakened or strengthened by the Hall effect. The quasi-Lagrangian nature of the moving mesh method used here automatically increases the resolution in the forming core, making it well suited for more realistic studies with non-constant magnetic diffusivities in the future.
△ Less
Submitted 27 September, 2023;
originally announced September 2023.
-
Magnetic field amplification in cosmological zoom simulations from dwarf galaxies to galaxy groups
Authors:
Ruediger Pakmor,
Rebekka Bieri,
Freeke van de Voort,
Maria Werhahn,
Azadeh Fattahi,
Thomas Guillet,
Christoph Pfrommer,
Volker Springel,
Rosie Y. Talbot
Abstract:
Magnetic fields are ubiquitous in the Universe. Recently, cosmological simulations of galaxies have successfully begun to incorporate magnetic fields and their evolution in galaxies and their haloes. However, so far they have mostly focused on Milky Way-like galaxies. Here we analyse a sample of high resolution cosmological zoom simulations of disc galaxies in haloes with mass $M_\mathrm{200c}$ fr…
▽ More
Magnetic fields are ubiquitous in the Universe. Recently, cosmological simulations of galaxies have successfully begun to incorporate magnetic fields and their evolution in galaxies and their haloes. However, so far they have mostly focused on Milky Way-like galaxies. Here we analyse a sample of high resolution cosmological zoom simulations of disc galaxies in haloes with mass $M_\mathrm{200c}$ from $10^{10}\,\mathrm{M}_\odot$ to $10^{13}\,\mathrm{M}_\odot$, simulated with the Auriga galaxy formation model. We show that with sufficient numerical resolution the magnetic field amplification and saturation is converged. The magnetic field strength reaches equipartition with turbulent energy density for galaxies in haloes with $M_\mathrm{200c}\gtrsim 10^{11.5}\,\mathrm{M_\odot}$. For galaxies in less massive haloes, the magnetic field strength saturates at a fraction of equipartition that decreases with decreasing halo mass. For our lowest mass haloes, the magnetic field saturates significantly below $10\%$ of equipartition. We quantify the resolution we need to obtain converged magnetic field strengths and discuss our resolution requirements also in the context of the IllustrisTNG cosmological box simulations. We show that, at $z=0$, rotation-dominated galaxies in our sample exhibit for the most part an ordered large scale magnetic field, with fewer field reversals in more massive galaxies. Finally, we compare the magnetic fields in our cosmological galaxies at $z=0$ with simulations of isolated galaxies in a collapsing halo setup. Our results pave the way for detailed studies of cosmic rays and other physical processes in similar cosmological galaxy simulations that crucially depend on the strength and structure of magnetic fields.
△ Less
Submitted 9 January, 2024; v1 submitted 22 September, 2023;
originally announced September 2023.
-
Evolution of cosmic filaments in the MTNG simulation
Authors:
Daniela Galárraga-Espinosa,
Corentin Cadiou,
Céline Gouin,
Simon D. M. White,
Volker Springel,
Rüdiger Pakmor,
Boryana Hadzhiyska,
Sownak Bose,
Fulvio Ferlito,
Lars Hernquist,
Rahul Kannan,
Monica Barrera,
Ana Maria Delgado,
César Hernández-Aguayo
Abstract:
We present a study of the evolution of cosmic filaments across redshift with an emphasis on some important properties: filament lengths, growth rates, and radial profiles of galaxy densities. Following an observation-driven approach, we build cosmic filament catalogues at z=0,1,2,3, and 4 from the galaxy distributions of the large hydro-dynamical run of the MilleniumTNG project. We employ the exte…
▽ More
We present a study of the evolution of cosmic filaments across redshift with an emphasis on some important properties: filament lengths, growth rates, and radial profiles of galaxy densities. Following an observation-driven approach, we build cosmic filament catalogues at z=0,1,2,3, and 4 from the galaxy distributions of the large hydro-dynamical run of the MilleniumTNG project. We employ the extensively used DisPerSE cosmic web finder code, for which we provide a user-friendly guide, including the details of a physics-driven calibration procedure, with the hope of helping future users. We perform the first statistical measurements of the evolution of connectivity in a large-scale simulation, finding that the connectivity of cosmic nodes (defined as the number of filaments attached) globally decreases from early to late times. The study of cosmic filaments in proper coordinates reveals that filaments grow in length and radial extent, as expected from large-scale structures in an expanding Universe. But the most interesting results arise once the Hubble flow is factored out. We find remarkably stable comoving filament length functions and over-density profiles, showing only little evolution of the total population of filaments in the past ~12.25 Gyrs. However, by tracking the spatial evolution of individual structures, we demonstrate that filaments of different lengths actually follow different evolutionary paths. While short filaments preferentially contract, long filaments expand along their longitudinal direction with growth rates that are the highest in the early, matter-dominated Universe. Filament diversity at fixed redshift is also shown by the different (~$5 σ$) density values between the shortest and longest filaments. Our results hint that cosmic filaments can be used as additional probes for dark energy, but further theoretical work is still needed.
△ Less
Submitted 12 January, 2024; v1 submitted 15 September, 2023;
originally announced September 2023.
-
The thesan project: public data release of radiation-hydrodynamic simulations matching reionization-era JWST observations
Authors:
Enrico Garaldi,
Rahul Kannan,
Aaron Smith,
Josh Borrow,
Mark Vogelsberger,
Rüdiger Pakmor,
Volker Springel,
Lars Hernquist,
Daniela Galárraga-Espinosa,
Jessica Y. -C. Yeh,
Xuejian Shen,
Clara Xu,
Meredith Neyer,
Benedetta Spina,
Mouza Almualla,
Yu Zhao
Abstract:
Cosmological simulations serve as invaluable tools for understanding the Universe. However, the technical complexity and substantial computational resources required to generate such simulations often limit their accessibility within the broader research community. Notable exceptions exist, but most are not suited for simultaneously studying the physics of galaxy formation and cosmic reionization…
▽ More
Cosmological simulations serve as invaluable tools for understanding the Universe. However, the technical complexity and substantial computational resources required to generate such simulations often limit their accessibility within the broader research community. Notable exceptions exist, but most are not suited for simultaneously studying the physics of galaxy formation and cosmic reionization during the first billion years of cosmic history. This is especially relevant now that a fleet of advanced observatories (e.g. James Webb Space Telescope, Nancy Grace Roman Space Telescope, SPHEREx, ELT, SKA) will soon provide an holistic picture of this defining epoch. To bridge this gap, we publicly release all simulation outputs and post-processing products generated within the THESAN simulation project at https://thesan-project.com. This project focuses on the $z \geq 5.5$ Universe, combining a radiation-hydrodynamics solver (AREPO-RT), a well-tested galaxy formation model (IllustrisTNG) and cosmic dust physics to provide a comprehensive view of the Epoch of Reionization. The THESAN suite includes 16 distinct simulations, each varying in volume, resolution, and underlying physical models. This paper outlines the unique features of these new simulations, the production and detailed format of the wide range of derived data products, and the process for data retrieval. Finally, as a case study, we compare our simulation data with a number of recent observations from the James Webb Space Telescope, affirming the accuracy and applicability of THESAN. The examples also serve as prototypes for how to utilise the released dataset to perform comparisons between predictions and observations.
△ Less
Submitted 21 March, 2024; v1 submitted 12 September, 2023;
originally announced September 2023.
-
Non-ideal magnetohydrodynamics on a moving mesh I: Ohmic and ambipolar diffusion
Authors:
Oliver Zier,
Volker Springel,
Alexander C. Mayer
Abstract:
Especially in cold and high-density regions, the assumptions of ideal magnetohydrodynamics (MHD) can break down, making first order non-ideal terms such as Ohmic and ambipolar diffusion as well as the Hall effect important. In this study we present a new numerical scheme for the first two resistive terms, which we implement in the moving-mesh code AREPO using the single-fluid approximation combine…
▽ More
Especially in cold and high-density regions, the assumptions of ideal magnetohydrodynamics (MHD) can break down, making first order non-ideal terms such as Ohmic and ambipolar diffusion as well as the Hall effect important. In this study we present a new numerical scheme for the first two resistive terms, which we implement in the moving-mesh code AREPO using the single-fluid approximation combined with a new gradient estimation technique based on a least-squares fit per interface. Through various test calculations including the diffusion of a magnetic peak, the structure of a magnetic C-shock, and the damping of an Alfvén wave, we show that we can achieve an accuracy comparable to the state-of-the-art code ATHENA++. We apply the scheme to the linear growth of the magnetorotational instability and find good agreement with the analytical growth rates. By simulating the collapse of a magnetised cloud with constant magnetic diffusion, we show that the new scheme is stable even for large density contrasts. Thanks to the Lagrangian nature of the moving mesh method the new scheme is thus well suited for intended future applications where a high resolution in the dense cores of collapsing protostellar clouds needs to be achieved. In a forthcoming work we will extend the scheme to the Hall effect.
△ Less
Submitted 21 July, 2023;
originally announced July 2023.
-
A reproduction of the Milky Way's Faraday rotation measure map in galaxy simulations from global to local scales
Authors:
Stefan Reissl,
Ralf S. Klessen,
Eric W. Pellegrini,
Daniel Rahner,
Rüdiger Pakmor,
Robert Grand,
Facundo Gomez,
Federico Marinacci,
Volker Springel
Abstract:
Magnetic fields are of critical importance for our understanding of the origin and long-term evolution of the Milky Way. This is due to their decisive role in the dynamical evolution of the interstellar medium (ISM) and their influence on the star-formation process. Faraday rotation measures (RM) along many different sightlines across the Galaxy are a primary means to infer the magnetic field topo…
▽ More
Magnetic fields are of critical importance for our understanding of the origin and long-term evolution of the Milky Way. This is due to their decisive role in the dynamical evolution of the interstellar medium (ISM) and their influence on the star-formation process. Faraday rotation measures (RM) along many different sightlines across the Galaxy are a primary means to infer the magnetic field topology and strength from observations. However, the interpretation of the data has been hampered by the failure of previous attempts to explain the observations in theoretical models and to synthesize a realistic multi-scale all-sky RM map. We here utilize a cosmological magnetohydrodynamic (MHD) simulation of the formation of the Milky Way, augment it with a novel star cluster population synthesis model for a more realistic structure of the local interstellar medium, and perform detailed polarized radiative transfer calculations on the resulting model. This yields a faithful first principles prediction of the Faraday sky as observed on Earth. The results reproduce the observations of the Galaxy not only on global scales, but also on local scales of individual star-forming clouds. They also imply that the Local Bubble containing our Sun dominates the RM signal over large regions of the sky. Modern cosmological MHD simulations of the Milky Way's formation, combined with a simple and plausible model for the fraction of free electrons in the ISM, explain the RM observations remarkably well, thus indicating the emergence of a firm theoretical understanding of the genesis of magnetic fields in our Universe across cosmic time.
△ Less
Submitted 11 July, 2023;
originally announced July 2023.
-
Where shadows lie: reconstruction of anisotropies in the neutrino sky
Authors:
Willem Elbers,
Carlos S. Frenk,
Adrian Jenkins,
Baojiu Li,
Silvia Pascoli,
Jens Jasche,
Guilhem Lavaux,
Volker Springel
Abstract:
The Cosmic Neutrino Background (CNB) encodes a wealth of information, but has not yet been observed directly. To determine the prospects of detection and to study its information content, we reconstruct the phase-space distribution of local relic neutrinos from the three-dimensional distribution of matter within 200 Mpc/h of the Milky Way. Our analysis relies on constrained realization simulations…
▽ More
The Cosmic Neutrino Background (CNB) encodes a wealth of information, but has not yet been observed directly. To determine the prospects of detection and to study its information content, we reconstruct the phase-space distribution of local relic neutrinos from the three-dimensional distribution of matter within 200 Mpc/h of the Milky Way. Our analysis relies on constrained realization simulations and forward modelling of the 2M++ galaxy catalogue. We find that the angular distribution of neutrinos is anti-correlated with the projected matter density, due to the capture and deflection of neutrinos by massive structures along the line of sight. Of relevance to tritium capture experiments, we find that the gravitational clustering effect of the large-scale structure on the local number density of neutrinos is more important than that of the Milky Way for neutrino masses less than 0.1 eV. Nevertheless, we predict that the density of relic neutrinos is close to the cosmic average, with a suppression or enhancement over the mean of (-0.3%, +7%, +27%) for masses of (0.01, 0.05, 0.1) eV. This implies no more than a marginal increase in the event rate for tritium capture experiments like PTOLEMY. We also predict that the CNB and CMB rest frames coincide for 0.01 eV neutrinos, but that neutrino velocities are significantly perturbed for masses larger than 0.05 eV. Regardless of mass, we find that the angle between the neutrino dipole and the ecliptic plane is small, implying a near-maximal annual modulation in the bulk velocity. Along with this paper, we publicly release our simulation data, comprising more than 100 simulations for six different neutrino masses.
△ Less
Submitted 28 December, 2023; v1 submitted 6 July, 2023;
originally announced July 2023.
-
Close Encounters of Star - Black Hole Binaries with Single Stars
Authors:
Taeho Ryu,
Selma de Mink,
Rob Farmer,
Ruediger Pakmor,
Rosalba Perna,
Volker Springel
Abstract:
Multi-body dynamical interactions of binaries with other objects are one of the main driving mechanisms for the evolution of star clusters. It is thus important to bring our understanding of three-body interactions beyond the commonly employed point-particle approximation. To this end we here investigate the hydrodynamics of three-body encounters between star-black hole (BH) binaries and single st…
▽ More
Multi-body dynamical interactions of binaries with other objects are one of the main driving mechanisms for the evolution of star clusters. It is thus important to bring our understanding of three-body interactions beyond the commonly employed point-particle approximation. To this end we here investigate the hydrodynamics of three-body encounters between star-black hole (BH) binaries and single stars, focusing on the identification of final outcomes and their long-term evolution and observational properties, using the moving-mesh hydrodynamics code AREPO. This type of encounters produces five types of outcomes: stellar disruption, stellar collision, weak perturbation of the original binary, binary member exchange, and triple formation. The two decisive parameters are the binary phase angle, which determines which two objects meet at the first closest approach, and the impact parameter, which sets the boundary between violent and non-violent interactions. When the impact parameter is smaller than the semimajor axis of the binary, tidal disruptions and star-BH collisions frequently occur when the BH and the incoming star first meet, while the two stars mostly merge when the two stars meet first instead. In both cases, the BHs accrete from an accretion disk at super-Eddington rates, possibly generating flares luminous enough to be observed. The stellar collision products either form a binary with the BH or remain unbound to the BH. Upon collision, the merged stars are hotter and larger than main sequence stars of the same mass at similar age. Even after recovering their thermal equilibrium state, stellar collision products, if isolated, would remain hotter and brighter than main sequence stars until becoming giants.
△ Less
Submitted 10 July, 2023; v1 submitted 6 July, 2023;
originally announced July 2023.
-
KETJU -- resolving small-scale supermassive black hole dynamics in GADGET-4
Authors:
Matias Mannerkoski,
Alexander Rawlings,
Peter H. Johansson,
Thorsten Naab,
Antti Rantala,
Volker Springel,
Dimitrios Irodotou,
Shihong Liao
Abstract:
We present the new public version of the KETJU supermassive black hole (SMBH) dynamics module, as implemented into GADGET-4. KETJU adds a small region around each SMBH where the dynamics of the SMBHs and stellar particles are integrated using an algorithmically regularised integrator instead of the leapfrog integrator with gravitational softening used by GADGET-4. This enables modelling SMBHs as p…
▽ More
We present the new public version of the KETJU supermassive black hole (SMBH) dynamics module, as implemented into GADGET-4. KETJU adds a small region around each SMBH where the dynamics of the SMBHs and stellar particles are integrated using an algorithmically regularised integrator instead of the leapfrog integrator with gravitational softening used by GADGET-4. This enables modelling SMBHs as point particles even during close interactions with stellar particles or other SMBHs, effectively removing the spatial resolution limitation caused by gravitational softening. KETJU also includes post-Newtonian corrections, which allows following the dynamics of SMBH binaries to sub-parsec scales and down to tens of Schwarzschild radii. Systems with multiple SMBHs are also supported, with the code also including the leading non-linear cross terms that appear in the post-Newtonian equations for such systems. We present tests of the code showing that it correctly captures, at sufficient mass resolution, the sinking driven by dynamical friction and binary hardening driven by stellar scattering. We also present an example application demonstrating how the code can be applied to study the dynamics of SMBHs in mergers of multiple galaxies and the effect they have on the properties of the surrounding galaxy. We expect that the presented KETJU SMBH dynamics module can also be straightforwardly incorporated into other codes similar to GADGET-4, which would allow coupling small-scale SMBH dynamics to the rich variety of galactic physics models that exist in the literature.
△ Less
Submitted 8 June, 2023;
originally announced June 2023.
-
The Structure and Dynamics of Massive High-$z$ Cosmic-Web Filaments: Three Radial Zones in Filament Cross-Sections
Authors:
Yue Samuel Lu,
Nir Mandelker,
S. Peng Oh,
Avishai Dekel,
Frank C. van den Bosch,
Volker Springel,
Daisuke Nagai,
Freeke van de Voort
Abstract:
We analyse the internal structure and dynamics of cosmic-web filaments that connect massive high-$z$ haloes. Our analysis is based on a high-resolution AREPO cosmological simulation zooming-in on a volume encompassing three ${\rm Mpc}$-scale filaments feeding three massive haloes of $\sim 10^{12}\,\text{M}_\odot$ at $z \sim 4$, embedded in a large-scale sheet. Each filament is surrounded by a cyli…
▽ More
We analyse the internal structure and dynamics of cosmic-web filaments that connect massive high-$z$ haloes. Our analysis is based on a high-resolution AREPO cosmological simulation zooming-in on a volume encompassing three ${\rm Mpc}$-scale filaments feeding three massive haloes of $\sim 10^{12}\,\text{M}_\odot$ at $z \sim 4$, embedded in a large-scale sheet. Each filament is surrounded by a cylindrical accretion shock of radius $r_{\rm shock} \sim 50 \,{\rm kpc}$. The post-shock gas is in virial equilibrium with the potential well set by an isothermal dark-matter filament. The filament line-mass is $\sim 9\times 10^8\,\text{M}_\odot\,{\rm kpc}^{-1}$, the gas fraction within $r_{\rm shock}$ is the universal baryon fraction, and the virial temperature is $\sim 7\times 10^5 {\rm K}$. In the outer ''thermal'' (T) zone, $r \geq 0.65 \, r_{\rm shock}$, inward gravity and ram-pressure forces are over-balanced by outwards thermal pressure forces, decelerating the inflowing gas expanding the shock outward. In the intermediate ''vortex'' (V) zone, $0.25 \leq r/ r_{\rm shock} \leq 0.65$, the velocity field is dominated by a quadrupolar vortex structure due to offset inflow along the sheet through the post-shock gas. The outwards force is dominated by centrifugal forces associated with these vortices, with additional contributions from global rotation and thermal pressure. The shear and turbulent forces associated with the vortices act inward. The inner ''stream'' (S) zone, $r < 0.25 \, r_{\rm shock}$, is a dense isothermal core, $T\sim 3 \times 10^4 \, {\rm K}$ and $n_{\rm H}\sim 0.01 \,{\rm cm^{-3}}$, defining the cold streams that feed galaxies. The core is formed by an isobaric cooling flow and is associated with a decrease in outwards forces, though it exhibits both inflows and outflows. [abridged]
△ Less
Submitted 20 November, 2023; v1 submitted 6 June, 2023;
originally announced June 2023.
-
Magnetic Fields in Cosmic Voids
Authors:
Agustín M. Rodríguez-Medrano,
Federico A. Stasyszyn,
Dante J. Paz,
Volker Springel
Abstract:
Magnetic fields are one of most concealed components of the universe. They are observed as part of the intergalactic medium and on galaxy cluster scales, however their origin and evolution is unclear. In this work we use the IllustrisTNG simulation to investigate the effects of magnetic fields in cosmic voids, the least dense regions of the universe. We find that, under the hypothesis of a uniform…
▽ More
Magnetic fields are one of most concealed components of the universe. They are observed as part of the intergalactic medium and on galaxy cluster scales, however their origin and evolution is unclear. In this work we use the IllustrisTNG simulation to investigate the effects of magnetic fields in cosmic voids, the least dense regions of the universe. We find that, under the hypothesis of a uniform primordial magnetic field, the voids still reflect the primordial properties of the fields. On the other hand, the galaxies in their interior acquire weaker magnetic fields than galaxies in denser environments.
△ Less
Submitted 2 May, 2023;
originally announced May 2023.
-
Interpreting Sunyaev-Zel'dovich observations with MillenniumTNG: Mass and environment scaling relations
Authors:
Boryana Hadzhiyska,
Simone Ferraro,
Rüdiger Pakmor,
Sownak Bose,
Ana Maria Delgado,
César Hernández-Aguayo,
Rahul Kannan,
Volker Springel,
Simon D. M. White,
Lars Hernquist
Abstract:
In the coming years, Sunyaev-Zel'dovich (SZ) measurements can dramatically improve our understanding of the Intergalactic Medium (IGM) and the role of feedback processes on galaxy formation, allowing us to calibrate important astrophysical systematics in cosmological constraints from weak lensing galaxy clustering surveys. However, the signal is only measured in a two-dimensional projection, and i…
▽ More
In the coming years, Sunyaev-Zel'dovich (SZ) measurements can dramatically improve our understanding of the Intergalactic Medium (IGM) and the role of feedback processes on galaxy formation, allowing us to calibrate important astrophysical systematics in cosmological constraints from weak lensing galaxy clustering surveys. However, the signal is only measured in a two-dimensional projection, and its correct interpretation relies on understanding the connection between observable quantities and the underlying intrinsic properties of the gas, in addition to the relation between the gas and the underlying matter distribution. One way to address these challenges is through the use of hydrodynamical simulations such as the high-resolution, large-volume MillenniumTNG suite. We find that measurements of the optical depth, $τ$, and the Compton-y parameter, $Y$, receive large line-of-sight contributions which can be removed effectively by applying a Compensated Aperture Photometry (CAP) filter. In contrast with other $τ$ probes (e.g., X-rays and Fast Radio Bursts), the kSZ-inferred $τ$ receives most of its signal from a confined cylindrical region around the halo due to the velocity decorrelation along the line-of-sight. Additionally, we perform fits to the $Y-M$ and $τ-M$ scaling relations and report best-fit parameters adopting the smoothly broken power law (SBPL) formalism. We note that subgrid physics modeling can broaden the error bar on these by 30\% for intermediate-mass halos ($\sim$$10^{13} \, {\rm M}_{\odot}$). The scatter of the scaling relations can be captured by an intrinsic dependence on concentration, and an extrinsic dependence on tidal shear. Finally, we comment on the effect of using galaxies rather than halos in real observations, which can bias the inferred SZ profiles by $\sim$20\% for $L_\ast$-galaxies.
△ Less
Submitted 1 May, 2023;
originally announced May 2023.
-
The MillenniumTNG Project: Intrinsic alignments of galaxies and halos
Authors:
Ana Maria Delgado,
Boryana Hadzhiyska,
Sownak Bose,
Volker Springel,
Lars Hernquist,
Monica Barrer,
Rüdiger Pakmor,
Fulvio Ferlito,
Rahul Kannan,
César Hernández-Aguayo,
Simon D. M. White,
Carlos Frenk
Abstract:
The intrinsic alignment (IA) of observed galaxy shapes with the underlying cosmic web is a source of contamination in weak lensing surveys. Sensitive methods to identify the IA signal will therefore need to be included in the upcoming weak lensing analysis pipelines. Hydrodynamical cosmological simulations allow us to directly measure the intrinsic ellipticities of galaxies and thus provide a powe…
▽ More
The intrinsic alignment (IA) of observed galaxy shapes with the underlying cosmic web is a source of contamination in weak lensing surveys. Sensitive methods to identify the IA signal will therefore need to be included in the upcoming weak lensing analysis pipelines. Hydrodynamical cosmological simulations allow us to directly measure the intrinsic ellipticities of galaxies and thus provide a powerful approach to predict and understand the IA signal. Here we employ the novel, large-volume hydrodynamical simulation MTNG740, a product of the MillenniumTNG (MTNG) project, to study the IA of galaxies. We measure the projected correlation functions between the intrinsic shape/shear of galaxies and various tracers of large-scale structure, $w_{+g},\ w_{+m},\ w_{++}$ over the radial range $r_{\rm p} \in [0.02 , 200]\,h^{-1}{\rm Mpc}$ and at redshifts $z=0.0$, $0.5$ and $1.0$. We detect significant signal-to-noise IA signals with the density field for both elliptical and spiral galaxies. We also find significant intrinsic shear-shear correlations for ellipticals. We further examine correlations of the intrinsic shape of galaxies with the local tidal field. Here we find a significant IA signal for elliptical galaxies assuming a linear model. We also detect a weak IA signal for spiral galaxies under a quadratic tidal torquing model. Lastly, we measure the alignment between central galaxies and their host dark-matter halos, finding small to moderate misalignments between their principal axes that decline with halo mass.
△ Less
Submitted 20 July, 2023; v1 submitted 24 April, 2023;
originally announced April 2023.
-
The MillenniumTNG Project: The impact of baryons and massive neutrinos on high-resolution weak gravitational lensing convergence maps
Authors:
Fulvio Ferlito,
Volker Springel,
Christopher T. Davies,
César Hernández-Aguayo,
Rüdiger Pakmor,
Monica Barrera,
Simon D. M. White,
Ana Maria Delgado,
Boryana Hadzhiyska,
Lars Hernquist,
Rahul Kannan,
Sownak Bose,
Carlos Frenk
Abstract:
We study weak gravitational lensing convergence maps produced from the MillenniumTNG (MTNG) simulations by direct projection of the mass distribution on the past backwards lightcone of a fiducial observer. We explore the lensing maps over a large dynamic range in simulation mass and angular resolution, allowing us to establish a clear assessment of numerical convergence. By comparing full physics…
▽ More
We study weak gravitational lensing convergence maps produced from the MillenniumTNG (MTNG) simulations by direct projection of the mass distribution on the past backwards lightcone of a fiducial observer. We explore the lensing maps over a large dynamic range in simulation mass and angular resolution, allowing us to establish a clear assessment of numerical convergence. By comparing full physics hydrodynamical simulations with corresponding dark-matter-only runs we quantify the impact of baryonic physics on the most important weak lensing statistics. Likewise, we predict the impact of massive neutrinos reliably far into the non-linear regime. We also demonstrate that the "fixed & paired" variance suppression technique increases the statistical robustness of the simulation predictions on large scales not only for time slices but also for continuously output lightcone data. We find that both baryonic and neutrino effects substantially impact weak lensing shear measurements, with the latter dominating over the former on large angular scales. Thus, both effects must explicitly be included to obtain sufficiently accurate predictions for stage IV lensing surveys. Reassuringly, our results agree accurately with other simulation results where available, supporting the promise of simulation modelling for precision cosmology far into the non-linear regime.
△ Less
Submitted 14 June, 2024; v1 submitted 24 April, 2023;
originally announced April 2023.
-
Close encounters of black hole - star binaries with stellar-mass black holes
Authors:
Taeho Ryu,
Ruggero Valli,
Rudiger Pakmor,
Rosalba Perna,
Selma E. de Mink,
Volker Springel
Abstract:
Dynamical interactions involving binaries play a crucial role in the evolution of star clusters and galaxies. We continue our investigation of the hydrodynamics of three-body encounters, focusing on binary black hole (BBH) formation, stellar disruption, and electromagnetic (EM) emission in dynamical interactions between a BH-star binary and a stellar-mass BH, using the moving-mesh hydrodynamics co…
▽ More
Dynamical interactions involving binaries play a crucial role in the evolution of star clusters and galaxies. We continue our investigation of the hydrodynamics of three-body encounters, focusing on binary black hole (BBH) formation, stellar disruption, and electromagnetic (EM) emission in dynamical interactions between a BH-star binary and a stellar-mass BH, using the moving-mesh hydrodynamics code {\small AREPO}. This type of encounters can be divided into two classes depending on whether the final outcome includes BBHs. This outcome is primarily determined by which two objects meet at the first closest approach. BBHs are more likely to form when the star and the incoming BH encounter first with an impact parameter smaller than the binary's semimajor axis. In this case, the star is frequently disrupted. On the other hand, when the two BHs encounter first, frequent consequences are an orbit perturbation of the original binary or a binary member exchange. For the parameters chosen in this study, BBH formation, accompanied by stellar disruption, happens in roughly 1 out of 4 encounters. The close correlation between BBH formation and stellar disruption has possible implications for EM counterparts at the binary's merger. The BH that disrupts the star is promptly surrounded by an optically and geometrically thick disk with accretion rates exceeding the Eddington limit. If the debris disk cools fast enough to become long-lived, EM counterparts can be produced at the time of the BBH merger.
△ Less
Submitted 10 July, 2023; v1 submitted 4 April, 2023;
originally announced April 2023.
-
Disk flaring with TNG50: diversity across Milky Way and M31 analogs
Authors:
Diego Sotillo-Ramos,
Martina Donnari,
Annalisa Pillepich,
Neige Frankel,
Dylan Nelson,
Volker Springel,
Lars Hernquist
Abstract:
We use the sample of 198 Milky Way (MW) and Andromeda (M31) analogs from TNG50 to quantify the level of disk flaring predicted by a modern, high-resolution cosmological hydrodynamical simulation. Disk flaring refers to the increase of vertical stellar disk height with galactocentric distance. The TNG50 galaxies are selected to have stellar disky morphology, a stellar mass in the range of…
▽ More
We use the sample of 198 Milky Way (MW) and Andromeda (M31) analogs from TNG50 to quantify the level of disk flaring predicted by a modern, high-resolution cosmological hydrodynamical simulation. Disk flaring refers to the increase of vertical stellar disk height with galactocentric distance. The TNG50 galaxies are selected to have stellar disky morphology, a stellar mass in the range of $M_* = 10^{10.5 - 11.2}~\rm{M_{\odot}}$, and a MW-like Mpc-scale environment at $z=0$. The stellar disks of such TNG50 MW/M31 analogs exhibit a wide diversity of structural properties, including a number of galaxies with disk scalelength and thin and thick disk scaleheights that are comparable to those measured or inferred for the Galaxy and Andromeda. With one set of physical ingredients, TNG50 returns a large variety of flaring flavours and amounts, also for mono-age stellar populations. With this paper, we hence propose a non-parametric characterization of flaring. The typical MW/M31 analogs exhibit disk scaleheights that are $1.5-2$ times larger in the outer than in the inner regions of the disk for both old and young stellar populations, but with a large galaxy-to-galaxy variation. Which stellar population flares more, and by how much, also varies from galaxy to galaxy. TNG50 de facto brackets existing observational constraints for the Galaxy and all previous numerical findings. A link between the amount of flaring and the $z=0$ global galaxy structural properties or merger history is complex. However, a connection between the scaleheights and the local stellar vertical kinematics and gravitational potential is clearly in place.
△ Less
Submitted 28 March, 2023;
originally announced March 2023.
-
Milky Way and Andromeda analogs from the TNG50 simulation
Authors:
Annalisa Pillepich,
Diego Sotillo-Ramos,
Rahul Ramesh,
Dylan Nelson,
Christoph Engler,
Vicente Rodriguez-Gomez,
Martin Fournier,
Martina Donnari,
Volker Springel,
Lars Hernquist
Abstract:
We present the properties of Milky Way- and Andromeda-like (MW/M31-like) galaxies simulated within TNG50, the highest-resolution run of the IllustrisTNG suite of $Λ$CDM magneto-hydrodynamical simulations. We introduce our fiducial selection for MW/M31 analogs, which we propose for direct usage as well as for reference in future analyses. TNG50 contains 198 MW/M31 analogs, i.e. galaxies with stella…
▽ More
We present the properties of Milky Way- and Andromeda-like (MW/M31-like) galaxies simulated within TNG50, the highest-resolution run of the IllustrisTNG suite of $Λ$CDM magneto-hydrodynamical simulations. We introduce our fiducial selection for MW/M31 analogs, which we propose for direct usage as well as for reference in future analyses. TNG50 contains 198 MW/M31 analogs, i.e. galaxies with stellar disky morphology, with a stellar mass in the range of $M_* = 10^{10.5 - 11.2}$ Msun, and within a MW-like Mpc-scale environment at z=0. These are resolved with baryonic (dark matter) mass resolution of $8.5\times10^4$ Msun ($4.5\times10^5$ Msun) and $\sim150$ pc of average spatial resolution in the star-forming regions: we therefore expand by many factors (2 orders of magnitude) the sample size of cosmologically-simulated analogs with similar ($\times 10$ better) numerical resolution. The majority of TNG50 MW/M31 analogs at $z=0$ exhibit a bar, 60 per cent are star-forming, the sample includes 3 Local Group (LG)-like systems, and a number of galaxies host one or more satellites as massive as e.g. the Magellanic Clouds. Even within such a relatively narrow selection, TNG50 reveals a great diversity in galaxy and halo properties, as well as in past histories. Within the TNG50 sample, it is possible to identify several simulated galaxies whose integral and structural properties are consistent, one or more at a time, with those measured for the Galaxy and Andromeda. With this paper, we document and release a series of broadly applicable data products that build upon the IllustrisTNG public release and aim to facilitate easy access and analysis by public users. These include datacubes across snapshots ($0 \le z \le 7$) for each TNG50 MW/M31-like galaxy, and a series of value-added catalogs that will be continually expanded to provide a convenient and up to date community resource.
△ Less
Submitted 28 March, 2023;
originally announced March 2023.
-
The mass accretion history of dark matter haloes down to Earth mass
Authors:
Yizhou Liu,
Liang Gao,
Sownak Bose,
Carlos S. Frenk,
Adrian Jenkins,
Volker Springel,
Jie Wang,
Simon D. M. White,
Haonan Zheng
Abstract:
We take advantage of the unprecedented dynamical range provided by the "Cosmic-Zoom" project to study the mass accretion history (MAH) of present-day dark matter haloes over the entire mass range present in the $Λ$CDM paradigm when the dark matter is made of weakly interacting massive particles of mass $100\ \mathrm{GeV}$. In particular, we complement previous studies by exploring the MAHs of halo…
▽ More
We take advantage of the unprecedented dynamical range provided by the "Cosmic-Zoom" project to study the mass accretion history (MAH) of present-day dark matter haloes over the entire mass range present in the $Λ$CDM paradigm when the dark matter is made of weakly interacting massive particles of mass $100\ \mathrm{GeV}$. In particular, we complement previous studies by exploring the MAHs of haloes with mass from $10^8\ h^{-1}\mathrm{M_{\odot}}$ down to Earth mass, $10^{-6}\ h^{-1}\mathrm{M_{\odot}}$. The formation redshift of low-mass haloes anti-correlates weakly with mass, peaking at $z=3$ for haloes of mass $10^{-4}\ h^{-1}\mathrm{M_{\odot}}$. Even lower masses are affected by the free-streaming cutoff in the primordial spectrum of density fluctuations and form at lower redshift. We compare MAHs in our simulations with predictions from two analytical models based on the extended Press-Schechter theory (EPS), and three empirical models derived by fitting and extrapolating either results from cosmological $N$-body simulations or Monte Carlo realizations of halo growth. All models fit our simulations reasonably well over the mass range for which they were calibrated. While the empirical models match better for more massive haloes, $M>10^{10}\ h^{-1}\mathrm{M_{\odot}}$, the analytical models do better when extrapolated down to Earth mass. At the higher masses, we explore the correlation between local environment density and MAH, finding that biases are relatively weak, with typical MAHs for haloes in extremely low-density and in typical regions differing by less than $20$ percent at high redshift. We conclude that EPS theory predicts the hierarchical build-up of dark matter haloes quite well over the entire mass range.
△ Less
Submitted 28 March, 2023;
originally announced March 2023.
-
Spectrally resolved cosmic rays -- III. Dynamical impact and properties of the circumgalactic medium
Authors:
Philipp Girichidis,
Maria Werhahn,
Christoph Pfrommer,
Rüdiger Pakmor,
Volker Springel
Abstract:
Cosmic rays (CRs) are dynamically important for the formation and evolution of galaxies by regulating star formation and by powering galactic outflows. However, to what extent CRs regulate galaxy formation depends on the coupling strength of CRs with the ambient plasma and the effective CR transport speed along the magnetic field. Moreover, both properties sensitively depend on the CR momentum, wh…
▽ More
Cosmic rays (CRs) are dynamically important for the formation and evolution of galaxies by regulating star formation and by powering galactic outflows. However, to what extent CRs regulate galaxy formation depends on the coupling strength of CRs with the ambient plasma and the effective CR transport speed along the magnetic field. Moreover, both properties sensitively depend on the CR momentum, which is largely unexplored in three-dimensional hydrodynamical simulations. We perform magneto-hydrodynamical simulations of entire galaxies with masses ranging from $10^{10}$ to $10^{12}\,\mathrm{M}_\odot$ and compare dynamically coupled CRs in the grey approximation with a spectrally resolved model that includes CR momenta from $0.1\,\mathrm{GeV}~c^{-1}$ to $100\,\mathrm{TeV}~c^{-1}$. We find that hadronic cooling of CRs dominates over Alfvén cooling, with the latter emulating CR losses as a result of streaming of CRs down their pressure gradient. While star formation rates and galaxy morphologies are only mildly affected by the spectral CR modelling, mass loading factors of galactic outflows can differ by up to a factor of four in dwarf galaxies. All simulated low-mass halos ($M=10^{10}$, $10^{11}$, and $3\times10^{11}\,\mathrm{M}_\odot$) drive strong outflows, where CR transport is temporally dominated by advection. In contrast, the Milky Way-mass galaxy with $M=10^{12}\,\mathrm{M}_\odot$ does not drive sustained outflows, so that CR transport is entirely dominated by diffusion. The effective energy weighted diffusion coefficients vary by two orders of magnitude from the canonical energy-weighted values of $\langle{D}\rangle_{e_\mathrm{cr}}\sim10^{28}\,\mathrm{cm^2\,s^{-1}}$ in the disc up to $3\times10^{29}\,\mathrm{cm^2\,s^{-1}}$ in the circumgalactic medium, where we observe substantial temperature and CR pressure differences between our grey and spectral CR models.
△ Less
Submitted 23 December, 2023; v1 submitted 6 March, 2023;
originally announced March 2023.
-
Exploring the diversity and similarity of radially anisotropic Milky Way-like stellar haloes: implications for disrupted dwarf galaxy searches
Authors:
Matthew D. A. Orkney,
Chervin F. P. Laporte,
Robert J. J. Grand,
Facundo A. Gómez,
Freeke van de Voort,
Azadeh Fattahi,
Federico Marinacci,
Rüdiger Pakmor,
Francesca Fragkoudi,
Volker Springel
Abstract:
We investigate the properties of mergers comparable to the Gaia-Sausage-Enceladus (GSE) using cosmological hydrodynamical simulations of Milky Way-like galaxies. The merger progenitors span an order of magnitude in their peak stellar mass ($3\times10^8<M_{\star}/\rm{M}_{\odot}<4\times10^9$) and include both rotation and pressure-supported galaxies ($0.10<D/T<0.77$). In a minority of cases, the GSE…
▽ More
We investigate the properties of mergers comparable to the Gaia-Sausage-Enceladus (GSE) using cosmological hydrodynamical simulations of Milky Way-like galaxies. The merger progenitors span an order of magnitude in their peak stellar mass ($3\times10^8<M_{\star}/\rm{M}_{\odot}<4\times10^9$) and include both rotation and pressure-supported galaxies ($0.10<D/T<0.77$). In a minority of cases, the GSE-like debris is comprised of stars from more than one merger progenitor. However, there is a close similarity in their chemodynamical properties and the triaxial shapes of their debris, and so it is not always possible to distinguish them. The merger progenitors host a variety of luminous satellites ($0-8$ with $M_{\star}>10^6\,\rm{M}_{\odot}$), but most of these do not follow the merger to low orbital energies. Between $0-1$ of these satellites may survive to $z=0$, but with no clear signatures of their past association. We show that the fraction of stars originating from GSE-like mergers is reduced for lower metallicities (reaching a minimum around $\text{[Fe/H]} = -2$), and also within $5\,$kpc of the galactic centre. Whilst these central regions are dominated by in-situ stars, the ex-situ fraction trends towards a 100 per cent asymptote when considering the most metal-poor stars ($\text{[Fe/H]}\ll-2.5$). Considering this, its near proximity, and its small volume on the sky, the Galactic centre lends itself as a prime environment in the search for the stars from the earliest galaxies, whilst avoiding contamination from GSE stars.
△ Less
Submitted 7 September, 2023; v1 submitted 3 March, 2023;
originally announced March 2023.
-
THESAN-HR: How does reionization impact early galaxy evolution?
Authors:
Josh Borrow,
Rahul Kannan,
Enrico Garaldi,
Aaron Smith,
Mark Vogelsberger,
Rüdiger Pakmor,
Volker Springel,
Lars Hernquist
Abstract:
Early galaxies were the radiation source for reionization, with the photoheating feedback from the reionization process expected to reduce the efficiency of star formation in low mass haloes. Hence, to fully understand reionization and galaxy formation, we must study their impact on each other. The THESAN project has so far aimed to study the impact of galaxy formation physics on reionization, but…
▽ More
Early galaxies were the radiation source for reionization, with the photoheating feedback from the reionization process expected to reduce the efficiency of star formation in low mass haloes. Hence, to fully understand reionization and galaxy formation, we must study their impact on each other. The THESAN project has so far aimed to study the impact of galaxy formation physics on reionization, but here we present the new THESAN simulations with a factor 50 higher resolution ($m_{\rm b} \approx 10^4$~M$_\odot$) that aim to self-consistently study the back-reaction of reionization on galaxies. By resolving haloes with virial temperatures $T_{\rm vir} < 10^4$~K, we are able to demonstrate that simplistic, spatially-uniform, reionization models are not sufficient to study early galaxy evolution. Comparing the self-consistent THESAN model (employing fully coupled radiation hydrodynamics) to a uniform UV background, we are able to show that galaxies in THESAN are predicted to be larger in physical extent (by a factor $\sim 2$), less metal enriched (by $\sim 0.2$~dex), and less abundant (by a factor $\sim 10$ at $M_{\rm 1500}~=~-10$) by $z=5$. We show that differences in star formation and enrichment patterns lead to significantly different predictions for star formation in low mass haloes, low-metallicity star formation, and even the occupation fraction of haloes. We posit that cosmological galaxy formation simulations aiming to study early galaxy formation $z \gtrsim 3$ must employ a spatially inhomogeneous UV background to accurately reproduce galaxy properties.
△ Less
Submitted 17 August, 2023; v1 submitted 6 December, 2022;
originally announced December 2022.
-
Gravito-turbulence in local disk simulations with an adaptive moving mesh
Authors:
Oliver Zier,
Volker Springel
Abstract:
Self-gravity plays an important role in the evolution of rotationally supported systems such as protoplanetary disks, accretion disks around black holes, or galactic disks, as it can both feed turbulence or lead to gravitational fragmentation. While such systems can be studied in the shearing box approximation with high local resolution, the large density contrasts that are possible in the case of…
▽ More
Self-gravity plays an important role in the evolution of rotationally supported systems such as protoplanetary disks, accretion disks around black holes, or galactic disks, as it can both feed turbulence or lead to gravitational fragmentation. While such systems can be studied in the shearing box approximation with high local resolution, the large density contrasts that are possible in the case of fragmentation still limit the utility of Eulerian codes with constant spatial resolution. In this paper, we present a novel self-gravity solver for the shearing box based on the TreePM method of the moving-mesh code AREPO. The spatial gravitational resolution is adaptive which is important to make full use of the quasi-Lagrangian hydrodynamical resolution of the code. We apply our new implementation to two- and three-dimensional, self-gravitating disks combined with a simple $β$-cooling prescription. For weak cooling we find a steady, gravito-turbulent state, while for strong cooling the formation of fragments is inevitable. To reach convergence for the critical cooling efficiency above which fragmentation occurs, we require a smoothing of the gravitational force in the two dimensional case that mimics the stratification of the three-dimensional simulations. The critical cooling efficiency we find, $β\approx 3$, as well as box-averaged quantities characterizing the gravito-turbulent state, agree well with various previous results in the literature. Interestingly, we observe stochastic fragmentation for $β> 3$, which slightly decreases the cooling efficiency required to observe fragmentation over the lifetime of a protoplanetary disk. The numerical method outlined here appears well suited to study the problem of galactic disks as well as magnetized, self-gravitating disks.
△ Less
Submitted 27 September, 2023; v1 submitted 5 December, 2022;
originally announced December 2022.
-
Euclid: Modelling massive neutrinos in cosmology -- a code comparison
Authors:
J. Adamek,
R. E. Angulo,
C. Arnold,
M. Baldi,
M. Biagetti,
B. Bose,
C. Carbone,
T. Castro,
J. Dakin,
K. Dolag,
W. Elbers,
C. Fidler,
C. Giocoli,
S. Hannestad,
F. Hassani,
C. Hernández-Aguayo,
K. Koyama,
B. Li,
R. Mauland,
P. Monaco,
C. Moretti,
D. F. Mota,
C. Partmann,
G. Parimbelli,
D. Potter
, et al. (111 additional authors not shown)
Abstract:
The measurement of the absolute neutrino mass scale from cosmological large-scale clustering data is one of the key science goals of the Euclid mission. Such a measurement relies on precise modelling of the impact of neutrinos on structure formation, which can be studied with $N$-body simulations. Here we present the results from a major code comparison effort to establish the maturity and reliabi…
▽ More
The measurement of the absolute neutrino mass scale from cosmological large-scale clustering data is one of the key science goals of the Euclid mission. Such a measurement relies on precise modelling of the impact of neutrinos on structure formation, which can be studied with $N$-body simulations. Here we present the results from a major code comparison effort to establish the maturity and reliability of numerical methods for treating massive neutrinos. The comparison includes eleven full $N$-body implementations (not all of them independent), two $N$-body schemes with approximate time integration, and four additional codes that directly predict or emulate the matter power spectrum. Using a common set of initial data we quantify the relative agreement on the nonlinear power spectrum of cold dark matter and baryons and, for the $N$-body codes, also the relative agreement on the bispectrum, halo mass function, and halo bias. We find that the different numerical implementations produce fully consistent results. We can therefore be confident that we can model the impact of massive neutrinos at the sub-percent level in the most common summary statistics. We also provide a code validation pipeline for future reference.
△ Less
Submitted 8 August, 2023; v1 submitted 22 November, 2022;
originally announced November 2022.
-
Active galactic nucleus jet feedback in hydrostatic halos
Authors:
Rainer Weinberger,
Kung-Yi Su,
Kristian Ehlert,
Christoph Pfrommer,
Lars Hernquist,
Greg L. Bryan,
Volker Springel,
Yuan Li,
Blakesley Burkhart,
Ena Choi,
Claude-André Faucher-Giguère
Abstract:
Feedback driven by jets from active galactic nuclei is believed to be responsible for reducing cooling flows in cool-core galaxy clusters. We use simulations to model feedback from hydrodynamic jets in isolated halos. While the jet propagation converges only after the diameter of the jet is well resolved, reliable predictions about the effects these jets have on the cooling time distribution funct…
▽ More
Feedback driven by jets from active galactic nuclei is believed to be responsible for reducing cooling flows in cool-core galaxy clusters. We use simulations to model feedback from hydrodynamic jets in isolated halos. While the jet propagation converges only after the diameter of the jet is well resolved, reliable predictions about the effects these jets have on the cooling time distribution function only require resolutions sufficient to keep the jet-inflated cavities stable. Comparing different model variations, as well as an independent jet model using a different hydrodynamics code, we show that the dominant uncertainties are the choices of jet properties within a given model. Independent of implementation, we find that light, thermal jets with low momentum flux tend to delay the onset of a cooling flow more efficiently on a $50$ Myr timescale than heavy, kinetic jets. The delay of the cooling flow originates from a displacement and boost in entropy of the central gas. If the jet kinetic luminosity depends on accretion rate, collimated, light, hydrodynamic jets are able to reduce cooling flows in halos, without a need for jet precession or wide opening angles. Comparing the jet feedback with a `kinetic wind' implementation shows that equal amounts of star formation rate reduction can be achieved by different interactions with the halo gas: the jet has a larger effect on the hot halo gas while leaving the denser, star forming phase in place, while the wind acts more locally on the star forming phase, which manifests itself in different time-variability properties.
△ Less
Submitted 5 May, 2023; v1 submitted 21 November, 2022;
originally announced November 2022.