Cosmological Zoom-in Simulations Research Articles

Constraining the physical properties of gas in high-z galaxies with far-infrared and submillimetre line ratios

Optical emission line diagnostics, which are a common tool for constraining the properties of the interstellar medium (ISM) of galaxies, become progressively inaccessible at higher redshifts for ground-based facilities. Far-infrared (FIR) emission lines, which are redshifted into atmospheric windows that are accessible for ground-based submillimetre facilities, could provide ISM diagnostics alternative to optical emission lines. We investigated FIR line ratios involving CII OIII OIII NII and NIII using synthetic emission lines applied to a high-resolution gas $= 883.4 M$_ odot $) cosmological zoom-in simulation, including radiative transfer post-processing with the code at z = 6.5. We find that the CII NII 122 ratio is sensitive to the temperature and density of photodissociation regions. It might therefore be a useful tool for tracing the properties of this gas phase in galaxies. We also find that NII NIII is a good tracer of the temperature and that OIII OIII 88 is a good tracer of the gas density of HII regions. Emission line ratios containing the OIII line are sensitive to high-velocity outflowing gas.

The AGORA High-resolution Galaxy Simulations Comparison Project. IV. Halo and Galaxy Mass Assembly in a Cosmological Zoom-in Simulation at z ≤ 2

In this fourth paper from the AGORA Collaboration, we study the evolution down to redshift z = 2 and below of a set of cosmological zoom-in simulations of a Milky Way mass galaxy by eight of the leading hydrodynamic simulation codes. We also compare this CosmoRun suite of simulations with dark matter-only simulations by the same eight codes. We analyze general properties of the halo and galaxy at z = 4 and 3, and before the last major merger, focusing on the formation of well-defined rotationally supported disks, the mass–metallicity relation, the specific star formation rate, the gas metallicity gradients, and the nonaxisymmetric structures in the stellar disks. Codes generally converge well to the stellar-to-halo mass ratios predicted by semianalytic models at z ∼ 2. We see that almost all the hydro codes develop rotationally supported structures at low redshifts. Most agree within 0.5 dex with the observed mass–metallicity relation at high and intermediate redshifts, and reproduce the gas metallicity gradients obtained from analytical models and low-redshift observations. We confirm that the intercode differences in the halo assembly history reported in the first paper of the collaboration also exist in CosmoRun, making the code-to-code comparison more difficult. We show that such differences are mainly due to variations in code-dependent parameters that control the time stepping strategy of the gravity solver. We find that variations in the early stellar feedback can also result in differences in the timing of the low-redshift mergers. All the simulation data down to z = 2 and the auxiliary data will be made publicly available.

Is GN-z11 powered by a super-Eddington massive black hole?

Context. Observations of z ∼ 6 quasars powered by supermassive black holes (SMBHs; MBH ∼ 108 − 10 M⊙) challenge our current understanding of early black hole (BH) formation and evolution. The advent of the James Webb Space Telescope (JWST) has enabled the study of massive BHs (MBHs; MBH ∼ 106 − 7 M⊙) up to z ∼ 11, thus bridging the properties of z ∼ 6 quasars to their ancestors. Aims. The JWST spectroscopic observations of GN-z11, a well-known z = 10.6 star-forming galaxy, have been interpreted with the presence of a super-Eddington (Eddington ratio ≡ λEdd ∼ 5.5) accreting MBH. To test this hypothesis, we used a zoom-in cosmological simulation of galaxy formation and BH co-evolution. Methods. We first tested the simulation results against the observed probability distribution function (PDF) of λEdd found in z ∼ 6 quasars. Then, in the simulation we selected the BHs that satisfy the following criteria: (a) 10 < z < 11, (b) MBH > 106 M⊙. Next, we applied the extreme value statistics to the PDF of λEdd resulting from the simulation. Results. We find that the probability of observing a z ∼ 10 − 11 MBH accreting with λEdd ∼ 5.5 in the volume surveyed by JWST is very low (< 0.2%). We compared our predictions with those in the literature, and discussed the main limitations of our work. Conclusions. Our simulation cannot explain the JWST observations of GN-z11. This might be due to: (i) poor resolution and statistics in simulations, (ii) simplistic sub-grid models (e.g. BH accretion and seeding), (iii) uncertainties in the data analysis and interpretation.

