Search

The Ly$\alpha$ forest, a series of HI absorption lines in the quasar spectra, is a powerful tool for probing the large-scale structure of the intergalactic medium. Its three-dimensional (3D) correlation and cross-correlations with quasars allow precise measurements of the baryon acoustic oscillation feature and redshift space distortions at redshifts $z>2$. Understanding small-scale astrophysical phenomena, such as star formation and feedback, is crucial for full-shape analyses. In this study, we measure the 3D auto-power spectrum of the Ly$\alpha$ forest and its cross-power spectrum with halos using hydrodynamic simulations from the GADGET3-OSAKA code, which includes models for star formation and supernova feedback. Across five astrophysical models, we find significant deviations from the Fiducial model, with $5-10\,\%$ differences for wavenumbers $k>2\,h\mathrm{Mpc}^{-1}$ in the Ly$\alpha$ auto-power spectrum. The Ly$\alpha\,\times\,$halo cross-power spectra show even larger deviations, exceeding $10\,\%$ in some cases. Using the fitting models of Arinyo-i-Prats et al. (2015) and Givans et al. (2022), we jointly fit the Ly$\alpha$ auto- and Ly$\alpha$ $\times$ halo cross-power spectra, and assess the accuracy of the estimated $f\sigma_8$ parameter by comparing it with the ground truth from the simulations, while varying the maximum wavenumber $k_\mathrm{max}$ and minimum halo mass $M_h$. Our results demonstrate that the extended model of Givans et al. (2022) is highly effective in reproducing $f\sigma_8$ at $k_\mathrm{max}\leq3.0\,h\mathrm{Mpc}^{-1}$ for $M_h>10^{10.5} M_\odot$, and remains robust against astrophysical uncertainties., Comment: 14 pages, 8 figures

The AGORA Cosmorun (arXiv:2106.09738) is a set of hydrodynamical cosmological zoom-in simulations carried out within the AGORA High-resolution Galaxy Simulations Comparison Project (arXiv:1308.2669,arXiv:1610.03066). These simulations show the formation and evolution of a Milky Way-sized galaxy using eight of the most widely used numerical codes in the community (Art-I, Enzo, Ramses, Changa, Gadget-3, Gear, Gizmo, and Arepo). In this short report, we describe the public release of the raw output data from all of these simulations at z = 8, 7, 6, 5, 4, 3, 2 (plus at z=1, 0 when available), and several metadata files containing the halo centers, virial quantities, and merger trees. The data from even thinner timesteps will be released as soon as the upcoming collaboration papers (VII-IX) are submitted and accepted., Comment: ArXiv Data Release from the AGORA Collaboration. arXiv admin note: text overlap with arXiv:2001.04354

We apply the convolutional neural networks (CNNs) to the mock 21cm maps from the post-reionization epoch to show that the $\Lambda$ cold dark matter and warm dark matter (WDM) model can be distinguished for WDM particle masses $m_{FD}<3$\,keV, under the assumption of thermal production of WDM following the Fermi-Dirac (FD) distribution. We demonstrate that the CNN is a potent tool in distinguishing the dark matter masses, highlighting its sensitivity to the subtle differences in the 21cm maps produced by varying dark matter masses. Furthermore, we extend our analysis to encompass different WDM production mechanisms, recognizing that the dark matter production mechanism in the early Universe is among the sources of the most significant uncertainty for the dark matter model building. In this work, given the mass of the dark matter, we discuss the feasibility of discriminating four different WDM models: Fermi-Dirac (FD) distribution model, neutrino minimal Standard Model ($\nu$MSM), Dodelson-Widrow (DW), and Shi-Fuller (SF) model. For instance, when the WDM mass is 2\,keV, we show that one can differentiate between CDM, FD, $\nu$MSM, and DW models while discerning between the DW and SF models turns out to be challenging. Our results reinforce the viability of the CNN as a robust analysis for 21cm maps and shed light on its potential to unravel the features associated with different dark matter production mechanisms., Comment: 16 pages, 13 figures, published in PRD

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 non-axisymmetric structures in the stellar disks. Codes generally converge well to the stellar-to-halo mass ratios predicted by semi-analytic models at $z\sim$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 MZR at high and intermediate redshifts, and reproduce the gas metallicity gradients obtained from analytical models and low-redshift observations. We confirm that the inter-code 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., Comment: Submitted to ApJ. Visit the AGORA Collaboration website (www.agorasimulations.org ) for more information

We analyze and compare the satellite halo populations at $z\sim2$ in the high-resolution cosmological zoom-in simulations of a $10^{12}\,{\rm M}_{\odot}$ target halo ($z=0$ mass) carried out on eight widely-used astrophysical simulation codes ({\sc Art-I}, {\sc Enzo}, {\sc Ramses}, {\sc Changa}, {\sc Gadget-3}, {\sc Gear}, {\sc Arepo-t}, and {\sc Gizmo}) for the {\it AGORA} High-resolution Galaxy Simulations Comparison Project. We use slightly different redshift epochs near $z=2$ for each code (hereafter ``$z\sim2$') 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 {\it CosmoRun} simulations and the DMO simulations, we discover that each {\it CosmoRun} halo tends to be less massive than its DMO counterpart. When we consider only the halos containing stellar particles at $z\sim2$, the number of satellite {\it galaxies} is significantly fewer than that of dark matter halos in all participating {\it 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 {\it AGORA} and the stellar feedback prescription commonly used in each code, with sufficient numerical resolution ($\lesssim100$ 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 {\it AGORA}, where outstanding problems in galaxy formation theory are studied simultaneously on multiple numerical platforms., Comment: Accepted for publication to the ApJ, 19 pages, 10 figures. Visit the AGORA Collaboration website (www.agorasimulations.org) for more information. For summary video, please see https://www.youtube.com/watch?v=HShRJYvPe1k

We analyze the circumgalactic medium (CGM) for eight commonly-used cosmological codes in the AGORA collaboration. The codes are calibrated to use identical initial conditions, cosmology, heating and cooling, and star formation thresholds, but each evolves with its own unique code architecture and stellar feedback implementation. Here, we analyze the results of these simulations in terms of the structure, composition, and phase dynamics of the CGM. We show properties such as metal distribution, ionization levels, and kinematics are effective tracers of the effects of the different code feedback and implementation methods, and as such they can be highly divergent between simulations. This is merely a fiducial set of models, against which we will in the future compare multiple feedback recipes for each code. Nevertheless, we find that the large parameter space these simulations establish can help disentangle the different variables that affect observable quantities in the CGM, e.g., showing that abundances for ions with higher ionization energy are more strongly determined by the simulation's metallicity, while abundances for ions with lower ionization energy are more strongly determined by the gas density and temperature., Comment: 32 pages, 17 figures. Published in ApJ