-
The MillenniumTNG Project: The hydrodynamical full physics simulation and a first look at its galaxy clusters
Authors:
Ruediger Pakmor,
Volker Springel,
Jonathan P. Coles,
Thomas Guillet,
Christoph Pfrommer,
Sownak Bose,
Monica Barrera,
Ana Maria Delgado,
Fulvio Ferlito,
Carlos Frenk,
Boryana Hadzhiyska,
César Hernández-Aguayo,
Lars Hernquist,
Rahul Kannan,
Simon D. M. White
Abstract:
Cosmological simulations are an important theoretical pillar for understanding nonlinear structure formation in our Universe and for relating it to observations on large scales. In several papers, we introduce our MillenniumTNG (MTNG) project that provides a comprehensive set of high-resolution, large volume simulations of cosmic structure formation aiming to better understand physical processes o…
▽ More
Cosmological simulations are an important theoretical pillar for understanding nonlinear structure formation in our Universe and for relating it to observations on large scales. In several papers, we introduce our MillenniumTNG (MTNG) project that provides a comprehensive set of high-resolution, large volume simulations of cosmic structure formation aiming to better understand physical processes on large scales and to help interpreting upcoming large-scale galaxy surveys. We here focus on the full physics box MTNG740 that computes a volume of $(740\,\mathrm{Mpc})^3$ with a baryonic mass resolution of $3.1\times~10^7\,\mathrm{M_\odot}$ using \textsc{arepo} with $80.6$~billion cells and the IllustrisTNG galaxy formation model. We verify that the galaxy properties produced by MTNG740 are consistent with the TNG simulations, including more recent observations. We focus on galaxy clusters and analyse cluster scaling relations and radial profiles. We show that both are broadly consistent with various observational constraints. We demonstrate that the SZ-signal on a deep lightcone is consistent with Planck limits. Finally, we compare MTNG740 clusters with galaxy clusters found in Planck and the SDSS-8 RedMaPPer richness catalogue in observational space, finding very good agreement as well. However, {\it simultaneously} matching cluster masses, richness, and Compton-$y$ requires us to assume that the SZ mass estimates for Planck clusters are underestimated by $0.2$~dex on average. Thanks to its unprecedented volume for a high-resolution hydrodynamical calculation, the MTNG740 simulation offers rich possibilities to study baryons in galaxies, galaxy clusters, and in large scale structure, and in particular their impact on upcoming large cosmological surveys.
△ Less
Submitted 6 December, 2022; v1 submitted 18 October, 2022;
originally announced October 2022.
-
The SATIN project I: Turbulent multi-phase ISM in Milky Way simulations with SNe feedback from stellar clusters
Authors:
Rebekka Bieri,
Thorsten Naab,
Sam Geen,
Jonathan P. Coles,
Rüdiger Pakmor,
Stefanie Walch
Abstract:
We introduce the star formation and Supernova (SN) feedback model of the SATIN (Simulating AGNs Through ISM with Non-Equilibrium Effects) project to simulate the evolution of the star forming multi-phase interstellar medium (ISM) of entire disk galaxies. This galaxy-wide implementation of a successful ISM feedback model naturally covers an order of magnitude in gas surface density, shear and radia…
▽ More
We introduce the star formation and Supernova (SN) feedback model of the SATIN (Simulating AGNs Through ISM with Non-Equilibrium Effects) project to simulate the evolution of the star forming multi-phase interstellar medium (ISM) of entire disk galaxies. This galaxy-wide implementation of a successful ISM feedback model naturally covers an order of magnitude in gas surface density, shear and radial motions. It is implemented in the adaptive mesh refinement code RAMSES at a peak resolution of 9 pc. New stars are represented by star cluster (sink) particles with individual SN delay times for massive stars. With SN feedback, cooling and gravity, the galactic ISM develops a realistic three-phase structure. The star formation rates naturally follow observed scaling relations for the local Milky Way gas surface density. SNe drive additional turbulence in the warm (300 K < $T$ < 10$^4$ K) gas and increase the kinetic energy of the cold gas, cooling out of the warm phase. The majority of the gas leaving the galactic ISM is warm and hot with mass loading factors of $3 \le η\le 10$. While the hot gas is leaving the system, the warm and cold gas falls back onto the disc in a galactic fountain flow.
△ Less
Submitted 14 September, 2022;
originally announced September 2022.
-
Bridging the Gap Between Simply Parametrized and Free-Form Pixelated Models of Galaxy Lenses: The Case of WFI 2033-4723 Quad
Authors:
Bernardo Barrera,
Liliya L. R. Williams,
Jonathan P. Coles,
Philipp Denzel
Abstract:
We study the radial and azimuthal mass distribution of the lensing galaxy in WFI2033-4723. Mindful of the fact that modeling results depend on modeling assumptions, we examine two very different recent models: simply parametrized (SP) models from the H0LiCOW collaboration, and pixelated free-form (FF) GLASS models. In addition, we fit our own models which are a compromise between the astrophysical…
▽ More
We study the radial and azimuthal mass distribution of the lensing galaxy in WFI2033-4723. Mindful of the fact that modeling results depend on modeling assumptions, we examine two very different recent models: simply parametrized (SP) models from the H0LiCOW collaboration, and pixelated free-form (FF) GLASS models. In addition, we fit our own models which are a compromise between the astrophysical grounding of SP, and the flexibility of FF approaches. Our models consist of two offset parametric mass components, and generate many solutions, all fitting the quasar point image data. Among other results, we show that to reproduce point image properties the lensing mass must be lopsided, but the origin of this asymmetry can reside in the main lens plane or along the line of sight. We also show that there is a degeneracy between the slope of the density profile and the magnitude of external shear, and that the models from various modeling approaches are connected not by the mass sheet degeneracy, but by a more generalized transformation. Finally, we discuss interpretation degeneracy which afflicts all mass modeling: inability to correctly assign mass to the main lensing galaxy vs. nearby galaxies or line of sight structures. While this may not be a problem for the determination of $H_0$, interpretation degeneracy may become a major issue for the detailed study of galaxy structure.
△ Less
Submitted 27 September, 2021; v1 submitted 9 August, 2021;
originally announced August 2021.
-
The lens SW05 J143454.4+522850: a fossil group at redshift 0.6?
Authors:
Philipp Denzel,
Onur Çatmabacak,
Jonathan P. Coles,
Claude Cornen,
Robert Feldmann,
Ignacio Ferreras,
Xanthe Gwyn Palmer,
Rafael Küng,
Dominik Leier,
Prasenjit Saha,
Aprajita Verma
Abstract:
Fossil groups are considered the end product of natural galaxy group evolution in which group members sink towards the centre of the gravitational potential due to dynamical friction, merging into a single, massive, and X-ray bright elliptical. Since gravitational lensing depends on the mass of a foreground object, its mass concentration, and distance to the observer, we can expect lensing effects…
▽ More
Fossil groups are considered the end product of natural galaxy group evolution in which group members sink towards the centre of the gravitational potential due to dynamical friction, merging into a single, massive, and X-ray bright elliptical. Since gravitational lensing depends on the mass of a foreground object, its mass concentration, and distance to the observer, we can expect lensing effects of such fossil groups to be particularly strong. This paper explores the exceptional system $\mathrm{J}143454.4+522850$. We combine gravitational lensing with stellar population-synthesis to separate the total mass of the lens into stars and dark matter. The enclosed mass profiles are contrasted with state-of-the-art galaxy formation simulations, to conclude that SW05 is likely a fossil group with a high stellar to dark matter mass fraction $0.027\pm0.003$ with respect to expectations from abundance matching $0.012\pm0.004$, indicative of a more efficient conversion of gas into stars in fossil groups.
△ Less
Submitted 7 April, 2021;
originally announced April 2021.
-
The Hubble constant from eight time-delay galaxy lenses
Authors:
Philipp Denzel,
Jonathan P. Coles,
Prasenjit Saha,
Liliya L. R. Williams
Abstract:
We present a determination of the Hubble constant from the joint, free-form analysis of 8 strongly, quadruply lensing systems. In the concordance cosmology, we find $H_0 = 71.8^{+3.9}_{-3.3}\,\mathrm{km}\,\mathrm{s}^{-1}\,\mathrm{Mpc}^{-1}$ with a precision of $4.97\%$. This is in agreement with the latest measurements from Supernovae Type Ia and Planck observations of the cosmic microwave backgro…
▽ More
We present a determination of the Hubble constant from the joint, free-form analysis of 8 strongly, quadruply lensing systems. In the concordance cosmology, we find $H_0 = 71.8^{+3.9}_{-3.3}\,\mathrm{km}\,\mathrm{s}^{-1}\,\mathrm{Mpc}^{-1}$ with a precision of $4.97\%$. This is in agreement with the latest measurements from Supernovae Type Ia and Planck observations of the cosmic microwave background. Our precision is lower compared to these and other recent time-delay cosmography determinations, because our modelling strategies reflect the systematic uncertainties of lensing degeneracies. We furthermore are able to find reasonable lensed image reconstructions by constraining to either value of $H_0$ from local and early Universe measurements. This leads us to conclude that current lensing constraints on $H_0$ are not strong enough to break the "Hubble tension" problem of cosmology.
△ Less
Submitted 19 February, 2021; v1 submitted 28 July, 2020;
originally announced July 2020.
-
Lessons from a blind study of simulated lenses: image reconstructions do not always reproduce true convergence
Authors:
Philipp Denzel,
Sampath Mukherjee,
Jonathan P. Coles,
Prasenjit Saha
Abstract:
In the coming years, strong gravitational lens discoveries are expected to increase in frequency by two orders of magnitude. Lens-modelling techniques are being developed to prepare for the coming massive influx of new lens data, and blind tests of lens reconstruction with simulated data are needed for validation. In this paper we present a systematic blind study of a sample of 15 simulated strong…
▽ More
In the coming years, strong gravitational lens discoveries are expected to increase in frequency by two orders of magnitude. Lens-modelling techniques are being developed to prepare for the coming massive influx of new lens data, and blind tests of lens reconstruction with simulated data are needed for validation. In this paper we present a systematic blind study of a sample of 15 simulated strong gravitational lenses from the EAGLE suite of hydrodynamic simulations. We model these lenses with a free-form technique and evaluate reconstructed mass distributions using criteria based on shape, orientation, and lensed image reconstruction. Especially useful is a lensing analogue of the Roche potential in binary star systems, which we call the $\textit{lensing Roche potential}$. This we introduce in order to factor out the well-known problem of steepness or mass-sheet degeneracy. Einstein radii are on average well recovered with a relative error of ${\sim}5\%$ for quads and ${\sim}25\%$ for doubles; the position angle of ellipticity is on average also reproduced well up to $\pm10^{\circ}$, but the reconstructed mass maps tend to be too round and too shallow. It is also easy to reproduce the lensed images, but optimising on this criterion does not guarantee better reconstruction of the mass distribution.
△ Less
Submitted 17 June, 2020; v1 submitted 22 October, 2019;
originally announced October 2019.
-
An Optimizing Symbolic Algebra Approach for Generating Fast Multipole Method Operators
Authors:
Jonathan P. Coles,
Rebekka Bieri
Abstract:
We have developed a symbolic algebra approach to automatically produce, verify, and optimize computer code for the Fast Multipole Method (FMM) operators. This approach allows for flexibility in choosing a basis set and kernel, and can generate computer code for any expansion order in multiple languages. The procedure is implemented in the publicly available Python program Mosaic. Optimizations per…
▽ More
We have developed a symbolic algebra approach to automatically produce, verify, and optimize computer code for the Fast Multipole Method (FMM) operators. This approach allows for flexibility in choosing a basis set and kernel, and can generate computer code for any expansion order in multiple languages. The procedure is implemented in the publicly available Python program Mosaic. Optimizations performed at the symbolic level through algebraic manipulations significantly reduce the number of mathematical operations compared with a straightforward implementation of the equations. We find that the optimizer is able to eliminate 20-80% of the floating-point operations and for the expansion orders $p \le 10$ it changes the observed scaling properties. We present our approach using three variants of the operators with the Cartesian basis set for the harmonic potential kernel $1/r$, including the use of totally symmetric and traceless multipole tensors.
△ Less
Submitted 28 May, 2020; v1 submitted 15 November, 2018;
originally announced November 2018.
-
Light versus dark in strong-lens galaxies: Dark matter haloes that are rounder than their stars
Authors:
Claudio Bruderer,
Justin I. Read,
Jonathan P. Coles,
Dominik Leier,
Emilio E. Falco,
Ignacio Ferreras,
Prasenjit Saha
Abstract:
We measure the projected density profile, shape and alignment of the stellar and dark matter mass distribution in 11 strong-lens galaxies. We find that the projected dark matter density profile - under the assumption of a Chabrier stellar initial mass function - shows significant variation from galaxy to galaxy. Those with an outermost image beyond $\sim 10$ kpc are very well fit by a projected NF…
▽ More
We measure the projected density profile, shape and alignment of the stellar and dark matter mass distribution in 11 strong-lens galaxies. We find that the projected dark matter density profile - under the assumption of a Chabrier stellar initial mass function - shows significant variation from galaxy to galaxy. Those with an outermost image beyond $\sim 10$ kpc are very well fit by a projected NFW profile; those with images within 10 kpc appear to be more concentrated than NFW, as expected if their dark haloes contract due to baryonic cooling. We find that over several half-light radii, the dark matter haloes of these lenses are rounder than their stellar mass distributions. While the haloes are never more elliptical than $e_{dm} = 0.2$, their stars can extend to $e_* > 0.2$. Galaxies with high dark matter ellipticity and weak external shear show strong alignment between light and dark; those with strong shear ($γ\gtrsim 0.1$) can be highly misaligned. This is reassuring since isolated misaligned galaxies are expected to be unstable. Our results provide a new constraint on galaxy formation models. For a given cosmology, these must explain the origin of both very round dark matter haloes and misaligned strong-lens systems.
△ Less
Submitted 3 November, 2015;
originally announced November 2015.
-
Gravitational Lens Recovery with GLASS: Measuring the mass profile and shape of a lens
Authors:
Jonathan P. Coles,
Justin I. Read,
Prasenjit Saha
Abstract:
We use a new non-parametric gravitational modelling tool -- \Glass{} -- to determine what quality of data (strong lensing, stellar kinematics, and/or stellar masses) are required to measure the circularly averaged mass profile of a lens and its shape. \Glass{} uses an under-constrained adaptive grid of mass pixels to model the lens, searching through thousands of models to marginalise over model u…
▽ More
We use a new non-parametric gravitational modelling tool -- \Glass{} -- to determine what quality of data (strong lensing, stellar kinematics, and/or stellar masses) are required to measure the circularly averaged mass profile of a lens and its shape. \Glass{} uses an under-constrained adaptive grid of mass pixels to model the lens, searching through thousands of models to marginalise over model uncertainties. Our key findings are as follows: (i) for pure lens data, multiple sources with wide redshift separation give the strongest constraints as this breaks the well-known mass-sheet or steepness degeneracy; (ii) a single quad with time delays also performs well, giving a good recovery of both the mass profile and its shape; (iii) stellar masses -- for lenses where the stars dominate the central potential -- can also break the steepness degeneracy, giving a recovery for doubles almost as good as having a quad with time delay data, or multiple source redshifts; (iv) stellar kinematics provide a robust measure of the mass at the half light radius of the stars $r_{1/2}$ that can also break the steepness degeneracy if the Einstein radius $r_E \neq r_{1/2}$; and (v) if $r_E \sim r_{1/2}$, then stellar kinematic data can be used to probe the stellar velocity anisotropy $β$ -- an interesting quantity in its own right. Where information on the mass distribution from lensing and/or other probes becomes redundant, this opens up the possibility of using strong lensing to constrain cosmological models.
△ Less
Submitted 29 August, 2014; v1 submitted 30 January, 2014;
originally announced January 2014.