Abstract

Context. The distribution of chemical elements in star-forming regions can store information on the chemical enrichment history of galaxies and particularly of recent events. Negative metallicity gradients are expected in galaxies forming inside-out. Azimuthal-averaged profiles are usually fit to the projected chemical distributions to quantify them. However, observations show that the metallicity profiles can be broken. Aims. We aim to study the diversity of metallicity profiles that can arise in the current cosmological context and compare them with available observations. Additionally, we seek to identify the physical processes responsible for breaks in metallicity profiles by using two galaxies as case studies. Methods. We analyzed central galaxies from the cosmological simulations of CIELO project, with stellar masses within the range of 108.5 to 1010.5 M? at z=0. A new algorithm, DB-A, was developed to fit multiple power laws to the metallicity profiles, enabling a flexible assessment of metallicity gradients in various galactic regions. The simulations include detailed modeling of gas components, metal-dependent cooling, star formation, and supernova feedback. Results. At z=0, we find a diversity of shapes, with inner and outer drops and rises, and there are a few galaxies with double breaks. Inner, outer, and middle gradients are in agreement with observations. We also find that using a single linear regression to fit gradients usually traces the middle gradient well. A detailed temporal analysis of the main galaxies of a Local Group analog revealed the occurrence of inner and outer breaks at all cosmic times, with the latter being the most common feature during the evolution of our case studies. Significant variability in the metallicity gradients was found at high redshift, transitioning to more gradual evolution at lower redshifts. Most of the inner breaks have enhanced oxygen abundances in the center, which are linked to gas accretion followed by efficient star formation. Inner breaks with diluted oxygen abundances are less common and are found in galaxies with disrupted gas distributions which are affected by feedback-driven ejection of enriched gas. Outer breaks with high abundances are linked to processes such as the re-accretion of enriched material, extended star formation, and enhanced gas mixing from the circumgalactic medium. Outer breaks with diluted metallicities in the outskirts are found mainly at high redshift and are associated with the accretion of metal-poor gas from cold flows. We also highlight and illustrate the complex interplay of these processes which act often together. © The Authors 2025.