We thank the referee for constructive criticism that improved the quality of the paper. MR thanks Astha Agarwal, Harley Katz, Taysun Kimm, Joakim Rosdahl, and Mahsa Sanati for useful discussions and comments during the preparation of this manuscript. MR acknowledges support from the Beecroft Fellowship funded by Adrian Beecroft. ET acknowledges the UKRI Science and Technology Facilities Council (STFC) for support (grant ST/V50712X/1). OA acknowledges support from the Knut and Alice Wallenberg Foundation, the Swedish Research Council (grant 2019-04659), the Swedish National Space Agency (SNSA Dnr 2023-00164), and the LMK Foundation. JIR would like to thank the STFC for support from grants ST/Y002865/1 and ST/Y002857/1. This project has received funding from the European Union's Horizon 2020 research and innovation programme under grant agreement no. 818085 GMGalaxies. This work was performed using the DiRAC Data Intensive service at Leicester, operated by the University of Leicester IT Services, which forms part of the STFC DiRAC HPC Facility (www.dirac.ac.uk). The equipment was funded by BEIS capital funding via STFC capital grants ST/K000373/1 and ST/R002363/1 and STFC DiRAC Operations grant ST/R001014/1. DiRAC is part of the National e-Infrastructure. The authors acknowledge the use of the UCL Grace High Performance Computing Facility, the Surrey Eureka supercomputer facility, and their associated support services. This work was partially supported by the UCL Cosmoparticle Initiative.