Open Access
Issue
A&A
Volume 668, December 2022
Article Number A191
Number of page(s) 14
Section Cosmology (including clusters of galaxies)
DOI https://doi.org/10.1051/0004-6361/202244334
Published online 20 December 2022

© The Authors 2022

Licence Creative CommonsOpen Access article, published by EDP Sciences, under the terms of the Creative Commons Attribution License (https://creativecommons.org/licenses/by/4.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

This article is published in open access under the Subscribe-to-Open model. Subscribe to A&A to support open access publication.

1. Introduction

In order to understand the chemical evolution of the entire Universe, it is vital that we understand the global mechanisms that have dominated the production of chemical elements since the Big Bang. Primordial nucleosynthesis is responsible for the synthesis of deuterium, 3He, 4He, and traces of 7Li, and it ceased after the first few minutes of the existence of the Universe. After that period, a global chemical enrichment would resume only when stellar nucleosynthesis started inside the nuclei of stars.

Two classes of stars are of primary interest for the chemical evolution scenario: primordial, metal-free stars, known as first stars or Population III (Pop III) stars, and second-generation, more enriched stars, known as Population II (Pop II) stars. It is believed that Pop III stars had very unusual proprieties, and despite intense observational efforts, these stars have not yet been observed, although a few candidates have been proposed (e.g., Kashikawa et al. 2012; Sobral et al. 2015; Vanzella et al. 2020). Thus, researchers have been combining efforts to build consistent modeling of these stars in recent decades (Heger & Woosley 2002, 2010; Schaerer 2002; Chieffi & Limongi 2004; Takahashi et al. 2018). Results indicate that they would have very large masses, between 100 and 200 M, or even higher, around 500 to 1000 M (Ohkubo et al. 2006), mainly due to the lack of metals in the gas, making cooling processes very inefficient (Galli & Palla 2013; Hirano & Yoshida 2013).

Also, the question about their initial mass function (IMF) and their role in the reionization of the Universe (Abel et al. 2000; Nakamura & Umemura 2001; Gnedin & Ostriker 1997; Cojazzi et al. 2000; Tumlinson & Shull 2000; Ciardi et al. 2001; Loeb & Barkana 2001; Venkatesan et al. 2003; Tumlinson et al. 2004; Barkana 2006), among other properties, have been discussed to better understand the physics of this high-mass stellar population.

In terms of chemical production, the above results indicate they were substantially important. For instance, Pop III stars with masses from 140 to 260 M produced huge amounts of metals. They ended their lives as pair-instability supernovae (PISNe), injecting highly enriched material back into the interstellar medium (ISM) and intergalactic medium (IGM), and leaving no remnants behind after a complete disruptive process (Heger & Woosley 2002; Takahashi et al. 2018).

After the first generation of Pop III stars started to die, injecting large amounts of metals into the ISM and IGM, Pop II stars started to form. With more and more enriched material available, cooling processes started to become more efficient, giving origin to less massive stars with physical proprieties close to those observed today, and also extensively modeled mainly according to their masses and metallicities (Chieffi & Limongi 2004; Kobayashi 2005; Campbell & Lattanzio 2008; Karakas 2010; Doherty et al. 2013, 2014).

Cosmological chemical evolution models have also been largely explored. There are several different models that seek to evaluate different aspects related to the chemical enrichment, such as the evolution of the mass-metallicity relation (Ma et al. 2016; Torrey et al. 2019), the establishment of a critical value for the metallicity of the Universe, enabling the transition from Pop III to Pop II stars (Bromm et al. 2001; Bromm & Loeb 2003; Fang & Cen 2004; Matteucci & Calura 2005; Santoro & Shull 2006; Tornatore et al. 2007; Maio et al. 2010; Schneider 2010), Hypernovae (HNe) feedback (Kobayashi et al. 2007), the role of galactic outflows (Davé & Oppenheimer 2007), the chemical properties of local galaxies based on their formation through the hierarchical model of structure formation (Calura & Menci 2009), the evolution of N abundance in the Universe and the reason for a large dispersion in observational data (Vangioni et al. 2018), and the influence of dark matter (DM) halos on the gas reservoir available for star formation (Lilly et al. 2013), among other examples of interesting contributions to the study of the Universe through its chemical enrichment.

The models can be generically classified into semi-analytic and hydrodynamic simulations. Nevertheless, there is increasing uncertainty connected to the chemical evolution as we move from local to cosmological scales, which are independent of the analytic or computational modeling. From the small scale represented by nuclear reaction rates and stellar masses to larger, galactic, and cosmological scales, there is a cumulative uncertainty since each scale carries its own set of considerations and uncertainties. For a specific discussion on this subject, see, in particular, Côté et al. (2016). These models help discuss general and particular aspects of metallicity evolution in cosmological terms, but several do not detail the contributions to the evolution of single elements. Moreover, comparing observations with high-redshift simulations is a challenge. Based upon the points mentioned above, we propose a semi-analytic model to investigate the contributions of Pop III and Pop II stars to the cosmological evolution of single elements across the redshift interval 0 ≤ z ≲ 20, not including the details of hydrodynamic simulations, and taking into account different perspectives to compare our results with observations in a range of different redshifts.

We start in Sect. 2, introducing and justifying the choices for the cosmological background, which is going to be the basis for the chemical evolution model: we describe the model developed by Pereira & Miranda (2010), and the incorporated changes in the scenario that allow for an adequate coupling of the star formation model with the equations of the chemical evolution of the Universe. We address the adapted cosmological model as Corazza, Miranda & Wuensche (hereinafter represented as CMW along with the text). In addition, the modifications of the model introduced in this work allow for a better adjustment of the cosmic star formation rate (CSFR) to the observational data available up to redshift ∼10, and satisfy all the points studied by Gribel et al. (2017) in their unified model connecting the CSFR with the local star formation.

We also adapted the chemical models developed over the past 40 years for the Galaxy (see, e.g., Tinsley & Larson 1978; Larson et al. 1980; Matteucci 2016). This adaptation allows us to build an adequate model for the chemical enrichment of the Universe. Implementing chemical yields for stars with masses between 0.85 and 260 M and metallicities from 0 up to Z = 0.02 Z allows us to provide, in Sect. 3, several results, data comparison, and discussions about the cosmic evolution of 11 chemical elements: iron (Fe), silicon (Si), zinc (Zn), nickel (Ni), phosphorus (P), magnesium (Mg), aluminum (Al), sulfur (S), carbon (C), nitrogen (N), and oxygen (O). We draw the conclusions in Sect. 4.

2. Methodology

In this section, we show how we obtain the CSFR from the process of large-scale structure formation and how we build a scenario that describes the cosmic chemical enrichment from redshift ∼20 to the present. The first DM halos decoupled from the Universe’s expansion, collapsing and virializing, probably between the end of recombination and redshift ∼20. The potential wells of these first halos generated the conditions for the baryonic matter to flow into these structures, agglomerate, and form the first stars. The characterization of the cosmological star formation and the consequent chemical enrichment of the Universe is, in this way, connected to the DM halo formation within a given mass range and as function of the redshift.

2.1. Cosmological scenario

Dark matter halos drag the baryonic matter into their interiors. We can describe this process through the adaptation of the formalism developed originally by Press & Schechter (1974), which allows us to directly estimate the fraction of baryons (fb) incorporated into the halos:

f b ( z ) = M min M max f ( M , z ) M d M 0 f ( M , z ) M d M , $$ \begin{aligned} f_{\rm b} (z) = \frac{\int _{M_{\rm min}}^{M_{\rm max}}\,f(M,z)\,M\,\mathrm{d}M}{\int _{0}^{\infty } f(M,z)\,M\,\mathrm{d}M}, \end{aligned} $$(1)

where

d f ( M , z ) = ρ m M d ln σ 1 d M f ST ( σ ) d M $$ \begin{aligned} \mathrm{d}f(M,z) = \frac{\rho _{\rm m}}{M}\frac{\mathrm{d}\ln {\sigma ^{-1}}}{\mathrm{d}M}f_{\rm ST}(\sigma )\, \mathrm{d}M \end{aligned} $$(2)

is the number of DM halos per comoving volume at a given redshift within the mass interval [M, M + dM], and ρm is the matter density of the Universe.

The halo mass function, fST(σ), proposed by Sheth & Tormen (1999), is:

f ST ( σ ) = 0.3222 2 a π δ c σ exp ( a δ c 2 2 σ 2 ) [ 1 + ( σ 2 a δ c 2 ) p ] , $$ \begin{aligned} f_{\rm ST}(\sigma ) = 0.3222 \sqrt{\frac{2a}{\pi }}\frac{\delta _{\rm c}}{\sigma } \exp \left(\frac{-a\delta _{\rm c}^{2}}{2\sigma ^{2}}\right) \left[1 + \left(\frac{\sigma ^{2}}{a\delta _{\rm c}^{2}}\right)^{p}\right], \end{aligned} $$(3)

with δc = 1.686, a = 0.707, p = 0.3, and σ(M, z) is the variance of the linear density field.

The fact that stars form only in suitably dense structures is parameterized in Eq. (1) by the threshold mass Mmin. We consider Mmin = 106M to be the minimum mass for the first star-forming halos to appear in hierarchical models. The upper limit Mmax can take values up to ≳1017M. This limit is set according to the mass scale of galaxy superclusters (Salvadori et al. 2007; Pereira & Miranda 2010), limiting the scale of the largest structures formed in the present Universe. In any case, the results have shown to be weakly dependent on the upper limit if Mmax >  1017M.

The function σ(M, z) in Eq. (3) can be determined from the power spectrum P(k) smoothed with a spherical top-hat filter function of radius R, which, on average, encloses a mass M (R = [3M/4πρ(z)]1/3). Thus,

σ 2 ( M , z ) = D 2 ( z ) 2 π 2 0 k 2 P ( k ) W 2 ( k , M ) d k , $$ \begin{aligned} \sigma ^{2}(M,z) = \frac{D^{2}(z)}{2\pi ^{2}} \int _{0}^{\infty }{k^{2}\, P(k)\, W^{2}(k,M)\,\mathrm{d}k}, \end{aligned} $$(4)

where W(k, M) is the top-hat filter in the k-space:

W ( k , M ) = 3 ( k R ) 3 [ sin ( k R ) k R cos ( k R ) ] . $$ \begin{aligned} W(k,M) = \frac{3}{(kR)^{3}}\,[\sin (kR)-k R\cos (kR)]. \end{aligned} $$(5)

The dependence on redshift comes from the growth factor D(z), that is, σ(M, z)=σ(M, 0)D(z). Here we use the analytical approach for D(z) as derived by Carroll et al. (1992).

The rate at which fluctuations grow on different scales depends on the interplay between self-gravitation, pressure support, and damping processes. All of these processes are part of the power spectrum given by P(k)∝knp (see, e.g., Gribel et al. 2017 for details).

From these equations, it is possible to determine how halos of different masses decouple from the Universe’s expansion and how baryonic matter is gradually incorporated into the center of the virialized halos. Structures more massive than ∼106M are formed at later times, as the redshift decreases. Thus, as more halos are formed, more baryonic matter flows into these structures, generating conditions for star formation. This allows us to define the baryon accretion rate as:

a b ( t ) = Ω 0 , b ρ 0 , c ( d t d z ) 1 | d f b d z | , $$ \begin{aligned} a_{\rm b}(t) = \Omega _{0,\mathrm{b} }\,\rho _{0,\mathrm{c} } \left(\frac{\mathrm{d}t}{\mathrm{d}z}\right)^{-1} \left| \frac{\mathrm{d}f_{\rm b}}{\mathrm{d}z} \right|, \end{aligned} $$(6)

where Ω0, b is the baryonic density parameter at z = 0, ρ 0 , c = 3 H 0 2 / 8 π G $ \rho_{0,\mathrm{c}} = 3H_{0}^{2}/8\pi G $ is the critical density of the Universe (H0 = 100 h km s−1 Mpc−1 is the value of the Hubble parameter at the current time), and

d t d z = 1 H 0 ( 1 + z ) Ω 0 , Λ + Ω 0 , m ( 1 + z ) 3 · $$ \begin{aligned} \frac{\mathrm{d}t}{\mathrm{d}z} = \frac{1}{H_{0}\, (1+z)\sqrt{\Omega _{0,\Lambda } + \Omega _{0,\mathrm{m} }\,(1+z)^{3}}}\cdot \end{aligned} $$(7)

The cosmological framework described through the set of equations presented above is similar to the one used by different authors (e.g., Daigne et al. 2006; Pereira & Miranda 2010; Tan et al. 2016; Gribel et al. 2017; Vangioni et al. 2018).

2.2. Cosmic star formation rate

Once we set the cosmological framework, it is possible to compute the CSFR by incorporating the IMF and the star formation rate (SFR). In particular, the number of stars formed per unit of mass (m), volume (V), and time (t) is given by

d 3 N ( m , V , t ) d m d V d t = φ ( m ) ψ ( t ) , $$ \begin{aligned} \frac{d^{3} N(m, V, t)}{\mathrm{d}m\mathrm{d}V\mathrm{d}t} = \varphi (m)\, \psi (t), \end{aligned} $$(8)

where ψ ( t ) ρ g α $ \psi(t) \propto \rho_{\mathrm{g}}^{\alpha} $ corresponds to the SFR (ρg is the gas density). It should be noted that ψ(t) follows the functional form known as Schmidt’s law (Schmidt 1959). On the other hand, the IMF is given by φ(m)∝m−(1 + x) and its functional form with x = 1.35 is called as Sapeter’s IMF (Salpeter 1959).

The IMF of the first stars is still an open question. We see that the Salpeter IMF favors the formation of low-mass stars, and various authors adopted it in their chemical evolution models (see, e.g., Calura & Matteucci 2004, 2006; Casey et al. 2012; Shu et al. 2016; Fraser et al. 2017; Vangioni et al. 2018) while some others (see, e.g., Nakamura & Umemura 2001; Schneider et al. 2006; Ma et al. 2017) also allow for the possibility of a top-heavy or bi-modal IMF.

In our study, we consider x = 1.35 as the reference value. However, to identify the influence of the IMF exponent on the chemical enrichment of the Universe, we also considered four other values, nominally, 0.85, 1.0, 1.7, and 1.85 allowing the formation of a higher (the first two values) or smaller (the last two values) number of high-mass stars when compared to the reference value 1.35. We also used α = 1 in agreement with Gribel et al. (2017), which shows that different properties from the star formation regions in the Galaxy, including the so-called Larson’s law, can be well reproduced with α = 1.

Therefore, Eq. (8) describes the number of stars formed within the DM halos that aggregate and concentrate baryons in their centers. A fraction of the mass in stars is ejected (through stellar winds and supernovae, for example) and returned to the ISM formed by these structures. The ejected mass fraction is given by:

d 2 M ej d V d t = m ( t ) m s ( m m r ) ψ ( t τ m ) φ ( m ) d m , $$ \begin{aligned} \frac{d^{2}M_{\rm ej}}{\mathrm{d}V\mathrm{d}t} = \int ^{m_{\rm s}}_{m{(t)}} (m-m_{\rm r})\, \psi (t-\tau _{\rm m})\,\varphi (m)\,\mathrm{d}m, \end{aligned} $$(9)

where m(t) is the stellar mass whose lifetime is equal to t, and mr represents the mass of the remnant, which depends both on the progenitor mass (m) and on the environment metallicity (Z). The star formation is taken at the time (t − τm), where τm, also a function of the metallicity (Z), is the lifetime of a star of mass m.

We used the results of Spera et al. (2015) to obtain the masses of the stellar remnants (mr) as functions of the metallicity and the initial stellar masses. The authors obtained their results from SEVN (population-synthesis code) coupled with the PARSEC code for stellar evolution tracks. In particular, in this work we use the fitting formulas presented in Appendix C of Spera et al. (2015).

Concerning the parameter τm, we use the metallicity-dependent formula given by Raiteri et al. (1996):

log τ m = a 0 ( Z ) + a 1 ( Z ) log ( M M ) + a 2 ( Z ) [ log ( M M ) ] 2 , $$ \begin{aligned} \log \,\tau _{\rm m} = a_{0}(Z) + a_{1}(Z)\log \left(\frac{M_{\star }}{M_{\odot }}\right) + a_{2}(Z) \left[\log \left(\frac{M_{\star }}{M_{\odot }}\right)\right]^{2}, \end{aligned} $$(10)

where τm is expressed in years, and the metallicity-dependent coefficients are (see Raiteri et al. 1996 for details):

a 0 ( Z ) = 10.13 + 0.07547 log Z 0.008084 ( log Z ) 2 , $$ \begin{aligned}&a_{0}(Z) = 10.13+0.07547\, \log Z - 0.008084\, (\log Z)^{2}, \end{aligned} $$(11)

a 1 ( Z ) = 4.424 0.7939 log Z 0.1187 ( log Z ) 2 , $$ \begin{aligned}&a_{1}(Z) = -4.424-0.7939\, \log Z - 0.1187\, (\log Z)^{2}, \end{aligned} $$(12)

a 2 ( Z ) = 1.262 + 0.3385 log Z + 0.05417 ( log Z ) 2 , $$ \begin{aligned}&a_{2}(Z) = 1.262+0.3385\, \log Z + 0.05417\, (\log Z)^{2}, \end{aligned} $$(13)

and Z is the absolute metallicity.

It is worth stressing that τm determined by Eq. (10) has an excellent agreement when compared to the stellar lifetimes presented in Table 2 of Ekström et al. (2008) for different values of mass and metallicity. In particular, the difference between the results for τm is lower than 5%, which has little effect on our results.

Following this formalism and combining the previous equations, we derive the equation that governs the total gas density ρg in the halos:

ρ ˙ g = d 2 M d V d t + d 2 M ej d V d t + a b ( t ) , $$ \begin{aligned} \dot{\rho }_{\rm g} = -\frac{d^{2}M_{\star }}{\mathrm{d}V\mathrm{d}t} + \frac{d^{2}M_{\rm ej}}{\mathrm{d}V\mathrm{d}t} + a_{\rm b}(t), \end{aligned} $$(14)

where the term ab(t) gives, for the halos, a matter of primordial composition. The system becomes closed without the term ab(t) in Eq. (14). Thus, this term corresponds to a primordial gas infall in the structures in formation. In other words, it describes the primordial baryonic matter that is captured by the potential wells generated by the halos.

On the other hand, the first term on the right side gives the mass of gas converted to stars per unit of volume and time. By Schmidt’s law, we have:

ψ ( t ) = d 2 M d V d t = k ρ g . $$ \begin{aligned} \psi (t) = \frac{d^{2}M_{\star }}{\mathrm{d}V\mathrm{d}t} = k\,\rho _{\rm g}. \end{aligned} $$(15)

It should be noted that the term k is the inverse of timescale for star formation, that is, k = 1/τs.

The total gas density can be calculated by numerical integration of Eq. (14), providing values for ρg(t) at each time t or redshift z as long as the τs parameter is set. The initial condition is zero gas density at t = 0 (z = 20) for solving Eq. (14). Moreover, there are some steps for obtaining the correct characterization of the function ρg. First of all, it is necessary to determine the function τs. This parameter is related to the CSFR via Schmidt’s Law, that is, ρ ˙ $ \dot\rho_{\star} $ is directly proportional to the gas density and inversely proportional to the characteristic timescale for the conversion of gas in stars. Second, if all the gas entering the system, plus the gas returning to the system through Eq. (9), is converted into stars, there is an overabundance of both stars and metals. Thus, ρ ˙ $ \dot\rho_{\star} $ must also be dependent on a parameter that measures the efficiency (⟨ε⟩) in which gas is converted into stars.

From the above considerations, we have:

ρ ˙ ( z ) = ε ρ g τ s · $$ \begin{aligned} \dot{\rho }_{\star } (z) = \langle \varepsilon _{\star }\rangle \frac{\rho _{\rm g}}{\tau _{\rm s}}\cdot \end{aligned} $$(16)

Once the CSFR is determined, the ab function is fixed by the structure formation scenario, and the IMF is determined by the choice of the x exponent, and then it will be possible to determine the function ρg by Eq. (14).

The aforementioned steps are essential to characterize the total gas density function. The calculation algorithm integrates the differential Eq. (14) through the sixth-order Runge-Kutta method. The differential equations for the various chemical elements and total metallicity of the Universe (Sect. 2.3) are solved by the same method.

2.2.1. Characterizing the functions ⟨ε⟩ and τs

The ⟨ε⟩ and τs functions work together to produce the CSFR with the best fit to the observational data. In particular, the cold gas used to form stars is given by:

ρ cold ( z ) = ε ρ g ( z ) , $$ \begin{aligned} \rho _{\rm cold}(z) = \langle \varepsilon _{\star }\rangle \,{\rho _{\rm g}(z)}, \end{aligned} $$(17)

where ⟨ε⟩ acts as efficiency for star formation. In principle, ρcold is composed of the sum of two components: molecular gas (ρH2) and atomic gas (ρHI). This allows us to rewrite τs as:

τ s ( z ) = ρ cold ρ ˙ = ρ H 2 ρ ˙ + ρ HI ρ ˙ = τ depl , H 2 + τ depl , HI , $$ \begin{aligned} \tau _{\rm s} (z) = \frac{\rho _{\rm cold}}{\dot{\rho }_{\star }} = \frac{\rho _{\rm H2}}{\dot{\rho }_{\star }} + \frac{\rho _{\rm HI}}{\dot{\rho }_{\star }} = \tau _{\rm depl,H2} + \tau _{\rm depl,HI}, \end{aligned} $$(18)

with τdepl, H2 and τdepl, HI representing, respectively, the depletion scales for molecular and atomic gases.

Equations (16)–(18) are solved together to characterize the CSFR and to determine the correct dependence of the ⟨ε⟩ and τs functions on the redshift. The constraints associated with these equations are the following: the best adjustment of the ρ ˙ ( z ) $ \dot\rho_{\star}(z) $ curve to the observational data available within the range [0−10] in redshift must be produced; ρ ˙ $ \dot\rho_{\star} $ should be normalize to return the value ∼0.016 M yr−1 Mpc−3 at z = 0, similar to the one determined by Madau & Dickinson (2014); the CSFR peak should be made at redshift z = 2. This value was chosen so that ρ ˙ $ \dot\rho_{\star} $ obtained here is in accordance with the peak of the CSFR used by Vangioni et al. (2018) and the one determined by Madau & Dickinson (2014). Furthermore, ⟨ε⟩∼0.01−0.02 should be produced at z = 0. This causes the value ⟨ε⟩ to be on the order of εff, the so-called SFR per free-fall time, inferred for the star-forming regions of the local Universe (see, e.g., Krumholz & McKee 2005; Gribel et al. 2017). Lastly, τs(z = 0)∼0.5−2.5 Gyr should be produced, similar to the value of τdep as inferred by Schinnerer et al. (2016) for the local Universe (z = 0).

Figure 1 shows the CSFR as a function of redshift and its behavior concerning the observational data. For comparison, we also present the CSFR used by Vangioni et al. (2018). Both our CSFR and the one used by Vangioni et al. (2018) reach a maximum value at z = 2, as mentioned above. Up to redshift ∼4, the two CSFRs have a very similar behavior. We fit both CSFRs to results from IR, UV, and GRB observations, as described in the figure caption.

thumbnail Fig. 1.

Solution for ρ ˙ ( z ) $ \dot\rho_{\star}(z) $ as derived in this work (red line) and, using black dashed line, the CSFR used by Vangioni et al. (2018) plotted for comparison. Left: evolution of the CSFR from the local Universe to z = 20. Right: same as for the left panel, but zooming into 0 ≤ z ≤ 10, allowing for a better visualization of the two CSFRs within the range of the available observational data. Data used in this figure: IR (Magnelli et al. 2011, 2013 – dark gray filled circles; Gruppioni et al. 2013 – dark orange open circles); UV (Wyder et al. 2005 – slate gray open triangle; Schiminovich et al. 2005 – blue crosses; Dahlen et al. 2007 – turquoise open squares; Reddy & Steidel 2009 – dark green filled squares; Robotham & Driver 2011 – chocolate cross; Cucciati et al. 2012 – green stars; Bouwens et al. 2012a, 2012b – teal filled triangles); and GRB (Kistler et al. 2009 – deep pink open triangles).

In Fig. 2 we present the behavior of the τs parameter for different values of the x exponent. Additionally, as determined by Péroux & Howk (2020) and Maio et al. (2022), the corresponding H2 depletion time closely follows the dynamical time (τdyn), taken to be 10% of Hubble time. The curves for τs and τdyn ∼ τdepl, H2 are shown in the left panel. It should be noted that the values for τs are above τdyn for higher redshifts, which means that the largest contribution to the characteristic timescale for star formation, in Eq. (18), is due to the atomic gas instead of molecular gas.

thumbnail Fig. 2.

Evolution of the characteristic timescale for star formation as a function of the redshift. Left: τs for different IMF exponent and the evolution of the dynamical time (tdyn), taken to be 10% of Hubble time. Right: our results compared with the timescales for the conversion of gas (filled areas) obtained by Péroux & Howk (2020) and Maio et al. (2022) for molecular and cold (HI + H2) gas.

At lower redshifts (z <  5), the curve τdyn gradually approaches τs, implying a greater contribution of molecular gas to star formation. We note that for the two IMFs we have τdyn >  τs, which happens at redshift ∼0.75 (∼0.38) to x = 0.85 (1.0). In these cases, the approximation τdyn ≃ thubble/10 is not adequate to describe the characteristic timescale for star formation. As a result, τdepl, H2 should be constant. For the other IMFs, τs is dominated by the molecular depletion time for z <  0.5.

The right panel of Fig. 2 shows our results compared with the timescales for the conversion of gas (filled areas) obtained by Péroux & Howk (2020) and Maio et al. (2022) for molecular and cold (HI + H2) gas. The behavior of τs between the two filled areas shows that, in our scenario, star formation is fueled by cold gas in atomic form for z >  6. When the curves describing τs approach the gray filled area, the contribution of molecular gas becomes gradually dominant to supply star formation.

As commented above, for IMF exponents 0.85 and 1.0, the approximation τdyn ≃ thubble/10 fails for z <  1 and τdyn should be constant in order to supply the star formation with gas in molecular form. Our model shows that for x in the range 1.35−1.85, about 70% to 90% of the star formation would come from molecular gas, for z ≲ 1, if τdepl, H2 = τdyn ≃ thubble/10.

The left panel of Fig. 3 shows the evolution of star formation efficiency (⟨ε⟩) with redshift. This parameter reaches values between ∼0.15−0.40, depending on the value of x for redshifts z >  3. Efficiency gradually decreases with redshift reaching values close to 0.01−0.02 at z = 0. In the right panel, we present the fraction of total cold gas mass (fgas) determined by the relation:

thumbnail Fig. 3.

Evolution of the fraction of gas converted into stars for different values of the IMF used in this work. Left: evolution of the average star formation efficiency (defined as ⟨ε⟩) as a function of z. Right: evolution of the fraction of total cold gas associated with star formation (fgas). The symbols show the measurements compiled as follows: circles by Tacconi et al. (2018), squares by Scoville et al. (2014), diamonds by Scoville et al. (2016), stars by Dessauges-Zavadsky et al. (2015), and triangles by Schinnerer et al. (2016). See also Fig. 13b of Hodge & da Cunha (2020).

f gas = M gas M gas + M , $$ \begin{aligned} f_{\rm gas} = \frac{M_{\rm gas}}{M_{\rm gas} + M_{\star }}, \end{aligned} $$(19)

where Mgas is the cold gas mass and M is the mass in stars.

Equation (19) permits a direct comparison with results obtained by other authors. In particular, Hodge & da Cunha (2020) present in their Fig. 13b measurements compiled for the cold gas fraction (H2 + HI) and compared with scaling relations derived in three different works (Scoville et al. 2017; Tacconi et al. 2018; Liu et al. 2019). Our curves for the cold gas fraction show good agreement with the scaling relation M gas M 0.65 $ M_{\mathrm{gas}} \sim M_{\star}^{0.65} $ derived by Tacconi et al. (2018).

In Table 1, we summarize the parameters used to obtain the CSFR – ρ ˙ ( z ) $ \dot\rho_{\star}(z) $. It depends on the cosmological parameters Ω0, m, Ω0, b, and Ω0, Λ, and the parameters related to the formation of large-scale structures of the Universe (σ8, np, and Mmin).

Table 1.

Cosmological and structure formation parameters used to obtain the CSFR.

The CSFR used by Vangioni et al. (2018) follows the expression originally derived by Springel & Hernquist (2003), but changing the values of the four free parameters to adjust ρ ˙ ( z ) $ \dot\rho_{\star}(z) $ to the most recent observational data. In this work, we include the available GRB data at high redshifts since the UV data will naturally suffer from the selection effect; only the brightest objects are observed.

Although our model is semi-analytic, by adding the redshift dependence to the functions ⟨ε⟩ and τs, it becomes possible to obtain ρ ˙ ( z ) $ \dot\rho_{\star}(z) $ with an adequate behavior within the redshift range where the CSFR data exist. In addition, the way we build ⟨ε⟩ and τs with the constraints that these functions must satisfy at z = 0, shows good agreement with the observational data. In particular, the ratio ⟨ε(z)⟩/⟨ε(z = 0)⟩ provided by our model is in good agreement with that obtained by Scoville et al. (2017) within the redshift range [0−3].

2.3. Chemical evolution scenario

The first chemical evolution models were developed for the framework of the Galaxy by Tinsley & Larson (1978), Larson et al. (1980), and later by Matteucci (2001). Their simple model of chemical evolution considers a closed-box evolving system with no inflows or outflows. Also, the IMF is constant in time, the chemical composition of the gas is primordial, and the mixing between the chemical products ejected by stars and the ISM is instantaneous.

We can adapt these concepts, which are the basis of the chemical evolution models of the Galaxy, straightforwardly. The main difference is that in the cosmological scenario, the halos continuously incorporate baryons (primordial gas) from the ambient (Universe). This is described by the function ab(t).

Once inside the halos, the gas is removed from the system to form stars at time t − τm. This is described using the CSFR ρ ˙ ( t τ m ) $ \dot\rho_{\star}(t-\tau_{\mathrm{m}}) $. Later, the gas returns to the system, at time t, when the stars die. A certain fraction of the gas used for the star formation in t − τm will be retained in the remnant population mr that forms in the time t. A new generation of stars will be formed at the instant t, removing gas from the system, and this processes is repeated in a cycle of continuous gas capture and chemical enrichment of the environment.

To determine the chemical enrichment of a given i-element, in addition to the functions ρ ˙ $ \dot\rho_{\star} $ and ab, we need to know how much mass of the i-element is returned when the star of mass m dies. This is described by the parameter PZi m that provides the “stellar yield” of the i-element.

Once all of these functions and parameters are characterized, we can write a differential equation for the mass density of the i-element as:

d ρ g i d t = m ( t ) m s [ ( m m r ) Z i ( t τ m ) + P Z i m ] ρ ˙ ( t τ m ) φ ( m ) d m Z i ρ ˙ ( t ) , $$ \begin{aligned} \frac{\mathrm{d}\rho _{\rm g\,i}}{\mathrm{d}t} =&\int _{m(t)}^{m_{\rm s}} \,[(m - m_{\rm r})\,Z_{\rm i}\,(t - \tau _{\rm m})\nonumber \\&+ P_{Z_{\rm i\,m}}]\, \dot{\rho }_{\star } (t - \tau _{\rm m})\varphi (m)\,\mathrm{d}m - Z_{\rm i}\,\dot{\rho }_{\star } (t), \end{aligned} $$(20)

where the term (m − mr) Zi(t − τm) accounts for the amount of i-element incorporated when the star was born, and which later returns to the ISM (we see that mrZi (t − τm) is the part of the i-element retained into the remnant). We note that the resulting ρg is dependent on ab(t) through Eq. (14), and thus the term Zi = ρg i/ρg takes into account ab(t) through ρg. The PZi m parameter is the mass produced of the i-element by a star of mass m. The term Z i ρ ˙ ( t ) $ Z_{\mathrm{i}}\,\dot\rho_{\star}(t) $ takes into account the removal of part of the i-element to form a new star generation.

Through the time integration of Eq. (20), we obtain the mass density ρg i of the i-element present in the gas contained within the halos. This allows us to determine quantities such as [Xi/H] as a function of redshift (or time) and to compare the results of our model with different observational data. It should be noted that Eq. (20) incorporates all the physics and constraints discussed in the previous sections.

In order to incorporate the contributions of particular stars, depending on their masses and metallicities, we selected stellar yields. These chemical yields are used to determine the elements that were ejected into the ISM at a given time by a star of a given mass and metallicity. They are calculated through detailed nucleosynthesis computational simulations, considering the main reactions that happen inside the stars. We consider the first stars to be zero-metallicity stars (Pop III); the subsequent more enriched Pop II stars are chosen within a range of different masses and metallicities. Tables 2 and 3 describe the stellar mass and metallicity ranges from where the chemical yields were chosen.

Table 2.

Masses selected for Pop III chemical yields.

Table 3.

Masses and metallicities selected for Pop II chemical yields.

Properly modeling chemical yields for Pop III stars is a challenging and complex task. For the range where they become PISNe, we chose to work with results from Heger & Woosley (2002), which are compatible with recent chemical yields calculated by Takahashi et al. (2018, hereafter TK18), which take into account rotating progenitors. The two models, HW02 and TK18, show no significant differences in the explosive yields for the elements chosen here, except for the large production of N in the TK18 nonmagnetic rotating models. The N behavior can be better understood in a detailed, recently developed model for the cosmological evolution of this element (Vangioni et al. 2018).

For Pop II stars, samples were chosen according to the best combination of mass and metallicity ranges, and also according to the stellar evolution models and parameters used to produce each sample; Karakas (2010) and Doherty et al. (2013, 2014) used the MONSTAR code for stellar evolution (Frost & Lattanzio 1996), OPAL opacities (Iglesias & Rogers 1996), and compatible mass-loss models (Reimers 1975; Vassiliadis & Wood 1993; Bloecker 1995).

The chemical evolution to be presented in the following sections makes use of the following definition:

[ X i / H ] = log 10 [ N ( X i ) / N ( H ) ] gas log 10 [ N ( X i ) / N ( H ) ] , $$ \begin{aligned}{\ }[\mathrm X_{i}/H ] = \log _{10}[\mathrm N(X_{i})/N(H) ]_{\rm gas} - \log _{10}[\mathrm N(X_{i})/N(H) ]_{\odot }, \end{aligned} $$(21)

where N(Xi) and N(H) are respectively the densities for the element Xi and for hydrogen.

2.3.1. Numerical technique to determine the elements ρg i and ρg

The solution for the set of differential equations for ρg and ρg i is obtained by a sixth-order Runge-Kutta algorithm. As an initial condition, we have zero values for the densities of gas and metals. To integrate the Pop III and Pop II chemical yields, we use a cubic spline interpolation algorithm for each dataset in Tables 2 and 3. For example, for the HW10 model (Pop III in Table 2) with a mass range of 10−100 M, we have 18 stellar mass values, nominally, 10, 12, 15, 17, 20, 22, 25, 27, 30, 33, 35, 40, 45, 50, 60, 75, 85, and 100 solar masses. Thus, the determination of the yields of a star with, for example, 70 M, occurs through the interpolation algorithm.

We see that there are no yields in the intervals ]3.0, 10[ and ]100, 140[ (Table 2). Thus, the interpolation algorithm does not include these open sets in the calculation. A similar implementation applies to Pop II yields (Table 3).

The total density of the gas depends on the term ab(t), which corresponds to the infall of primordial gas, basically hydrogen and helium. To check the consistency of our results, we determined the maximum numerical discrepancy (D) between ρg and the sum over all chemical species (hydrogen, helium, and metals), which is:

D = ρ g ρ H ρ He ρ g i ρ g · $$ \begin{aligned} D = \frac{\rho _{\rm g}-\rho _{\rm H} - \rho _{\rm He}-\sum \,\rho _{\rm g\,i}}{\rho _{\rm g}}\cdot \end{aligned} $$(22)

The maximum numerical discrepancy in the various time steps is |D| = 4.4 × 10−7.

2.4. Constructing the models

We start the calculation with all the elements for the cosmological scenario and chemical evolution set, and explore two scenarios. The first (model A) considers the chemical evolution of the Universe starting with the Pop III (masses and yields as shown in Table 2) stars. Once the metallicity of the Universe reaches Z = 10−6Z, no more Pop III stars can be formed. As a consequence, the Pop II stellar branch with Z = 10−6Z is born and evolves (masses and yields as shown in Table 3). This second step finishes when the metallicity of the Universe reaches Z = 10−4Z and, as a consequence, no more Pop II stars can be formed within the branch Z = 10−6Z. This process repeats every time the metallicity of the Universe crosses the limits indicated in Tables 2 and 3. It should be noted that the stars within the low-metallicity branches cannot form anymore as a consequence of the increase in the chemical enrichment of the Universe. However, stars with low mass born within the lowest-metal branch can still be alive now. Thus, they can coexist with stars of much higher metallicity during part of their lives.

For the second scenario (model B), we consider the chemical evolution only with Pop II stars. In this case, the first-star generation (Z = 0) of the Universe was composed of stars with masses and chemical yields similar to the Pop II stars of the branch Z = 10−6Z studied by Chieffi & Limongi (2004). The following steps are similar to those described for model A. Both scenarios are generated with the CSFR described in Sect. 2.1.

An important observation is that models A and B have different normalizations for the IMF. For model A, the normalization is obtained in the range 0.1−260 M, while model B (Pop II only) is normalized in the interval 0.1−120 M. This choice is due to the nonexistence of the stellar branch 140−260 M in the Pop II models discussed in the literature.

We describe the evolution of models A and B in Sects. 2.4.1 and 2.4.2.

2.4.1. Model A

Model A runs the following steps for the entire calculation. First, the total metallicity of the Universe (Ztotal) is used as a guide for the chemical yields of the different classes (or branches) presented in Tables 2 and 3. Assuming that the first stars, formed from pristine (H and He only) gas, started to die and enrich the ISM at redshift z = 20, Ztotal provides values for the evolution of the production of all elements heavier than He for the entire redshift interval. This parameter is then used as a “switch” between different metallicities, removing the chemical yields of a given metallicity and successively introducing those of the higher metallicity classes according to Tables 2 and 3, and as discussed above.

Abundances of individual chemical elements are then computed for the entire redshift range. Ztotal starts at zero, producing Pop III stars. The higher-mass Pop III stars (∼260 M) start dying first, throwing metals into the ISM and enriching the medium around them. Pop III stars continue to die and enrich the ISM until the medium reaches 10−6Z. At that point, new Pop III stars stop forming. It is important to note that Pop III stars with masses ≲0.9 M have longer lifetimes and should be still in their main-sequence phase today, just following the increase in the total metallicity of the Universe, and will participate in subsequent steps of the enrichment of the Universe when they leave the main sequence, along with the contribution of new, higher-metallicity stars formed later than that population, at much lower redshifts.

Once the metallicity of the ISM and/or IGM reaches 10−6Z, new stars with this metallicity signature start being formed. As they die, the model starts processing the yields from this class of stars, and the same process repeats. The metallicity from the medium increases as the higher-mass stars die first, while lower-mass stars live longer. Even with stars with a higher metallicity (10−4Z, for example) starting to form, the lower-mass ones will continue their lives unaffected by the external increase in metallicity.

The process continues as the Universe progresses toward the present metallicity, as described in Tables 2 and 3. It is important to emphasize that each metallicity class has its own stellar population clock triggered when the Universe’s metallicity crosses the various thresholds indicated in Tables 2 and 3. Thus, at certain intervals of time, the chemical enrichment of the Universe takes place by the joint action of stars of different metallicity classes.

2.4.2. Model B

Model B uses only Pop II yields. The chemical enrichment of the Universe starts from Ztotal = 0 at redshift 20. For the second model, we consider that the stellar branch Z = 10−6Z also represents the metal-free stars. Thus, chemical enrichment will occur through this class (or branch) until reaching Z = 10−4, when then the next metallicity class (with Z = 10−3Z) will assume the control of chemical enrichment. The following branches will take the enrichment command at the points indicated in Table 3.

Chemical elements analyzed in this work were selected considering the availability of chemical yields in the literature and observational data available for chemical abundances in damped Lyman-α systems (DLAs) and globular clusters (GCs), as described below (Sect. 3).

3. Results and discussion

The results from models A and B for [Ztotal/H] are presented in Fig. 4. The results are obtained from the CMW-CSFR with five different IMF values and integrated for all ranges from Pop III to Pop II stars, for model A, and for all ranges of Pop II stars, for model B. We take into account the progressive enrichment of the Universe, and consequently the transition between Pop III stars and the next, more metal-rich Pop II generations until Z = 0.02 Z.

thumbnail Fig. 4.

Evolution of metallicity for different values of x. Left: model A. Right: model B. Colors represent different values of x for the IMF.

For model A, at redshift z = 20, the first stars formed from metal-free gas start to die, and the chemical enrichment is very fast. For x = 0.85 and x = 1.00, the pristine Universe leaves from Z = 0 to reach Z = 10−6Z in less than ∼4 × 105 yr, given the higher number of high-mass stars that would form in this scenario. For x = 1.70, the same metallicity is reached ∼30 Myr after the death of the first Pop III star, while for x = 1.85, it would take ∼70 Myr for the same process to occur. The mean behavior is described by x = 1.35, where the Universe would reach Z = 10−6Z in ∼3 Myr. For model B, the same process takes from ∼2 Myr up to 85 Myr, depending on the IMF.

This rapid chemical enrichment in the initial phase can be explained by the metal production of Pop III–PISNe, which characterizes a chemical “flood” in the high-redshift Universe, in the case of model A. For model B, the chemical enrichment occurs mainly through stars with masses ∼30−35 M, and the condition Z = 10−6Z is reached eight times slower than when considering higher-mass Pop III stars.

Except for N, all other elements are mainly produced by PISNe in the Pop III era. According to Abia et al. (2001), the metallicity observed at high redshifts can be easily obtained from stellar pregalactic (Pop III) nucleosynthesis by postulating that only ∼10−2 of the total pristine gas is converted into stars. Considering that the star formation efficiency of the CMW scenario is ∼0.3 in the redshift range [5−20], which is ∼30 times larger than the value estimated by Abia et al. (2001), we verify that adding Pop III stars in the CMW scenario for the CSFR quickly floods the primordial Universe with metals.

Other evidence that PISNe are very efficient in enriching the ISM comes from the work of Matteucci & Calura (2005), where they show that only 110 to 115 PISNe would be needed to enrich a cubic megaparsec of the IGM to Z = 10−4Z (with the index of the IMF varying between 1.35 and 0.5).

In our model the rate for PISNe can be calculated using the relation:

R PISNe = M ˙ ( t ) M PISNe × 140 260 ϕ ( m ) m d m , $$ \begin{aligned} R_{\rm PISNe} = \frac{\dot{M}_{\star }(t)}{\langle M_{\rm PISNe}\rangle } \times \int _{140}^{260} \phi (m)\, m\, \mathrm{d}m, \end{aligned} $$(23)

where ⟨MPISNe⟩ is the average mass of the stars that ended their lives as PISNe, and (t) is the CSFR in M yr−1 when Z = 10−6Z.

If we consider the average mass of PISNe as ∼200 M, then RPISNe ∼ 6 × 10−5 yr−1 or 1 PISNe every ∼16 000 yr. The number of PISNe needed to enrich the Universe from Z = 0 to Z = 10−6Z is ∼175 for our models.

Similar to Matteucci & Calura (2005), our work shows the importance of PISNe for rapidly enriching the ISM. Matteucci & Calura (2005) find that ∼110−115 would be needed to enrich a cubic megaparsec of the IGM to ∼10−4Z. Our results show that ∼175 PISNe are needed to enrich the medium to ∼10−6Z. However, our results involve a volume equivalent to ∼105 Mpc3. Thus, the numerical comparison is not so direct between the two works.

In order to adequately address the question of the rapid contamination of the early Universe, we propose comparing the results with abundances from old GCs. Such GCs (with ages close to the age of the Universe) present an opportunity to explore the chemical and physical conditions of the earliest star-forming environments in the Universe (Dotter et al. 2011; Frebel & Norris 2015), in other words, they should present a metallicity value very similar to the Universe’s mean metallicity at the time they were formed.

Analysis of old GCs’ data can provide information about the age and metallicity for the entire cluster, enabling better estimations than for isolated, metal-poor stars, for example. Figure 5 shows the behavior of Ztotal for the two models compared with observations from GCs. Model A accounts for the majority of observations, regardless of the IMF, while model B is unable to fit the metal abundances for Ztotal, even for the lower values of x. Although Pop II stars with low metallicity and large masses do indeed play an essential role in the first steps of cosmic enrichment, it is clear that their contribution alone is insufficient to allow the ISM to maintain efficient star formation in order to reach the observed metal abundances along the cosmic history.

thumbnail Fig. 5.

Comparison between models A (left) and B (right) with data from GCs from Frebel et al. (2007), Bond et al. (2013), Cowan et al. (2002), Cayrel et al. (2001), Dotter et al. (2011), and Wagner-Kaiser et al. (2017).

Another analysis is performed for the interval z = [0−6], where we compare the results with dust-corrected abundances from DLAs (Fig. 6). The DLAs provide the most accurate measurements of chemical abundances on the gas-phase for the high-redshift Universe (Wolfe et al. 2005), and abundances can be determined with errors ≤0.1 dex (Vladilo 2002). The DLAs are also the perfect site for the initial stages of gas cooling and star formation (Maio & Tescari 2015); they dominate the neutral gas content of the Universe in the redshift interval z = [0−5], and therefore are the most crucial neutral gas reservoir for star formation.

thumbnail Fig. 6.

Comparison of models A (left) and B (right) with data from dust-corrected DLAs (gray crosses) from De Cia et al. (2018). We include an α-enhancement correction, with [Z/H]=[Fe/H]+0.3 dex, as suggested by Rafelski et al. (2012).

In this redshift interval, regardless of the observations, we expect to see an increase in [Ztotal/H] with decreasing redshift, with the total metallicity reaching values close to ∼0 (solar) near redshift z ∼ 0. This behavior is consistent with observations presented in similar contexts (Fynbo et al. 2006; Davé & Oppenheimer 2007; Kobayashi et al. 2007; Rollinde et al. 2009; Vangioni et al. 2018). Also, as pointed out by Calura & Matteucci (2004), main metal production in spirals and irregulars is always increasing with time.

It is possible to observe that, for model B, for x <  1.35, the abundances reach a maximum at redshift z ∼ 4 and start decreasing toward z = 0, while for x ≥ 1.35, metallicities tend to rise with decreasing redshift. Nevertheless, all the B models remain between −1.0 and −1.5 dex lower than the expected value for z ∼ 0.

We also note that Pop II models with x <  1.35 produce more high-mass stars returning more metals through the CL04 channel. For these IMFs, a bottleneck occurs when the metallicity of the system reaches 8 × 10−3Z since there is no contribution to chemical enrichment through the CL04 channel. This explains the flat behavior verified in Fig. 5.

On the other hand, when taking into account Pop III stars (model A), metallicities increase continuously as redshift approaches z ∼ 0. For x = 0.85, x = 1.00, and x = 1.35, models reach [Ztotal/H] close to 0, as expected, while for x = 1.70 and 1.85, the total metallicity is underestimated by approximately 0.25 to 0.30 dex.

When comparing results with DLA observations, there are two main problems that are relevant to the interpretation of our results. The first is the high dispersion between points relative to the same (or very close) redshift. There are a variety of models that investigate dispersion in DLAs (see, e.g., Dvorkin et al. 2015), and some authors agree that it happens due to peculiar nucleosynthetic signatures from each system and also due to different star formation histories (Centurion et al. 1998; Pettini et al. 2000; Dessauges-Zavadsky et al. 2002; De Cia et al. 2016), which leads to the production of different amounts of each chemical element.

According to De Cia et al. (2016), regardless of the star formation history, the availability of refractory metals in the ISM is a crucial driver of dust production, and DLA galaxies may have a wide range of star formation histories, which in principle are also different from those of the Galaxy (De Cia et al. 2016). Therefore, we plot the mean value between the models with different IMFs, suggesting that results for x = 0.85 represent the upper limit (due to the favorable formation of high-mass stars in this model) and x = 1.85 as the lower limit (due to the favorable formation of lower-mass stars).

The second problem in comparing results with DLA observations relates to dust depletion. Some chemical elements react with different species, forming molecular compounds that can get trapped on the surface of dust grains and cannot be detected by the observations of abundances in the gas phase, that is to say, abundances would look lower than their actual values. The majority of the results indicate that the behavior of dust depletion on DLAs is complex and varies from system to system (Vladilo et al. 2011; De Cia et al. 2016, 2018).

Therefore, we compare our results with dust-corrected DLA metallicities from De Cia et al. (2018). The author shows that, when including dust corrections, the average DLA metallicities are between 0.4 and 0.5 dex higher than without corrections. Where the author provides values for [Fe/H], we include an α-enhancement correction, with [Z/H] = [Fe/H] + 0.3 dex as suggested by Rafelski et al. (2012).

Either way, we reinforce that we aim to demonstrate that Pop III stars are required to represent mean cosmic abundances, which can be straightforwardly observed in Figs. 5 and 6. The impact of different depletion-corrected methodologies, fitting of the data, and the use of different chemical yields should be addressed in detail in the future.

The knees that appear for model B in Figs. 5 and 6 are associated with the IMF exponents 0.85 and 1.00. The main reason for that comes from the nonexistence of stars enriching the medium with masses above 9 M for metallicity classes 4 × 10−3 and 8 × 10−3. We see that these x exponents are responsible for forming a higher number of high-mass stars when compared to the other IMFs. This introduces a flat behavior for the [Z/H] relation for these IMF exponents.

Although there are uncertainties about the mass spectrum of Pop III stars, our results show the importance of this stellar population for reproducing the observational data. There is no good agreement for the metallicity observed for the Universe only with Pop II stars. We see that Pop III stars rapidly raise the metallicity of the Universe between redshifts ∼15−20, mainly due to the HW02 branch. They reinforce chemical enrichment in the range ∼5−15 through the HW10 model and complement the interval 0−5 through the CL08 channel. The chemical avalanche produced by Pop III stars at high and moderate redshifts acts as a booster so that Pop II stars can add their contribution, through the different branches of Pop II, to the chemical enrichment of the Universe.

3.1. Properties of yields

In this section, we present the cosmic chemical evolution for 11 chemical elements for models A and B, compared with data from DLAs taken from the literature (Fig. 7), and briefly discuss the main subjects regarding each of the elements. Observational data for other elements are not dust-corrected due to the lack of sufficient data points with enough information for correction (such as [Zn/Fe] or [Si/Fe]).

thumbnail Fig. 7.

Chemical evolution for Fe, Si, Zn, Ni, P, Mg, Al, S, C, N, and O since the first stars started to die (z = 20) until z = 0. Model A starts with zero-metallicity stars, and as the Universe gets enriched, subsequent Pop II stars with increasing metallicity start to appear, until reaching Z = 2 × 10−2Z, according to the model described in Sect. 2.4. It is possible to observe the chemical avalanche in the early Universe given by the high production of metals from Pop III stars. As discussed in the text, model B considers only Pop II stars. The gray crosses represent data from Pettini et al. (1997, 2000, 2008), Centurion et al. (1998), Vladilo (1998), Dessauges-Zavadsky et al. (2002, 2003), Nissen et al. (2004), Akerman et al. (2005), Kulkarni et al. (2007), Prochaska et al. (2007), Penprase et al. (2010), Cooke et al. (2011), Kulkarni et al. (2012), Rafelski et al. (2012), Jorgenson et al. (2013), Neeleman et al. (2013), Zafar et al. (2014), and De Cia et al. (2016). All abundances were rescaled to solar values from Asplund et al. (2009).

Iron and silicon. Fe and Si could be altered by depletion (De Cia et al. 2013). Observational abundances could increase by ∼0.5 dex if depletion is considered in the comparison of the model with DLA data. Details about the methodology used to correct Fe depletion can be checked in the work by Vladilo (2002). As for Si, Prochaska & Wolfe (2002) show that although it is a refractory element, its depletion is not strong enough to significantly alter the abundances of DLA systems. Vladilo et al. (2011) show that Si depletion is mild in the ISM and it is expected to be weaker in most DLA systems. The depletion of Fe and Mg are measured for comparison, and it is found that the mean depletion of Si is almost as high as that of Fe, despite Fe being much more depleted than Si in the galactic ISM. They also explain that Si depletion in DLA systems does not correlate with metallicity, unlike Fe, whose depletion rises along with metallicity increase.

Zinc. Zn is produced mainly in HNe explosions characterized by a more significant production of Zn, Co, V, and Ti than normal SNe (Nomoto et al. 2006). Stars with 500−1000 M produce high amounts of Zn compared to O, C, and other metals (Ohkubo et al. 2006). Kobayashi et al. (2006) suggest that HNe can enhance the production of Zn, and that Zn is considered to be undepleted in DLAs.

Nickel and phosphorus. The lack of Ni observations in DLAs poses a challenge in analyzing this element; nevertheless, SNe Ia produce between 4 × 10−3 and 1.4 × 10−2M of Ni, and from 8.5 × 10−5 to 4 × 10−4M of P, depending on the specific model (Nomoto et al. 1997), and it is important to observe the outcome of these types of stars in the present model.

Magnesium. Mg is a refractory element, and its depleting effect must be considered. A challenge that arises in Mg determination comes from the saturation of the doublet used for its characterization, leaving only one possible line to provide Mg abundance. Given the problems related to its determination, current observations should be confirmed by additional Mg measurements to conclude if it could have a particular nucleosynthesis effect in DLA systems (Dessauges-Zavadsky et al. 2002).

Aluminum. Al has the strongest metal line transition observed in DLAs, the Al II λ 1670 line (Prochaska & Wolfe 2002). However, in the majority of systems, the line is heavily saturated and, together with the blending of lines and blending with the Ly-α forest, determining Al abundances can be a real challenge (Dessauges-Zavadsky et al. 2003).

Sulfur. S is considered non-refractory by some authors (Prochaska & Wolfe 2002), but there is still discussion about its actual behavior and if it could be used as a parameter for measuring depletion (Jenkins 2009; De Cia et al. 2016).

Carbon, nitrogen, and oxygen. There is an excess in the abundance of C, N, and O. Regarding depletion effects, O is only mildly refractory according to observations of DLAs and is not highly affected by depletion, although it is challenging to observe once the majority of the lines fall into the Ly-α forest and tend to be saturated (Prochaska & Wolfe 2002). On the other hand, C is considered mildly refractory (Prochaska & Wolfe 2002). Once it is a major constituent of interstellar dust (Henry et al. 2000; Jenkins 2009), a substantial part of C could be trapped on dust grains. Also, the lack of observations of C in DLAs is a problem (Jenkins 2009). N does not exhibit progressively stronger depletions (Jenkins 2009) and appears to be better represented by model B, that is, by the behavior of Pop II stars.

There are, however, other physical processes that participate in the C, N, and O production dynamics, which could be interfering with the results. Jenkins (2009) shows that, depending on the case, the consumption of O for producing oxides and silicates is not consistent with results for differential depletion for this element. The lack of O in gas-phase observations is much higher than what is needed for producing these silicates and oxides, and it is very hard to correlate the lack of O in the ISM with models of interstellar grain production. The author suggests that the formation of compounds involving elements such as H or C could play an important role in taking these elements from the ISM. Therefore, although cooling processes considerably demand C, N, and O for gas cooling and fragmentation, the grain formation processes do not entirely solve the problem for all three of these overabundant elements. An interesting result from Ioppolo et al. (2008) suggests that O is incorporated in the form of amorphous H2O ice on the grain surfaces. Work recently developed by Loeb et al. (2016) suggests that there is a possibility that carbon-enhanced metal-poor (CEMP) stars from the second generation of stars could host or have hosted planetary systems in their habitable zones. The planets would likely have a major C component in their composition. Also, the degree of C enhancement in CEMP stars has been shown to notably increase as a function of decreasing metallicity (Carollo et al. 2012), that is, the C enhancement in this type of star is likely much higher in the primordial Universe. Sonnentrucker et al. (2010) also show that the abundance of water vapor in gas clouds in the Galaxy holds ∼0.1% of the available O.

Lastly, the DLA data have scattering larger than that produced by our models. This result can be associated with several effects. For example, DLA data depend on the characteristics of the host galaxies. As our model is semi-analytic, it cannot resolve individual galaxies. In addition, the CSFR obtained in our models ends up being a weighted average over the halo masses through the formulation presented in Sect. 2.1. On the other hand, diffusion is a physical process that can transport metals away from their production sites, which could introduce scattering into our models. Aggregate diffusion is a task to be explored in future work.

3.2. The effect of feedback on our results

The feedback effects configure in a complex task, and they are beyond the scope of this work. However, through our modeling, we can make some inferences about the possible impacts on our results. In a recent article, Lancaster et al. (2021) presented simulations in which one of the main characteristics of feedback effects is the reduction of star formation efficiency. The authors’ simulations involved clouds with characteristic radii of 2.5 to 20 pc. The efficiency reduction achieved in these simulations ranged from 3% to 47% over the no-feedback efficiency.

In our case, the reduction of the star formation efficiency implies a decrease in the τs scale to keep the CSFR adjusted to the observational data. In this case, some chemical elements such as Zn, P, and S would deviate more from the DLA data, especially with the IMF exponent x ≥ 1.70. A possible way to solve this problem would be the inclusion of other channels for the production of chemical elements, such as the inclusion of hypernova yields and SNe Ia.

4. Summary and conclusions

The main goal of this work was to investigate cosmic chemical enrichment through the evolution of chemical elements in the redshift interval 0 ≤ z ≤ 20, as well as the contributions of Pop III and Pop II stars to the cosmic enrichment of the Universe. It was achieved by building a cosmic chemical evolution model that couples a semi-analytic cosmological model, which computes the CSFR, to chemical evolution models for the galactic framework. We computed the evolution of production of Fe, Si, Zn, Ni, P, Mg, Al, S, C, N, and O and compared our results with observational data taken from DLAs in the redshift interval [0−6] and with GCs.

Our main results show that we can consistently model the evolution of cosmic abundances in the Universe using a semi-analytic approach. Also, the “chemical avalanche” on the primordial Universe, which quickly enriches the medium and provides conditions for Pop II stars to appear, is consistent with the literature on Pop III stars’ behavior and chemical evolution models (Heger & Woosley 2002, 2010; Takahashi et al. 2018).

Regarding the behavior of Pop III and Pop II stars separately, the main difference appears in the behavior of abundances toward z = 0. At the same time, our model considering regular intermediate and high-mass Pop II stars (model B) shows a decrease in the abundances (except for N and Ztotal), while the model including very massive Pop III stars (model A) reproduces increasing abundances as redshift decreases, which is consistent with observations and similar models in the literature. Model A also offers a better fit of Ztotal to GC data than model B.

We conclude stating that model A, where the inclusion of Pop III stars appears as the main difference, presents a very good description of mean chemical values across the studied redshift range and the key behavior for the evolution of cosmic abundances in the Universe. Our main results are summarized below:

  • The chemical enrichment process in the early Universe occurs very quickly regardless of the stellar population. The pristine Universe reaches Z = 10−6Z in ∼3.0 Myr for the model with both Pop II and Pop III stars and with IMF 1.35, and ∼25 Myr for the model with only Pop II stars (with x = 1.35). However, when considering only high-mass Pop II stars, the metals are quickly consumed, and the scenario cannot represent chemical abundances at lower redshifts.

  • Abundances from GCs for Ztotal are consistently represented by the model with both Pop II and Pop III stars, while the model without Pop III stars is unable to represent observational data, regardless of the IMF.

  • Abundances from DLAs for Ztotal are consistently represented by our model with Pop III and Pop II stars. When comparing the model with abundances corrected for dust depletion and alpha enhancement, the observations show proper accordance with the model considering both Pop II and Pop III stars, while the model with only Pop II stars cannot account for the behavior of metals toward z = 0.

  • Regarding the modeling for other elements, there are a few deviations in the results when comparing the models with data from DLAs. However, the combination of mechanisms needed to improve the results is self-completing and can be easily understood, such as the absence of some mechanisms SNe Ia, HNe, dust depletion affecting observational data, and the combination of yields from Pop II stars. HNe and maybe a higher-mass branch of stars ∼500−1000 M, (Ohkubo et al. 2006) should improve results for Zn, P, and Ni without raising O and C. In principle, these mechanisms are all consistent with each other and will be studied in a subsequent work.

  • The reason for the overabundances of C, N, and O shown in our results remains an open question. New observations focusing on depletion processes in the ISM could explain the overabundances found in the present work (and/or the lack of these elements in the ISM).

Altogether, our results indicate that the evolution of chemical abundances in the cosmological framework can be consistently tracked. Our most important result shows that Pop III stars’ contribution to the Universe’s chemical history should be better understood, and observational campaigns with instruments capable of actually identifying these objects should be seriously considered and implemented. Pop III observations are a long-awaited result, and a firm detection will shed new light on the cosmic history in earlier times. The other questions raised in this paper are being studied and will be the subject of forthcoming works.

Acknowledgments

We thank the Brazilian agency FAPESP for support under the thematic project 2014/11156-4. L.C.C. would like to thank the Coordenação de Aperfeiçoamento de Pessoal de Nível Superior (CAPES) – Finance code 001 – for a graduate research fellowship. O.D.M. and C.A.W. thank CNPq for partial financial support (grants 303350/2015-6 and 313597/2014-6, respectively).

References

  1. Abel, T., Bryan, G. L., & Norman, M. L. 2000, ApJ, 540, 39 [NASA ADS] [CrossRef] [Google Scholar]
  2. Abia, C., Dominguez, I., Straniero, O., et al. 2001, ApJ, 557, 126 [NASA ADS] [CrossRef] [Google Scholar]
  3. Akerman, C. J., Ellison, S. L., Pettini, M., & Steidel, C. C. 2005, A&A, 440, 499 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  4. Asplund, M., Grevesse, N., Sauval, A. J., & Scott, P. 2009, ARA&A, 47, 481 [NASA ADS] [CrossRef] [Google Scholar]
  5. Barkana, R. 2006, Science, 313, 931 [NASA ADS] [CrossRef] [Google Scholar]
  6. Bloecker, T. 1995, A&A, 299, 755 [NASA ADS] [Google Scholar]
  7. Bond, H. E., Nelan, E. P., VandenBerg, D. A., Schaefer, G. H., & Harmer, D. 2013, ApJ, 765, L12 [NASA ADS] [CrossRef] [Google Scholar]
  8. Bouwens, R. J., Illingworth, G. D., Oesch, P. A., et al. 2012a, ApJ, 752, L5 [Google Scholar]
  9. Bouwens, R. J., Illingworth, G. D., Oesch, P. A., et al. 2012b, ApJ, 754, 83 [Google Scholar]
  10. Bromm, V., & Loeb, A. 2003, Nature, 425, 812 [NASA ADS] [CrossRef] [Google Scholar]
  11. Bromm, V., Ferrara, A., Coppi, P. S., & Larson, R. B. 2001, MNRAS, 328, 969 [NASA ADS] [CrossRef] [Google Scholar]
  12. Calura, F., & Matteucci, F. 2004, MNRAS, 350, 351 [NASA ADS] [CrossRef] [Google Scholar]
  13. Calura, F., & Matteucci, F. 2006, in Origin of Matter and Evolution of Galaxies, eds. S. Kubono, W. Aoki, T. Kajino, T. Motobayashi, & K. Nomoto, AIP Conf. Ser., 847, 371 [NASA ADS] [Google Scholar]
  14. Calura, F., & Menci, N. 2009, MNRAS, 400, 1347 [NASA ADS] [CrossRef] [Google Scholar]
  15. Campbell, S. W., & Lattanzio, J. C. 2008, A&A, 490, 769 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  16. Carollo, D., Beers, T. C., Bovy, J., et al. 2012, ApJ, 744, 195 [NASA ADS] [CrossRef] [Google Scholar]
  17. Carroll, S. M., Press, W. H., & Turner, E. L. 1992, ARA&A, 30, 499 [NASA ADS] [CrossRef] [Google Scholar]
  18. Casey, C. M., Berta, S., Béthermin, M., et al. 2012, ApJ, 761, 139 [Google Scholar]
  19. Cayrel, R., Hill, V., Beers, T. C., et al. 2001, Nature, 409, 691 [CrossRef] [PubMed] [Google Scholar]
  20. Centurion, M., Bonifacio, P., Molaro, P., & Vladilo, G. 1998, ApJ, 509, 620 [NASA ADS] [CrossRef] [Google Scholar]
  21. Chieffi, A., & Limongi, M. 2004, ApJ, 608, 405 [NASA ADS] [CrossRef] [Google Scholar]
  22. Ciardi, B., Ferrara, A., Marri, S., & Raimondo, G. 2001, MNRAS, 324, 381 [NASA ADS] [CrossRef] [Google Scholar]
  23. Cojazzi, P., Bressan, A., Lucchin, F., Pantano, O., & Chavez, M. 2000, MNRAS, 315, L51 [NASA ADS] [CrossRef] [Google Scholar]
  24. Cooke, R., Pettini, M., Steidel, C. C., Rudie, G. C., & Nissen, P. E. 2011, MNRAS, 417, 1534 [CrossRef] [Google Scholar]
  25. Côté, B., Ritter, C., O’Shea, B. W., et al. 2016, ApJ, 824, 82 [CrossRef] [Google Scholar]
  26. Cowan, J. J., Sneden, C., Burles, S., et al. 2002, ApJ, 572, 861 [NASA ADS] [CrossRef] [Google Scholar]
  27. Cucciati, O., Tresse, L., Ilbert, O., et al. 2012, A&A, 539, A31 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  28. Dahlen, T., Mobasher, B., Dickinson, M., et al. 2007, ApJ, 654, 172 [Google Scholar]
  29. Daigne, F., Olive, K. A., Silk, J., Stoehr, F., & Vangioni, E. 2006, ApJ, 647, L773 [NASA ADS] [CrossRef] [Google Scholar]
  30. Davé, R., & Oppenheimer, B. D. 2007, MNRAS, 374, 427 [CrossRef] [Google Scholar]
  31. De Cia, A., Ledoux, C., Savaglio, S., Schady, P., & Vreeswijk, P. M. 2013, A&A, 560, A88 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  32. De Cia, A., Ledoux, C., Mattsson, L., et al. 2016, A&A, 596, A97 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  33. De Cia, A., Ledoux, C., Petitjean, P., & Savaglio, S. 2018, A&A, 611, A76 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  34. Dessauges-Zavadsky, M., Prochaska, J. X., & D’Odorico, S. 2002, A&A, 391, 801 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  35. Dessauges-Zavadsky, M., Peroux, C., Kim, T.-S., D’Odorico, S., & McMahon, R. G. 2003, MNRAS, 345, 447 [NASA ADS] [CrossRef] [Google Scholar]
  36. Dessauges-Zavadsky, M., Zamojski, M., Schaerer, D., et al. 2015, A&A, 577, A50 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  37. Doherty, C. L., Gil-Pons, P., Lau, H. H. B., Lattanzio, J. C., & Siess, L. 2013, MNRAS, 437, 195 [CrossRef] [Google Scholar]
  38. Doherty, C. L., Gil-Pons, P., Lau, H. H. B., et al. 2014, MNRAS, 441, 582 [NASA ADS] [CrossRef] [Google Scholar]
  39. Dotter, A., Sarajedini, A., & Anderson, J. 2011, ApJ, 738, 74 [NASA ADS] [CrossRef] [Google Scholar]
  40. Dvorkin, I., Silk, J., Vangioni, E., Petitjean, P., & Olive, K. A. 2015, MNRAS, 452, L36 [NASA ADS] [CrossRef] [Google Scholar]
  41. Ekström, S., Meynet, G., Maeder, A., & Barblan, F. 2008, A&A, 478, 467 [Google Scholar]
  42. Fang, T., & Cen, R. 2004, ApJ, 616, L87 [NASA ADS] [CrossRef] [Google Scholar]
  43. Fraser, M., Casey, A. R., Gilmore, G., Heger, A., & Chan, C. 2017, MNRAS, 468, 418 [NASA ADS] [CrossRef] [Google Scholar]
  44. Frebel, A., & Norris, J. E. 2015, ARA&A, 53, 631 [NASA ADS] [CrossRef] [Google Scholar]
  45. Frebel, A., Christlieb, N., Norris, J. E., et al. 2007, ApJ, 660, L117 [NASA ADS] [CrossRef] [Google Scholar]
  46. Frost, C. A., & Lattanzio, J. C. 1996, ApJ, 473, 383 [NASA ADS] [CrossRef] [Google Scholar]
  47. Fynbo, J. P. U., Starling, R. L. C., Ledoux, C., et al. 2006, A&A, 451, L47 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  48. Galli, D., & Palla, F. 2013, ARA&A, 51, 163 [NASA ADS] [CrossRef] [Google Scholar]
  49. Gnedin, N. Y., & Ostriker, J. P. 1997, ApJ, 486, 581 [Google Scholar]
  50. Gribel, C., Miranda, O. D., & Vilas-Boas, J. W. 2017, ApJ, 849, 108 [NASA ADS] [CrossRef] [Google Scholar]
  51. Gruppioni, C., Pozzi, F., Rodighiero, G., et al. 2013, MNRAS, 432, 23 [Google Scholar]
  52. Heger, A., & Woosley, S. E. 2002, ApJ, 567, 532 [Google Scholar]
  53. Heger, A., & Woosley, S. E. 2010, ApJ, 724, 341 [Google Scholar]
  54. Henry, R. B. C., Edmunds, M. G., & Koppen, J. 2000, ApJ, 541, 660 [NASA ADS] [CrossRef] [Google Scholar]
  55. Hirano, S., & Yoshida, N. 2013, ApJ, 763, 52 [NASA ADS] [CrossRef] [Google Scholar]
  56. Hodge, J. A., & da Cunha, E. 2020, R. Soc. Open Sci., 7, 200556 [NASA ADS] [CrossRef] [Google Scholar]
  57. Iglesias, C. A., & Rogers, F. J. 1996, ApJ, 464, 943 [NASA ADS] [CrossRef] [Google Scholar]
  58. Ioppolo, S., Cuppen, H. M., Romanzin, C., van Dishoeck, E. F., & Linnartz, H. 2008, ApJ, 686, 1474 [NASA ADS] [CrossRef] [Google Scholar]
  59. Jenkins, E. B. 2009, ApJ, 700, 1299 [Google Scholar]
  60. Jorgenson, R. A., Murphy, M. T., & Thompson, R. 2013, MNRAS, 435, 482 [NASA ADS] [CrossRef] [Google Scholar]
  61. Karakas, A. I. 2010, MNRAS, 403, 1413 [NASA ADS] [CrossRef] [Google Scholar]
  62. Kashikawa, N., Nagao, T., Toshikawa, J., et al. 2012, ApJ, 761, 85 [NASA ADS] [CrossRef] [Google Scholar]
  63. Kistler, M. D., Yüksel, H., Beacom, J. F., Hopkins, A. M., & Wyithe, J. S. B. 2009, ApJ, 705, L104 [Google Scholar]
  64. Kobayashi, C. 2005, in From Lithium to Uranium: Elemental Tracers of Early Cosmic Evolution, eds. V. Hill, P. Francois, & F. Primas, 228, 315 [NASA ADS] [Google Scholar]
  65. Kobayashi, C., Umeda, H., Nomoto, K., Tominaga, N., & Ohkubo, T. 2006, ApJ, 653, 1145 [NASA ADS] [CrossRef] [Google Scholar]
  66. Kobayashi, C., Springel, V., & White, S. D. M. 2007, MNRAS, 376, 1465 [Google Scholar]
  67. Krumholz, M. R., & McKee, C. F. 2005, ApJ, 630, 250 [Google Scholar]
  68. Kulkarni, V. P., Khare, P., Péroux, C., et al. 2007, ApJ, 661, 88 [NASA ADS] [CrossRef] [Google Scholar]
  69. Kulkarni, V. P., Meiring, J., Som, D., et al. 2012, ApJ, 749, 176 [NASA ADS] [CrossRef] [Google Scholar]
  70. Lancaster, L., Ostriker, E. C., Kim, J.-G., & Kim, C.-G. 2021, ApJ, 922, L3 [NASA ADS] [CrossRef] [Google Scholar]
  71. Larson, R. B., Tinsley, B. M., & Caldwell, C. N. 1980, ApJ, 237, 692 [Google Scholar]
  72. Lilly, S. J., Carollo, C. M., Pipino, A., Renzini, A., & Peng, Y. 2013, ApJ, 772, 119 [NASA ADS] [CrossRef] [Google Scholar]
  73. Liu, D., Schinnerer, E., Groves, B., et al. 2019, ApJ, 887, 235 [Google Scholar]
  74. Loeb, A., & Barkana, R. 2001, ARA&A, 39, 19 [NASA ADS] [CrossRef] [Google Scholar]
  75. Loeb, A., Batista, R. A., & Sloan, D. 2016, JCAP, 2016, 040 [CrossRef] [Google Scholar]
  76. Ma, X., Hopkins, P. F., Faucher-Giguère, C.-A., et al. 2016, MNRAS, 456, 2140 [NASA ADS] [CrossRef] [Google Scholar]
  77. Ma, Q., Maio, U., Ciardi, B., & Salvaterra, R. 2017, MNRAS, 466, 1140 [NASA ADS] [CrossRef] [Google Scholar]
  78. Madau, P., & Dickinson, M. 2014, ARA&A, 52, 415 [Google Scholar]
  79. Magnelli, B., Elbaz, D., Chary, R. R., et al. 2011, A&A, 528, A35 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  80. Magnelli, B., Popesso, P., Berta, S., et al. 2013, A&A, 553, A132 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  81. Maio, U., & Tescari, E. 2015, MNRAS, 453, 3799 [NASA ADS] [CrossRef] [Google Scholar]
  82. Maio, U., Péroux, C., & Ciardi, B. 2022, A&A, 657, A47 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  83. Maio, U., Ciardi, B., Dolag, K., Tornatore, L., & Khochfar, S. 2010, MNRAS, 407, 1003 [CrossRef] [Google Scholar]
  84. Matteucci, F. 2001, The Chemical Evolution of the Galaxy (Springer), 253 [Google Scholar]
  85. Matteucci, F. 2016, J. Phys. Conf. Ser., 703, 012004 [NASA ADS] [CrossRef] [Google Scholar]
  86. Matteucci, F., & Calura, F. 2005, MNRAS, 360, 447 [NASA ADS] [CrossRef] [Google Scholar]
  87. Nakamura, F., & Umemura, M. 2001, ApJ, 548, 19 [NASA ADS] [CrossRef] [Google Scholar]
  88. Neeleman, M., Wolfe, A. M., Prochaska, J. X., & Rafelski, M. 2013, ApJ, 769, 54 [NASA ADS] [CrossRef] [Google Scholar]
  89. Nissen, P. E., Chen, Y. Q., Asplund, M., & Pettini, M. 2004, A&A, 415, 993 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  90. Nomoto, K., Iwamoto, K., Nakasato, N., et al. 1997, Nucl. Phys. A, 621, 467 [NASA ADS] [CrossRef] [Google Scholar]
  91. Nomoto, K., Tominaga, N., Umeda, H., Kobayashi, C., & Maeda, K. 2006, Nucl. Phys. A, 777, 424 [CrossRef] [Google Scholar]
  92. Ohkubo, T., Umeda, H., Maeda, K., et al. 2006, ApJ, 645, 1352 [NASA ADS] [CrossRef] [Google Scholar]
  93. Penprase, B. E., Prochaska, J. X., Sargent, W. L. W., Toro-Martinez, I., & Beeler, D. J. 2010, ApJ, 721, 1 [NASA ADS] [CrossRef] [Google Scholar]
  94. Pereira, E. S., & Miranda, O. D. 2010, MNRAS, 401, 1924 [NASA ADS] [CrossRef] [Google Scholar]
  95. Péroux, C., & Howk, J. C. 2020, ARA&A, 58, 363 [CrossRef] [Google Scholar]
  96. Pettini, M., King, D. L., Smith, L. J., & Hunstead, R. W. 1997, ApJ, 478, 536 [NASA ADS] [CrossRef] [Google Scholar]
  97. Pettini, M., Ellison, S. L., Steidel, C. C., Shapley, A. E., & Bowen, D. V. 2000, ApJ, 532, 65 [NASA ADS] [CrossRef] [Google Scholar]
  98. Pettini, M., Zych, B. J., Steidel, C. C., & Chaffee, F. H. 2008, MNRAS, 385, 2011 [NASA ADS] [CrossRef] [Google Scholar]
  99. Press, W. H., & Schechter, P. 1974, ApJ, 187, 425 [Google Scholar]
  100. Prochaska, J. X., & Wolfe, A. M. 2002, ApJ, 566, 68 [NASA ADS] [CrossRef] [Google Scholar]
  101. Prochaska, J. X., Wolfe, A. M., Howk, J. C., et al. 2007, ApJS, 171, 29 [NASA ADS] [CrossRef] [Google Scholar]
  102. Rafelski, M., Wolfe, A. M., Prochaska, J. X., Neeleman, M., & Mendez, A. J. 2012, ApJ, 755, 89 [NASA ADS] [CrossRef] [Google Scholar]
  103. Raiteri, C. M., Villata, M., & Navarro, J. F. 1996, A&A, 315, 105 [NASA ADS] [Google Scholar]
  104. Reddy, N. A., & Steidel, C. C. 2009, ApJ, 692, 778 [Google Scholar]
  105. Reimers, D. 1975, Mem. Soc. R. Sci. Liège, 8, 369 [Google Scholar]
  106. Robotham, A. S. G., & Driver, S. P. 2011, MNRAS, 413, 2570 [Google Scholar]
  107. Rollinde, E., Vangioni, E., Maurin, D., et al. 2009, MNRAS, 398, 1782 [NASA ADS] [CrossRef] [Google Scholar]
  108. Salpeter, E. E. 1959, ApJ, 129, 608 [NASA ADS] [CrossRef] [Google Scholar]
  109. Salvadori, S., Schneider, R., & Ferrara, A. 2007, MNRAS, 381, 647 [NASA ADS] [CrossRef] [Google Scholar]
  110. Santoro, F., & Shull, J. M. 2006, ApJ, 643, 26 [NASA ADS] [CrossRef] [Google Scholar]
  111. Schaerer, D. 2002, A&A, 382, 28 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  112. Schiminovich, D., Ilbert, O., Arnouts, S., et al. 2005, ApJ, 619, L47 [Google Scholar]
  113. Schinnerer, E., Groves, B., Sargent, M. T., et al. 2016, ApJ, 833, 112 [CrossRef] [Google Scholar]
  114. Schmidt, M. 1959, ApJ, 129, 243 [NASA ADS] [CrossRef] [Google Scholar]
  115. Schneider, R. 2010, in First Stars and Galaxies: Challenges for the Next Decade, eds. D. J. Whalen, V. Bromm, & N. Yoshida, AIP Conf. Ser., 1294, 102 [NASA ADS] [Google Scholar]
  116. Schneider, R., Salvaterra, R., Ferrara, A., & Ciardi, B. 2006, MNRAS, 369, 825 [NASA ADS] [CrossRef] [Google Scholar]
  117. Scoville, N., Aussel, H., Sheth, K., et al. 2014, ApJ, 783, 84 [CrossRef] [Google Scholar]
  118. Scoville, N., Sheth, K., Aussel, H., et al. 2016, ApJ, 820, 83 [Google Scholar]
  119. Scoville, N., Lee, N., Vanden Bout, P., et al. 2017, ApJ, 837, 150 [Google Scholar]
  120. Sheth, R. K., & Tormen, G. 1999, MNRAS, 308, 119 [Google Scholar]
  121. Shu, X. W., Elbaz, D., Bourne, N., et al. 2016, ApJS, 222, 4 [NASA ADS] [CrossRef] [Google Scholar]
  122. Sobral, D., Matthee, J., Darvish, B., et al. 2015, ApJ, 808, 139 [Google Scholar]
  123. Sonnentrucker, P., Neufeld, D. A., Phillips, T. G., et al. 2010, A&A, 521, L12 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  124. Spera, M., Mapelli, M., & Bressan, A. 2015, MNRAS, 451, 4086 [NASA ADS] [CrossRef] [Google Scholar]
  125. Springel, V., & Hernquist, L. 2003, MNRAS, 339, 312 [NASA ADS] [CrossRef] [Google Scholar]
  126. Tacconi, L. J., Genzel, R., Saintonge, A., et al. 2018, ApJ, 853, 179 [Google Scholar]
  127. Takahashi, K., Yoshida, T., & Umeda, H. 2018, ApJ, 857, 111 [NASA ADS] [CrossRef] [Google Scholar]
  128. Tan, W.-W., Wang, F. Y., & Cheng, K. S. 2016, ApJ, 829, 29 [NASA ADS] [CrossRef] [Google Scholar]
  129. Tinsley, B. M., & Larson, R. B. 1978, ApJ, 221, 554 [Google Scholar]
  130. Tornatore, L., Ferrara, A., & Schneider, R. 2007, MNRAS, 382, 945 [NASA ADS] [CrossRef] [Google Scholar]
  131. Torrey, P., Vogelsberger, M., Marinacci, F., et al. 2019, MNRAS, 484, 5587 [NASA ADS] [Google Scholar]
  132. Tumlinson, J., & Shull, J. M. 2000, ApJ, 528, L65 [NASA ADS] [CrossRef] [Google Scholar]
  133. Tumlinson, J., Venkatesan, A., & Shull, J. M. 2004, ApJ, 612, 602 [NASA ADS] [CrossRef] [Google Scholar]
  134. Vangioni, E., Dvorkin, I., Olive, K. A., et al. 2018, MNRAS, 477, 56 [NASA ADS] [CrossRef] [Google Scholar]
  135. Vanzella, E., Meneghetti, M., Caminha, G. B., et al. 2020, MNRAS, 494, L81 [Google Scholar]
  136. Vassiliadis, E., & Wood, P. R. 1993, ApJ, 413, 641 [Google Scholar]
  137. Venkatesan, A., Tumlinson, J., & Shull, J. M. 2003, ApJ, 584, 621 [NASA ADS] [CrossRef] [Google Scholar]
  138. Vladilo, G. 1998, ApJ, 493, 583 [NASA ADS] [CrossRef] [Google Scholar]
  139. Vladilo, G. 2002, A&A, 391, 407 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  140. Vladilo, G., Abate, C., Yin, J., Cescutti, G., & Matteucci, F. 2011, A&A, 530, A33 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  141. Wagner-Kaiser, R., Mackey, D., Sarajedini, A., et al. 2017, MNRAS, 471, 3347 [Google Scholar]
  142. Wolfe, A. M., Gawiser, E., & Prochaska, J. X. 2005, ARA&A, 43, 861 [NASA ADS] [CrossRef] [Google Scholar]
  143. Wyder, T. K., Treyer, M. A., Milliard, B., et al. 2005, ApJ, 619, L15 [Google Scholar]
  144. Zafar, T., Centurión, M., Péroux, C., et al. 2014, MNRAS, 444, 744 [NASA ADS] [CrossRef] [Google Scholar]

All Tables

Table 1.

Cosmological and structure formation parameters used to obtain the CSFR.

Table 2.

Masses selected for Pop III chemical yields.

Table 3.

Masses and metallicities selected for Pop II chemical yields.

All Figures

thumbnail Fig. 1.

Solution for ρ ˙ ( z ) $ \dot\rho_{\star}(z) $ as derived in this work (red line) and, using black dashed line, the CSFR used by Vangioni et al. (2018) plotted for comparison. Left: evolution of the CSFR from the local Universe to z = 20. Right: same as for the left panel, but zooming into 0 ≤ z ≤ 10, allowing for a better visualization of the two CSFRs within the range of the available observational data. Data used in this figure: IR (Magnelli et al. 2011, 2013 – dark gray filled circles; Gruppioni et al. 2013 – dark orange open circles); UV (Wyder et al. 2005 – slate gray open triangle; Schiminovich et al. 2005 – blue crosses; Dahlen et al. 2007 – turquoise open squares; Reddy & Steidel 2009 – dark green filled squares; Robotham & Driver 2011 – chocolate cross; Cucciati et al. 2012 – green stars; Bouwens et al. 2012a, 2012b – teal filled triangles); and GRB (Kistler et al. 2009 – deep pink open triangles).

In the text
thumbnail Fig. 2.

Evolution of the characteristic timescale for star formation as a function of the redshift. Left: τs for different IMF exponent and the evolution of the dynamical time (tdyn), taken to be 10% of Hubble time. Right: our results compared with the timescales for the conversion of gas (filled areas) obtained by Péroux & Howk (2020) and Maio et al. (2022) for molecular and cold (HI + H2) gas.

In the text
thumbnail Fig. 3.

Evolution of the fraction of gas converted into stars for different values of the IMF used in this work. Left: evolution of the average star formation efficiency (defined as ⟨ε⟩) as a function of z. Right: evolution of the fraction of total cold gas associated with star formation (fgas). The symbols show the measurements compiled as follows: circles by Tacconi et al. (2018), squares by Scoville et al. (2014), diamonds by Scoville et al. (2016), stars by Dessauges-Zavadsky et al. (2015), and triangles by Schinnerer et al. (2016). See also Fig. 13b of Hodge & da Cunha (2020).

In the text
thumbnail Fig. 4.

Evolution of metallicity for different values of x. Left: model A. Right: model B. Colors represent different values of x for the IMF.

In the text
thumbnail Fig. 6.

Comparison of models A (left) and B (right) with data from dust-corrected DLAs (gray crosses) from De Cia et al. (2018). We include an α-enhancement correction, with [Z/H]=[Fe/H]+0.3 dex, as suggested by Rafelski et al. (2012).

In the text
thumbnail Fig. 7.

Chemical evolution for Fe, Si, Zn, Ni, P, Mg, Al, S, C, N, and O since the first stars started to die (z = 20) until z = 0. Model A starts with zero-metallicity stars, and as the Universe gets enriched, subsequent Pop II stars with increasing metallicity start to appear, until reaching Z = 2 × 10−2Z, according to the model described in Sect. 2.4. It is possible to observe the chemical avalanche in the early Universe given by the high production of metals from Pop III stars. As discussed in the text, model B considers only Pop II stars. The gray crosses represent data from Pettini et al. (1997, 2000, 2008), Centurion et al. (1998), Vladilo (1998), Dessauges-Zavadsky et al. (2002, 2003), Nissen et al. (2004), Akerman et al. (2005), Kulkarni et al. (2007), Prochaska et al. (2007), Penprase et al. (2010), Cooke et al. (2011), Kulkarni et al. (2012), Rafelski et al. (2012), Jorgenson et al. (2013), Neeleman et al. (2013), Zafar et al. (2014), and De Cia et al. (2016). All abundances were rescaled to solar values from Asplund et al. (2009).

In the text

Current usage metrics show cumulative count of Article Views (full-text article views including HTML views, PDF and ePub downloads, according to the available data) and Abstracts Views on Vision4Press platform.

Data correspond to usage on the plateform after 2015. The current usage metrics is available 48-96 hours after online publication and is updated daily on week days.

Initial download of the metrics may take a while.