Cosmological Simulations Research Articles

Sensitivity Analysis of Simulation-Based Inference for Galaxy Clustering

Abstract Simulation-based inference (SBI) is a promising approach to leverage high fidelity cosmological simulations and extract information from the non-Gaussian, non-linear scales that cannot be modeled analytically. However, scaling SBI to the next generation of cosmological surveys faces the computational challenge of requiring a large number of accurate simulations over a wide range of cosmologies, while simultaneously encompassing large cosmological volumes at high resolution. This challenge can potentially be mitigated by balancing the accuracy and computational cost for different components of the the forward model while ensuring robust inference. To guide our steps in this, we perform a sensitivity analysis of SBI for galaxy clustering on various components of the cosmological simulations: gravity model, halo-finder and the galaxy-halo distribution models (halo-occupation distribution, HOD). We infer the σ8 and Ωm using galaxy power spectrum multipoles and the bispectrum monopole assuming a galaxy number density expected from the luminous red galaxies observed using the Dark Energy Spectroscopy Instrument (DESI). We find that SBI is insensitive to changing gravity model between N-body simulations and particle mesh (PM) simulations. However, changing the halo-finder from friends-of-friends (FoF) to Rockstar can lead to biased estimate of σ8 based on the bispectrum. For galaxy models, training SBI on more complex HOD leads to consistent inference for less complex HOD models, but SBI trained on simpler HOD models fails when applied to analyze data from a more complex HOD model. Based on our results, we discuss the outlook on cosmological simulations with a focus on applying SBI approaches to future galaxy surveys.

Stellar stripping efficiencies of satellites in numerical simulations: the effect of resolution, satellite properties and numerical disruption

Abstract The stellar stripping of satellites in cluster haloes is understood to play an important role in the production of intracluster light. Increasingly, cosmological simulations have been utilised to investigate its origin and assembly. However, such simulations typically model individual galaxies at relatively coarse resolutions, raising concerns about their accuracy. Although there is a growing literature on the importance of numerical resolution for the accurate recovery of the mass loss rates of dark matter (DM) haloes, there has been no comparable investigation into the numerical resolution required to accurately recover stellar mass loss rates in galaxy clusters. Using N-body simulations of satellite galaxies orbiting in a cluster halo represented by a static external potential, we conduct a set of convergence tests in order to explore the role of numerical resolution and force softening length on stellar stripping efficiency. We consider a number of orbital configurations, satellite masses and satellite morphologies. We find that stellar mass resolution is of minor importance relative to DM resolution. Resolving the central regions of satellite DM halos is critical to accurately recover stellar mass loss rates. Poorly resolved DM haloes develop cored inner profiles and, if this core is of comparable size to the stellar component of the satellite galaxy, this leads to significant over-stripping. To prevent this, relatively high DM mass resolutions of around MDM ∼ 106 M⊙, better than those achieved by many contemporary cosmological simulations, are necessary.

Cold Gas Subgrid Model (CGSM): A two-fluid framework for modeling unresolved cold gas in galaxy simulations

Abstract The cold (∼104 K) component of the circumgalactic medium (CGM) accounts for a significant fraction of all galactic baryons. However, using current galaxy-scale simulations to determine the origin and evolution of cold CGM gas poses a significant challenge, since it is computationally infeasible to directly simulate a galactic halo alongside the sub-pc scales that are crucial for understanding the interactions between cold CGM gas and the surrounding ‘hot’ medium. In this work, we introduce a new approach: the Cold Gas Subgrid Model (CGSM), which models unresolved cold gas as a second fluid in addition to the standard ‘normal’ gas fluid. The CGSM tracks the total mass density and bulk momentum of unresolved cold gas, deriving the properties of its unresolved cloudlets from the resolved gas phase. The interactions between the subgrid cold fluid and the resolved fluid are modeled by prescriptions from high-resolution simulations of ‘cloud crushing’ and thermal instability. Through a series of idealized tests, we demonstrate the CGSM’s ability to overcome the resolution limitations of traditional hydrodynamics simulations, successfully capturing the correct cold gas mass, its spatial distribution, and the timescales for cloud destruction and growth. We discuss the implications of using this model in cosmological simulations to more accurately represent the microphysics that govern the galactic baryon cycle.

The Stellar Bar–Dark Matter Halo Connection in the TNG50 Simulations

