Researcher profile

Klaus Dolag

Klaus Dolag contributes to research discovery and scholarly infrastructure.

ResearcherAffiliation not importedOpen to collaborate

Trust snapshot

Quick read

Trust 21 - EmergingVerification L1Unclaimed author
27works
0followers
12topics
4close collaborators

Actions

Decide how to stay connected

Follow researcher0

Identity and collaboration

How to connect with this researcher

Claiming links this public author record to a researcher profile and unlocks direct collaboration workflows.

Log in to claim

Direct collaboration

Open a focused conversation when the fit is right

Claim this author entity first to unlock direct invitations.

Research graph

See the researcher in context

Open full explorer

Inspect adjacent work, topics, institutions and collaborators without jumping out to a separate graph page.

Building this graph slice

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

27 published item(s)

preprint2025arXiv

Cosmological Inference with Cosmic Voids and Neural Network Emulators

Cosmic Voids are a promising probe of cosmology for spectroscopic galaxy surveys due to their unique response to cosmological parameters. Their combination with other probes promises to break parameter degeneracies. Due to simplifying assumptions, analytical models for void statistics are only representative of a subset of the full void population. We present a set of neural-based emulators for void summary statistics of watershed voids, which retain more information about the full void population than simplified analytical models. We build emulators for the void size function and void density profiles traced by the halo number density using the Quijote suite of simulations for a broad range of the $Λ\mathrm{CDM}$ parameter space. The emulators replace the computation of these statistics from computationally expensive cosmological simulations. We demonstrate the cosmological constraining power of voids using our emulators, which offer orders-of-magnitude acceleration in parameter estimation, capture more cosmological information compared to analytic models, and produce more realistic posteriors compared to Fisher forecasts. We find that the parameters $Ω_m$ and $σ_8$ in this Quijote setup can be recovered to $14.4\%$ and $8.4\%$ accuracy respectively using void density profiles; including the additional information in the void size function improves the accuracy on $σ_8$ to $6.8\%$. We demonstrate the robustness of our approach to two important variables in the underlying simulations, the resolution, and the inclusion of baryons. We find that our pipeline is robust to variations in resolution, and we show that the posteriors derived from the emulated void statistics are unaffected by the inclusion of baryons with the Magneticum hydrodynamic simulations. This opens up the possibility of a baryon-independent probe of the large-scale structure.

preprint2023arXiv

KiDS-1000 Cosmology: Constraints from density split statistics

Context. Weak lensing and clustering statistics beyond two-point functions can capture non-Gaussian information about the matter density field, thereby improving the constraints on cosmological parameters relative to the mainstream methods based on correlation functions and power spectra. Aims. This paper presents a cosmological analysis of the fourth data release of the Kilo Degree Survey based on the density split statistics, which measures the mean shear profiles around regions classified according to foreground densities. The latter is constructed from a bright galaxy sample, which we further split into red and blue samples, allowing us to probe their respective connection to the underlying dark matter density. Methods. We use the state-of-the-art model of the density splitting statistics and validate its robustness against mock data infused with known systematic effects such as intrinsic galaxy alignment and baryonic feedback. Results. After marginalising over the photometric redshift uncertainty and the residual shear calibration bias, we measure for the full KiDS-bright sample a structure growth parameter of $S_8 = σ_8 \sqrt{Ω_\mathrm{m}/0.3} = 0.74^{+0.03}_{-0.02}$ that is competitive to and consistent with two-point cosmic shear results, a matter density of $Ω_\mathrm{m} = 0.28 \pm 0.02$, and a constant galaxy bias of $b = 1.32^{+0.12}_{-0.10}$.

preprint2022arXiv

\textsc{The Three Hundred} project: The \textsc{Gizmo-Simba} run

We introduce \textsc{Gizmo-Simba}, a new suite of galaxy cluster simulations within \textsc{The Three Hundred} project. \textsc{The Three Hundred} consists of zoom re-simulations of 324 clusters with $M_{200}\gtrsim 10^{14.8}M_\odot$ drawn from the MultiDark-Planck $N$-body simulation, run using several hydrodynamic and semi-analytic codes. The \textsc{Gizmo-Simba} suite adds a state-of-the-art galaxy formation model based on the highly successful {\sc Simba} simulation, mildly re-calibrated to match $z=0$ cluster stellar properties. Comparing to \textsc{The Three Hundred} zooms run with \textsc{Gadget-X}, we find intrinsic differences in the evolution of the stellar and gas mass fractions, BCG ages, and galaxy colour-magnitude diagrams, with \textsc{Gizmo-Simba} generally providing a good match to available data at $z \approx 0$. \textsc{Gizmo-Simba}'s unique black hole growth and feedback model yields agreement with the observed BH scaling relations at the intermediate-mass range and predicts a slightly different slope at high masses where few observations currently lie. \textsc{Gizmo-Simba} provides a new and novel platform to elucidate the co-evolution of galaxies, gas, and black holes within the densest cosmic environments.

