The gas, metal and dust evolution in low-metallicity local and high-redshift galaxies

The chemical enrichment in the interstellar medium (ISM) of galaxies is regulated by several physical processes: stellar evolution, grain formation and destruction, galactic inflows and outflows. Understanding such processes is essential to follow the chemical enrichment of galaxies through the cosmic epochs, and to interpret the observations. Despite the importance of such topics, the efficiency of the different processes driving the evolution of baryons in galaxies, remain controversial. We revise the current description of metal and dust evolution in local low-metallicity dwarf galaxies and we develop a description for Lyman Break Galaxies. Our main goal is to reproduce i) the peak in the mass of dust over the mass of stars (sMdust) observed within few hundred Myrs; ii) the decrease of the sMdust at later time. The spectral energy distribution of the galaxies is fitted with the"Code Investigating GALaxies Emission"(CIGALE), through which the stellar and dust masses, and the star formation rate are estimated. For some of the dwarf galaxies, the metal and gas content are also available. We run different calculations of chemical evolution in galaxies, and we fit the observed properties through the model predictions. We show that i) a top-heavy initial mass function that favours massive stars and a dust condensation fraction for Type II Supernovae (SNe II) of 50% or more help to reproduce the peak of sMdust observed after 100 Myrs since the beginning of the cycle; ii) galactic outflows play a crucial role in reproducing the decline in sMdust with age, and they are more efficient than grain destruction from SNe II; iii) a star formation efficiency (mass of gas converted into stars) of few per cent is required to explain the metallicity of local dwarf galaxies; iv) dust growth in the ISM is not necessary to reproduce the sMdust and, if present, its effect is erased by galactic outflows.


