We study the effect of the gas accretion rate ($\dot M_{\rm accr}$) on the radial gas metallicity profile (RMP) of galaxies using the EAGLE cosmological hydrodynamic simulations, focusing on central galaxies of stellar mass $M_\star \gtrsim 10^9 \, {\rm M_\odot}$ at $z \le 1$. We find clear relations between $\dot M_{\rm accr}$ and the slope of the RMP (measured within an effective radius), where higher $\dot M_{\rm accr}$ are associated with more negative slopes. The slope of the RMPs depends more strongly on $\dot M_{\rm accr}$ than on stellar mass, star formation rate or gas fraction, suggesting $\dot M_{\rm accr}$ to be a more fundamental driver of the RMP slope of galaxies. We find that eliminating the dependence on stellar mass is essential for pinning down the properties that shape the slope of the RMP. Although $\dot M_{\rm accr}$ is the main property modulating the slope of the RMP, we find that it causes other correlations that are more easily testable observationally: at fixed stellar mass, galaxies with more negative RMP slopes tend to have higher gas fractions and SFRs, while galaxies with lower gas fractions and SFRs tend to have flatter metallicity profiles within an effective radius.

We present results from the semi-analytic model of galaxy formation SAG applied on the MultiDark simulation MDPL2. SAG features an updated supernova (SN) feedback scheme and a robust modelling of the environmental effects on satellite galaxies. This incorporates a gradual starvation of the hot gas halo driven by the action of ram pressure stripping (RPS), that can affect the cold gas disc, and tidal stripping (TS), which can act on all baryonic components. Galaxy orbits of orphan satellites are integrated providing adequate positions and velocities for the estimation of RPS and TS. The star formation history and stellar mass assembly of galaxies are sensitive to the redshift dependence implemented in the SN feedback model. We discuss a variant of our model that allows to reconcile the predicted star formation rate density at $z \gtrsim 3$ with the observed one, at the expense of an excess in the faint end of the stellar mass function at $z=2$. The fractions of passive galaxies as a function of stellar mass, halo mass and the halo-centric distances are consistent with observational measurements. The model also reproduces the evolution of the main sequence of star forming central and satellite galaxies. The similarity between them is a result of the gradual starvation of the hot gas halo suffered by satellites, in which RPS plays a dominant role. RPS of the cold gas does not affect the fraction of quenched satellites but it contributes to reach the right atomic hydrogen gas content for more massive satellites ($M_{\star}\gtrsim 10^{10}\,{\rm M}_{\odot}$).

We use the semi-analytic model (SAM) of galaxy formation and evolution SAG coupled with the MULTIDARK simulation MDPL2 to study the evolution of the stellar mass-gas metallicity relation of galaxies (MZR). We test several implementations of the dependence of the mass loading due to supernovae (SNe). We find that no evolution in the normalization of the MZR is obtained unless we introduce an explicit scaling of the reheated and ejected mass with redshift as $(1+z)^\beta$. The latter is in agreement with results from the FIRE simulations, and it should encompass small scale properties of the interstellar medium varying over time, which are not captured in SAMs, as well as other energy sources in addition to SNe. Increasing $\beta$ leads to stronger evolution of the MZR normalization; $\beta = 1.9$ reproduces the observed MZR in the range $0 < z < 3.5$. A stronger redshift dependence of outflows reduces the levels of star formation at earlier epochs with the consequent decrease of metal production. This leads to a slower increase of the gas metallicity compared to the stellar mass build-up. The cold gas can be contaminated either by receiving a direct injection of the material recycled by stellar winds and SNe or by gas cooling. The relative role of each process for a given stellar mass depends on the criterion adopted to regulate the fate of the recycled material. However, modifying the metal loading of the outflows has mild impact on the zero-point evolution and does not affect our conclusions.