preprint2022arXiv

Cosmological simulations with rare and frequent dark matter self-interactions

Dark matter (DM) with self-interactions is a promising solution for the small-scale problems of the standard cosmological model. Here we perform the first cosmological simulation of frequent DM self-interactions, corresponding to small-angle DM scatterings. The focus of our analysis lies in finding and understanding differences to the traditionally assumed rare DM (large-angle) self scatterings. For this purpose, we compute the distribution of DM densities, the matter power spectrum, the two-point correlation function and the halo and subhalo mass functions. Furthermore, we investigate the density profiles of the DM haloes and their shapes. We find that overall large-angle and small-angle scatterings behave fairly similarly with a few exceptions. In particular, the number of satellites is considerably suppressed for frequent compared to rare self-interactions with the same cross-section. Overall we observe that while differences between the two cases may be difficult to establish using a single measure, the degeneracy may be broken through a combination of multiple ones. For instance, the combination of satellite counts with halo density or shape profiles could allow discriminating between rare and frequent self-interactions. As a by-product of our analysis, we provide - for the first time - upper limits on the cross-section for frequent self-interactions.

preprint2022arXiv

Decomposition of galactic X-ray emission with Phox: Contributions from hot gas and XRBs

We provide a numerical framework with which spatially and spectrally accurate representations of X-ray binary populations can be studied from hydrodynamical cosmological simulations. We construct average spectra accounting for a hot gas component and verify the emergence of observed scaling relations between galaxy wide X-ray luminosity ($L_{X}$) and stellar mass ($M_{\star}$) as well as star-formation rate (SFR). Using simulated galaxy halos extracted from the $(48\,h^{-1} \mathrm{cMpc})^3$ volume of the Magneticum Pathfinder cosmological simulations at $z = 0.07$ we generate mock spectra with the X-ray photon-simulator Phox. We extend the Phox code to account for the stellar component in the simulation and study the resulting contribution in composite galactic spectra. Average X-ray luminosity functions are perfectly reproduced up to the one-photon luminosity limit. Comparing our resulting $L_{X}-\mathrm{SFR}-M_{\star}$ relation for X-ray binaries with recent observations of field galaxies in the Virgo galaxy cluster we find significant overlap. Invoking a metallicity dependent model for high-mass X-ray binaries yields an anti-correlation between mass-weighted stellar metallicity and SFR normalized luminosity. The spatial distribution of high-mass X-ray binaries coincides with star-formation regions of simulated galaxies while low-mass X-ray binaries follow the stellar mass surface density. X-ray binary emission is the dominant contribution in the 2-10 keV band in the absence of an actively accreting central super-massive black hole with 50% contribution in the 0.5-2 keV band rivaling the hot gas component. Our modelling remains consistent with observations despite uncertainties connected to our approach. The predictive power and easily extendable framework hold great value for future investigations of galactic X-ray spectra.

preprint2022arXiv

Galaxies in the central regions of simulated galaxy clusters

In this paper, we assess the impact of numerical resolution and of the implementation of energy input from AGN feedback models on the inner structure of cluster sub-haloes in hydrodynamic simulations. We compare several zoom-in re-simulations of a sub-sample of the cluster-sized haloes studied in Meneghetti et al. (2020), obtained by varying mass resolution, softening length and AGN energy feedback scheme. We study the impact of these different setups on the subhalo abundances, their radial distribution, their density and mass profiles and the relation between the maximum circular velocity, which is a proxy for subhalo compactness. Regardless of the adopted numerical resolution and feedback model, subhaloes with masses Msub < 1e11Msun/h, the most relevant mass-range for galaxy-galaxy strong lensing, have maximum circular velocities ~30% smaller than those measured from strong lensing observations of Bergamini et al. (2019). We also find that simulations with less effective AGN energy feedback produce massive subhaloes (Msub> 1e11 Msun/h ) with higher maximum circular velocity and that their Vmax - Msub relation approaches the observed one. However the stellar-mass number count of these objects exceeds the one found in observations and we find that the compactness of these simulated subhaloes is the result of an extremely over-efficient star formation in their cores, also leading to larger-than-observed subhalo stellar mass. We conclude that simulations are unable to simultaneously reproduce the observed stellar masses and compactness (or maximum circular velocities) of cluster galaxies. Thus, the discrepancy between theory and observations that emerged from the analysis of Meneghetti et al. (2020) persists. It remains an open question as to whether such a discrepancy reflects limitations of the current implementation of galaxy formation models or the LCDM paradigm.

