Local variations of the radial metallicity gradient in a simulated NIHAO-UHD Milky Way analogue and their implications for (extra-)galactic studies

Radial metallicity gradients are fundamental to understanding galaxy formation and evolution. In our high-resolution simulation of a NIHAO-UHD Milky Way analogue, we analyze the linearity, scatter, spatial coherence, and age-related variations of metallicity gradients using young stars and gas. While a global linear model generally captures the gradient, it ever so slightly overestimates metallicity in the inner galaxy and underestimates it in the outer regions of our simulated galaxy. Both a quadratic model, showing an initially steeper gradient that smoothly flattens outward, and a piecewise linear model with a break radius around 9.3-11.5~kpc (2.4-3.0 effective radii) fit the data equally better. The spread of [Fe/H] of young stars in the simulation increases by tenfold from the innermost to the outer galaxy at a radius of 20 kpc. We find that stars born at similar times along radial spirals drive this spread in the outer galaxy, with a chemical under- and over-enhancement of up to 0.1 dex at leading and trailing regions of such spirals, respectively. This localised chemical variance highlights the need to examine radial and azimuthal selection effects for both Galactic and extragalactic observational studies. The arguably idealised but volume-complete simulations suggest that future studies should not only test linear and piecewise linear gradients, but also non-linear functions such as quadratic ones to test for a smooth gradient rather than one with a break radius. Either finding would help to determine the importance of different enrichment or mixing pathways and thus our understanding of galaxy formation and evolution scenarios.