Stellar bars in disk galaxies grow as stars in near-circular orbits lose angular momentum to their environments, including their dark matter (DM) halo, and transform into elongated bar orbits. This angular momentum exchange during galaxy evolution hints at a connection between bar properties and the DM halo spin λ, the dimensionless form of DM angular momentum. We investigate the connection between halo spin λ and galaxy properties in the presence/absence of stellar bars, using the cosmological magnetohydrodynamic TNG50 simulations at multiple redshifts (0 < z r < 1). We determine the bar strength (or bar amplitude, A 2/A 0), using Fourier decomposition of the face-on stellar density distribution. We determine the halo spin for barred and unbarred galaxies (0 < A 2/A 0 < 0.7) in the center of the DM halo, close to the galaxy’s stellar disk. At z r = 0, there is an anticorrelation between halo spin and bar strength. Strongly barred galaxies (A 2/A 0 > 0.4) reside in DM halos with low spin and low specific angular momentum at their centers. In contrast, unbarred/weakly barred galaxies (A 2/A 0 < 0.2) exist in halos with higher central spin and higher specific angular momentum. The anticorrelation is due to the barred galaxies’ higher DM mass and lower angular momentum than the unbarred galaxies at z r = 0, as a result of galaxy evolution. At high redshifts (z r = 1), all galaxies have higher halo spin compared to those at lower redshifts (z r = 0), with a weak anticorrelation for galaxies having A 2/A 0 > 0.2. The formation of DM bars in strongly barred systems highlights how angular momentum transfer to the halo can influence its central spin.

Osaka Feedback Model. III. Cosmological Simulation CROCODILE

We introduce our new cosmological simulation data set CROCODILE, executed using the GADGET4-Osaka smoothed particle hydrodynamics code. This simulation incorporates an updated supernova (SN) feedback model of Y. Oku et al. and an active galactic nuclei (AGN) feedback model. A key innovation in our SN feedback model is the integration of a metallicity- and redshift-dependent, top-heavy initial mass function. Our SN model introduces a new consideration that results in an order of magnitude difference in the energy injection rate per unit stellar mass formed at high redshift. The CROCODILE data set is comprehensive, encompassing a variety of runs with diverse feedback parameters. This allows for an in-depth exploration of the relative impacts of different feedback processes in galactic evolution. Our initial comparisons with observational data, spanning the galaxy stellar mass function, the star formation main sequence, and the mass–metallicity relation, show promising agreement, especially for the Fiducial run. These results establish a solid foundation for our future work. We find that SN feedback is a key driver in the chemical enrichment of the intergalactic medium (IGM). Additionally, the AGN feedback creates metal-rich, bipolar outflows that extend and enrich the circumgalactic medium and IGM over a few Mpc scales.

Detecting Detached Black Hole Binaries through Photometric Variability

Understanding the connection between the properties of black holes (BHs) and their progenitors is interesting in many branches of astrophysics. Discovering BHs in detached orbits with luminous companions (LCs) promises to help establish this connection since the LC and BH progenitor are expected to have the same metallicity and formation time. We explore the possibility of detecting BH–LC binaries in detached orbits using photometric variations of the LC flux, induced by tidal ellipsoidal variation, relativistic beaming, and self-lensing. We create realistic present-day populations of detached BH–LC binaries in the Milky Way (MW) using binary population synthesis where we adopt observationally motivated initial stellar and binary properties, star formation history, and the present-day distribution of these sources in the MW based on detailed cosmological simulations. We test detectability of these sources via photometric variability by Gaia and TESS missions by incorporating their respective detailed detection biases as well as interstellar extinction. We find that Gaia is expected to resolve ∼300–1000 (∼700–1500) detached BH–LC binaries with signal-to-noise ratio (SNR) ≥ 10 (1) depending on the photometric precision and details of supernova physics. Similarly, the number of resolved BH–LC binaries with TESS is ∼50–200 (∼140–350). We find that 136−15+15 BH–LC binaries would be common between Gaia and TESS. Moreover, ∼60–70 (∼50–200) BH–LC binaries identifiable using photometry with SNR ≥ 10 may also be resolved using Gaia’s radial velocity (astrometry).

Toward Implementation of the Pressure-regulated, Feedback-modulated Model of Star Formation in Cosmological Simulations: Methods and Application to TNG