preprint2022arXiv

Planes of Satellite Galaxies in the Magneticum Pathfinder Simulations

Planes of satellites are observed around many galaxies. However, these observations are still considered a point of tension for the $Λ$CDM paradigm. We use the fully hydrodynamical cosmological $Λ$CDM state-of-the-art simulation Magneticum Pathfinder to investigate the existence of such planes over a large range of haloes, from Milky Way to galaxy cluster masses. To this end, we develop the Momentum in Thinnest Plane (MTP) method to identify planes and quantify the properties of their constituent satellites, considering both position and momentum. We find that thin planes ($20\%\cdot R_\mathrm{halo}$) containing at least $50\%$ of the total number of satellites can be found in almost all systems. In Milky Way mass-like systems, around 86\% of such planes are even aligned in momentum ($90\%$ of the total satellite momentum), where the fraction is smaller if more satellites are required to be inside the plane. We further find a mass dependency, with more massive systems exhibiting systematically thicker planes. This may point towards the change from continuous accretion of small objects along filaments and sheets for less massive haloes to the accretion of large objects (e.g., major mergers) dominating the growth of more massive haloes. There is no correlation between the existence of a plane and the main galaxy&#39;s morphology. Finally, we find a clear preference for the minor axes of the satellite planes and the host galaxy to be aligned, in agreement with recent observations.

preprint2022arXiv

PSI: Constructing ad-hoc Simplices to Interpolate High-Dimensional Unstructured Data

Interpolating unstructured data using barycentric coordinates becomes infeasible at high dimensions due to the prohibitive memory requirements of building a Delaunay triangulation. We present a new algorithm to construct ad-hoc simplices that are empirically guaranteed to contain the target coordinates, based on a nearest neighbor heuristic and an iterative dimensionality reduction through projection. We use these simplices to interpolate the astrophysical cooling function $Λ$ and show that this new approach produces good results with just a fraction of the previously required memory.

preprint2022arXiv

Robust field-level inference with dark matter halos

We train graph neural networks on halo catalogues from Gadget N-body simulations to perform field-level likelihood-free inference of cosmological parameters. The catalogues contain $\lesssim$5,000 halos with masses $\gtrsim 10^{10}~h^{-1}M_\odot$ in a periodic volume of $(25~h^{-1}{\rm Mpc})^3$; every halo in the catalogue is characterized by several properties such as position, mass, velocity, concentration, and maximum circular velocity. Our models, built to be permutationally, translationally, and rotationally invariant, do not impose a minimum scale on which to extract information and are able to infer the values of $Ω_{\rm m}$ and $σ_8$ with a mean relative error of $\sim6\%$, when using positions plus velocities and positions plus masses, respectively. More importantly, we find that our models are very robust: they can infer the value of $Ω_{\rm m}$ and $σ_8$ when tested using halo catalogues from thousands of N-body simulations run with five different N-body codes: Abacus, CUBEP$^3$M, Enzo, PKDGrav3, and Ramses. Surprisingly, the model trained to infer $Ω_{\rm m}$ also works when tested on thousands of state-of-the-art CAMELS hydrodynamic simulations run with four different codes and subgrid physics implementations. Using halo properties such as concentration and maximum circular velocity allow our models to extract more information, at the expense of breaking the robustness of the models. This may happen because the different N-body codes are not converged on the relevant scales corresponding to these parameters.

preprint2022arXiv

The eROSITA view of the Abell 3391/95 field: The Northern Clump. The largest infalling structure in the longest known gas filament observed with eROSITA, XMM-Newton, and Chandra

SRG/eROSITA PV observations revealed the A3391/95 cluster system and the Northern Clump (MCXC J0621.7-5242 galaxy cluster) are aligning along a cosmic filament in soft X-rays, similarly to what has been seen in simulations before. We aim to understand the dynamical state of the Northern Clump as it enters the atmosphere ($3\times R_{200}$) of A3391. We analyzed joint eROSITA, XMM-Newton, and Chandra observations to probe the morphological, thermal, and chemical properties of the Northern Clump from its center out to a radius of 988 kpc ($R_{200}$). We utilized the ASKAP/EMU radio data, DECam optical image, and Planck y-map to study the influence of the WAT radio source on the Northern Clump central ICM. From the Magneticum simulation, we identified an analog of the A3391/95 system along with an infalling group resembling the Northern Clump. The Northern Clump is a WCC cluster centered on a WAT radio galaxy. The gas temperature over $0.2-0.5R_{500}$ is $k_BT_{500}=1.99\pm0.04$ keV. We employed the $M-T$ scaling relation and obtained a mass estimate of $M_{500}=(7.68\pm0.43)\times10^{13}M_{\odot}$ and $R_{500}=(636\pm12)$ kpc. Its atmosphere has a boxy shape and deviates from spherical symmetry. We identify a southern surface brightness edge, likely caused by subsonic motion relative to the filament gas. At $\sim\! R_{500}$, the southern atmosphere appears to be 42% hotter than its northern atmosphere. We detect a downstream tail pointing toward the north with a projected length of $\sim318$ kpc, plausibly the result of ram pressure stripping. The analog group in the Magneticum simulation is experiencing changes in its gas properties and a shift between the position of the halo center and that of the bound gas while approaching the main cluster pair.

