POSYDON Version 2: Population Synthesis with Detailed Binary-Evolution Simulations across a Cosmological Range of Metallicities

Whether considering rare astrophysical events on cosmological scales or unresolved stellar populations, accurate models must account for the integrated contribution from the entire history of star formation upon which that population is built. Here, we describe the second version of POSYDON, an open-source binary population synthesis code based on extensive grids of detailed binary evolution models computed using the MESA code, which follows both stars' structures as a binary system evolves through its complete evolution from the zero-age main sequence, through multiple phases of mass transfer and supernovae, to their death as compact objects. To generate synthetic binary populations, POSYDON uses advanced methods to interpolate between our large, densely spaced grids of simulated binaries. In our updated version of POSYDON, we account for the evolution of stellar binaries across a cosmological range of metallicities, extending from $10^{-4}\,Z_{\odot}$ to $2\,Z_{\odot}$, including grids specifically focused on the Small and Large Magellanic Clouds ($0.2\,Z_{\odot}$ and $0.45\,Z_{\odot}$). In addition to describing our model grids and detailing our methodology, we outline several improvements to POSYDON. These include the incorporation of single stars in stellar populations, a treatment for stellar mergers, and a careful modeling of "reverse-mass transferring" binaries, in which an once-accreting star later becomes a donor star. Our simulations are focused on binaries with at least one high-mass component, such as those that host neutron stars and black holes, and we provide post-processing methods to account for the cosmological evolution of metallicity and star formation as well as rate calculations for gravitational wave events, gamma-ray bursts, and other transients.