Great Balls of FIRE

Context. Despite the nearly hundred gravitational-wave detections reported by the LIGO-Virgo-KAGRA Collaboration, the question of the cosmological origin of merging binary black holes (BBHs) remains open. The two main formation channels generally considered are from isolated field binaries or via dynamical assembly in dense star clusters. Aims. Here we focus on understanding the dynamical formation of merging BBHs within massive clusters in galaxies of different masses. Methods. To this end, we applied a new framework to consistently model the formation and evolution of massive star clusters in zoom-in cosmological simulations of galaxies. Each simulation, taken from the FIRE project, provides a realistic star formation environment, with a unique star formation history, that hosts realistic giant molecular clouds that constitute the birthplace of star clusters. Combined with the code for star cluster evolution CMC, we are able to produce populations of dynamically formed merging BBHs across cosmic time in different environments. Results. As the most massive star clusters preferentially form in dense massive clouds of gas, we find that, despite their low metallicities favouring the creation of black holes, low-mass galaxies contain few massive clusters and therefore make a limited contribution to the global production of dynamically formed merging BBHs. Furthermore, we find that massive clusters can host hierarchical BBH mergers with clear, identifiable physical properties. Looking at the evolution of the BBH merger rate in different galaxies, we find strong correlations between BBH mergers and the most extreme episodes of star formation. Finally, we discuss the implications for future LIGO-Virgo-KAGRA gravitational wave observations.

Sustained super-Eddington accretion in high-redshift quasars

Observations of z ≳ 6 quasars provide information on the early evolution of the most massive black holes (MBHs) and galaxies. Current observations, able to trace both gas and stellar properties, reveal a population of MBHs that is significantly more massive than expected from the local MBH-stellar mass relation. The population lies on, but mostly above, the relation observed in the nearby Universe. This suggests that these objects grew very rapidly. To explain their presence when the Universe was less than 1 Gyr old and to assess the physical conditions for their rapid growth, we explored whether episodes of accretion above the Eddington limit can occur across cosmic epochs. By employing state-of-the-art high-resolution cosmological zoom-in simulations of a z ∼ 7 quasar, where different accretion regimes are included consistently, together with their associated radiative and kinetic feedback, we show that super-Eddington phases can be sustained for relatively long timescales (tens of millions of years). This allows the MBH to rapidly grow by up to three orders of magnitude, depending on the strength of the kinetic feedback. We also show by means of a semianalytic calculation that the MBH spin remains moderate and does not take on extremely high values during the super-Eddington phases. This results in a lower feedback efficiency, which may allow the rapid growth required to explain over-massive high-redshift MBHs.

The High-Redshift Gas-Phase Mass–Metallicity Relation in FIRE-2

The unprecedented infrared spectroscopic capabilities of JWST have provided high-quality interstellar medium metallicity measurements and enabled characterization of the gas-phase mass–metallicity relation (MZR) for galaxies at z ≳ 5 for the first time. We analyze the gas-phase MZR and its evolution in a high-redshift suite of FIRE-2 cosmological zoom-in simulations at z = 5–12 and for stellar masses M * ∼ 106–1010 M ⊙. These simulations implement a multichannel stellar feedback model and produce broadly realistic galaxy properties, including when evolved to z = 0. The simulations predict very weak redshift evolution of the MZR over the redshift range studied, with the normalization of the MZR increasing by less than 0.01 dex as redshift decreases from z = 12 to z = 5. The median MZR in the simulations is well approximated as a constant power-law relation across this redshift range given by log(Z/Z⊙)=0.37log(M*/M⊙)−4.3 . We find good agreement between our best-fit model and recent observations made by JWST at high redshift. The weak evolution of the MZR at z > 5 contrasts with the evolution at z ≲ 3, where increasing normalization of the MZR with decreasing redshift is observed and predicted by most models. The FIRE-2 simulations predict increasing scatter in the gas-phase MZR with decreasing stellar mass, in qualitative agreement with some observations.

Radio relics in massive galaxy cluster mergers in the TNG-Cluster simulation