preprint2022arXiv

Unequal-mass mergers of dark matter haloes with rare and frequent self-interactions

Dark matter (DM) self-interactions have been proposed to solve problems on small length scales within the standard cold DM cosmology. Here, we investigate the effects of DM self-interactions in merging systems of galaxies and galaxy clusters with equal and unequal mass ratios. We perform N-body DM-only simulations of idealized setups to study the effects of DM self-interactions that are elastic and velocity-independent. We go beyond the commonly adopted assumption of large-angle (rare) DM scatterings, paying attention to the impact of small-angle (frequent) scatterings on astrophysical observables and related quantities. Specifically, we focus on DM-galaxy offsets, galaxy--galaxy distances, halo shapes, morphology, and the phase--space distribution. Moreover, we compare two methods to identify peaks: one based on the gravitational potential and one based on isodensity contours. We find that the results are sensitive to the peak finding method, which poses a challenge for the analysis of merging systems in simulations and observations, especially for minor mergers. Large DM-galaxy offsets can occur in minor mergers, especially with frequent self-interactions. The subhalo tends to dissolve quickly for these cases. While clusters in late merger phases lead to potentially large differences between rare and frequent scatterings, we believe that these differences are non-trivial to extract from observations. We therefore study the galaxy/star populations which remain distinct even after the DM haloes have coalesced. We find that these collisionless tracers behave differently for rare and frequent scatterings, potentially giving a handle to learn about the micro-physics of DM.

preprint2021arXiv

A New View of Observed Galaxies through 3D Modelling and Visualisation

Observational astronomers survey the sky in great detail to gain a better understanding of many types of astronomical phenomena. In particular, the formation and evolution of galaxies, including our own, is a wide field of research. Three dimensional (spatial 3D) scientific visualisation is typically limited to simulated galaxies, due to the inherently two dimensional spatial resolution of Earth-based observations. However, with appropriate means of reconstruction, such visualisation can also be used to bring out the inherent 3D structure that exists in 2D observations of known galaxies, providing new views of these galaxies and visually illustrating the spatial relationships within galaxy groups that are not obvious in 2D. We present a novel approach to reconstruct and visualise 3D representations of nearby galaxies based on observational data using the scientific visualisation software Splotch. We apply our approach to a case study of the nearby barred spiral galaxy known as M83, presenting a new perspective of the M83 local group and highlighting the similarities between our reconstructed views of M83 and other known galaxies of similar inclinations.

preprint2021arXiv

Effects of boosting on extragalactic components: methods and statistical studies

In this work we examine the impact of our motion with respect to the CMB rest frame on statistics of CMB maps by examining the one-, two-, three- and four- point statistics of simulated maps of the CMB and Sunyaev-Zeldovich (SZ) effects. We validate boosting codes by comparing their outcomes for temperature and polarization power spectra up to $\ell \simeq 6000$. We derive and validate a new analytical formula for the computation of the boosted power spectrum of a signal with a generic frequency dependence. As an example we show how this increases the boosting correction to the power spectrum of CMB intensity measurements by $\sim 30\%$ at 150 GHz. We examine the effect of boosting on thermal and kinetic SZ power spectra from semianalytical and hydrodynamical simulations; the boosting correction is generally small for both simulations, except when considering frequencies near the tSZ null. For the non-Gaussian statistics, in general we find that boosting has no impact with two exceptions. We find that, whilst the statistics of the CMB convergence field are unaffected, quadratic estimators that are used to measure this field can become biased at the $O(1)\%$ level by boosting effects. We present a simple modification to the standard estimators that removes this bias. Second, bispectrum estimators can receive a systematic bias from the Doppler induced quadrupole when there is anisotropy in the sky -- in practice this anisotropy comes from masking and inhomogenous noise. This effect is unobservable and already removed by existing analysis methods.

preprint2021arXiv

Galaxy Velocity Bias in Cosmological Simulations: Towards Percent-level Calibration

