The metallicity gradients of star-forming regions store information of the assembly history of galaxies

Abstract: The variations in metallicity and spatial patterns within star-forming regions of galaxies result from diverse physical processes unfolding throughout their evolutionary history, with a particular emphasis in recent events. Analysing MaNGA and \textsc{eagle} galaxies, we discovered an additional dependence of the mass-metallicity relation (MZR) on metallicity gradients ($\nabla_{{\rm (O/H)}}$). Two regimes emerged for low and high stellar mass galaxies, distinctly separated at approximately ${\rm M_{\star}} >10^{9.75}$. Low-mass galaxies with strong positive $\nabla_{{\rm (O/H)}}$ appear less enriched than the MZR median, while those with strong negative gradients are consistently more enriched in both simulated and observed samples. Interestingly, low-mass galaxies with strong negative $\nabla_{{\rm (O/H)}}$ exhibit high star-forming activity, regardless of stellar surface density or $\nabla_{{\rm (O/H)}}$. In contrast, a discrepancy arises for massive galaxies between MaNGA and \textsc{eagle} datasets. The latter exhibit a notable anticorrelation between specific star formation rate and stellar surface density, independent of $\nabla_{{\rm (O/H)}}$, while MaNGA galaxies show this trend mainly for strong positive $\nabla_{{\rm (O/H)}}$. Further investigation indicates that galaxies with strong negative gradients tend to host smaller central black holes in observed datasets, a trend not replicated in simulations. These findings suggest disparities in metallicity recycling and mixing history between observations and simulations, particularly in massive galaxies with varying metallicity gradients. These distinctions could contribute to a more comprehensive understanding of the underlying physics.