Radio relics are diffuse synchrotron sources in the outskirts of merging galaxy clusters energized by the merger shocks. In this paper, we present an overview of the radio relics in massive cluster mergers identified in the new TNG-Cluster simulation. This is a suite of magnetohydrodynamical cosmological zoom-in simulations of 352 massive galaxy clusters with M500c = 1014.0 − 15.3 M⊙ sampled from a 1 Gpc-sized cosmological box. The simulations were performed using the moving-mesh code AREPO with the galaxy formation model and high numerical resolution consistent with the TNG300 run of the IllustrisTNG series. We post-processed the shock properties obtained from the on-the-fly shock finder to estimate the diffuse radio emission generated by cosmological shockwaves for a total of ∼300 radio relics at redshift z = 0 − 1. TNG-Cluster returned a variety of radio relics with diverse morphologies, encompassing classical examples of double radio relics, single relics, and “inverted” radio relics that are convex to the cluster center. Moreover, the simulated radio relics reproduced both the abundance and statistical relations of observed relics. We find that extremely large radio relics (> 2 Mpc) are predominantly produced in massive cluster mergers with M500c ≳ 8 × 1014 M⊙. This underscores the significance of simulating massive mergers to study giant radio relics similar to those found in observations. We released a library of radio relics from the TNG-Cluster simulation, which will serve as a crucial reference for upcoming next-generation surveys.

Toward Cosmological Simulations of the Magnetized Intracluster Medium with Resolved Coulomb Collision Scale

We present the first results of one extremely high-resolution, nonradiative magnetohydrodynamical cosmological zoom-in simulation of a massive cluster with a virial mass of M vir = 2.0 × 1015 solar masses. We adopt a mass resolution of 4 × 105 M ⊙ with a maximum spatial resolution of around 250 pc in the central regions of the cluster. We follow the detailed amplification process in a resolved small-scale turbulent dynamo in the intracluster medium (ICM) with strong exponential growth until redshift 4, after which the field grows weakly in the adiabatic compression limit until redshift 2. The energy in the field is slightly reduced as the system approaches redshift zero in agreement with adiabatic decompression. The field structure is highly turbulent in the center and shows field reversals on a length scale of a few tens of kiloparsecs and an anticorrelation between the radial and angular field components in the central region that is ordered by small-scale turbulent dynamo action. The large-scale field on megaparsec scales is almost isotropic, indicating that the structure formation process in massive galaxy cluster formation suppresses any memory of both the initial field configuration and the amplified morphology via the turbulent dynamo. We demonstrate that extremely high-resolution simulations of the magnetized ICM are within reach that can simultaneously resolve the small-scale magnetic field structure, which is of major importance for the injection of and transport of cosmic rays in the ICM. This work is a major cornerstone for follow-up studies with an on-the-fly treatment of cosmic rays to model in detail electron-synchrotron and gamma-ray emissions.

Varying primordial state fractions in exo- and endothermic SIDM simulations of Milky Way-mass haloes

Abstract Self-interacting dark matter (SIDM) is increasingly studied as a potential solution to small-scale discrepancies between simulations of cold dark matter (CDM) and observations. We examine a physically motivated two-state SIDM model with both elastic and inelastic scatterings. In particular, endothermic, exothermic, and elastic scattering have equal transfer cross sections at high relative velocities (vrel ≳ 400 km s−1). In a suite of cosmological zoom-in simulation of Milky Way-size haloes, we vary the primordial state fractions to understand the impact of inelastic dark matter self-interactions on halo structure and evolution. In particular, we test how the initial conditions impact the present-day properties of dark matter haloes. Depending on the primordial state fraction, scattering reactions will be dominated by either exothermic or endothermic effects for high and low initial excited state fractions respectively. We find that increasing the initial excited fraction reduces the mass of the main halo, as well as the number of subhaloes on all mass scales. The main haloes are cored, with lower inner densities and higher outer densities compared with CDM. Additionally, we find that the shape of the main halo becomes more spherical the higher the initial excited state fraction is. Finally, we show that the number of satellites steadily decreases with initial excited state fraction across all satellite masses.

Dissipative Dark Substructure: The Consequences of Atomic Dark Matter on Milky Way Analog Subhalos