Introduction
Understanding the metals and dust cycle in the interstellar medium (ISM) of galaxies is essential to study their properties and evolution, and to interpret different observations of local and very distant galaxies. Indeed, dust shapes the spectra of galaxies by absorbing the stellar radiation in the ultraviolet (UV) and visible wavelengths and re-emitting it in the infrared (IR) bands, while the metals in the gas phase are responsible for various emission and absorption lines. Different works have been devoted to the study of the dust content in galaxies in the Local Universe (Lisenfeld & Ferrara 1998;Dwek 1998;Hirashita 1999;da Cunha et al. 2010;Zhukovska 2014;Schneider et al. 2016;Rémy-Ruyer et al. 2015;Popping et al. 2017;Gioannini et al. 2017;De Vis et al. 2017;McKinnon et al. 2018;Ginolfi et al. 2018;De Vis et al. 2019) as well as in a variety of high-redshift galaxies (Graziani et al. 2019;Burgarella et al. 2020), including quasars at z > 6 ( Valiante et al. 2014;Calura et al. 2014;Mancini et al. 2015;Aoyama et al. 2017;Wang et al. 2017) and sub-millimiter galaxies (Rowlands et al. 2014;Cai et al. 2020).
Dust is condensed in Type II supernovae (SNe II) remnants and in the dense outflows of low-mass stars at the end of their evolution, during the thermally pulsing asymptotic giant branch (TP-AGB) phase. Massive stars evolve on short timescales (< 30 Myrs), while AGB stars need a longer time in order to start to enrich the ISM. Furthermore, massive and low-mass stars produce different types of elements and dust. Massive stars are factories of oxygen together with silicate and dust (Limongi & Chieffi 2018;Marassi et al. 2019), while low-mass stars can release large amount of carbon already at very low metallicities (Nanni et al. 2013(Nanni et al. , 2014. SN Ia, originating from exploding Article number, page 1 of 28 arXiv:2006.15146v1 [astro-ph.GA] 26 Jun 2020 A&A proofs: manuscript no. aanda white dwarfs in binary systems, enrich the ISM predominantly with iron on a timescale comparable with the evolution of low mass stars that, contrary to massive stars, end their evolution as white dwarfs. The role of SNe as dust producers is still debated. On the one hand, it has been argued that shock waves propagating after their explosion heat up the gas and can severely reduce the amount of dust condensed (Nozawa et al. 2011;Bocchio et al. 2016;Bianchi & Schneider 2007). On the other hand, for SNe II, Gall et al. (2014) have found evidence of the presence of 0.1-0.5 M of large grains in SN2010jl, able to survive to the passage of the reverse shock. In addition, recent observations from Matsuura et al. (2019) of SN1987 have shown that dust might re-condense in the cooling gas that experienced the passage of the forward shock or survive to the passage of the reverse shock. In the former case, this might indicate that dust can be efficiently reformed after being destroyed by shocks. The presence of dust has been also suggested for Type Ia SNe (Gomez et al. 2012).
Around TP-AGB stars, dust can be condensed in considerable amount (Ventura et al. 2012;Nanni et al. 2013Nanni et al. , 2014, but since they evolve on longer time-scales their contribution may be irrelevant at early epochs. In different investigations, dust growth in the ISM of galaxies is needed for explaining the large mass of dust observed that may be difficult to reproduce from stellar sources only (Asano et al. 2013;Zhukovska 2014;Michałowski 2015;Mancini et al. 2015;Leśniewska & Michałowski 2019). However, dust accretion in the ISM may be difficult to explain from the micro physical point of view (Ferrara et al. 2016). Icy mantles can contribute in depleting elements from the gas phase, but they are expected to evaporate in presence of a strong radiation field (Ferrara et al. 2016). Furthermore, the probability that two silicon-bearing molecules can encounter in the grain icy mantle is very low (Ceccarelli et al. 2018). Some laboratory experiments showed that dust could accrete in low-temperature environments by addition of molecules on the surface of dust grains (Rouillé et al. 2015). This process might be prevented in space by the formation of icy mantles on the grain surface (Ceccarelli et al. 2018), and grains may grow when the icy mantle evaporates (Rouillé et al. 2020). Theoretical studies and experiments also suggest charged grains and molecules could favour the process of grain growth (Hollenstein 2000;Bleecker et al. 2006). In the astrophysical context, such a possibility has been investigated in the conditions typical of the Milky Way (Zhukovska et al. 2016) and further studies are needed. Grain growth also appears to be more efficient in galaxies which large stellar masses and related and to their environment (e.g. more efficient grain growth fro larger gas densities) (Mancini et al. 2015;Schneider et al. 2016;Popping et al. 2017;Graziani et al. 2019), or enhanced if the ISM is turbulent (Mattsson et al. 2019;Mattsson 2020). Because of these controversial results, it is important to critically examine dust growth in comparison with observational data.
From the Atacama Large Millimiter Array (ALMA) detection, in Burgarella et al. (2020) we have proposed a template for the dust emission of Lyman Break Galaxies (LBGs) in the end of reionization that might be the first dust grains in the Universe. Along with rest-frame UV and optical data, this template is used to fit and model the multi-wavelength emission of these objects. From the derived physical parameters, we build a chemical evolutionary model that explains the (sub-)mm detection (or non-detection) of high-redshift LBGs in the end of reionization without any dust growth in the ISM. This model successfully explains different diagnostic diagrams (e.g. the evolution of the dust content in these galaxies). The non-detection is explained by dust destruction from shocks produced by SNe plus removal in the circum-and intergalactic media by outflows.
In this paper, we interpret the observations of low-metallicity dwarf galaxies (Rémy-Ruyer et al. 2015) in the Local Universe from the Dwarf Galaxy Survey (DGS Madden et al. 2013), together with the LBGs which share similar properties in terms of star formation activity, metallicity and stellar mass. The spectral energy distribution (SED) of these galaxies is fitted as discussed in (Burgarella et al. 2020). In particular, we aim to simultaneously reproduce the available constraints for such galaxies: total dust mass, star formation rate (SFR), metallicity, dust-to-gas ratio, gas fraction, age and, for few local galaxies, circumgalactic dust. By doing this, we identify the critical processes that drive the evolution of baryons in these galaxies: star formation efficiency, metal and dust enrichment from stellar sources (SNe II, TP-AGB stars and SNe Ia), dust destruction from SNe, galactic outflows, and the role of dust accretion in the ISM.
The paper is organised as follows. In Section 2 we describe the two samples of galaxies. In Section 3 we explain the prescriptions adopted for following the gas and dust evolution in the ISM of galaxies. In Section 4 we present the method for fitting the SED of galaxies and for comparing the model predictions with the properties derived. In Section 5 the results are provided, while Section 6 is devoted to discuss the results of this work in comparison with the ones in the literature. Our conclusions are provided in Section 7.

Selected sample
In this work, we analyse two data samples of DGS galaxies and LBGs for which the SED fitting has been performed in Burgarella et al. (2020) by means of the code cigale (Boquien et al. 2019). In cigale the synthetic spectra are calculated from the energy balance between the photons absorbed and re-emitted by dust. Emission and absorption lines of gas are also included. The physical properties of the galaxies are estimated through the analysis of the likelihood distribution. The DGS galaxies form a low-metallicity (12 + log(O/H) < 8.5) sub-sample of the galaxies presented in Rémy-Ruyer et al. (2015) for which the photometric data from the infrared to the sub-mm range are available. The SED is then built by adding the photometry at shorter wavelengths (UV, B and R) taken from the NASA Extragalactic Database (NED). We take care of using the most recent data corresponding to the integrated photometry of the galaxy. Employing data at short wavelengths allows us to better constrain the stellar mass and the SFR of the selected galaxies. LBGs were selected in the rest-frame UV and observed by ALMA. Burgarella et al. (2020) used the detections to build an infrared SED template that is later used to re-fit all the objects individually. The SED fitting has been performed employing a Chabrier initial mass function (IMF) in cigale with stellar mass in the range 0.1 ≤ M * ≤ 100 M (Bruzual & Charlot 2003). A top-heavy IMF is not available in the current version of the cigale and will be a future improvement in the code. For the galaxies selected, the code cigale provides a good fit between the synthetic and observed photometry (reduced χ 2 ≤ 5) for 31 out of 48 DGS galaxies and for 18 LBGs which represents the two samples studied here.

Model for chemical evolution of galaxies
In order to model the chemical evolution of the ISM of galaxies, we need to understand how the gas, metals and dust evolve with time, given the IMF of the stars, the star formation history (SFH) of the galaxy, and the theoretical yields from the stars.

Basic assumptions and parameters
The stars in the galaxies considered are assumed to form with a certain SFR. Consistently with the SED fitting performed in cigale, the functional form of the SFR used through this work is the "delayed" one: where the SFR is given in M yr −1 , τ is a characteristic time-scale and C SFR a constant defined in such a way that the final mass of stars is normalised to one solar mass: where M ev is the total mass released into the ISM from the evolving stars and t end = 13 Gyrs. In this work we tested two different types of IMF for the stars: the Chabrier IMF (Chabrier 2003) and a "top-heavy" IMF defined as: with α = 1, 1.35, 1.5. The constant C IMF is such that the integral of the IMF is normalised to one solar mass: where M min and M max are the minimum and maximum masses within which the theoretical metal yields are computed. Different metal yields available in the literature (characterised by diverse values of M min and M max ) are tested in this work (see Section 5.2.6). A top-heavy IMF implies a higher number of SNe II: the lower is the value of α the larger is the number of SNe II.

Gas and metals
Given the gas and metal yields for stars, the IMF, the SFR and outflows, we are able to describe the evolution of the mass of gas and chemical element, i, in the ISM. Under the aforementioned assumptions, the evolution of the mass of gas and metals is: where the first term is the integral between the minimum (M L ) and the maximum (M U ) mass of stars having life-time t that are enriching the ISM at time t with a mass gas fraction or metal fraction f g and f i,g , respectively. The stars enriching the ISM at a certain time were formed at t − t and have metallicity Z(t − t ). The fractions of gas and metal ejected into the ISM depend on the metallicity and on the stellar mass M * . The contribution to the ISM enrichment is weighed for the IMF. The second term of the equations is the astration of gas and metals due to formation of stars. The chemical enrichment of the ISM and astraction due to star formation (given by the first and second terms of Eqs. 5, 6) are computed by means of the omega code (Côté et al. 2017;Ritter et al. 2018, One-zone Model for the Evolution of GAlaxies). Starting from these calculations, the overall evolution of gas and metals is followed by an external routine. At each time-step the galactic outflow is assumed to be proportional to the SFR through the "mass-loading factor" (ML), which assumes the outflow to be due to the stellar feedback in the ISM (Murray et al. 2005). A corresponding mass of gas and metals is removed from the ISM. The code allows the possibility of introducing a galactic inflow of pristine gas provided by the last term of Eq. 5. Such a term is proportional to ML and regulated by the ratio between the mass of inflowing and outflowing gas (I/O). The total initial mass of baryons in the galaxy, M bar , is another input quantity of our model. This quantity is entirely composed by gas at the beginning of the simulation, and it is partially converted into stars according to Eq. 1. We set M bar to a multiple of the stellar mass (that we call M gas ) which is normalised to 1 M at 13 Gyrs, and we add to this quantity the mass of gas ejected by all the stars during their evolution, M ev in Eq. 2. The value of M ev is evaluated with the omega code.

Dust
We assume that the dust formed in supernovae remnants, around evolved low-mass stars and in the ISM is composed predominantly by silicates (olivine and pyroxene), amorphous carbon dust (characterised by different fractions of sp 2 and sp 3 bonds) and metallic iron. Carbon is also partially locked in hydrocarbons, such as in polycyclic aromatic hydrocarbons. Iron dust can be partially included in silicates or being in different chemical forms (e.g. iron oxides, iron carbides, iron sulphides, pyrrhotite). Throughout this work we consider iron-free silicates (olivine and pyroxene), amorphous carbon dust and metallic iron. Since we are interested in estimating the total mass of dust rather than its specific chemical composition, the presence of dust species of chemical composition different from the one here considered (such as iron-rich silicates, or other dust species that includes carbon or iron atoms) are not expected to significantly modify the predicted mass of dust with respect to our estimates.
Starting from the calculations from omega, we then follow the evolution of the dust content in the ISM, by assuming at each time step that a fraction of the available metals is condensed into dust grains. In addition, we take into account different physical processes that change the mass of dust as a function of time: outflows, dust destruction from SN shock waves that propagate in the ISM and dust growth in the ISM. Consistently with the calculations of gas and metal evolution, we assume that the inflow material is composed by pristine gas that does not change the total content of dust. For each dust species, the variation of the dust content in the galaxy will be given by: where the first term of the equation takes into account the metal enrichment from the composite stellar population and the astraction computed with omega and f i,d is the condensation frac-Article number, page 3 of 28 A&A proofs: manuscript no. aanda tion of dust 1 . Different condensation fractions for SNe and TP-AGB stars are assumed. The quantity dM i,d /dt in the first term is the maximum mass of each dust species that can be formed and it is computed by dividing dM i /dt in Eq. 6 by the mass of each element and by multiplying for the mass of the nominal monomer of dust given its stoichiometric formula. A fraction of dust is lost due to outflows and it is provided by the second term in the right hand of Eq. 7, where δ i,d (t) is the dust-to-gas ratio at time "t". The terms dM d,i,destr /dt is the variation of the mass of dust due to the dust destruction from SN shocks and dM d,i,gr /dt is dust growth in the ISM. We discuss the different mechanisms separately in the next sections.

Dust condensation fractions for TP-AGB stars and SNe
For TP-AGB stars, the chemistry of dust is mainly determined by the number of carbon atoms in the stellar atmosphere over the number of oxygen ones (C/O ratio). For C/O < 1 all the carbon atoms are locked in the very stable carbon monoxide molecules (CO) and the remaining oxygen is present in the atmosphere to form molecules and dust. The star is classified as oxygen-rich. For C/O > 1 an excess of carbon is present in the atmosphere from which molecules and dust are formed. The star is classified as carbon-rich. Normally, the value of the C/O is less than one, but during the TP-AGB evolution values of C/O > 1 are attained as a consequence of third dredge-up events that follow the thermal pulse (Herwig 2005). Around oxygen-rich TP-AGB stars, a certain fraction of silicon, oxygen and magnesium are used to build silicate dust (mainly olivine and pyroxene) and, in smaller amounts, metal oxides, while during their carbon-rich phase, the main dust species produced is carbon dust. Other species such as silicon carbide and metallic iron are also produced but in smaller amounts (Höfner & Olofsson 2018). The chemistry of the dust formed around SNe II is instead more uncertain, since it depends on the mixing between the C and O layers which determines the amount of such species locked into CO molecules (Nozawa et al. 2003). We here consider two different sets of condensation fractions for the different species formed around TP-AGB stars or in SN remnants. Around low-mass stars, the condensation fraction of silicate and carbon dust can reach up to 50-60% (or slightly more) during the superwind phase when the star loses most of its mass (Ventura et al. 2012;Nanni et al. 2013Nanni et al. , 2014, while for SN remnants the dust-to-metal mass fraction is estimated to be between 30 and 60% without considering the effect of dust destruction from the reverse shock and/or possible dust reformation (Marassi et al. 2019).
Following Nanni et al. (2013), for each time-step we assume the following condensation fractions for the different species for TP-AGB stars: f py = 0.3 for pyroxene, f ol = 0.3 for olivine, f car = 0.5 for carbon dust and f ir = 0.01 for metallic iron. Different condensation fractions for the dust species formed in SNe II are considered. Given the uncertainties in the dust condensation fraction for SNe II, we tested a range of values between 0.1 and close to 1 for all the dust species. In case the condensation fraction is close to the maximum, we set the condensation fraction of carbon equal to 0.5, since some should be available in the atmosphere to form CO molecules (Spyromilio et al. 2001).
Despite being a simplified approach, since dust yields are not consistently computed in the circumstellar envelopes of TP-AGB stars and in SN remnants for any of the metal yields considered in this work, such a choice allows us to test the results for different choices of the metal yields available in the literature.

Dust destruction from SN shock waves
We adopt a macroscopic description for dust destruction due to the passage of SN forward shocks in the ISM. We assume that all the dust species are destroyed with the same efficiency by such a process. The time-scale which regulates the destruction from SN shock waves is given by: where M g is the gas mass in the ISM which evolves according to equation 5. The quantity M swept , given in M , is a model parameter that provides how much of the ISM mass is swept by each SN event, and is the destruction efficiency. The quantity R SN is instead the SN rate that depends on the SFH and the IMF. For each dust species i the destruction term is therefore given by: In this work we consider two values of M swept = 1000, 6800 M , while = 0.1 in all the cases (McKee 1989). M swept = 6800 M represents the standard case in the literature (e.g. Dwek et al. 2007), while the assumption M swept = 1000 M has been also considered on the basis of the calculations performed by Nozawa et al. (2006).

Dust growth in the ISM
We here analyse the efficiency of the dust growth process in the standard framework where molecules are added from the gas phase to the grain, forming the bulk of silicates (olivine or pyroxene), amorphous carbon or metallic iron. No icy mantle formation is considered in this work. Different parameters are adopted in order to model such a process. Due to the lack of detailed laboratory measurements of grain growth at low temperature, many of these parameters are uncertain. We treat such a process following grain growth and the variation of the size of each dust species as a function of time, as calculated in Nanni et al. (2013Nanni et al. ( , 2014. We assume accretion to start on pre-existing grains (seed nuclei) on which other molecules from the gas phase are added to form the bulk of the grain. Such grains can be reprocessed in the ISM by SN shocks or aggregate in dense clouds modifying their size distribution. For each of the dust species the overall reaction for the formation of the grain needs to be chosen a priori. The reactions assumed in this work are provided in Table 1. Carbon dust is formed from acetylenic radicals (C 2 H (0<=y<=4) ) on the grain surface. We here only consider the formation of carbon through the addition of C 2 H 2 since we are interested in estimating the total amount of dust formed in the ISM. The exact pathway of grain formation is uncertain and such an investigation is beyond the scope of this paper. The possible destruction of dust due to the sputtering of the grains by H 2 molecules is neglected in this work, since this process is efficient above ≈ 1000 K (Gail & Sedlmayr 1999;Nanni et al. 2013). The variation of the grain size is regulated by the equation: where we assume the initial size of dust grains to be 0.01 µm (Hirashita & Kuo 2011) which is derived from the size distribution of grains in our Galaxy (Mathis et al. 1977). If smaller grains are initially available in the ISM dust growth in the ISM may occur more rapidly (Hirashita 2012). The quantity V 0,i is the volume of the nominal monomer of each dust species and the growth rate, J gr i , is instead evaluated as: where s j and s i are the stoichiometric coefficient in the reaction for the species in the gas phase and for the monomer of dust, respectively. The quantity α i is a unitless parameter known as "sticking coefficient" and represents the probability that a species in the gas phase sticks on the grain surface. The value of the sticking coefficient varies from 0 to 1. On the basis of recent experiments at low temperature, the reaction between SiO molecules and the grain surface is expected to have zero energy barrier that corresponds to a sticking coefficient equal to 1 (Rouillé et al. 2015). We here assume α i = 1 for all the dust species. Such an assumption, however, represents an upper limit of the efficiency of grain accretion. The quantity n j is the number density of the species j in the gas phase involved in the formation of the dust grain, computed as: where j is the abundance of the element j and we adopt n H = 10 3 cm −3 (Hirashita 2000(Hirashita , 2012. For larger values of the gas density grain growth would proceed more rapidly. At each timestep the quantity j in the gas phase is computed by subtracting the fraction of each element condensed into dust grains from the total available. The incorporation of CO molecules in grains is not considered. Carbon monoxide is a very stable molecule that does not contribute to the accretion of dust grains. The fraction of carbon monoxide in the gas phase is estimated to be around 40% of the total amount of carbon (Lodders 2010;Agúndez & Wakelam 2013). The abundance of SiO, Mg, H 2 O, C 2 H 2 and Fe in the gas phase is therefore estimated respectively as the abundance of Si, Mg, O, C and Fe neither locked in dust grains nor in CO molecules. The quantity v th,j , is its thermal velocity, given by: where k B is the Boltzmann constant, T gas = 25 K the gas temperature and m j the mass of the gas species. From the grain size it is possible to directly compute the total amount of dust produced in the dense phase of the ISM.
where ρ i is the mass density of the dust species i, n s,i is the number of seed nuclei. We estimate the total number of seed particles n s,i by dividing the mass of each dust species by the mass of a dust grain which depends on its size, as done in Asano et al. (2013). With this choice we implicitly assume that all the dust grains already present in galaxies can act as seed particle on which other molecules can accrete. The number of seed nuclei is expected to change as the ISM is progressively enriched by newly formed dust grains of a certain typical size and will depend on the different physical processes that modify the number of grains (e.g. shattering, coagulation). The larger is the number of seed nuclei the faster dust growth in the ISM is. If star formation mainly occurs at the beginning of the cycle (which is the case for the galaxies considered here), a further increase of the mass of dust in the ISM is mainly due to grain growth in the ISM. In this regime, the number of seed nuclei can be assumed to be roughly constant (if new seeds are not formed in the ISM) and equal to the total mass of a certain dust species divided by the current mass of the grain for each dust species, i as computed in Asano et al. (2013).

Method
The observed photometry for galaxies is fitted with the code cigale. The features of the code are explained in Boquien et al. (2019) while the assumptions adopted for computing the synthetic spectra are discussed in Burgarella et al. (2020). The SED fitting procedure of the sources described in Section 2 allows us to derive the best value for the parameter τ in Eq. 1 for each of the galaxies, as well as the total mass of stars and dust, and the SFR. Furthermore, throughout this paper we use the specific mass of dust (sMdust), i.e. the mass of dust divided by the mass of stars, specific SFR (sSFR), i.e. the SFR divided by the mass of stars, age. For DGS galaxies other properties are available from the literature: metallicity, gas fraction, dust-to-gas ratio, circumgalactic dust fraction. The properties derived are compared with the output of a set of chemical evolutionary tracks obtained by varying the input parameters described in Section 3. The different combinations of the parameters adopted in Section 3 (IMF, condensation fraction, stellar mass range of the SN progenitors, efficiency of the galactic outflow, SN destruction, and initial gas content in the ISM) are listed in Table 2. For all the models (one at each time-step) inside the different chemical evolution tracks, we find the residual between the properties of the galaxies predicted by the chemical evolution calculation and the ones derived from the SED fitting of cigale: where k represents the property of the galaxy. For the error on the age we adopted the maximum value between the uncertainty given by cigale and 10% of the age. The metallicity and the gas mass are taken from Rémy-Ruyer et al. (2013, 2014 for DGS galaxies. We compute the total mass of gas by selecting Article number, page 5 of 28 A&A proofs: manuscript no. aanda from Rémy-Ruyer et al. (2014) the atomic and molecular hydrogen mass. The total hydrogen gas is then multiplied by the mean atomic weight of the galaxy (µ gal ). The molecular hydrogen mass is derived from a metallicity-dependent conversion factor between the CO and H 2 masses (M H2,z ). All the data adopted are provided in Table A.1. In case only upper limits are available, we adopt the maximum possible value. The error associated to the gas fraction and to the dust-to-gas is assumed to be 100% given the uncertainties affecting the conversion between CO and H 2 abundances. For LBGs, we do not have direct measurements of the metallicity and of the gas content, and therefore we do not use such quantities for fitting the galaxies. The estimates of the circumgalactic dust fraction are taken from McCormick et al. (2018) for three of the DGS galaxies in the sample that are He2-10, NGC 1569, NGC 5253.
From the residual R gal we calculate the probability density distribution (in analogy with the χ 2 distribution), p = exp (−R 2 gal /2). Each of the properties of the galaxies is then estimated through the likelihood analysis by computing its average value and standard deviation by using as weights the probability density distribution among all the calculated models. The predicted values in different diagnostic diagrams (e.g. sMdust vs sSFR) are compared with the ones derived from the SED fitting with cigale or taken from the literature. This procedure allows to constrain the properties of individual galaxies.

Properties of galaxies derived from the SED fitting
We here discuss the main properties derived for the galaxies considered from the SED fitting procedure. The SED fitting performed for the DGS galaxies in this work in which the infrared data from Rémy-Ruyer et al. (2015) are combined with UV and optical photometry (see Section 2) allow us to better constrain the stellar mass and SFR of these galaxies with respect to the ones in Rémy-Ruyer et al. (2015) where only infrared data have been employed.
The average value of τ is estimated to be of 83 Myrs for LBGs, while for DGS galaxies, we selected two representative values to run the chemical evolutionary models which are of τ = 83 Myrs for those objects with τ < 100 Myrs (5 galaxies), and the average value of τ = 300 Myrs for the others 26 galaxies. LBGs are characterised by stellar population no older than ≈ 700 Myrs, while stellar populations up to ≈ 800 Myrs are derived for DGS. As shown in Fig. 1 the sSFR of DGS galaxies is between ≈ 10 −11 and ≈ 10 −7 yrs −1 , while for LBGs the range is narrower, between ≈ 10 −9 and ≈ 10 −7 yrs −1 . The sMdust is between ≈ 10 −4 and ≈ 10 −2 for DGS galaxies, while LBGs attain larger values up to ≈ 6 × 10 −2 . In the plot shown in Fig. 1 the evolutionary time increases from right to left. This represents a relevant diagnostic diagram largely employed throughout this work and in the literature. The main features for both the DGS galaxies and LBGs in the sMdust against sSFR (and age) shown in Fig. 1, is i) a peak at the beginning of the cycle; ii) a later decrease. LBGs attain a value of sMdust larger than DGS galaxies, while the decline of sMdust is faster. The large error bars in of sMdust obtained for some of the galaxies are due to the non detection for these objects of the dust continuum emission for which only upper limits are available.
For DGS galaxies we have additional information on the metallicity and on the gas content (See Rémy-Ruyer et al. 2013, 2014, and references therein). The metallicity values are 7.1 12 + log(O/H) 8.4, the gas fraction 0.4 f gas 1 and the dust-to-gas ratio 10 −5 D/G 10 −3 , where this latter quantity has been estimated by considering the dust mass derived in the analysis here presented. All the properties of DGS galaxies derived from the SED fitting performed with cigale are provided in Table A.1 in the Appendix together with all plots showing the best fitting spectra. For LBGs properties and SED fitting we refer to the work of Burgarella et al. (2020).
In Fig. 2 we show the comparison between the mass of dust, stellar content and SFR derived in Burgarella et al. (2020) and the ones in Rémy-Ruyer et al. (2015). We find that the stellar mass and SFR derived by Rémy-Ruyer et al. (2015) are systematically larger than the one derived in this work, while the mass of dust, constrained by the infrared emission is comparable. In Fig. 3 we compare the sMdust against the sSFR of the work by Rémy-Ruyer et al. (2015) with the ones derived from the SED fitting with cigale for the same galaxies. Our data points appear to be shifted at higher sMdust and a few of them present high sSFR. We will discuss in the following sections the consequences for our analysis. The trend shown has been derived through the SED fitting performed with a Chabrier IMF since a top-heavy one is not yet available in the SED fitting code cigale. By employing a top-heavy IMF in the SED fitting we do not expect significant variations in the stellar and dust masses derived, while the predicted age of the stellar population will tend to be older with lower values of sSFR (Pforr et al. 2012). Since the values of sMdust would remain approximately the same, we expect that the trend between sMdust and sSFR will still be present, but shifted to lower sSFR.

Reproducing the gas, metal and dust content of galaxies
We here present the results obtained from several tests from which we derive the suitable parameters to perform the system- Table 2. List of parameters adopted in the simulations of metal and dust evolution described in Section 3. First tests are run in order to select the reference parameters adopted to run systematic calculations. The stellar mass produced after 13 Gyrs is always normalised to 1 M . Different theoretical metal yields are tested in Section 5.2.6.

First tests τ [Myrs]
83, 300 M stars [M ] 1 M gas [M ] (2,10,20,100) ×M star M bar [M ] M gas + M ev M swept [M ] 1000, 6800 ML (0, 0.5, 0.6, 0.65, 0.8)×M gas .01, f car = 0.5 Dust growth in the ISM NO atic calculations and to build a grid of models through which we fit the properties of galaxies derived from cigale. The parameters adopted for these calculations and for computing the grid are listed in Table 2. We use different diagnostic diagrams such as sMdust, metallicity, gas fraction as a function of the sSFR (or age) to investigate the metal and dust enrichment in galaxies, as well as the gas and dust removal from their ISM.

Initial mass function and dust condensation fraction for SNe II
The first constrain that we are able to provide to our models is on the IMF, since not all the choices are able to reproduce the largest values of sMdust at large sSFR (corresponding to the beginning of the cycle). In the upper panel of Fig. 4 we show an example of calculation adopting the Chabrier IMF and the topheavy IMF with α = 1.35. The cases with α = 1, 1.5 yield results close to this latter case. We focus our investigation to models with τ = 300 Myrs, M gas = 100 × M stars and M swept = 1000 M . Such a choice of parameters allows to minimise the effect of dust destruction in the ISM (see Eq. 8) that may prevent the models to attain the largest values of sMdust at the beginning of the baryon cycle. For each of the models shown the efficiency of the outflow reasonably reproduces the observations. Similar trends are recovered by employing τ = 83 Myrs. The observed values of sMdust at the beginning of the cycle are not well reproduced by a Chabrier IMF with 50% condensation fraction for SNe. This holds both in case dust growth in the ISM is considered or neglected. The comparison improves by increasing the condensation fraction close to 100% for SNe which represents an extreme case, and still does not reproduce the sMdust for most of the galaxies at the beginning of the cycle. The observations are better reproduced if a top-heavy IMF is assumed. For such choice of the IMF the ISM is rapidly enriched with metals and dust: a large value of sMdust (≈ 10 −2 ) is attained within 100 Myrs (≈ sS FR ≈ 10 −8 yr −1 ) by assuming 50% condensation fraction for SNe. In the lower panel of Fig. 4 we also show the effect of changing condensation fraction for SNe II. Even when dust growth in the ISM is included, a low condensation fraction of dust in SNe (≈ 10%) prevents the models from reproducing the largest values of sMdust attained around sS FR ≈ 10 −8 yr −1 . Therefore, a larger amount of dust needs to be condensed (25-50%). In the same panel we also show the case with condensation fraction for SNe II close to the maximum. Also for this extreme value of the IMF and of the condensation fraction of dust the largest values of sMdust for some LBGs is not attained. This discrepancy between theoretical calculations and observations may depend on several factors, as for example, the chemical composition of dust assumed to derive the mass of dust from the far-infrared emission, which is affected by uncertainties of a factor of 10 (Ysard et al. 2019).
In summary: a top-heavy IMF helps to reproduce the largest sMdust derived from observations. With such a choice a dust condensation fraction around 50% is required. In this work, we consider condensation fractions between 25% and ≈100% for the systematic calculations.

Gas reservoir and star formation efficiency
For DGS galaxies it is possible to constrain the star formation efficiency (i.e. the mass of gas converted into stars) by studying the trend between their metallicity and the sSFR (or age). The comparison between model predictions and observations is shown in Fig. 5 where the metallicity (log(O/H) + 12) is plotted as a function of the sSFR. As representative case, we study the behaviour of models with τ = 300 Myrs, M swept = 1000 M , condensation fraction of dust for SNe II of 50%, a topheavy IMF with α = 1.35, different mass-loading factors and M gas = 2 − 100 × M star . The initial mass of baryons is composed only by gas at the beginning that will then form stars according the star formation law define by Eq. 1. The quantity (log(O/H) + 12) is the abundance of oxygen in the gas phase, which is obtained in the models by subtracting the oxygen condensed into dust grains (silicates) and in CO molecules from the total. From this plot we conclude that the initial mass of gas in the galaxy should be between 10 and 100 times the final stellar mass. For lower mass of gas, the ISM is enriched too quickly with respect to the observations. In absence of outflow, the metallicity increases for a given choice of the mass of baryons, until a plateau is reached. If some galactic outflow is present, the final metallicity attained is larger and it rises more quickly. This is due to the fact that the outflow remove simultaneously both metal and hydrogen from the ISM, however, since the mass of hydrogen is larger than the mass of oxygen, hydrogen is removed in larger amounts, and the ratio decreases with time (see also Eq. 6). Different combinations of initial baryon mass and outflow are able to cover the observed range of metallicity.
In summary: an initial gas reservoir between ≈10 and ≈100 times larger than the final stellar mass is needed to reproduce the observed metallicity of DGS galaxies. Different choices of the mass-loading factor allow us to cover the observed range of metallicity values. This constraint is used to set the input param-eters for the evolution of LBGs for which the measurements of the metallicity are not available.

Dust destruction by SNe, star formation and removal from galactic outflows
After a peak in sMdust, a decline is observed for decreasing sSFR or increasing age, as shown in Fig. 1. In general, this effect can be due to combination of (i) star formation which increases the total stellar mass, (ii) dust astraction by star formation, (iii) dust destruction from SN shocks, (iv) galactic outflows. In the two panels of Fig. 6 Fig. 6, it is clear that if outflow is not included, only in the case with the lowest mass of gas and no grain growth in the ISM, the combination of dust destruction from SNe and astraction due to star formation partially reduces sMdust. In any case, all the models without outflow remain too flat for decreasing values of the sSFR. This trend indicates that dust astraction and destruction from SNe are not sufficient to decrease sMdust for decreasing sSFR, even for the maximum efficiency of dust destruction from SNe (M swept = 6800 M ). Indeed, since for most of DGS galaxies a large amount gas is required in order to be able to reproduce their metallicity trend, dust destruction by SNe is not efficient by construction (see Eq. 9). By employing different values of ML combined with M gas we are able to reproduce different sMdust in the galaxies. Therefore, we conclude that galactic outflows are required in order to reproduce the observed trend between sMdust and the sSFR for sS FR < 10 −8 yr −1 for both DGS galaxies and LBGs. We also notice that despite the fact that a combination with a mass of baryons equal to 100 the final mass of the stars without galactic outflow will be suitable to reproduce the observed metallicity of galaxies with log(O/H) + 12 < 7.8 and sS FR < 10 −9 yrs −1 the corresponding sMdust which is less than few 10 −3 at those sSFR is not reproduced. Galactic outflows are also needed in order to reproduce the observed gas fraction as a function of the sSFR for DGS galaxies. In Fig. 7 we show the observations for DGS galaxies overplotted with the evolutionary models characterised by different M gas and outflow efficiencies. In the case without outflow, only astraction is at work in reducing the amount of available gas. As expected, gas consumption by star formation only mildly affects the gas fraction in case M gas = 10 × M stars , while the effect is negligible for larger gas content. Galactic outflows are therefore necessary to efficiently decrease gas fractions for decreasing sSFR.
In summary: galactic outflow is an essential feature to reproduce the decline in sMdust for sS FR 10 −8 yr −1 observed for both DGS galaxies and LBGs, and for reproducing the trend between the observed gas fraction and the sSFR in DGS galaxies. We therefore considered different efficiency for this process in the systematic calculations.

Efficiency of dust growth in the ISM
We here discuss the efficiency of dust growth in the framework introduced in Section 3.3.3 in order to assess the relevance of this process in the galaxies under study. By looking at Figs. 4 and 6 it is possible to compare the efficiency of dust growth for different physical assumptions adopted in the simulations. In all the models shown, the chemical evolutionary tracks computed by including dust growth overlap with the ones without dust growth between sSFR of ≈ 10 −8 yr −1 and 10 −7 yr −1 which corresponds to the early dust enrichment from SNe II. Around sSFR of ≈ 10 −8 yr −1 the mass of dust tends to increase due to the effect of grain growth in the ISM. When the outflow is not included, the theoretical curves rise before reaching a plateau which correspond to the maximum possible condensation of dust from the available metals (Fig. 6). This trend is qualitatively similar to the ones obtained by Asano et al. (2013). In the models without outflow, the increase of the dust in the theoretical tracks due to grain growth does not reproduce the observed decrease in sMdust, as already noticed for the models without grain growth. The models with and without dust growth become more similar for increasing efficiency of the outflow, which is an essential feature to reproduce the observed trends. In this case, the effect of grain growth in the ISM tends to be cancelled by dust removal. Furthermore, the effect of dust growth combined with a more efficient outflow is degenerate with the models without grain growth and less outflow.
From Fig. 4 we also notice that for low dust condensation fractions of SNe II and/or an IMF different from the top-heavy one, the galaxies with the largest sMdust observed at the beginning of the cycle are not reproduced even if grain growth in the ISM is included in the calculations. This might indicate that such a process is not the main one responsible for dust production at the beginning of the baryon cycle, unless we assume that grain growth occurs at larger densities (> 10 3 cm −3 ) or/and the dust grains initially present in the ISM are small (< 10 −6 cm).
In summary: dust growth in the ISM is not dominant compared to outflow. Thus, it is not necessary to include it in order to reproduce the trend between sMdust and the sSFR for both DGS galaxies and LBGs. If present, the effect of such a process is not evident because combined with the one of galactic outflow.

Galactic inflows
We here explore the possible effect of introducing galactic inflows in our models in addition to the outflow required to reproduce different observations. In the two panels of Fig. 8 we show as example the effect introduced by including different amount of inflow besides galactic outflow in the sMdust against sSFR and in the metallicity against sSFR plots for a few selected models with τ = 300 Myrs, M gas = 50 × M stars , M swept = 6800 M and different choices of the galactic inflows and outflows. Since the dust and gas need to be removed from the galaxy I/O in Eq. 5 needs to be < 1. We select as test cases I/O = 0.2, 0.5 combined with efficient outflow (ML = 0.6, 0.8 × M gas ). As shown in Fig. 8, in those models in which galactic inflows is include an efficient outflow is needed in order to remove the dust in the galaxies and to reproduce the observations. The difficulty in removing the dust in the ISM in presence of an inflow is caused by the fact that the dust is diluted in the ISM. This process decreases the dust fraction δ i,d in Eq. 7 and dust is removed less efficiently from the galaxy. From the lower panel in the same figure, it is instead possible to appreciate the slower metal enrichment due to the inflow of pristine gas, and the corresponding decrease of the fraction of metals ejected in the outflow (due to the decrease of f i,g in Eq. 6). Furthermore, the various combinations of galactic inflow and outflow introduce some degeneracy in the models with the ones discussed in the previous sections. We chose not to introduce any further degeneracy in our systematic calculations and to assume that all the gas is already present in the galaxy before the beginning of star formation.

Metals from different theoretical data sets and the contribution of different sources
We here explore the effect on the chemical enrichment by employing different sets of theoretical metal yields taken from the literature, combined with the different time-scales of the delayed SFH and top-heavy IMF. We consider 7 different theoretical sets for SNe II, 3 for TP-AGB stars and 2 for Pop III stars (characterised by zero metallicity). For pop III stars we assume the same IMF as for the other stars. The metal yields for SNe II by Kobayashi et al. (2006) and Nomoto et al. (2013) Table 2.
Starting from these metal yields we calculate the maximum possible amount of dust produced considering silicates, carbon and iron dust. For simplicity, we consider the case in which silicate dust is composed entirely by pyroxene (MgSiO 3 ). However, the mass of silicates derived if only olivine (Mg 2 SiO 4 ) is considered does not change considerably. For silicates, the maximum possible mass of dust is estimated at each time-step by evaluating the least abundant of all the elements that form this dust species. The final mass of dust is then obtained by estimating the maximum number of monomer and mass of dust that can be formed from the stoichiometric formula of pyroxene.
A representative example of the results for SNe II Fig. 9 where the maximum possible sMdust of silicate, carbon and iron dust is plotted a function of time. A top-heavy IMF with α = 1.5 and τ = 300 Myrs have been selected. The trend obtained by employing different data sets is qualitatively similar for different IMF and τ. The spread in the maximum sMdust spans a factor of about 5 for all the dust species considered. For all the data sets but Limongi & Chieffi (2018), the final mass of silicate dust is larger than carbon by a factor between 1.2 and 3.6. For the three data sets by Limongi & Chieffi (2018) the mass of silicate dust is down to one third the mass of carbon dust. The maximum iron mass is around 10% and 15% of the total dust mass for Kobayashi et al. (2006), Limongi & Chieffi (2018) and for Nomoto et al. (2013), respectively. For the two data sets by Ritter et al. (2018) the mass of iron dust is between 30-35% of the total. The sets of models which produce the largest maximum mass of dust are from Kobayashi et al. (2006) and Limongi & Chieffi (2018) with rotational velocity R=150, 300 km s −1 . The three data sets yield very similar values of sMdust from SNe II close to the expected observed peak of dust around 100-200 Myrs. The most remarkable difference among the dust produced by these three data sets is the predicted chemistry of dust which would be dominated by silicates for Kobayashi et al. (2006) and by Limongi & Chieffi (2018). For Kobayashi et al. (2006) the maximum amount of silicates produced is three times larger than carbon dust. The differences between the dust chemistry of the ISM obtained by adopting different theoretical yields might be compared and constrained with future observations. For these three data sets we also check the amount of oxygen released in the ISM, since oxygen is the most abundant metal and it traces the metallicity of galaxies. The predicted oxygen abundances are compared with observations. We find that the values are always comparable even for different τ and for the different top-heavy IMF. Furthermore, since reproducing the large mass of observed dust in the Early Universe is challenging, we expect that the comparison with observations would worsen by employing the theoretical data sets yielding a lower amount of metals and dust with respect to Kobayashi et al. (2006) and Limongi & Chieffi (2018).
To test the effect on the yields from TP-AGB stars we consider different combinations of theoretical yields for SNe II and TP-AGB stars. The chemical enrichment from TP-AGB stars is affected by the choice of the SN theoretical yields that change the metallicity of low-mass stars as a function of time and therefore the overall metal and dust masses released in the ISM. The combinations considered are provided in Table 2. The results are shown in Fig. 10 where we plot the maximum possible sMdust for SNe II, of TP-AGB stars and of SNe Ia. For all the combinations of yields, we find that the contribution from TP-AGB stars and SNe Ia is negligible, due to the top-heavy IMF. Furthermore, in case a larger amount of dust would be released by TP-AGB stars or SNe Ia, it would start to be relevant only after a few hundreds Myrs, when the dust is expected to start to decrease in the observed galaxies. Comparable amount of sMdust could be obtained for TP-AGB stars and SN remnants by assuming drastic low condensation efficiency in SN remnants, that is not our preferred scenario, as discussed in the previous sections. The contribution from TP-AGB stars is more relevant in case the selected IMF is the Chabrier one. However, this scenario is never able to reproduce the observed trends for the galaxies considered in this work.
We additionally tested two different theoretical yields for Pop III stars (Heger & Woosley 2010;Nomoto et al. 2013). We find that the contribution to the metal enrichment of Pop III stars is negligible with respect to the total.
In summary: the metal yields selected throughout this paper for the systematic calculations are from Kobayashi et al. (2006) since it allows for the most favourable conditions for attaining the largest values of sMdust already after 100 Myrs since the beginning of the baryon cycle. The choice of the metal yields for Pop III, TP-AGB stars and Type Ia SNe do not largely affect the results, since these stellar sources provide only a minor contribution to the dust enrichment. We therefore arbitrarily select the theoretical yields by Heger & Woosley (2010) for Pop III stars (at the ages appropriate for the sample), Cristallo et al. (2015) for TP-AGB stars and Iwamoto et al. (1999) for Type Ia SNe.

Input quantities for the systematic calculations
On the basis of the previous tests we are able to define the characteristics required to systematically compute the chemical evolutionary models in order to perform the best fit between model predictions and observations that we here recall: metal yields favouring a fast metal enrichment from SNe II (e.g. Kobayashi et al. 2006); top-heavy IMF; total mass of baryons ≥ 10 times the final stellar mass; efficient galactic outflow; condensation fraction for dust produced in SNe II ≥ 25 %; no necessity of dust growth in the ISM.
The input quantities adopted to compute the grid of models are summarised in Table 2 ("Systematic calculations").

Characterisation of the individual galaxies
Following the procedure described in Section 4, we estimate the properties of DGS galaxies and LBGs, such as the chemical composition of their dust, the fraction between the gas and dust ejected in the circumgalactic medium and the fraction of dust destroyed by SN shocks. In Figs. from 12 to 15 we show the different properties of galaxies derived from the SED fitting or taken from the liteature, i.e. sMdust, metallicity, gas fraction and dust-to-gas ratio, as a function of the sSFR and age together with the corresponding distribution obtained from the fit of individual objects through the models in the chemical evolutionary tracks. In Figs. 13-15 we also show the metallicity, gas fraction, and dust-to-gas ratio as a function of the sSFR and time for LBGs, for which reliable constraints are not available yet. We additionally show the dust-to-gas ratio as a function of the metallicity in Fig. 16. For the DGS galaxies such trend is compared with the available estimates, while for LBGs only the values derived from the chemical evolution models are shown. In Fig. 17 we show and the total mass of dust inside the galaxy versus the mass of stars.
The best fit between the properties of individual galaxies and model predictions are provided in Table 3. The complete version of the table is provided in electronic form.
From the performed analysis we derive the following trends: -Overall properties. The models developed here are in general able to fit the main properties of the observed DGS galaxies and LBGs. The stellar production efficiency (defined as the mass of stars over the total mass of baryons initially available) is usually of few % for both DGS galaxies and LBGs. A strong outflow is required in order to remove the mass of dust and gas present in the galaxy and to reproduce the observation both for DGS galaxies and LBGs. Depending on the specific galaxy under consideration, the amount of gas removed from the ISM ranges between 6 or 12% for DGS galaxies and LBGs respectively and almost the totality. The exponent of the IMF in Eq. 3 is α ≈ 1 for the majority of the DGS galaxies, while this value is not well constrained for LBGs for which the results are more degenerate. This is not surprising since for LBGs less information are available. The set of observations for DGS galaxies is reproduced by a condensation fraction of SNe II which is between 25% and 64% with an average value of (35±12)% considering the entire sample of DGS, while for LBGs a condensation fraction above 40% is required with an average value of 56±17% if all the LBGs are considered. Also in this case, a larger uncertainty on the values of the condensation fraction is found for LBGs. We also keep in mind that there is some degeneracy between the metal yields provided in the literature and the assumed dust condensation fraction. The dust-to-gas ratio as a function of the metallicity is fairly well reproduced for DGS galaxies (Fig. 16), even though the observed distribution of data points is more tilted than the one we recovered. This discrepancy between our model predictions and observations might be partially due to the large uncertainties affecting the determination of the gas mass from the observations, especially in low-metallicity galaxies. The trend between the total mass of dust and stellar mass is naturally reproduced (Fig. 17). -Dust chemical composition and condensation fraction.
The predicted chemical composition of dust changes as a function of time. This is determined by the IMF and by the condensation fractions assumed for the different species. The dust chemical composition is dominated by silicates which are always between 65-80% of the total mass fraction. Carbon dust mass is instead a factor between 0.6 and 1.6 the mass of iron. We are aware of the fact that such chemical Table 3.
Results of the sample of DGS galaxies and LBGs that have been individually fitted through the models provided in Table 2. For each of galaxy analysed we provide the averaged residual, the predicted mass of the different dust components (silicates, carbon and iron) normalised for the final stellar mass, the dust condensation fraction for SNe II, the initial mass of baryons normalised for the final stellar mass, the slope of the top-heavy IMF, the M swept in M , the fraction of the dust in the circumgalactic medium over the total ( f dust,out ), the fraction of gas in the circumgalactic medium over the total ( f gas,out ) and the ratio between the dust removed by the galactic outflow over the total removed and destroyed ( f dust,out/(out+SN) ). A&A proofs: manuscript no. aanda composition of dust is expected to be dependent on the metal yields selected for SNe II. -Efficiency of the outflow and of SN shocks for dust removal. We estimate the amount of dust removed from the ISM through galactic winds and destroyed by SN shocks for each of the galaxy considered. Depending on the galaxy, between ≈20% and ≈90% of the dust disappearing from the ISM is removed by galactic outflows rather than destroyed by SN shocks. For three of the DGS galaxies in our sample NGC 1569, He2-10, NGC 5253) estimates of circumgalactic dust are available from McCormick et al. (2018). We use such a piece of information for fitting the galaxies against the models, the observed and predicted distributions are shown in Fig. 11. For the three galaxies considered, our approach provides a satisfactory fit of the observed circumgalactic dust fraction. -Comparison between DGS galaxies and LBGs. The metallicity values and gas fractions predicted by our models for LBGs are comparable with the ones derived for DGS, however we predict a faster metal enrichment for the LBGs than for DGS galaxies (lower panel of Fig. 13) and the gas is ejected in a shorter time-scale for LBGs (≈ 400 Myrs) than for DGS galaxies (≈ 1Gyr). From Fig. 15 it is possible to notice that the typical dust-to-gas ratio of LBGs is comparable with the upper limit observed for DGS galaxies as consequence of the large condensation fraction required to explain the sMdust of these galaxies.

Discussion
In the present work we show that a top-heavy IMF and high condensation fraction in SN remnants are required in order to have a fast increase of the dust content in the ISM of galaxies within 100 Myrs since the beginning of the baryon cycle. This finding is in agreement with other works in the literature in which the large amount of dust in Quasars at z > 6 is explained by adopting similar assumptions (Gall et al. 2011a,b). This result is also in agreement with the work from Liu & Hirashita (2019) who find that the most of the ALMA observations of galaxies at z > 7 can be explain by assuming a dust condensation fraction of ≈ 0.5 for stellar sources. In the local Universe, top-heavy IMF have been adopted to explain some observational properties of ultra compact dwarf galaxies , ultra faint dwarf galaxies (Geha et al. 2013;McWilliam et al. 2013) and Galactic globular clusters (Marks et al. 2012).
On the other hand, in many works in the literature dust growth in the ISM has been considered of fundamental importance in order to explain the amount of dust observed in local and high-redshift galaxies (Asano et al. 2013;Zhukovska 2014;Michałowski 2015;Aoyama et al. 2017;Ginolfi et al. 2018;Leśniewska & Michałowski 2019). Popping et al. (2017) have predicted an efficient dust growth in the ISM for galaxies with stellar masses log(M * /M ) > 7, while dust enrichment from stellar sources appears to be dominant for galaxies characterised by lower masses. Mancini et al. (2015) and Graziani et al. (2019) have found a similar trend but for larger stellar masses of log(M * /M ) > 9 and 8.5 < log(M * /M ) < 9.5, respectively. The dependence between the efficiency of grain growth in the ISM and the stellar mass of the galaxies found in these works is ascribed to the presence of more metals in more massive galaxies. Schneider et al. (2016) have also pointed out the dependence of grain growth efficiency with the density of the gas in the galaxy. da Cunha et al. (2010) use the model from Calura et al. (2008) for dwarf galaxies which includes grain growth in the ISM, and in which a Salpeter IMF is adopted for these galaxies. However, the galaxies at larger sMdust and sSFR are not reproduced by their models (see Fig. 10, left panel of da Cunha et al. (2010)). This result is not in contradiction with our analysis, since we need to assume a top-heavy IMF in order to reproduce the larger values of sMdust and this option has not been explored in their analysis.
In the analysis performed by Rémy-Ruyer et al. (2015) the models developed by Asano et al. (2013) are employed. As shown in Fig. 11 of Rémy-Ruyer et al. (2015), the sMdust of galaxies with sSFR between 10 −9 and 10 −8 is not reproduced. Following the analysis by Rémy-Ruyer et al. (2015), De Vis et al. (2017 reinterpreted the properties of the same sample of galaxies by studying their chemical evolution under different scenarios (e.g. with and without dust growth in the ISM). They have obtained their best fit models and the properties of DGS derived by Rémy-Ruyer et al. (2015) by employing a Chabrier IMF and assuming that a) a star formation history characterised by several bursts as proposed by Zhukovska (2014), b) dust can grow in the ISM but at lower rates with less efficient star formation, and c) a stronger inflow and outflow in DGS galaxies than in less star-forming galaxies. Despite the ability of such a model to reproduce the overall observed trends of different samples of local galaxies, the DGS characterised by the largest sSFR remain challenging to explain (see the top-left panel of Fig. C1 in De Vis et al. (2017)). The tension between this kind of models and observations is even more critical in our analysis since our SED fitting yields values of the sMdust at earlier epochs (larger sSFR) that are even larger than the ones obtained by Rémy-Ruyer et al. (2015) and employed by De Vis et al. (2017Vis et al. ( , 2019, as discussed in Section 5.1. Such a difference implies that larger values of sMdust should be achieved in DGS galaxies earlier than in the aforementioned works. In our work, this effect is obtained by adopting a top-heavy IMF which increases the amount of metals available to form dust already at early epochs coupled with a condensation fraction > 25% for SNe II. In such a framework, dust growth in the ISM is not needed to produce the total amount of dust observed. Moreover, we show that this process would become efficient when the galactic outflow is expected to start to reduce the amount of gas and dust in the galaxies. Therefore, even if present, the outflow would cancel the effect of dust growth in the ISM. The trend between the dust-to-gas ratio and the metallicity is well reproduced for DGS galaxies by Zhukovska (2014) who has included dust growth in the ISM and very low condensation fraction for silicates in her calculations ( f sil = 10 −3 ). The same observed trend is not well reproduced by assuming a more efficient condensation in SN remnants and no dust growth in the ISM, similarly to what we derive in Fig. 16. Contrary to our work and to the analysis by De Vis et al. (2017Vis et al. ( , 2019 a closed-box model is employed in Zhukovska (2014) together with a variety of SFH for the galaxies. However, the comparison with the observations is limited to the trend between the dust-to-gas and metallicity. In our analysis, the low condensation fraction from SNe II required to reproduce this trend does not reproduce the large value of sMdust at the beginning of the cycle. The discrepancy between our model predictions and observations might be partially due to the large uncertanties in the determination of the mass of gas in low metallicity galaxies. Ginolfi et al. (2018) have also proposed dust growth in the ISM of dwarf galaxies to explain the observed amount of dust. However, the investigation is limited to a Salpeter IMF, and the possibility of producing a larger amount of dust by employing a top-heavy IMF is not explored in their work.
From the point of view of the micro physics of grain accretion, Ferrara et al. (2016) has indicated that dust accretion in the diffuse ISM is difficult to obtain due to the low gas densities that imply a low accretion rate, even taking into account that the energy barrier between SiO molecules and silicate grain surface has been experimentally found to be zero even at low temperatures (Rouillé et al. 2015). One way of increasing the depletion of the gas elements in the ISM and the size of dust grains is by forming icy mantles around the bulk of the grain in dense molecular clouds, where the density is higher than in the diffuse medium. The icy mantle formed around dust grains can prevent the process of grain accretion by creating an energy barrier between molecules in the gas phase and the bulk of the grain. Icy mantles can however also partly evaporate when stars are formed. Such a process might be particularly efficient if the SFR is relatively high as it is the case for DGS galaxies and LBGs. Ceccarelli et al. (2018) have further developed the idea proposed by Ferrara et al. (2016) by studying how the chemistry of the icy mantles changes as a function of different parameters, including the metallicity. The aim of such a work, is to estimate the probability for Si-bearing species to encounter in the icy mantle and to form clusters that can build up silicate grains. The authors have found a very low probability for cluster formation due to this process, especially at the low metallicities that characterise the galaxies considered in this paper. Conversely, the results of a recent experiment suggests that silicate and carbon dust may be formed on the grain surface following the evaporation of icy mantles (Rouillé et al. 2020). Other laboratory experiments and theoretical works suggest that a non-negligible amount of dust can be formed in low-temperature dusty plasma (Bleecker et al. 2006;Hollenstein 2000). However, under which conditions such a process is relevant (e.g. intensity of the radiation field providing different ionisation fractions for the elements) has only been investigated for the typical conditions of the Milky Way (Zhukovska et al. 2016).
The inability of models to reproduce the largest sMdust observed may suggest that dust condensation in SN remnants might be the preferred scenario to reproduce the observations discussed in this paper.
However, how much dust is produced and destroyed in the SN remnants is still very debated. A large amount of dust (0.1 − 0.5 M ) is formed in situ in SN remnants (Matsuura et al. 2011), but it is still unclear how much dust is destroyed by the reverse shock or reformed. Matsuura et al. (2019) have interpreted the increase in the 31.5 µm photometry in SN1987A obtained from SOFIA FORCAST 11.1, 19.7 in 2016 with respect to 10 years earlier as a possible indication of dust re-condensation in the forward shock. Indeed, the mass of dust needed to explain the emission of this dust component at a temperature of ∼ 85 K is more than 10 times larger than the one estimated 10 years before. This might indicate that dust could be reformed after being destroyed by shocks. Such an interpretation should be confirmed by the future time of the observed spectra. Indeed, another possible explanation for the observed emission is that the decrease in the density of the expanding ejecta has allowed X-rays to efficiently heat the freshly formed dust. The re-formation of dust in the post-shocked regions has also found theoretical confirmation (Sarangi et al. 2018). From the observational point of view Gall et al. (2014) have found that large dust grains are rapidly formed in SN remnants. Such grains are large enough to be able to survive to the reverse shock. From the theoretical side, different authors have studied the formation and dust survival in SN ejecta. On one hand, Bianchi & Schneider (2007); Bocchio et al. (2016) have estimated that only between 2-20% of the dust mass observed in SN remnants will survive to the reverse shock, and that the mass of dust observed around SN1987A has not been yet destroyed by the passage of the reverse shock. On the other hand, different calculations suggest that a lager fraction (42-98%) of the newly produced dust can survive to the reverse shock in clumpy ejecta (Biscaro & Cherchneff 2016). The survival of dust grains can also depend on the their size, on the energetic of the explosion, on the gas density of the ISM and on the thickness of the hydrogen envelope (Nozawa et al. 2007).
In our work we also show that another condition to reproduce the metallicity trend, is related to the initial gas content in the galaxies that has to be at least more than 10 times the final mass of stars. In this context, galactic outflows are needed to efficiently remove gas and dust from the galaxy. Such a result is observationally supported for LBGs in which outflows are commonly observed (Shapley et al. 2003;Pettini et al. 2002;Gallerani et al. 2018), while star formation is rather inefficient. As far as DGS galaxies are concerned, circumgalactic dust has been detected in different galaxies included in our sample (e.g.

Conclusions
In this paper we develop and revise the existing prescriptions for dust evolution in galaxies in order to interpret the observations of DGS galaxies in the local Universe and of LBGs at high-redshift.
Different investigations in the literature show the inadequacy of the current models for reproducing the largest values of the sMdust of mass observed at the largest sSFR for DGS galaxies. The inability of the current models for reproducing the observations, is even more severe in our investigation, in which the sMdust and sSFR are even larger than the previous estimates in the literature (Rémy-Ruyer et al. 2015). This indicates that the stateof-the-art framework adopted for interpreting the dust content in these galaxies need to be revised. In particular, the dust prescriptions that assume low condensation efficiency for SN and dust growth in the ISM are not able to reproduce the largest sMdust at early epochs estimated for DGS galaxies and LBGs. For LBGs we develop a new description for metal and dust evolution that is in part based on the results obtained for DGS that are their nearby counterparts.
For both samples of galaxies, we reproduce the observations by adopting a top-heavy IMF which favour the fast enrichment of the ISM through SNe II and dust condensation fraction > 25% for DGS galaxies and > 40% for LBGs. In this context, the dust enrichment from TP-AGB stars and from Type Ia SNe plays a minor role, since a top-heavy IMF is not favourable for the formation of massive stars. Galactic outflows are essential to reproduce the decline of the sMdust (and of the gas fraction for DGS galaxies) as a function of sSFR. Dust destruction by SN usually plays a secondary role with respect to dust removal from the galactic outflow. A typical star formation efficiency (the mass of gas converted into stars) of few per cent is instead required to reproduce the trends with the metallicity.
Despite grain growth is often considered to be a fundamental process to explain the amount of dust observed, we find that the contribution of such a process starts to be relevant when the outflow start to remove dust from the ISM, and therefore its pres-        LBGs (black tringles) overplotted with models characterised by different choices the mass-loading factor and inflows as indicated in the legend. For all the models we select a value of τ = 300 Myrs in Eq. 1, α = 1.35 for the top-heavy IMF in Eq. 3 and M swept = 6800 M in Eq. 8, M gas = 50 × M stars and dust condensation fraction for SNe II of 50%. Lower panel: metallicity vs sSFR for DGS galaxies overplotted with the same models as in the upper panel.      Lower panel: the dust-to-gas ratio as a function of age. The same colourcoding as in Fig. 12 Fig. 16. The observed and predicted dust-to-gas ratio as a function of the metallicity. The same colour-coding as in Fig. 12 is adopted. The predicted distribution for LBGs are here shown for comparison.