Galaxy cluster masses, rich with cosmological information, can be estimated from internal dark matter (DM) velocity dispersions, which in turn can be observationally inferred from satellite galaxy velocities. However, galaxies are biased tracers of the DM, and the bias can vary over host halo and galaxy properties as well as time. We precisely calibrate the velocity bias, b_v -- defined as the ratio of galaxy and DM velocity dispersions -- as a function of redshift, host halo mass, and galaxy stellar mass threshold (Mstarsat), for massive halos (M200c > 1e13.5 msun) from five cosmological simulations: IllustrisTNG, Magneticum, Bahamas + Macsis, The Three Hundred Project, and MultiDark Planck-2. We first compare scaling relations for galaxy and DM velocity dispersion across simulations; the former is estimated using a new ensemble velocity likelihood method that is unbiased for low galaxy counts per halo, while the latter uses a local linear regression. The simulations show consistent trends of b_v increasing with M200c and decreasing with redshift and Mstarsat. The ensemble-estimated theoretical uncertainty in b_v is 2-3% but becomes percent-level when considering only the three highest resolution simulations. We update the mass-richness normalization previously estimated by Farahi et al. (2016) for an SDSS redMaPPer cluster sample. The improved accuracy of our b_v estimates reduces the mass normalization uncertainty from 22% to 8%, demonstrating that dynamical estimation techniques can be competitive with weak lensing in calibrating population mean masses. We discuss necessary steps for further improving this precision. Our estimates for b_v(M200c, Mstarsat, z) are made publicly available.

preprint2021arXiv

On the small scale turbulent dynamo in the intracluster medium: A comparison to dynamo theory

We present non-radiative, cosmological zoom-simulations of galaxy cluster formation with magnetic fields and (anisotropic) thermal conduction of one very massive galaxy cluster with a mass at redshift zero that corresponds to $M_\mathrm{vir} \sim 2 \times 10^{15} M_{\odot}$. We run the cluster on three resolution levels (1X, 10X, 25X), starting with an effective mass resolution of $2 \times 10^8M_{\odot}$, subsequently increasing the particle number to reach $4 \times 10^6M_{\odot}$. The maximum spatial resolution obtained in the simulations is limited by the gravitational softening reaching $ε=1.0$ kpc at the highest resolution level, allowing to resolve the hierarchical assembly of the structures in very fine detail. All simulations presented, have been carried out with the SPMHD-code Gadget-3 with a heavily updated SPMHD prescription. The primary focus is to investigate magnetic field amplification in the Intracluster Medium (ICM). We show that the main amplification mechanism is the small scale-turbulent-dynamo in the limit of reconnection diffusion. In our two highest resolution models we start to resolve the magnetic field amplification driven by this process and we explicitly quantify this with the magnetic power-spectra and the magnetic tension that limits the bending of the magnetic field lines consistent with dynamo theory. Furthermore, we investigate the $\nabla \cdot \mathbf{B}=0$ constraint within our simulations and show that we achieve comparable results to state-of-the-art AMR or moving-mesh techniques, used in codes such as Enzo and Arepo. Our results show for the first time in a fully cosmological simulation of a galaxy cluster that dynamo action can be resolved in the framework of a modern Lagrangian magnetohydrodynamic (MHD) method, a study that is currently missing in the literature.

preprint2021arXiv

The eROSITA view of the Abell 3391/95 field: a case study from the Magneticum cosmological simulation

We investigated the origin and gas properties of a simulated galaxy cluster pair, to connect simulation predictions to the SRG/eROSITA X-ray observations of the Abell 3391/95 field. The simulated system has been extracted from the (352 cMpc/h)^3 volume of the Magneticum Pathfinder cosmological simulations at z=0.07. We tracked back in time the main progenitors of the pair clusters and surrounding groups to study the assembly history of the system and its evolution. Similarly to the observed A3391/95 system, the simulated pair is embedded in a complex network of gas filaments, with structures aligned over more than 20 projected Mpc and the whole region collapsing towards the central overdense node. The spheres of influence (3*R200) of the two main clusters already overlap at z=0.07, but their virial boundaries are still physically separated. The diffuse gas located in the interconnecting bridge closely reflects the WHIM, with typical temperature of ~1 keV and overdensity $\sim 100$, with respect to the mean baryon density of the Universe, and lower enrichment level compared to the ICM in clusters. We find that most of the bridge gas collapsed from directions roughly orthogonal to the intra-cluster gas accretion directions, and its origin is mostly unrelated to the two cluster progenitors. We find clear signatures in the surrounding groups of infall motion towards the pair, such as significant radial velocities and slowdown of gas compared to dark matter. These findings further support the picture of the Northern Clump (MCXC J0621.7-5242) cluster infalling along a cosmic gas filament towards Abell 3391, possibly merging with it. We conclude that, in such a configuration, the pair clusters of the A3391/95-like system are in a pre-merger phase, and did not interact yet. The diffuse gas in the interconnecting bridge is mostly warm filament gas, rather than tidally-stripped cluster gas.