Dissipative Dark Substructure: The Consequences of Atomic Dark Matter on Milky Way Analog Subhalos

Forecasts for Galaxy Formation and Dark Matter Constraints from Dwarf Galaxy Surveys

The abundance of faint dwarf galaxies is determined by the underlying population of low-mass dark matter (DM) halos and the efficiency of galaxy formation in these systems. Here, we quantify potential galaxy formation and DM constraints from future dwarf satellite galaxy surveys. We generate satellite populations using a suite of Milky Way (MW)–mass cosmological zoom-in simulations and an empirical galaxy–halo connection model, and assess sensitivity to galaxy formation and DM signals when marginalizing over galaxy–halo connection uncertainties. We find that a survey of all satellites around one MW-mass host can constrain a galaxy formation cutoff at peak virial masses of M50=108M⊙ at the 1σ level; however, a tail toward low M50 prevents a 2σ measurement. In this scenario, combining hosts with differing bright satellite abundances significantly reduces uncertainties on M50 at the 1σ level, but the 2σ tail toward low M50 persists. We project that observations of one (two) complete satellite populations can constrain warm DM models with m WDM ≈ 10 keV (20 keV). Subhalo mass function (SHMF) suppression can be constrained to ≈70%, 60%, and 50% that in cold dark matter (CDM) at peak virial masses of 108, 109, and 1010 M ⊙, respectively; SHMF enhancement constraints are weaker (≈20, 4, and 2 times that in CDM, respectively) due to galaxy–halo connection degeneracies. These results motivate searches for faint dwarf galaxies beyond the MW and indicate that ongoing missions like Euclid and upcoming facilities including the Vera C. Rubin Observatory and Nancy Grace Roman Space Telescope will probe new galaxy formation and DM physics.

Dissipative Dark Matter on FIRE. II. Observational Signatures and Constraints from Local Dwarf Galaxies

We analyze the first cosmological baryonic zoom-in simulations of galaxies in dissipative self-interacting dark matter (dSIDM). The simulations utilize the FIRE-2 galaxy formation physics with the inclusion of dissipative dark matter self-interactions modeled as a constant fractional energy dissipation (f diss = 0.75). In this paper, we examine the properties of dwarf galaxies with M * ∼ 105–109 M ⊙ in both isolation and within Milky Way–mass hosts. For isolated dwarfs, we find more compact galaxy sizes and promotion of disk formation in dSIDM with (σ/m) ≤ 1 cm2 g−1. On the contrary, models with (σ/m) = 10 cm2 g−1 produce puffier stellar distributions that are in tension with the observed size–mass relation. In addition, owing to the steeper central density profiles, the subkiloparsec circular velocities of isolated dwarfs when (σ/m) ≥ 0.1 cm2 g−1 are enhanced by about a factor of 2, which are still consistent with the kinematic measurements of Local Group dwarfs but in tension with the H i rotation curves of more massive field dwarfs. Meanwhile, for satellites of Milky Way–mass hosts, the median circular velocity profiles are marginally affected by dSIDM physics, but dSIDM may help promote the structural diversity of dwarf satellites. The number of satellites is slightly enhanced in dSIDM, but the differences are small compared with the large host-to-host variations. In conclusion, the dSIDM models with (σ/m) ≳ 0.1 cm2 g−1, f diss = 0.75 are in tension in massive dwarfs (M halo ∼ 1011 M ⊙) due to circular velocity constraints. However, models with lower effective cross sections (at this halo mass/velocity scale) are still viable and can produce nontrivial observable signatures.

The origin of cold gas in the circumgalactic medium