Traditional star formation subgrid models implemented in cosmological galaxy formation simulations, such as that of V. Springel & L. Hernquist (hereafter SH03), employ adjustable parameters to satisfy constraints measured in the local Universe. In recent years, however, theory and spatially resolved simulations of the turbulent, multiphase, star-forming interstellar medium (ISM) have begun to produce new first-principles models, which when fully developed can replace traditional subgrid prescriptions. This approach has advantages of being physically motivated and predictive rather than empirically tuned, and allowing for varying environmental conditions rather than being tied to local-Universe conditions. As a prototype of this new approach, by combining calibrations from the TIGRESS numerical framework with the pressure-regulated feedback-modulated (PRFM) theory, simple formulae can be obtained for both the gas depletion time and an effective equation of state. Considering galaxies in TNG50, we compare the “native” simulation outputs with postprocessed predictions from PRFM. At TNG50 resolution, the total midplane pressure is nearly equal to the total ISM weight, indicating that galaxies in TNG50 are close to satisfying vertical equilibrium. The measured gas scale height is also close to theoretical equilibrium predictions. The slopes of the effective equations of states are similar, but with effective velocity dispersion normalization from SH03 slightly larger than that from current TIGRESS simulations. Because of this and the decrease in PRFM feedback yield at high pressure, the PRFM model predicts shorter gas depletion times than the SH03 model at high densities and redshift. Our results represent a first step toward implementing new, numerically calibrated subgrid algorithms in cosmological galaxy formation simulations.

CHEX-MATE: The intracluster medium entropy distribution in the gravity-dominated regime