preprint2021arXiv

The Three Hundred Project: The gas disruption of infalling objects in cluster environments

We analyse the gas content evolution of infalling haloes in cluster environments from THE THREE HUNDRED project, a collection of 324 numerically modelled galaxy clusters. The haloes in our sample were selected within $5R_{200}$ of the main cluster halo at $z=0$ and have total halo mass $M_{200}\geq10^{11} h^{-1} M_{\odot}$. We track their main progenitors and study their gas evolution since their crossing into the infall region, which we define as $1-4R_{200}$. Studying the radial trends of our populations using both the full phase space information and a line-of-sight projection, we confirm the Arthur et al. (2019) result and identify a characteristic radius around $1.7R_{200}$ in 3D and at $R_{200}$ in projection at which infalling haloes lose nearly all of the gas prior their infall. Splitting the trends by subhalo status we show that subhaloes residing in group-mass and low-mass host haloes in the infall region follow similar radial gas-loss trends as their hosts, whereas subhaloes of cluster-mass host haloes are stripped of their gas much further out. Our results show that infalling objects suffer significant gaseous disruption that correlates with time-since-infall, cluster-centric distance and host mass, and that the gaseous disruption they experience is a combination of subhalo pre-processing and object gas depletion at a radius which behaves like an accretion shock.

preprint2020arXiv

Characterizing the Uncertainty in Cluster Magnetic Fields derived from Rotation Measures

Magnetic fields play vital roles in intracluster media (ICMs), but estimating their strengths and distributions from observations is a major challenge. Faraday rotation measures (RMs) are widely applied to this task, so it is critical to understand inherent uncertainties in RM analysis. In this paper, we seek to characterize those uncertainties given the types of information available today, independent of the specific technique used. We conduct synthetic RM observations through the ICM of a galaxy cluster drawn from an MHD cosmological simulation in which the magnetic field is known. We analyze the synthetic RM observations using an analytical formalism based on commonly used model assumptions allowing us to relate model physical variables to outcome uncertainties. Despite the simplicity of some assumptions, and unknown physical parameters, we are able to extract an approximate magnitude of the central magnetic field within an apparently irreducible uncertain factor approximately 3. Principal, largely irreducible, uncertainties come from the unknown depth along the line of sight of embedded polarized sources, the lack of robust coherence lengths from area-constrained polarization sampling, and the unknown scaling between ICM electron density and magnetic field strength. The RM-estimated central magnetic field strengths span more than an order of magnitude including the full range of synthetic experiments.

preprint2020arXiv

Collision of Merger and Accretion Shocks: Formation of Mpc-scale Contact Discontinuity in the Perseus Cluster

Two Mpc-size contact discontinuities have recently been identified in the XMM-Newton and Suzaku X-ray observations of the outskirts of the Perseus cluster (Walker et al. 2020). These structures have been tentatively interpreted as &#34;sloshing cold fronts&#34;, which are customarily associated with differential motions of the cluster gas, perturbed by a merger. In this study we consider an alternative scenario, namely, that the most prominent discontinuity near the cluster virial radius is the result of the collision between the accretion shock and the &#34;runaway&#34; merger shock. We also discuss the possible origin of the second discontinuity at ~1.2 Mpc.

preprint2020arXiv

Cosmology dependence of galaxy cluster scaling relations

The abundance of galaxy clusters as a function of mass and redshift is a well known powerful cosmological probe, which relies on underlying modelling assumptions on the mass-observable relations (MOR). Some of the MOR parameters can be constrained directly from multi-wavelength observations, as the normalization at some reference cosmology, the mass-slope, the redshift evolution and the intrinsic scatter. However, the cosmology dependence of MORs cannot be tested with multi-wavelength observations alone. We use {\tt Magneticum} simulations to explore the cosmology dependence of galaxy cluster scaling relations. We run fifteen hydro-dynamical cosmological simulations varying $Ω_m$, $Ω_b$, $h_0$ and $σ_8$ (around a reference cosmological model). The MORs considered are gas mass, baryonic mass, gas temperature, $Y$ and velocity dispersion as a function of virial mass. We verify that the mass and redshift slopes and the intrinsic scatter of the MORs are nearly independent of cosmology with variations significantly smaller than current observational uncertainties. We show that the gas mass and baryonic mass sensitively depends only on the baryon fraction, velocity dispersion and gas temperature on $h_0$, and $Y$ on both baryon fraction and $h_0$. We investigate the cosmological implications of our MOR parameterization on a mock catalog created for an idealized eROSITA-like experiment. We show that our parametrization introduces a strong degeneracy between the cosmological parameters and the normalization of the MOR. Finally, the parameter constraints derived at different overdensity ($Δ_{500c}$), for X-ray bolometric gas luminosity, and for different subgrid physics prescriptions are shown in the appendix.