Context. The presence of cold gas (T ≲ 104 K) in the circumgalactic medium (CGM) of galaxies has been confirmed in observations and in high-resolution simulations, but its origin is still a puzzle. Possible mechanisms are cold accretion from the intergalactic medium (IGM), clumps embedded in outflows and transported from the disk, and gas detaching from the hot CGM phase via thermal instabilities. Aims. In this work we characterize the history of cold CGM gas in order to identify the dominant origin channels at different evolutionary stages of the main galaxy. Methods. To this end, we tracked gas particles in different snapshots of the smoothed particle hydrodynamics (SPH) cosmological zoom-in simulation Eris2k. We performed a backward tracking of cold gas, starting from different redshifts until we could identify one of the followings origins for the particle: cold inflow, ejection from the disk, cooling down in situ, or stripping from a satellite. We also performed a forward tracking of gas in different components of the galaxy (such as the disk and outflows). Results. We find a clear transition between two epochs. For z > 2, most cold gas (up to 80%) in the CGM comes from cold accretion streams as the galaxy is accreting in the cold mode from the IGM. At lower z, gas either cools down in situ after several recycles (with 10–20% of the gas cooling in outflow), or it is ejected directly from the disk (up to 30%). Outflows have a major contribution to the cold CGM gas budget at z < 1, with almost 50% of the hot gas cooling in outflow. Finally, we discuss possible mechanisms for CGM cooling, showing that the thermally unstable gas with tcool/tff < 1 (precipitation-regulated feedback) is abundant up to r ∼ 100 kpc and cooling times are shorter than 50 Myr for densities n > 10−2 cm−3.

Gas and Star Formation in Satellites of Milky Way Analogs

We have imaged the entirety of eight (plus one partial) Milky Way (MW)–like satellite systems, a total of 42 (45) satellites, from the Satellites Around Galactic Analogs II catalog in both Hα and H i with the Canada–France–Hawaii Telescope and the Jansky Very Large Array. In these eight systems we have identified four cases where a satellite appears to be currently undergoing ram pressure stripping (RPS) as its H i gas collides with the circumgalactic medium (CGM) of its host. We also see a clear suppression of gas fraction (M HI/M *) with decreasing (projected) satellite–host separation—to our knowledge, the first time this has been observed in a sample of MW-like systems. Comparisons to the Auriga, A Project Of Simulating The Local Environment, and TNG50 cosmological zoom-in simulations show consistent global behavior, but they systematically underpredict gas fractions across all satellites by roughly 0.5 dex. Using a simplistic RPS model, we estimate the average peak CGM density that satellites in these systems have encountered to be logρcgm/gcm−3≈−27.3 . Furthermore, we see tentative evidence that these satellites are following a specific star formation rate to gas fraction relation that is distinct from field galaxies. Finally, we detect one new gas-rich satellite in the UGC 903 system with an optical size and surface brightness meeting the standard criteria to be considered an ultra-diffuse galaxy.

Extragalactic Magnetism with SOFIA (SALSA Legacy Program). VII. A Tomographic View of Far-infrared and Radio Polarimetric Observations through MHD Simulations of Galaxies

The structure of magnetic fields in galaxies remains poorly constrained, despite the importance of magnetism in the evolution of galaxies. Radio synchrotron and far-infrared (FIR) polarization and polarimetric observations are the best methods to measure galactic scale properties of magnetic fields in galaxies beyond the Milky Way. We use synthetic polarimetric observations of a simulated galaxy to identify and quantify the regions, scales, and interstellar medium (ISM) phases probed at FIR and radio wavelengths. Our studied suite of magnetohydrodynamical cosmological zoom-in simulations features high-resolutions (10 pc full-cell size) and multiple magnetization models. Our synthetic observations have a striking resemblance to those of observed galaxies. We find that the total and polarized radio emission extends to approximately double the altitude above the galactic disk (half-intensity disk thickness of h I radio ∼ h PI radio = 0.23 ± 0.03 kpc) relative to the total FIR and polarized emission that are concentrated in the disk midplane (h I FIR ∼ h PI FIR = 0.11 ± 0.01 kpc). Radio emission traces magnetic fields at scales of ≳300 pc, whereas FIR emission probes magnetic fields at the smallest scales of our simulations. These scales are comparable to our spatial resolution and well below the spatial resolution (<300 pc) of existing FIR polarimetric measurements. Finally, we confirm that synchrotron emission traces a combination of the warm neutral and cold neutral gas phases, whereas FIR emission follows the densest gas in the cold neutral phase in the simulation. These results are independent of the ISM magnetic field strength. The complementarity we measure between radio and FIR wavelengths motivates future multiwavelength polarimetric observations to advance our knowledge of extragalactic magnetism.

