2023MNRAS.524.3252H

The Lilly-Madau plot is commonly interpreted as the history of the cosmic star formation of the Universe by showing the co-moving star formation rate density (SFRD) over cosmic time. Therefore, the Lilly-Madau plot is sensitive not only to the star formation history (SFH) but also to the number density of galaxies. Assessing the Catalogue of Neighbouring Galaxies, we reconstruct the SFHs of galaxies located in the Local Cosmological Volume (LV) based on delayed-τ and power-law SFH models. Galaxies with stellar masses of $M_{*} \gtrsim 10^{10}\, {\rm M_{\odot }}$ typically evolve according to the delayed-τ model by having first increasing followed by exponentially declining SFRs, while the majority of less massive star-forming galaxies have an almost constant or increasing SFH. Deducing the cosmic SFRD evolution of the LV reveals that the SFHs of local galaxies are inconsistent with the Lilly-Madau plot. The SFRDs of the LV are significantly lower at redshifts of $z \lesssim 3$ underestimating the Lilly-Madau peak at z = 1.86 by a factor of 2.16 ± 0.32 (delayed-τ model) and 5.90 ± 0.88 (power-law model). Assuming the delayed-τ model for galaxies with $M_{*} \ge 10^{10}\, {\rm M_{\odot }}$ and a power-law model for less massive galaxies, the SFRD is 2.22 ± 0.33 lower than measured at z = 1.86. This inconsistency between the evolution of the local and global SFRD has cosmological implications. Since the Lilly-Madau plot also constrains the cosmological matter field, the near-constancy of SFHs of LV galaxies could imply that the peak of the Lilly-Madau plot at z = 1.86 is the imprint of an ≈5 co-moving Gpc-scale inhomogeneity.