Welcome to the new version of CaltechAUTHORS. Login is currently restricted to library staff. If you notice any issues, please email coda@library.caltech.edu
Published August 20, 2021 | Submitted + Published
Journal Article Open

AGORA High-resolution Galaxy Simulations Comparison Project. III. Cosmological Zoom-in Simulation of a Milky Way–mass Halo

Abstract

We present a suite of high-resolution cosmological zoom-in simulations to z = 4 of a 10¹² M_⊙ halo at z = 0, obtained using seven contemporary astrophysical simulation codes (Art-I, Enzo, Ramses, Changa, Gadget-3, Gear, and Gizmo) widely used in the numerical galaxy formation community. The physics prescriptions for gas cooling and heating and star formation are the same as the ones used in our previous Assembling Galaxies of Resolved Anatomy (AGORA) disk comparison but now account for the effects of cosmological processes such as the expansion of the universe, intergalactic gas inflow, and the cosmic ultraviolet background radiation emitted by massive stars and quasars. In this work, we introduce the most careful comparison yet of galaxy formation simulations run by different code groups, together with a series of four calibration steps each of which is designed to reduce the number of tunable simulation parameters adopted in the final run. In the first two steps, we methodically calibrate the gas physics, such as cooling and heating, in simulations without star formation. In the third step, we seek agreement on the total stellar mass produced with the common star formation prescription used in the AGORA disk comparison, in stellar-feedback-free simulations. In the last calibration step, we activate stellar feedback, where each code group is asked to set the feedback prescription to as close to the most widely used one in its code community as possible, while aiming for convergence in the stellar mass at z = 4 to the values predicted by semiempirical models. After all the participating code groups successfully complete the calibration steps, we achieve a suite of cosmological simulations with similar mass assembly histories down to z = 4. With numerical accuracy that resolves the internal structure of a target halo (≲100 physical pc at z = 4), we find that the codes overall agree well with one another, e.g., in gas and stellar properties, but also show differences, e.g., in circumgalactic medium (CGM) properties. We argue that, if adequately tested in accordance with our proposed calibration steps and common parameters, high-resolution cosmological zoom-in simulations can have robust and reproducible results. New code groups are invited to join and enrich this comparison by generating equivalent models or to test the code's compatibility on their own, by adopting the common initial conditions, the common easy-to-implement physics package, and the proposed calibration steps. Further analyses of the zoom-in simulations presented here will be presented in forthcoming reports from the AGORA Collaboration, including studies of the CGM, simulations by additional codes, and results at lower redshift.

Additional Information

© 2021. The American Astronomical Society. Received 2021 March 30; revised 2021 June 2; accepted 2021 June 4; published 2021 August 17. We thank all of our colleagues who participate in the AGORA Project for their collaborative spirit, which has allowed the AGORA Collaboration to remain strong as a platform to foster and launch multiple science-oriented comparison efforts. We thank Aldo Rodríguez-Puebla for sharing results from abundance matching semiempirical models, and Volker Springel for providing the original versions of Gadget-3 to be used in the AGORA Project. We also thank the anonymous referee for insightful comments and suggestions. This research used resources of the National Energy Research Scientific Computing Center, a user facility supported by the Office of Science of the U.S. Department of Energy under contract No. DE-AC02-05CH11231. S.R.-F. acknowledges support from a Spanish postdoctoral fellowship, under grant No. 2017-T2/TIC-5592. His work has been supported by the Madrid Government (Comunidad de Madrid–Spain) under the Multiannual Agreement with Complutense University in the line Program to Stimulate Research for Young Doctors in the context of V PRICIT (Regional Programme of Research and Technological Innovation). He also acknowledges financial support from the Spanish Ministry of Economy and Competitiveness under grant Nos. AYA2016-75808-R, AYA2017-90589-REDT, and S2018/NMT-429, and from CAM-UCM under grant No. PR65/19-22462. J.K. acknowledges support from the Samsung Science and Technology Foundation under project No. SSTF-BA1802-04. His work was also supported by the National Institute of Supercomputing and Networking, Korea Institute of Science and Technology Information, with supercomputing resources, including technical support (grants KSC-2018-CRE-0052 and KSC-2019-CRE-0163). K.N. acknowledges support from MEXT/JSPS KAKENHI grant Nos. JP17H01111, 19H05810, and 20H00180, as well as travel support from Kavli IPMU, World Premier Research Center Initiative, where part of this work was conducted. A.L. acknowledges funding by the MIUR under the grant PRIN 2017-MB8AEZ. D.C. is a Ramon Cajal Researcher and is supported by Ministerio de Ciencia, Innovación y Universidades (FEDER) under research grant PGC2018-094975-C21. H.V. acknowledges support from PAPIIT of Universidad Nacional Autónoma de México (UNAM) under grant No. IN101918 and also from Centro Nacional de Supercomputo (CNS-IPICYT-CONACYT). The Art-I simulations were performed on the Brigit/Eolo cluster at Centro de Proceso de Datos, Universidad Complutense de Madrid, and on the Stócatl supercomputer at Instituto de Astronomía de la UNAM. The Ramses simulations were performed on the Miztli supercomputer at LANCAD, UNAM, within the research project LANCAD-UNAM-DGTIC-151 and on Laboratorio Nacional de Supercómputo del Sureste, CONACYT. The Changa simulations were performed on the Atócatl supercomputer at Instituto de Astronomía de la UNAM, and on the Extreme Science and Engineering Discovery Environment (XSEDE) allocations TG-AST20020 and TG-MCA94P018. XSEDE is supported by the National Science Foundation grant ACI-1053575. The Gadget3-Osaka simulations and analyses were performed on the XC50 systems at the Center for Computational Astrophysics of the National Astronomical Observatory of Japan, on Octopus at the Cybermedia Center of Osaka University, and on Oakforest-PACS at the University of Tokyo as part of the HPCI System Research Project (hp190050 and hp200041). The publicly available Enzo and yt codes used in this work are the products of collaborative efforts by many independent scientists from numerous institutions around the world. Their commitment to open science has helped make this work possible.

Attached Files

Published - Roca-Fabrega_2021_ApJ_917_64.pdf

Submitted - 2106-09738.pdf

Files

2106-09738.pdf
Files (35.9 MB)
Name Size Download all
md5:eb6168b6d9ddfd1e41e623c047f3bfea
28.7 MB Preview Download
md5:41f346dec5fe8d4ca3d2f65e2ce71b39
7.2 MB Preview Download

Additional details

Created:
October 4, 2023
Modified:
October 24, 2023