The impact of cosmic rays on the interstellar medium and galactic outflows of Milky Way analogues

Abstract During the last decade, cosmological simulations have managed to reproduce realistic and morphologically diverse galaxies, spanning the Hubble sequence. Central to this success was a phenomenological calibration of the few included feedback processes, whilst glossing over higher complexity baryonic physics. This approach diminishes the predictive power of such simulations, preventing to further our understanding of galaxy formation. To tackle this fundamental issue, we investigate the impact of cosmic rays (CRs) and magnetic fields on the interstellar medium (ISM) and the launching of outflows in a cosmological zoom-in simulation of a Milky Way-like galaxy. We find that including CRs decreases the stellar mass of the galaxy by a factor of 10 at high redshift and ∼4 at cosmic noon, leading to a stellar mass to halo mass ratio in good agreement with abundance matching models. Such decrease is caused by two effects: i) a reduction of cold, high-density, star-forming gas, and ii) a larger fraction of SN events exploding at lower densities, where they have a higher impact. SN-injected CRs produce enhanced, multi-phase galactic outflows, which are accelerated by CR pressure gradients in the circumgalactic medium of the galaxy. While the mass budget of these outflows is dominated by the warm ionised gas, warm neutral and cold gas phases contribute significantly at high redshifts. Importantly, our work shows that future JWST observations of galaxies and their multi-phase outflows across cosmic time have the ability to constrain the role of CRs in regulating star formation.

Dense stellar clump formation driven by strong quasar winds in the FIRE cosmological hydrodynamic simulations

ABSTRACT We investigate the formation of dense stellar clumps in a suite of high-resolution cosmological zoom-in simulations of a massive, star-forming galaxy at z ∼ 2 under the presence of strong quasar winds. Our simulations include multiphase ISM physics from the Feedback In Realistic Environments (FIRE) project and a novel implementation of hyper-refined accretion disc winds. We show that powerful quasar winds can have a global negative impact on galaxy growth while in the strongest cases triggering the formation of an off-centre clump with stellar mass ${\rm M}_{\star }\sim 10^{7}\, {\rm M}_{\odot }$, effective radius ${\rm R}_{\rm 1/2\, \rm Clump}\sim 20\, {\rm pc}$, and surface density $\Sigma _{\star } \sim 10^{4}\, {\rm M}_{\odot }\, {\rm pc}^{-2}$. The clump progenitor gas cloud is originally not star-forming, but strong ram pressure gradients driven by the quasar winds (orders of magnitude stronger than experienced in the absence of winds) lead to rapid compression and subsequent conversion of gas into stars at densities much higher than the average density of star-forming gas. The AGN-triggered star-forming clump reaches ${\rm SFR} \sim 50\, {\rm M}_{\odot }\, {\rm yr}^{-1}$ and $\Sigma _{\rm SFR} \sim 10^{4}\, {\rm M}_{\odot }\, {\rm yr}^{-1}\, {\rm kpc}^{-2}$, converting most of the progenitor gas cloud into stars in ∼2 Myr, significantly faster than its initial free-fall time and with stellar feedback unable to stop star formation. In contrast, the same gas cloud in the absence of quasar winds forms stars over a much longer period of time (∼35 Myr), at lower densities, and losing spatial coherency. The presence of young, ultra-dense, gravitationally bound stellar clumps in recently quenched galaxies could thus indicate local positive feedback acting alongside the strong negative impact of powerful quasar winds, providing a plausible formation scenario for globular clusters.

Zooming in on the circumgalactic medium with GIBLE

We used a cosmological zoom-in simulation of a Milky Way-like galaxy to study and quantify the topology of magnetic field lines around cold gas clouds in the circumgalactic medium (CGM). This simulation is a new addition to Project GIBLE, a suite of cosmological magnetohydrodynamic simulations of galaxy formation with preferential super-Lagrangian refinement in the CGM, reaching an unprecedented CGM gas mass resolution of ∼225 M⊙. To maximize statistics and resolution, we focused on a sample of ∼200 clouds with masses of ∼106 M⊙. The topology of magnetic field lines around clouds is diverse, from threading to draping, and there is large variation in the magnetic curvature (κ) within cloud-background interfaces. We typically find little variation of κ between upstream and downstream cloud faces, implying that strongly draped configurations are rare. In addition, κ correlates strongly with multiple properties of the interface and the ambient background, including cloud overdensity and relative velocity, suggesting that cloud properties impact the topology of interface magnetic fields.