preprint2020arXiv

Encounters of Merger and Accretion Shocks in Galaxy Clusters and their Effects on Intracluster Medium

Several types/classes of shocks naturally arise during formation and evolution of galaxy clusters. One such class is represented by accretion shocks, associated with deceleration of infalling baryons. Such shocks, characterized by a very high Mach number, are present even in 1D models of cluster evolution. Another class is composed of &#34;runaway merger shocks&#34;, which appear when a merger shock, driven by a sufficiently massive infalling subcluster, propagates away from the main-cluster center. We argue that, when the merger shock overtakes the accretion shock, a new long-living shock is formed that propagates to large distances from the main cluster (well beyond its virial radius) affecting the cold gas around the cluster. We refer to these structures as Merger-accelerated Accretion shocks (MA-shocks) in this paper. We show examples of such MA-shocks in 1D and 3D simulations and discuss their characteristic properties. In particular, (1) MA-shocks shape the boundary separating the hot intracluster medium (ICM) from the unshocked gas, giving this boundary a &#34;flower-like&#34; morphology. In 3D, MA-shocks occupy space between the dense accreting filaments. (2) Evolution of MA-shocks highly depends on the Mach number of the runaway merger shock and the mass accretion rate parameter of the cluster. (3) MA-shocks may lead to the misalignment of the ICM boundary and the splashback radius.

preprint2020arXiv

First Simulations of Axion Minicluster Halos

We study the gravitational collapse of axion dark matter fluctuations in the post-inflationary scenario, so-called axion miniclusters, with N-body simulations. Largely confirming theoretical expectations, overdensities begin to collapse in the radiation-dominated epoch and form an early distribution of miniclusters with masses up to $10^{-12}\,M_\odot$. After matter-radiation equality, ongoing mergers give rise to a steep power-law distribution of minicluster halo masses. The density profiles of well-resolved halos are NFW-like to good approximation. The fraction of axion dark matter in these bound structures is $\sim 0.75$ at redshift $z=100$.

preprint2020arXiv

Gadget3 on GPUs with OpenACC

We present preliminary results of a GPU porting of all main Gadget3 modules (gravity computation, SPH density computation, SPH hydrodynamic force, and thermal conduction) using OpenACC directives. Here we assign one GPU to each MPI rank and exploit both the host and accellerator capabilities by overlapping computations on the CPUs and GPUs: while GPUs asynchronously compute interactions between particles within their MPI ranks, CPUs perform tree-walks and MPI communications of neighbouring particles. We profile various portions of the code to understand the origin of our speedup, where we find that a peak speedup is not achieved because of time-steps with few active particles. We run a hydrodynamic cosmological simulation from the Magneticum project, with $2\cdot10^{7}$ particles, where we find a final total speedup of $\approx 2.$ We also present the results of an encouraging scaling test of a preliminary gravity-only OpenACC porting, run in the context of the EuroHack17 event, where the prototype of the porting proved to keep a constant speedup up to $1024$ GPUs.

preprint2020arXiv

Kinematics of Simulated Galaxies II: Probing the Stellar Kinematics of Galaxies out to Large Radii

We investigate the stellar kinematics of a sample of galaxies extracted from the hydrodynamic cosmological Magneticum Pathfinder simulations out to $5$ half-mass radii. We construct differential radial stellar spin profiles quantified by the observationally widely used $λ_\mathrm{R}$ and the closely related $(V/σ)$ parameters. We find three characteristic profile shapes: profiles exhibiting a (i) peak within $2.5$ half-mass radii and a subsequent decrease (ii) continuous increase that plateaus at larger radii typically with a high amplitude (iii) completely flat behaviour typically with low amplitude, in agreement with observations. This shows that the kinematic state of the stellar component can vary significantly with radius, suggesting a distinct interplay between in-situ star formation and ex-situ accretion of stars. Following the evolution of our sample through time, we provide evidence that the accretion history of galaxies with decreasing profiles is dominated by the anisotropic accretion of low mass satellites that get disrupted beyond $ \sim 2.0$ half-mass radii, building up a stellar halo with non-ordered motion while maintaining the central rotation already present at $z=2$. In fact, at $z=2$ decreasing profiles are the predominant profile class. Hence, we can predict a distinct formation pathway for galaxies with a decreasing profile and show that the centre resembles an old embedded disk. Furthermore, we show that the radius of the kinematic transition provides a good estimation for the transition radius from in-situ stars in the centre to accreted stars in the halo.