We characterise the intracluster gas entropy profiles of 32 very high-mass ($M_ $ M$_ odot Planck SZ-detected galaxy clusters (HIGHMz), selected from the CHEX-MATE sample, allowing us to study the intracluster medium (ICM) entropy distribution in a regime where non-gravitational effects are expected to be minimised. Using XMM- Newton measurements, we determined the entropy profiles up to $ $ for all objects. We assessed the relative role of gas density and temperature measurements on the uncertainty in entropy reconstruction, showing that in the outer regions the largest contribution comes from the temperature. The scaled profiles exhibit a large dispersion in the central regions, but converge rapidly to the value expected from simple gravitational collapse beyond the core regions. We quantified the correlation between the ICM morphological parameters and scaled entropy as a function of radius, showing that centrally peaked objects have low central entropy, while morphologically disturbed objects have high central entropy. We compared the scaled HIGHMz entropy profiles to results from other observational samples, finding differences in normalisation, which appear linked to the average mass of the samples in question. Combining HIGHMz with other samples, we found that a weaker mass dependence than self-similar in the scaling ($A_m -0.25$) allows us to minimise the dispersion in the radial range $ for clusters spanning over a decade in mass. The deviation from self-similar predictions is radially dependent and is more pronounced at small and intermediate radii than at $R_ $. We also investigated the distribution of central entropy $K_0$, finding no evidence for bimodality in the data and outer slope alpha , which peaks at $ 1.1$ with tails at both low and high alpha that correlate with dynamical state. Using weak-lensing masses for half of the sample, we found an indication for a small suppression of the scatter ($ beyond the core when using masses derived from $Y_X$ in the rescaling. Finally, we compared our results to recent cosmological numerical simulations from HE T HREE H UNDRED and MACSIS, finding good agreement with the observational data in this mass regime. These results provide a robust observational benchmark in the gravity-dominated regime, and will serve as a future reference for samples at lower masses, higher redshifts, and for ongoing work using cosmological numerical simulations.

How Do Primordial Black Holes Change the Halo Mass Function and Structure?

We examine the effects of massive primordial black holes (PBHs) on cosmic structure formation, employing both a semianalytical approach and cosmological simulations. Our simulations incorporate PBHs with a monochromatic mass distribution centered around 106 M ⊙, constituting a fraction of 10−2 to 10−4 of the dark matter (DM) in the Universe, with the remainder being collisionless particle DM. Additionally, we conduct a ΛCDM simulation for comparative analysis with runs that include PBHs. At smaller scales, halos containing PBHs exhibit similar density and velocity dispersion profiles to those without PBHs. Conversely, at larger scales, PBHs can expedite the formation of massive halos and reside at their centers owing to the “seed effect.” To analyze the relative distribution of PBH host halos compared to non-PBH halos, we apply nearest neighbor statistics. Our results suggest that PBH host halos, through gravitational influence, significantly impact the structure formation process, compared to the ΛCDM case, by attracting and engulfing nearby newly formed minihalos. Should PBHs constitute a fraction of DM significantly larger than ∼10−3, almost all newly formed halos will be absorbed by PBH-seeded halos. Consequently, our simulations predict a bimodal feature in the halo mass function, with most of the massive halos containing at least one PBH at their core and the rest being less massive non-PBH halos.

Learning the Universe: GalactISM Simulations of Resolved Star Formation and Galactic Outflows across Main-sequence and Quenched Galactic Environments

We present a suite of six high-resolution chemodynamical simulations of isolated galaxies, spanning observed disk-dominated environments on the star-forming main sequence, as well as quenched, bulge-dominated environments. We compare and contrast the physics driving star formation and stellar feedback among the galaxies, with a view to modeling these processes in cosmological simulations. We find that the mass loading of galactic outflows is coupled to the clustering of supernova explosions, which varies strongly with the rate of galactic rotation Ω = v circ/R via the Toomre length, leading to smoother gas disks in the bulge-dominated galaxies. This sets an equation of state in the star-forming gas that also varies strongly with Ω, so that the bulge-dominated galaxies have higher midplane densities, lower velocity dispersions, and higher molecular gas fractions than their main-sequence counterparts. The star formation rate in five out of six galaxies is independent of Ω and is consistent with regulation by the midplane gas pressure alone. In the sixth galaxy, which has the most centrally concentrated bulge and thus the highest Ω, we reproduce dynamical suppression of the star formation efficiency in agreement with observations. This produces a transition away from pressure-regulated star formation.

Magnetic fields in the outskirts of PSZ2 G096.88+24.18 from a depolarization analysis of radio relics

Context. Radio relics are diffuse, non-thermal radio sources present in a number of merging galaxy clusters. They are characterized by elongated arc-like shapes and highly polarized emission (up to ∼60%) at gigahertz frequencies, and are expected to trace shock waves in the cluster outskirts induced by galaxy cluster mergers. Their polarized emission can be used to study the magnetic field properties of the host cluster. Aims. In this paper, we investigate the polarization properties of the double radio relics in PSZ2 G096.88+24.18 using the rotation measure (RM) synthesis, and try to constrain the characteristics of the magnetic field that reproduce the observed depolarization as a function of resolution (beam depolarization). Our aim is to understand the nature of the low polarization fraction that characterizes the southern relic with respect to the northern relic. Methods. We present new 1–2 GHz Karl G. Jansky Very Large Array (VLA) observations in multiple configurations. We derived the RM and polarization of the two relics by applying the RM synthesis technique, and thus solved for bandwidth depolarization in the wide observing bandwidth. To study the effect of beam depolarization, we degraded the image resolution and studied the decreasing trend of polarization fraction with increasing beam size. Finally, we performed 3D magnetic field simulations using multiple models for the magnetic field power spectrum over a wide range of scales, in order to constrain the characteristics of the cluster magnetic field that can reproduce the observed beam depolarization trend. Results. Using RM synthesis, we obtained a polarization fraction of (18.6 ± 0.3)% for the northern relic and (14.6 ± 0.1)% for the southern one. Having corrected for bandwidth depolarization, and after noticing the absence of relevant complex Faraday spectrum, we inferred that the nature of the depolarization for the southern relic is external, and possibly related to the turbulent gas distribution within the cluster, or to the complex spatial structure of the relic. The best-fit magnetic field power spectrum, which reproduces the observed depolarization trend for the southern relic, was obtained for a turbulent magnetic field model, described by a power spectrum derived from cosmological simulations, and defined within the scales of Λmin = 35 kpc and Λmax = 400 kpc. This yields an average magnetic field of the cluster within 1 Mpc3 volume of ∼2 μG.

Simba-EoR: Early galaxy formation in the Simba simulation including a new sub-grid interstellar medium model

Abstract We update the dust model present within the Simba galaxy simulations with a self-consistent framework for the co-evolution of dust and molecular hydrogen populations in the interstellar medium, and use this to explore z ≥ 6 galaxy evolution. In addition to tracking the evolution of dust and molecular hydrogen abundances, our model fully integrates these species into the Simba simulation, explicitly modelling their impact on physical processes such as star formation and cooling through the inclusion of a novel two-phase sub-grid model for interstellar gas. Running two cosmological simulations down to z ∼ 6 we find that our Simba-EoR model displays a generally tighter concordance with observational data than fiducial Simba. Additionally we observe that our Simba-EoR models increase star formation activity at early epochs, producing larger dust-to-gas ratios consequently. Finally, we discover a significant population of hot dust at ∼100 K, aligning with contemporaneous observations of high-redshift dusty galaxies, alongside the large ∼20 K population typically identified.

Cosmological Radiation Simulation Analysis: Taking Arepo-RT as an Example

With the fast-growing progress in computational physics, astrophysical simulations are made possible and continually upgraded with novel physical models. Contemporarily, many modern simulations feature the capability of cosmological radiation. In order to collate the recent improvements in this field, the numerical solver Arepo-RT is discussed as an important example. This study reorganizes and interprets the main physical models and numerical techniques used in Arepo-RT and similar solvers, with some details simplified to suit common readers. This paper describes simply and briefly the finite volume method, the radiative transfer model, the radiation source model, the Godunov-type schemes and the time integration, all of which are frequently used in modern radiation simulations. Apart from the realizations of the solvers, the results and limitations are also outlined in this paper, including the justification of reduced speed of light approximation, the results compared with JWST, and the issue of missing photons. These results shed light on guiding further improvements for Arepo-RT.

Comparing Performance in the Codes AREPO, GIZMO and GAMER-2

Modern astrophysical simulation is the main aspect of human getting insight from the cosmos. Numerical simulations became irreplaceable tools for generating various phases of the universe, the formation of galaxies, and other astrophysical phenomena. As the advances of computational power, validity and diversity of astrophysical modeling have been growing. In addition to the development of the processor performance, the diversity of the software community is also rising therefore, research between various simulation codes is worthwhile. It allowed researchers to know the nature of different simulation techniques linked to their basic algorithms. This study summarizes the basics of common simulation methods used by researchers including N-body Simulations, Smoothed Particle Hydrodynamics, Monte Carlo Methods, Finite Difference and Finite Element Methods. Furthermore, this research analyzes the basics of AREPO, GIZMO, and GAMER-2 codes, and compared the performance of these codes. Last, according to the analysis, this study discusses and demonstrates the advantages and limitations of these codes. These results shed light on guiding further exploration of cosmology simulations.

Recovering chemical bimodalities in observed edge-on stellar disks: Insights from AURIGA simulations

The well-known bimodal distribution of Milky Way disk stars in the [α/Fe]–metallicity plane is often used to define thick and thin disks. In external edge-on galaxies, there have been attempts to identify this type of bimodality using integral field spectroscopy (IFS) data. However, for unresolved stellar populations, observations only contain integrated information, making these studies challenging. We assessed the ability to recover chemical bimodalities in IFS observations of edge-on galaxies, using 24 Milky Way-mass galaxies from the AURIGA zoom-in cosmological simulations. We first analyzed the distribution of single stellar particles in the [Mg/Fe]–[Fe/H] plane, finding that bimodality is frequent but not ubiquitous and often unclear. Then we produced mock IFS [Mg/Fe] and [Fe/H] maps of galaxies seen edge on, and considered integrated stellar-population properties (projected and spatially binned). We investigated how the distribution of stars in the [Mg/Fe]–[Fe/H] plane is affected by edge-on projection and spatial binning. Bimodality is preserved, while distributions change their shapes. Naturally, broad distributions of individual star particles are narrowed into smaller [Mg/Fe] and [Fe/H] ranges for spatial bins. We observe continuous distributions from high [Mg/Fe] and low [Fe/H], to lower [Mg/Fe] values and higher [Fe/H]. Despite being continuous, these distributions are bimodal in most cases. The overlap in [Fe/H] is small, and different [Mg/Fe] components show up as peaks instead of sequences (even when the latter are present for individual particles). The larger the spatial bins, the narrower the [Mg/Fe]–[Fe/H] distribution. This narrowing helps amplify the density of different [Mg/Fe] peaks, often leading to a clearer bimodality in mock IFS observations than for original star particles. We also assessed the correspondence of chemical bimodalities with the distinction between geometric thick and thin disks. Their individual particles have different distributions, but mostly overlap in [Mg/Fe] and [Fe/H]. However, integrated properties of geometric thick and thin disks in mock maps do mostly segregate into different regions of the [Mg/Fe]–[Fe/H] plane. In bimodal distributions, they correspond to the two distinct peaks. Our results show that this approach can be used for bimodality studies in future IFS observations of edge-on external galaxies.

Gas-phase metallicity gradients in galaxies at z ∼ 6–8

The study of gas-phase metallicity and its spatial distribution at high redshift is crucial to understand the processes that shaped the growth and evolution of galaxies in the early Universe. Here we study the spatially resolved metallicity in three systems at z ∼ 6 − 8, namely A2744-YD4, BDF-3299, and COSMOS24108, with JWST NIRSpec IFU low-resolution (R ∼ 100) spectroscopic observations. These are among the highest-z sources in which metallicity gradients have been probed so far. Each of these systems hosts several spatial components in the process of merging within a few kiloparsecs, identified from the rest-frame UV and optical stellar continuum and ionised gas emission line maps. The sources have heterogeneous properties, with stellar masses log(M*/M⊙) ∼7.6–9.3, star formation rates (SFRs) ∼1–15 M⊙ yr−1, and gas-phase metallicities 12+log(O/H) ∼7.7–8.3, which exhibit a large scatter within each system. Their properties are generally consistent with those of the highest-redshift samples to date (z ∼ 3 − 10), though the sources in A2744-YD4 and COSMOS24108 are at the high end of the mass-metallicity relation (MZR) defined by the z ∼ 3 − 10 sources. Moreover, the targets in this work follow the predicted slope of the MZR at z ∼ 6 − 8 from most cosmological simulations. The gas-phase metallicity gradients are consistent with being flat in the main sources of each system. Flat metallicity gradients are thought to arise from gas mixing processes on galaxy scales, such as mergers or galactic outflows and supernova winds driven by intense stellar feedback, which wash out any gradient formed in the galaxy. The existence of flat gradients at z ∼ 6 − 8 sets also important constraints on future cosmological simulations and chemical evolution models, whose predictions on the cosmic evolution of metallicity gradients often differ significantly, especially at high redshift, but are mostly limited to z ≲ 3 so far.

Topological Approach to Void Finding Applied to the SDSS Galaxy Map

The structure of the low redshift Universe is dominated by a multiscale void distribution delineated by filaments and walls of galaxies. The characteristics of voids, such as morphology, average density profile, and correlation function, can be used as cosmological probes. However, their physical properties are difficult to infer due to shot noise and the general lack of tracer particles used to define them. In this work, we construct a robust, topology-based void-finding algorithm that utilizes Persistent Homology to detect persistent features in the data. We apply this approach to a volume-limited subsample of galaxies in the SDSS I/II Main Galaxy catalog with the r-band absolute magnitude brighter than M r = −20.19, and a set of mock catalogs constructed using the Horizon Run 4 cosmological N-body simulation. We measure the size distribution of voids, their averaged radial profile, sphericity, and the centroid nearest neighbor separation, using conservative values for the threshold and persistence. We find 32 topologically robust voids in the SDSS data over the redshift range 0.02 ≤ z ≤ 0.116, with effective radii in the range 21−56 h −1 Mpc. The median nearest neighbor void separation is found to be ∼57 h −1 Mpc, and the median radial void profile is consistent with the expected shape from the mock data.

On the Corotation of Milky Way Satellites: LMC-mass Satellites Induce Apparent Motions in Outer Halo Tracers

Abstract Understanding the physical mechanism behind the formation of a corotating thin plane of satellite galaxies, like the one observed around the Milky Way (MW), has been challenging. The perturbations induced by a massive satellite galaxy, like the Large Magellanic Cloud (LMC), provide valuable insight into this problem. The LMC induces an apparent corotating motion in the outer halo by displacing the inner regions of the halo with respect to the outer halo. Using the Latte suite of Feedback In Realistic Environments cosmological simulations of MW-mass galaxies, we confirm that the apparent motion of the outer halo induced by the infall of a massive satellite changes the observed distribution of orbital poles of outer-halo tracers, including satellites. We quantify the changes in the distribution of orbital poles using the two-point angular correlation function and find that all satellites induce changes. However, the most massive satellites with pericentric passages between ≈30 and 100 kpc induce the largest changes. The best LMC-like satellite analog shows the largest change in orbital pole distribution. The dispersion of orbital poles decreases by 20° during the first two pericentric passages. Even when excluding the satellites brought in with the LMC-like satellite, there is clustering of orbital poles. These results suggest that in the MW, the recent pericentric passage of the LMC should have changed the observed distribution of orbital poles of all other satellites. Therefore, studies of kinematically coherent planes of satellites that seek to place the MW in a cosmological context should account for the existence of a massive satellite like the LMC.

Untangling stellar components of galaxies: Evaluation of dynamical decomposition methods in simulated galaxies with GalaxyChop

Galaxy formation is intrinsically connected to the distinct evolutionary processes of disk and spheroidal systems, which are the fundamental stellar components of galaxies. Understanding the mutual, dynamical interplay and co-evolution of these components requires a detailed dynamical analysis to allow for disentanglement of these systems. We introduce JEHistogram, a new method for the dynamical decomposition of simulated galaxies into disk and spheroidal stellar components that utilizes the angular momentum and energy of star particles. We evaluate its performance against five previously established methods using a sample of equilibrium galaxies with stellar masses in the range $10^ M_ gal /M_ $. Our assessment involves several metrics, including the completeness and purity of stellar particle classification, scale lengths, mass density profiles, velocity dispersion, and rotational velocity profiles. While all methods approximate the properties of the original components, such as mass fractions and density or velocity profiles, JEHistogram demonstrates a better accuracy, particularly in the inner regions of galaxies where component overlap complicates separation. Additionally, we apply JEHistogram to a Milky Way-like galaxy from the IllustrisTNG cosmological simulations, showcasing its capability to derive properties such as the size, mass, velocity, color, and age of dynamically defined disk and spheroidal components. All the dynamical decomposition methods we analyzed are publicly accessible through the Python package GalaxyChop

The glow of axion quark nugget dark matter

Context. The existence of axion quark nuggets is a potential consequence of the axion field, which provides a possible solution to the charge-conjugation parity violation in quantum chromodynamics. In addition to explaining the cosmological discrepancy of matter-antimatter asymmetry and a visible-to-dark-matter ratio of Ωdark/Ωvisible ≃ 5, these composite compact objects are expected to represent a potentially ubiquitous electromagnetic background radiation by interacting with ordinary baryonic matter. We conducted an in-depth analysis of axion quark nugget-baryonic matter interactions in the environment of the intracluster medium in the constrained cosmological Simulation of the LOcal Web (SLOW). Aims. Here, we aim to provide upper limit predictions on electromagnetic counterparts of axion quark nuggets in the environment of galaxy clusters by inferring their thermal and non-thermal emission spectrum originating from axion quark nugget-cluster gas interactions. Methods. We analyzed the emission of axion quark nuggets in a large sample of 161 simulated galaxy clusters using the SLOW simulation. These clusters are divided into a sub-sample of 150 galaxy clusters, ordered in five mass bins ranging from 0.8 to 31.7 × 1014 M⊙, along with 11 cross-identified galaxy clusters from observations. We investigated dark matter-baryonic matter interactions in galaxy clusters in their present stage at the redshift of z = 0 by assuming all dark matter consists of axion quark nuggets. The resulting electromagnetic signatures were compared to thermal Bremsstrahlung and non-thermal cosmic ray (CR) synchrotron emission in each galaxy cluster. We further investigated individual frequency bands imitating the observable range of the WMAP, Planck, Euclid, and XRISM telescopes for the most promising cross-identified galaxy clusters hosting detectable signatures of axion quark nugget emission. Results. We observed a positive excess in the low- and high-energy frequency windows, where thermal and non-thermal axion quark nugget emission can significantly contribute to (or even outshine) the emission of the intracluster medium (ICM) in frequencies up to νT ≲ 3842.19 GHz and νT ϵ [3.97, 10.99] × 1010GHz, respectively. Emission signatures of axion quark nuggets are found to be observable if CR synchrotron emission of individual clusters is sufficiently low. The degeneracy in the parameters contributing to an emission excess makes it challenging to offer predictions with respect to pinpointing specific regions of a positive axion quark nugget excess; however, a general increase in the total galaxy cluster emission is expected based on this dark matter model. Axion quark nuggets constitute an increment of 4.80% of the total galaxy cluster emission in the low-energy regime of νT ≲ 3842.19 GHz for a selection of cross-identified galaxy clusters. We propose that the Fornax and Virgo clusters represent the most promising candidates in the search for axion quark nugget emission signatures. Conclusions. The results from our simulations point towards the possibility of detecting an axion quark nugget excess in galaxy clusters in observations if their signatures can be sufficiently disentangled from the ICM radiation. While this model proposes a promising explanation for the composition of dark matter, with the potential to have this outcome verified by observations, we propose further changes that are aimed at refining our methods. Our ultimate goal is to identify the extracted electromagnetic counterparts of axion quark nuggets with even greater precision in the near future.