The AGORA High-resolution Galaxy Simulations Comparison Project. V. Satellite Galaxy Populations in a Cosmological Zoom-in Simulation of a Milky Way–Mass Halo

We analyze and compare the satellite halo populations at z ∼ 2 in the high-resolution cosmological zoom-in simulations of a 1012 M ⊙ target halo (z = 0 mass) carried out on eight widely used astrophysical simulation codes (Art-I, Enzo, Ramses, Changa, Gadget-3, Gear, Arepo-t, and Gizmo) for the AGORA High-resolution Galaxy Simulations Comparison Project. We use slightly different redshift epochs near z = 2 for each code (hereafter “z ∼ 2”) at which the eight simulations are in the same stage in the target halo’s merger history. After identifying the matched pairs of halos between the CosmoRun simulations and the DMO simulations, we discover that each CosmoRun halo tends to be less massive than its DMO counterpart. When we consider only the halos containing stellar particles at z ∼ 2, the number of satellite galaxies is significantly fewer than that of dark matter halos in all participating AGORA simulations and is comparable to the number of present-day satellites near the Milky Way or M31. The so-called “missing satellite problem” is fully resolved across all participating codes simply by implementing the common baryonic physics adopted in AGORA and the stellar feedback prescription commonly used in each code, with sufficient numerical resolution (≲100 proper pc at z = 2). We also compare other properties such as the stellar mass–halo mass relation and the mass–metallicity relation. Our work highlights the value of comparison studies such as AGORA, where outstanding problems in galaxy formation theory are studied simultaneously on multiple numerical platforms.

Cosmic evolution of black hole spin and galaxy orientations: Clues from the NewHorizon and Galactica simulations

Black holes (BHs) are ubiquitous components of the center of most galaxies. In addition to their mass, the BH spin, through its amplitude and orientation, is a key factor in the galaxy formation process, as it controls the radiative efficiency of the accretion disk and relativistic jets. Using the recent cosmological high-resolution zoom-in simulations and in which the evolution of the BH spin is followed on the fly, we have tracked the cosmic history of a hundred BHs with a mass greater than $2 For each of them, we have studied the variations of the three-dimensional angle (Psi ) subtended between the BH spins and the angular momentum vectors of their host galaxies (estimated from the stellar component). The analysis of the individual evolution of the most massive BHs suggests that they are generally passing by three different regimes. First, for a short period after their birth, low-mass BHs BH &lt;$3times 10$^4$M$_ are rapidly spun up by gas accretion and their spin tends to be aligned with their host galaxy spin. Then follows a second phase in which the accretion of gas onto low-mass BHs BH is quite chaotic and inefficient, reflecting the complex and disturbed morphologies of forming proto-galaxies at high redshifts. The variations of Psi are rather erratic during this phase and are mainly driven by the rapid changes of the direction of the galaxy angular momentum. Then in a third and long phase, BHs are generally well settled in the center of galaxies around which the gas accretion becomes much more coherent BH &gt;$10$^5$ M$_ In this case, the BH spins tend to be well aligned with the angular momentum of their host galaxy and this configuration is generally stable even though BH merger episodes can temporally induce misalignment. We even find a few cases of BH-galaxy spin anti-alignment that lasts for a long time in which the gas component is counter-rotating with respect to the stellar component. We have also derived the distributions of cos(Psi ) at different redshifts and found that BHs and galaxy spins are generally aligned. Our analysis suggests that the fraction of BH-galaxy pairs with low Psi values reaches maximum at zsim 4-3 and then decreases until zsim 1.5 due to the high BH-merger rate. Afterward, it remains almost constant probably due to the fact that BH mergers becomes rare, except for a slight increase at late times. Finally, based on a Monte Carlo method, we also predict statistics for the 2-d projected spin-orbit angles lambda . In particular, the distribution of lambda traces the alignment tendency well in the three-dimensional analysis. Such predictions provide an interesting background for future observational analyses.