preprint2020arXiv

Stellar Property Statistics of Massive Halos from Cosmological Hydrodynamics Simulations: Common Kernel Shapes

We study stellar property statistics, including satellite galaxy occupation, of massive halo populations realized by three cosmological hydrodynamics simulations: BAHAMAS + MACSIS, TNG300 of the IllustrisTNG suite, and Magneticum Pathfinder. The simulations incorporate independent sub-grid methods for astrophysical processes with spatial resolutions ranging from $1.5$ to $6$ kpc, and each generates samples of $1000$ or more halos with $M_{\rm halo}> 10^{13.5} M_{\odot}$ at redshift $z=0$. Applying localized, linear regression (LLR), we extract halo mass-conditioned statistics (normalizations, slopes, and intrinsic covariance) for a three-element stellar property vector consisting of: i) $N_{sat}$, the number of satellite galaxies with stellar mass, $M_{\star, \rm sat} > 10^{10} M_{\odot}$ within radius $R_{200c}$ of the halo; ii) $M_{\star,\rm tot}$, the total stellar mass within that radius, and; iii) $M_{\star,\rm BCG}$, the gravitationally-bound stellar mass of the central galaxy within a $100 \, \rm kpc$ radius. Scaling parameters for the three properties with halo mass show mild differences among the simulations, in part due to numerical resolution, but there is qualitative agreement on property correlations, with halos having smaller than average central galaxies tending to also have smaller total stellar mass and a larger number of satellite galaxies. Marginalizing over total halo mass, we find the satellite galaxy kernel, $p(\ln N_{sat}\,|\,M_{\rm halo},z)$ to be consistently skewed left, with skewness parameter $γ= -0.91 \pm 0.02$, while that of $\ln M_{\star,\rm tot}$ is closer to log-normal, in all three simulations. The highest resolution simulations find $γ\simeq -0.8$ for the $z=0$ shape of $p(\ln M_{\star,\rm BCG}\,|\,M_{\rm halo},z)$ and also that the fractional scatter in total stellar mass is below $10\%$ in halos more massive than $10^{14.3} M_{\odot}$.

preprint2019arXiv

On the origin of magnetic driven winds and the structure of the galactic dynamo in isolated galaxies

We investigate the build-up of the galactic dynamo and subsequently the origin of a magnetic driven outflow. We use a setup of an isolated disc galaxy with a realistic circum-galactic medium (CGM). We find good agreement of the galactic dynamo with theoretical and observational predictions from the radial and toroidal components of the magnetic field as function of radius and disc scale height. We find several field reversals indicating dipole structure at early times and quadrupole structure at late times. Together with the magnetic pitch angle and the dynamo control parameters $R_α$, $R_ω$ and $D$ we present strong evidence for an $α^2$-$Ω$ dynamo. The formation of a bar in the centre leads to further amplification of the magnetic field via adiabatic compression which subsequently drives an outflow. Due to the Parker-Instability the magnetic field lines rise to the edge of the disc, break out and expand freely in the CGM driven by the magnetic pressure. Finally, we investigate the correlation between magnetic field and star formation rate. Globally, we find that the magnetic field is increasing as function of the star formation rate surface density with a slope between $0.3$ and $0.45$ in good agreement with predictions from theory and observations. Locally, we find that the magnetic field can decrease while star formation increases. We find that this effect is correlated with the diffusion of magnetic field from the spiral arms to the inter-arm regions which we explicitly include by solving the induction equation and accounting for non-linear terms.

preprint2017arXiv

Simulation-based marginal likelihood for cluster strong lensing cosmology

Comparisons between observed and predicted strong lensing properties of galaxy clusters have been routinely used to claim either tension or consistency with $Λ$CDM cosmology. However, standard approaches to such cosmological tests are unable to quantify the preference for one cosmology over another. We advocate approximating the relevant Bayes factor using a marginal likelihood that is based on the following summary statistic: the posterior probability distribution function for the parameters of the scaling relation between Einstein radii and cluster mass, $α$ and $β$. We demonstrate, for the first time, a method of estimating the marginal likelihood using the X-ray selected $z>0.5$ MACS clusters as a case in point and employing both N-body and hydrodynamic simulations of clusters. We investigate the uncertainty in this estimate and consequential ability to compare competing cosmologies, that arises from incomplete descriptions of baryonic processes, discrepancies in cluster selection criteria, redshift distribution, and dynamical state. The relation between triaxial cluster masses at various overdensities provide a promising alternative to the strong lensing test.