# Observational and theoretical constraints on the formation and early evolution of the first dust grains in galaxies at $5 < z < 10$

D. Burgarella, A. Nanni, H. Hirashita, P. Theulé, A. K. Inoue, T. T. Takeuchi

► To cite this version:

D. Burgarella, A. Nanni, H. Hirashita, P. Theulé, A. K. Inoue, et al.. Observational and theoretical constraints on the formation and early evolution of the first dust grains in galaxies at $5 < z < 10$. Astronomy and Astrophysics - A&A, 2020, 637, pp.A32. 10.1051/0004-6361/201937143. hal-02903859

HAL Id: hal-02903859

https://hal.science/hal-02903859

Submitted on 21 Jul 2020

HAL is a multi-disciplinary open access archive for the deposit and dissemination of scientific research documents, whether they are published or not. The documents may come from teaching and research institutions in France or abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est destinée au dépôt et à la diffusion de documents scientifiques de niveau recherche, publiés ou non, émanant des établissements d'enseignement et de recherche français ou étrangers, des laboratoires publics ou privés.

<figure><img src="image_1.png" /></figure>

Distributed under a Creative Commons Attribution 4.0 International License
# Observational and theoretical constraints on the formation and early evolution of the first dust grains in galaxies at $5 < z < 10$

D. Burgarella<fnref n="1" />, A. Nanni<fnref n="1" />, H. Hirashita<fnref n="2" />, P. Theulé<fnref n="1" />, A. K. Inoue<fnref n="3" />,<fnref n="4" />, and T. T. Takeuchi<fnref n="5" />

Received 19 November 2019 / Accepted 3 February 2020

## ABSTRACT

*Context.* The first generation of stars were born a few hundred million years after the big bang. These stars synthesise elements heavier than H and He, which are later expelled into the interstellar medium, initiating the rise of metals. Within this enriched medium, the first dust grains were formed. This event is cosmologically crucial for molecule formation, as dust plays a major role by cooling low-metallicity star-forming clouds, which can fragment to create lower mass stars. Collecting information on these first dust grains is difficult because of the negative alliance of large distances and low dust masses.

*Aims.* We aim to combine the observational information from galaxies at redshifts $5 \le z \le 10$ to constrain their dust emission and theoretically understand the first evolutionary phases of the dust cycle.

*Methods.* Spectral energy distributions (SEDs) are fitted with CIGALE and the physical parameters and their evolution are modelled. From this SED fitting, we built a dust-emission template for this population of galaxies in the reionisation epoch.

*Results.* Our new models explain why some early galaxies are observed and others are not. We follow in time the formation of the first grains by supernovae later destroyed by other supernova blasts and expelled in the circumgalactic and intergalactic media.

*Conclusions.* We find evidence for the first dust grains formed in the universe. But above all, this work underlines the need to collect more data and to develop new facilities to further constrain the dust cycle in galaxies in the reionisation epoch.

**Key words.** galaxies: formation – Galaxy: evolution – galaxies: high-redshift – dark ages, reionization, first stars – dust, extinction – infrared: galaxies

## 1. Introduction

Understanding the characteristics of the dust cycle is still an issue. This statement is even truer for Lyman break galaxies at all redshifts (LBGs, e.g. Burgarella et al. 2007) and at redshifts $z > 5$ (Bouwens et al. 2011; Finkelstein et al. 2010) at the epoch of reionization (EoR). This is due to the severe limits due to sub-millimetre (submm) and millimetre (mm) observations but also because we have a poor idea on the chemical conditions, the dust grain characteristics and even the stellar populations.

Fundamentally, as in the local universe, dust grains in high-redshift (Hi-z) galaxies absorb far ultraviolet (FUV) photons from young and massive stars. By doing so, these grains warm up to dust temperatures of a few tens of degrees and emit far-infrared (FIR) photons.

This paper re-evaluates the roles of the relevant physical processes in the formation and early evolution of dust grains in these LBGs: formation in supernovae (SNe) explosion and outflows from lower mass stars (e.g. Todini & Ferrara 2001; Matsuura et al. 2015; Ventura et al. 2012; Nanni et al. 2013, 2014; Marassi et al. 2019), destruction by SNe shocks, growth by accretion (e.g. Hirashita & Omukai 2009; Dwek & Cherchneff 2011; Asano et al. 2013), and galactic outflows (e.g. Jones et al.

2018; Ohyama et al. 2019). Within the context of this paper, two main questions can be asked:

*What is the dust cycle in high-redshift objects?*. Two types of stellar sources form and eject dust grains in the interstellar medium (ISM) of galaxies: stars with an initial mass in the $1–8 M_\odot$ range during the asymptotic giant branch (AGB) phase, and stars with an initial mass in the $8–40 M_\odot$ range during the core-collapse SNe phase. Dust grains also form seeds that grow in the ISM (e.g. Dwek 1998; Draine 2009; Jones & Nuth 2011; Asano et al. 2013). The AGB dominates dust production in galaxies old enough to allow low-mass stars to evolve to the AGB phase (Valiante et al. 2009). The AGB’s contribution to the global dust production in metal-poor environments reaches at most 30% after about 2 Gyr (Dell’Agli et al. 2019). This is confirmed in other works (e.g. Dwek & Cherchneff 2011; Dwek et al. 2011). However, SNe blasts produce a shockwave able to destroy dust grains in the ISM (e.g. Slavin et al. 2015; Matsuura et al. 2019; Nozawa et al. 2003; Dwek & Cherchneff 2011).

There are numerous works (e.g. review by Rupke 2018) showing that galactic winds are ubiquitous at low redshifts. These multi-phase winds carry neutral and ionised gas and dust grains outside the galaxies. For instance, in $M_{82}$, dust is found well beyond the radius where gas is thought to be

<footnote n="1">¹ Aix Marseille Univ, CNRS, CNES, LAM, Marseille, France
e-mail: denis.burgarella@lam.fr</footnote>

<footnote n="2">² Institute of Astronomy and Astrophysics, Academia Sinica, Astronomy-Mathematics Building, AS/NTU, No. 1, Sec. 4, Roosevelt Road, Taipei 10617, Taiwan</footnote>

<footnote n="3">³ Department of Physics, School of Advanced Science and Engineering, Faculty of Science and Engineering, Waseda University, 3-4-1 Okubo, Shinjuku, Tokyo 169-8555, Japan</footnote>

<footnote n="4">⁴ Waseda Research Institute for Science and Engineering, Faculty of Science and Engineering, Waseda University, 3-4-1 Okubo, Shinjuku Tokyo 169-8555, Japan</footnote>

<footnote n="5">⁵ Division of Particle and Astrophysical Science, Nagoya University, Furo-cho, Chikusa-ku, Nagoya 464-8602, Japan</footnote>
outflowing from the galaxy. This suggests that dust escapes from the gravitational potential of galaxies into the intergalactic space (e.g. Engelbracht et al. 2006; Yoshida et al. 2019). Dust is also observed in the outflows from other local star-forming galaxies and strong evidence for dust in the outflowing neutral clouds is also found at $z = 3$ (Shapley et al. 2003). Very recently, in a very detailed analysis, Jones et al. (2018) found that their measurements favour a picture in which the majority of heavy elements are ejected in a predominantly low-ionisation outflow, which regulates galactic chemical evolution.

The origin and the evolution of dust at high redshift has become a prominent issue since the early works of Watson et al. (2015). At $z>5$, large amounts of dust (~$10^8 M_\odot$) are found in galaxies (e.g. Riechers et al. 2013; Mortlock et al. 2011; Venemans et al. 2012; Michałowski 2015) or quasars (e.g. Bertoldi et al. 2003; Priddey et al. 2003; Robson et al. 2004; Beelen et al. 2006; Michałowski et al. 2010; Cullen et al. 2017; Leśniewska & Michałowski 2019). Recently, we started to collect new and interesting data on the galaxy host of γ-ray bursts that will be very useful for the study of the properties of dust and the comparison of dust extinction and dust attenuation laws in these objects (e.g. Stratta et al. 2011; Hjorth et al. 2013; Bolmer et al. 2018).

**How can we explain the rest-frame UV colour properties and detectability of Hi-z LBGs?** Most Hi-z LBGs observed in FIR are not detected (e.g. Capak et al. 2015; Bouwens et al. 2016; Hirashita et al. 2017). Explanations are proposed to understand the non-detections. Ferrara et al. (2017) suggested that these objects can contain a large molecular gas mass. While dust in the diffuse ISM attains relatively high temperatures ($T \approx 70–80 K$), dust embedded in dense gas remains cold ($T \approx 30–40 K$). Faisst et al. (2017) suggested that radiation pressure causes a spatial offset between dust clouds and young star-forming regions within the lifetime of O/B stars. These offsets modify the radiation balance and create viewing-angle effects that can change UV colours at fixed IRX (IRX = log ($L_{dust}/L_{FUV}$)).

We use a seven-year WMAP7 cosmology (Komatsu et al. 2011) as adopted by Astropy, of $\Omega_M = 0.273$, $\Omega_\Lambda = 0.727$ and $H_0 = 70.4$. We assume a Chabrier initial mass function (IMF, Chabrier 2003) for the SED fitting because top-heavy IMF ($m^{-\alpha}$ with $\alpha = 1.0, 1.35, 1.5$, Dabringhausen et al. 2010; Gall et al. 2011) single stellar population (SSPs) models are not available in CIGALE. The Chabrier IMF for the SED fitting is fixed in the one assumed by Bruzual & Charlot (2003), that is, the lower and upper mass cut-offs $m_{Low} = 0.1 M_\odot$ and $m_{Up} = 100 M_\odot$. Chabrier and top-heavy IMFs are assumed for the dust modelling lower and upper mass cut-offs $m_{Low} = 1.3 M_\odot$ and $m_{Up} = 40 M_\odot$ (with Kobayashi yields). In agreement with Sect. 5, with the exception of the pop III stars used between 10 and 30 $M_\odot$.

## 2. Building spectral energy distributions

The Hi-z LBG sample is inhomogeneous and collected from several papers as listed in Table 1. It corresponds to an FUV selection. The criteria are as follows: we selected all LBGs to have spectroscopic redshifts in the $5 \le z \le 10$ range for which a mm detection was attempted (successful or not) in the rest-frame FIR (Table 1). Practically speaking, this means that we focused on LBGs observed in ALMA's bands 6 and 7. This redshift range (at EoR) provides constraints in rest-frame FIR, meaning, observed in the mm range. As an additional criterion, we only use LBGs with at least five data points in the rest-frame UV and optical.

We collected photometric data to build SEDs and derive stellar masses ($M_{star}$), dust masses ($M_{dust}$), star formation rates (SFRs), FUV and luminosities ($L_{FUV}$ and $L_{dust}$), the age of the dominant stellar population ($A_{ge_{main}}$), and the UV slope $\beta$ (Calzetti et al. 1994). The set of data is consistently fitted with the same code, CIGALE (Burgarella et al. 2005; Noll et al. 2009; Boquien et al. 2019), and using the same dust models (Draine & Li 2007; Draine et al. 2014). The CIGALE modules and input parameters used for all the fits are described in Table 2. To model the dust emission, we also tested a modified, single dust temperature grey body, plus a mid-IR power law that approximates the hot-dust emission from the AGN (Active Galactic Nuclei) heating or clumpy hot star bursting regions (Casey 2012). The impact is negligible, and we base this paper on the former to constrain our models. The dust models are calibrated with a sub-sample of dwarf (compact with sizes ranging from 0.08' to 1.2', Madden et al. 2013) low-metallicity galaxies (Low-zZ) with $12 + \log(O/H) \le 8.4$ in the Dwarf Galaxy Survey (DGS, Madden et al. 2013; Rémy-Ruyer et al. 2015). The DGS was originally selected from the Herschel Dwarf Galaxy Survey. It is a compilation of a SPIRE SAG2 guaranteed time (GT) key program and SPIRE GT2 observations, using the PACS and SPIRE instruments onboard the Herschel Space Observatory to obtain 55–500 µm photometry and spectroscopy of 48 dwarf galaxies in our local Universe, chosen to cover a broad range of physical conditions. We add data from the NASA Extragalactic Database (NED) in the UV (GALEX) and in the optical (B and R bands). When selecting the FUV, NUV, B, and R data in NED, we kept the data that provide total fluxes for our galaxies, and, when several total fluxes were available, we favoured new or newly reprocessed ones. If these criteria are not respected, no data were added for the corresponding objects, which in the end meant that we did not keep them in the Low-zZ sample.

The IR information on these Hi-z LBGs is limited. Our approach relies on the assumption that the physical conditions for dust grains are shared by these LBGs, which means that these LBGs present the same FIR SED shape within uncertainties. Here, we used luminosity-weighted SEDs normalised at $\lambda = 200 \mu m$, where the dust is optically thin (e.g. Casey 2012). Our LBG sample presents different FIR-to-FUV ratios, and the dust attenuations are different. We leverage the redshift range to combine and build the IR SEDs of our ALMA-detected LBGs, and thus gain access to better spectral information for these objects (Fig. 1). However, we have no observational information on the dust emission, at wavelengths shorter than the peak, that is, near- and mid-IR. From the assumption that these Hi-z LBGs are metal-poor (Castellano et al. 2014; Yuan et al. 2019), we use the minimum value for $q_{PAH}$ (fraction of the dust mass in the PAH) because $q_{PAH}$ is correlated with the metallicity (e.g. Ciesla et al. 2014). From our SEDs, $\gamma$, the fraction of the dust heated by starlight above the lower cutoff $U_{min}$ (the minimum of the distribution of starlight intensity relative to the local interstellar radiation field), in Draine and Li's models cannot be constrained. We assume that the best value for the DGS sample is valid for the Hi-z LBGs (but the impact of this parameter on our result is limited to ±8%, estimated using the range found by Álvarez-Márquez et al. (2019) on a sample of 22000 LBGs at $z \sim 3.0$). In addition to the dust continuum information, we used the emission in the rest-frame FUV and optical range to model the stellar populations. Whenever available, [CII]158 µm and [OIII]88.3 µm were used to constrain the SED fitting. The nebular module in CIGALE uses nebular templates based on Inoue (2011), which were generated using CLOUDY 08.00 (Ferland et al. 1998). Each line has a Gaussian shape with
<table><caption>**Table 1.** Origins of data used in this paper.</caption><thead><tr><th>Source</th><th>Number of objects</th><th>Selection</th><th>Redshift</th><th>Note</th></tr></thead><tbody><tr><td>Bouwens et al. (2016)</td><td>3</td><td>UV</td><td>5 &le; z &le; 10</td><td>* ALMA 6 upper limits<br/>* S/N<sub>HST</sub> &gt; 1.5<br/>* LBGs with &gt;5 data points in UV-optical only</td></tr><tr><td>Capak et al. (2015) and Faisst et al. (2017)</td><td>9 (HZ1-HZ4 and HZ6-HZ10)</td><td>UV</td><td>z ~ 5.6</td><td>* [CII]158 &mu;m for all Hi-z LBGs<br/>* ALMA 7 detections: HZ4, HZ6 (3) HZ9 and HZ10 (5)<br/>* ALMA 7 upper limits for the others<br/>* HZ5 detected in Chandra and not included in the sample<br/>* Additional data from Pavesi et al. (2016)</td></tr><tr><td>Scoville et al. (2016)</td><td>1 (566428)</td><td>UV</td><td>z = 5.89</td><td>* ALMA 6 detection<br/>* [CII]158 &mu;m measurement</td></tr><tr><td>Willott et al. (2015)</td><td>2 (CLM1 and WMH5)</td><td>UV</td><td>z ~ 6.1</td><td>* ALMA 6 detections<br/>* [CII]158 &mu;m measurement</td></tr><tr><td>Aravena et al. (2016)</td><td>2 (ID27, ID31)</td><td>UV</td><td>z ~ 7.5</td><td>* ALMA 6 detections<br/>* [CII]158 &mu;m measurement</td></tr><tr><td>Hashimoto et al. (2018)</td><td>1 (MACSJD1149)</td><td>UV</td><td>z = 9.1</td><td>* ALMA 7 upper limit<br/>* [OIII]88 &mu;m measurement</td></tr><tr><td>Rémy-Ruyer et al. (2015)</td><td>31</td><td>DGS and Z &lt; 0.1 Z<sub>&odot;</sub></td><td>z ~ 0</td><td>* GALEX FUV and NUV data from NED<br/>* B and R data from NED<br/>* IR data (including Herschel) from Rémy-Ruyer et al. (2013) and Rémy-Ruyer et al. (2015)</td></tr></tbody></table>

**Notes.** The Hi-z LBG sample contains 18 objects, while the Low-zZ sample contains 31 objects for which the final fits are adequate, i.e. $\chi^2_{reduced} \le 5$. 9 Hi-z LBGs are detected in ALMA in Band 6 or Band 7, they provide 15 detected measurements, but only 11 measurements with $S/N_{submm} > 3$ are used to build the IR template.

a user-defined line width to take rotation into account, which can be especially important for narrowband filters and high-redshift objects due to the apparent line broadening with redshift in the observed frame. This normalised nebular emission-line spectrum is rescaled to the appropriate flux level by multiplying by the ionising photon luminosity that was computed along with the composite stellar population.

### 3. Methodology to derive the physical parameters

The SED analysis is performed in five phases.

**Phase 0.** A necessary preliminary phase was to perform a lot of initial fits on individual Low-zZ and Hi-z LBGs to understand the limits of the physical parameters to be explored (star formation history (SFH), dust attenuation, AGN fraction, etc.) in the CIGALE fits. Because of the exploration nature of these fits, these initial fits are not reported in Table 2. In order to perform this exploration of the parameter space, we first used different types of SFHs available in CIGALE: a periodic SFH, that is the same star formation event regularly occurring during the galaxy lifetime and a single or double exponential SFH. None of them provided an improvement of the fits. We also attempted to modify the metallicity of the SSPs but we did not manage to get any

strong constraint. Finally, we tested various dust emission models within the set available into CIGALE. But the limited number of FIR data points for individual galaxies prevented us to carry on this study and led us to the template approach described below.

**Phase 1.** CIGALE fits individual Hi-z LBGs with a large distribution of priors for the dust parameters and SFH (delayed<fnref n="1" />, delayed + burst, and constant SFHs are assumed) with average reduced $\chi^2$ < 5. On the other hand, constant SFHs provide results very close to the delayed ones. This was expected, given the ages of the stellar populations and the age of the universe. Fits with a delayed + burst SFH are not as good. As none of the stellar models available in CIGALE feature a top-heavy IMF, we used a Chabrier IMF. Using a top-heavy IMF in the SED fitting instead of a Chabrier one could change the results of the SED fitting in the UV-optical range, but less so in the IR range. Future development of CIGALE will address this issue. Similarly, stellar models that include binaries (e.g. BPASS Stanway & Eldridge 2018) would change the UV-optical emission of these galaxies. We will address this possibility in a future paper. These initial fits provide an estimated rest-frame flux<sub>200&mu;m</sub> to normalise the SEDs (with FIR detections) at $\lambda = 200&mu;m$

<footnote n="1">¹ SFR = $\frac{1}{\tau^2}e^{-\frac{1}{\tau}}$.</footnote>
<table><caption>**Table 2.** CIGALE modules and input parameters used for all the fits.</caption><thead><tr><th>Parameters</th><th>Symbol</th><th>Range (Ph. 1, Hi-z LBGs)</th><th>Range (Ph. 2, Low-zZ)</th><th>Range (Ph. 3, Hi-z LBGs)</th><th>Range (Ph. 4, Hi-z LBGs)</th></tr></thead><tbody><tr><td>Target sample</td><td></td><td>Individual Hi-z LBGs</td><td>Individual Low-zZ</td><td>Combined Hi-z LBGs</td><td>Individual Hi-z LBGs</td></tr><tr><td>e-folding timescale of the delayed SFH</td><td>τ<sub>main</sub> [Myr]</td><td>Delayed SFH and recent burst</td><td>25, 50, 100, 250, 500, 1000</td><td>25, 50, 100, 250, 500, 1000</td><td>25, 50, 100, 250, 500, 1000</td></tr><tr><td>Age of the main population</td><td>Age<sub>main</sub>[Myr]</td><td>61 log values in [2–1000]</td><td>61 log values in [2–1500]</td><td>17 log values in [1, 800]</td><td>101 log values in [50, 1260]</td></tr><tr><td>Burst</td><td>f<sub>burst</sub></td><td>No burst</td><td>No burst</td><td>No burst</td><td>No burst</td></tr><tr><td>SSP</td><td></td><td>SSP</td><td></td><td></td><td></td></tr><tr><td>Initial mass function</td><td>IMF</td><td>BC03</td><td>BC03</td><td>BC03</td><td>BC03</td></tr><tr><td>Metallicity</td><td>Z</td><td>Chabrier</td><td>Chabrier</td><td>Chabrier</td><td>Chabrier</td></tr><tr><td>Ionisation parameter</td><td>logU</td><td>-1.5</td><td>-1.5</td><td>-1.5</td><td>-1.5</td></tr><tr><td>Line width [km s<sup>-1</sup>]</td><td>–</td><td>100</td><td>100</td><td>100</td><td>100</td></tr><tr><td>Colour excess for both the old and young stellar populations</td><td>E_BV_lines</td><td>61 log values in [0.001, 0.3]</td><td>61 log values in [0.001, 1.0]</td><td>11 log values in [1e-3, 1e0]</td><td>101 log values in [1e-3, 1e0]</td></tr><tr><td>Bump amplitude</td><td>uv_bump_amplitude</td><td>0.0</td><td>0.0</td><td>0.0</td><td>0.0</td></tr><tr><td>Power law slope</td><td>power law slope</td><td>-0.7, -0.35, 0.0, 0.35, 0.7</td><td>0.0</td><td>0.0</td><td>0.0</td></tr><tr><td>Mass fraction of PAH</td><td>q<sub>PAH</sub></td><td>0.47,1.12, 1.77, 2.50, 3.19</td><td>0.47,1.12, 1.77, 2.50, 3.19</td><td>0.47</td><td>0.47</td></tr><tr><td>Minimum radiation field</td><td>U<sub>min</sub></td><td>0.1, 0.5, 1.0, 5.0, 10.0, 50.0</td><td>0.1, 0.5, 1.0, 5.0, 10.0, 50.0</td><td>4.0</td><td>4.0</td></tr><tr><td>Power law slope dU/dM ≈ U<sup>α</sup></td><td>α</td><td>2.0, 2.2, 2.4, 2.6, 2.8, 3.0</td><td>2.0, 2.2, 2.4, 2.6, 2.8, 3.0</td><td>2.2</td><td>2.2</td></tr><tr><td>Dust fraction in PDRs</td><td>γ</td><td>0.00, 0.01, 0.1, 0.2, 0.3, 0.4, 0.5, 0.6, 0.7, 0.8, 0.9, 1.0</td><td>0.00, 0.01, 0.1, 0.2, 0.3, 0.4, 0.5, 0.6, 0.7, 0.8, 0.9, 1.0</td><td>0.75</td><td>0.75</td></tr></tbody></table>

**Notes.** Note that we did not include any AGN emission.

**References.** BC03: Bruzual & Charlot (2003); Chabrier IMF: Chabrier (2003).

where the dust is optically thin (e.g. Casey 2012). Detected LBGs at $5 \le z \le 10$ made it possible to sample the spectrum of dust emission (Fig. 1). Since we do not have data beyond the rest-frame wavelength $\lambda = 200 \mu m$, normalizing at any longer wavelength would not be meaningful. Although only the Rayleigh-Jeans range is covered, dust temperatures in the range 40 K to 70 K fit the data better.

**Phase 2.** Local Low-zZ galaxies are often said to be analogues to early galaxies (Madden et al. 2013; Capak et al. 2015; Hou et al. 2019). We compiled a sample of Low-zZ galaxies from Rémy-Ruyer et al. (2013, 2015). We kept DGS galaxies with $12 + \log(O/H) < 8.4$. With the same assumptions on the priors as for the LBGs, these galaxies were fitted with CIGALE. We checked that the main dust physical parameters found for the Low-zZ sample were in agreement with the ones from the Hi-z LBG sample (which confirms the partial analogy), and we leveraged the former fits to set the value of the $\gamma$ parameter in the models from Draine & Li (2007) and Draine et al. (2014). We also checked that the value assumed for $q_{PAH}$ is consistent with our assumption of low metallicity (Fig. 2).

**Phase 3.** The previous phases suggested that a common template could be used to fit most of the objects in our samples. We interpret this possible generalization of the template as a similarity in the physical conditions in the Hi-z LBG sample. Note however, that there does not seem to be an universality as some objects present different dust emission (for instance SBS0335-052) and it would interesting to explore the origin of the difference). In this phase, $q_{PAH}$ and $\gamma$ are fixed to

the value given in Table 2, and we fitted the IR-combined SEDs of the Hi-z LBGs (Fig. 3). Phase 3 provides physical parameters for the average combined SEDs of the Hi-z LBGs used to build the Hi-z IR template. The fit with a modified black body gives: $\beta_{submm} = 1.51$ (slope on rest-frame sub-mm side of the dust emission) and $T_{dust} = 56.6$ K (Fig. 4).

**Phase 4.** We re-fitted all the Hi-z LBGs over the FUV-to-FIR range with this Hi-z IR template (Tables 3 and 4). As a safety check, we also fit the Low-zZ galaxies with the same IR template. The quality of the fit is shown in the Annex.

In conclusion, after studying the individual dust emission of the Hi-z LBG sample, we found some similarities with the local Low-zZ sample. Furthermore, we find that most of the Hi-z LBG sample show a similar dust emission. We have characterised this dust emission using a modified blackbody and a single model from the library of Draine & Li (2007), Draine et al. (2014). We use the latter in the rest of the paper as it presents a more physical understanding of the physical conditions in the ISM of galaxies in the early universe. The above methodology was successful in fitting the majority of the Hi-z and Low-zZ samples.

## 4. Results

From Table 3 and 4, we compute the specific SFR (sSFR = SFR/$M_{star}$), the specific dust mass ($sM_{dust} = M_{dust}/M_{star}$), and $IRX = L_{dust}/L_{FUV}$. In Figs. 5 and 6, we define the dust-formation rate diagram (DFRD): $sM_{dust}$ vs. sSFR. The x-axis provides information on the stellar population and the y-axis on dust
<figure><img src="image_2.png" /></figure>

**Fig. 1.** Dust emission of the detected Hi-z LBGs with $S/N_{submm} > 3$ (11 data points) are normalised at $\lambda_{rest-frame} = 200~\mu$m and combined to form a single Hi-z LBG template. For each detected LBG in our sample, we use the number of $S/N_{submm} > 3$ ALMA detections, which are validated in agreement with our criteria, as described in Table 1. The redshift range ($0.5 \le z \le 10$) makes it possible to cover a wide wavelength range and explore dust emission. We over-plot colour-coded curves for SEDs based on modified blackbodies in the range $30 \le T_{dust} [K] \le 85$.

<figure><img src="image_3.png" /></figure>

**Fig. 2.** Distribution of $\gamma_{Low-zZ}$ (left panel) and $q_{PAH\_Low-zZ}$ (right panel) dust parameters from Draine & Li's models for the Low-zZ galaxy sample. The $\gamma_{Low-zZ}$ distribution peaks at $0.75 \pm 0.07$ and most of the $q_{PAH} = 0.47$, which is the lowest value in Draine & Li's models. $\gamma_{Low-zZ}$ is assumed to be valid for the Hi-z LBG sample, and $q_{PAH\_Low-zZ}$ is fully consistent with our assumed $q_{PAH\_Hi-z}$.

grains. For some of the galaxies, the fit is not good, for example, SBS 035-052 (see Sect. 4). The two possible origins for these bad fits are the quality of the data, and/or the fact that our template does not match the dust properties of these objects. These galaxies are not used in our analysis.

The physical parameters derived for the Hi-z LBG and the Low-zZ DGS samples by CIGALE are listed in Tables 3, 4, and C.1.

By definition, CIGALE checks that the observed SED is consistent with the assumed SFH, and we tried several SFHs before selecting the delayed one. A first result from the SED fitting, both for the Hiz-LBGs and the Low-zZ DGS objects is that both $\tau_{main}$ and the range of ages are quite short (Tables 3, 4 and C.1) with ages in the range of 100–500 Myrs for the Hiz-LBGs, and 35–1300 Myrs for the Low-zZ sample. From this result, two important points should be noted: (i) we know that some young, and even very young, stellar populations are present in some of these objects. The most extreme one is the galaxy at $z = 9.1$ (Hashimoto et al. 2018) and observational evidences ([OIII]88.3 $\mu$m line, for instance) suggest the presence of a very young stellar population (a few Myrs). However, this paper deals with dust characteristics and any dust produced/destroyed over such short timescales do not impact on the dust budget. This is why we assume a delayed (no final bursts) star formation history in the SED fitting; (ii) this age range is in agreement with the properties of this sample as they are dwarf, low-metallicity, star-forming galaxies (Maiolino et al. 2008; Madden et al. 2013; Cormier et al. 2019). This is the reason why they are often cited as being high-redshift LBG analogues.

<figure><img src="image_4.png" /></figure>

<figure><img src="image_5.png" /></figure>

**Fig. 3.** Fit of IR template SED built from individual ones normalised at 200 $\mu$m and combined. We used Draine and Li's models in the top panel and a modified black body Casey (2012) in the bottom panel.

<figure><img src="image_6.png" /></figure>

**Fig. 4.** Probability distribution functions for the dust temperature and the emissivity derived from the fit of the Hi-z LBGs with CIGALE and a modified black body Casey (2012).

Figure 5 is colour-coded by $age/\tau_{main}$ of the stellar population, assuming a delayed SFH. Both panels of Fig. 5 exhibit an age-evolutionary sequence with younger objects at higher sSFR. To the right-hand side of both panels, $sM_{dust}$ presents (or reaches) an apparent maximum at sSFR $\sim 10^{-8} \text{ yr}^{-1}$ with a decline at lower sSFRs. Figure 6 (top and central panels) is colour-coded by IRX. Galaxies have low IRX to the right-hand side, meaning at large sSFRs, but also to the left-hand side, meaning at low sSFRs. In between, objects have larger IRX values.

Combining the information from Figs. 5 and 6, a qualitative scenario seems to emerge. We first focus on the Hi-z LBG sample, which is the core objective of this work: Hi-z LBGs with
<table><caption>**Table 3.** Parameters derived from fitting the SEDs of the Hi-z LBG sample.</caption><thead><tr><th>Id</th><th>χ<sup>2</sup><sub>ν</sub></th><th>sSFR(yr<sup>-1</sup>)</th><th>sSFR err[yr<sup>-1</sup>]</th><th>sM<sub>dust</sub></th><th>sM<sub>dust_err</sub></th><th>IRX</th><th>IRX_err</th><th>A<sub>FUV</sub></th><th>A<sub>FUV_err</sub></th><th>β<sub>cal94</sub></th><th>Age[Myr]</th><th>Age_err[Myr]</th><th>Age/τ</th><th>Age/τ_err</th><th>12+log<sub>10</sub> (O/H) <sup>(a)</sup></th></tr></thead><tbody><tr><td>ID27</td><td>0.442</td><td>7.126 × 10<sup>-8</sup></td><td>7.126 × 10<sup>-8</sup></td><td>0.026</td><td>0.026</td><td>0.277</td><td>0.277</td><td>0.105</td><td>0.123</td><td>-2.616</td><td>0.060</td><td>34.498</td><td>43.313</td><td>0.443</td><td>0.428</td><td>7.627</td></tr><tr><td>ID31</td><td>0.155</td><td>3.701 × 10<sup>-8</sup></td><td>3.699 × 10<sup>-8</sup></td><td>0.009</td><td>0.009</td><td>0.145</td><td>0.145</td><td>0.053</td><td>0.068</td><td>-2.613</td><td>0.052</td><td>53.981</td><td>56.521</td><td>0.651</td><td>0.534</td><td>7.540</td></tr><tr><td>lbg10</td><td>0.272</td><td>2.120 × 10<sup>-8</sup></td><td>1.998 × 10<sup>-8</sup></td><td>0.046</td><td>0.037</td><td>4.597</td><td>2.823</td><td>1.423</td><td>0.622</td><td>-2.007</td><td>0.107</td><td>106.721</td><td>96.855</td><td>1.724</td><td>1.822</td><td>8.003</td></tr><tr><td>HZ4</td><td>0.545</td><td>3.158 × 10<sup>-9</sup></td><td>9.228 × 10<sup>-10</sup></td><td>0.006</td><td>0.002</td><td>0.905</td><td>0.295</td><td>0.538</td><td>0.149</td><td>-2.200</td><td>0.058</td><td>302.940</td><td>137.273</td><td>2.763</td><td>1.018</td><td>7.803</td></tr><tr><td>HZ9</td><td>3.260</td><td>7.867 × 10<sup>-9</sup></td><td>2.008 × 10<sup>-9</sup></td><td>0.029</td><td>0.007</td><td>4.936</td><td>0.638</td><td>1.662</td><td>0.138</td><td>-1.848</td><td>0.073</td><td>101.509</td><td>42.765</td><td>2.874</td><td>0.751</td><td>8.011</td></tr><tr><td>lbg53</td><td>1.134</td><td>7.293 × 10<sup>-9</sup></td><td>6.141 × 10<sup>-9</sup></td><td>0.015</td><td>0.014</td><td>1.285</td><td>1.058</td><td>0.627</td><td>x 0.429</td><td>-2.226</td><td>0.081</td><td>134.700</td><td>83.369</td><td>2.908</td><td>1.879</td><td>7.846</td></tr><tr><td>MACSJD1149</td><td>0.470</td><td>2.969 × 10<sup>-9</sup></td><td>1.059 × 10<sup>-9</sup></td><td>0.010</td><td>0.004</td><td>3.002</td><td>0.925</td><td>1.229</td><td>0.251</td><td>-1.900</td><td>0.090</td><td>278.095</td><td>126.282</td><td>3.067</td><td>1.047</td><td>7.950</td></tr><tr><td>566428</td><td>0.309</td><td>2.358 × 10<sup>-9</sup></td><td>1.599 × 10<sup>-9</sup></td><td>0.007</td><td>0.005</td><td>1.543</td><td>1.127</td><td>0.734</td><td>0.444</td><td>-2.038</td><td>0.106</td><td>299.734</td><td>181.600</td><td>3.751</td><td>1.754</td><td>7.868</td></tr><tr><td>HZ2</td><td>2.282</td><td>8.208 × 10<sup>-10</sup></td><td>2.117 × 10<sup>-10</sup></td><td>2.326 × 10<sup>-4</sup></td><td>2.320 × 10<sup>-4</sup></td><td>0.094</td><td>0.094</td><td>0.041</td><td>0.051</td><td>-2.220</td><td>0.056</td><td>691.970</td><td>217.765</td><td>3.836</td><td>0.778</td><td>7.495</td></tr><tr><td>HZ10</td><td>3.652</td><td>2.309 × 10<sup>-9</sup></td><td>3.945 × 10<sup>-10</sup></td><td>0.016</td><td>0.002</td><td>10.397</td><td>1.120</td><td>2.347</td><td>0.140</td><td>-1.373</td><td>0.081</td><td>119.493</td><td>6.843</td><td>4.779</td><td>0.273</td><td>8.103</td></tr><tr><td>HZ8</td><td>1.735</td><td>9.080 × 10<sup>-10</sup></td><td>2.268 × 10<sup>-10</sup></td><td>3.180 × 10<sup>-4</sup></td><td>3.180 × 10<sup>-4</sup></td><td>0.065</td><td>0.065</td><td>0.042</td><td>0.047</td><td>-2.196</td><td>0.037</td><td>290.782</td><td>100.363</td><td>5.028</td><td>0.604</td><td>7.462</td></tr><tr><td>HZ7</td><td>0.593</td><td>7.686 × 10<sup>-10</sup></td><td>3.740 × 10<sup>-10</sup></td><td>8.429 × 10<sup>-4</sup></td><td>8.429 × 10<sup>-4</sup></td><td>0.200</td><td>0.200</td><td>0.130</td><td>0.130</td><td>-2.129</td><td>0.054</td><td>303.331</td><td>159.922</td><td>5.415</td><td>1.157</td><td>7.609</td></tr><tr><td>HZ1</td><td>0.484</td><td>4.871 × 10<sup>-10</sup></td><td>1.471 × 10<sup>-10</sup></td><td>2.367 × 10<sup>-4</sup></td><td>2.366 × 10<sup>-4</sup></td><td>0.089</td><td>0.089</td><td>0.053</td><td>0.060</td><td>-2.072</td><td>0.049</td><td>440.325</td><td>152.345</td><td>5.421</td><td>0.667</td><td>7.496</td></tr><tr><td>lbg48</td><td>0.216</td><td>9.947 × 10<sup>-10</sup></td><td>9.946 × 10<sup>-10</sup></td><td>0.004</td><td>0.004</td><td>1.499</td><td>1.499</td><td>0.477</td><td>0.521</td><td>-1.904</td><td>0.084</td><td>335.492</td><td>150.453</td><td>5.825</td><td>1.978</td><td>7.825</td></tr><tr><td>WMH5</td><td>1.951</td><td>8.477 × 10<sup>-10</sup></td><td>1.493 × 10<sup>-10</sup></td><td>0.006</td><td>0.001</td><td>2.084</td><td>3.030</td><td>1.024</td><td>0.121</td><td>-1.787</td><td>0.070</td><td>154.089</td><td>22.361</td><td>5.988</td><td>0.295</td><td>7.906</td></tr><tr><td>HZ3</td><td>0.837</td><td>5.131 × 10<sup>-10</sup></td><td>3.539 × 10<sup>-10</sup></td><td>0.001</td><td>0.001</td><td>0.305</td><td>0.305</td><td>0.168</td><td>0.184</td><td>-2.030</td><td>0.078</td><td>289.824</td><td>151.571</td><td>6.251</td><td>1.424</td><td>7.646</td></tr><tr><td>CLM1</td><td>0.990</td><td>3.938 × 10<sup>-10</sup></td><td>1.150 × 10<sup>-10</sup></td><td>1.166 × 10<sup>-4</sup></td><td>7.915 × 10<sup>-5</sup></td><td>0.029</td><td>0.018</td><td>0.023</td><td>0.015</td><td>-2.079</td><td>0.062</td><td>224.548</td><td>70.546</td><td>6.628</td><td>0.358</td><td>7.380</td></tr></tbody></table>

Notes. We note that the parameters for which the uncertainty (e.g. param_err) is ≥ parameter (i.e. param) should be considered as upper limits. <sup>(a)</sup>Estimated from the relation described in the text. However, this method being very indirect (see text), these values should only be considered as indicative.

<table><caption>**Table 4.** Luminosity, mass and SFR for the Hi-z LBG sample from the CIGALE SED fitting.</caption><thead><tr><th>Id</th><th>χ<sup>2</sup><sub>ν</sub></th><th>L<sub>dust</sub>[L<sub>☉</sub>]</th><th>L<sub>dust_err</sub>[L<sub>☉</sub>]</th><th>L<sub>FUV</sub>[L<sub>☉</sub>]</th><th>L<sub>FUV_err</sub>[L<sub>☉</sub>]</th><th>M<sub>dust</sub>[M<sub>☉</sub>]</th><th>M<sub>dust_err</sub>[M<sub>☉</sub>]</th><th>M<sub>star</sub>[M<sub>☉</sub>]</th><th>M<sub>star_err</sub>[M<sub>☉</sub>]</th><th>SFR[M<sub>☉</sub>yr<sup>-1</sup>]</th><th>SFR_err[M<sub>☉</sub>yr<sup>-1</sup>]</th></tr></thead><tbody><tr><td>ID27</td><td>0.442</td><td>6.352 × 10<sup>9</sup></td><td>8.123 × 10<sup>9</sup></td><td>3.759 × 10<sup>10</sup></td><td>1.880 × 10<sup>9</sup></td><td>2.335 × 10<sup>6</sup></td><td>2.986 × 10<sup>6</sup></td><td>1.549 × 10<sup>8</sup></td><td>1.725 × 10<sup>8</sup></td><td>9.705</td><td>2.255</td></tr><tr><td>ID31</td><td>0.155</td><td>4.827 × 10<sup>9</sup></td><td>6.545 × 10<sup>9</sup></td><td>6.133 × 10<sup>10</sup></td><td>3.066 × 10<sup>9</sup></td><td>1.775 × 10<sup>9</sup></td><td>2.406 × 10<sup>6</sup></td><td>3.342 × 10<sup>8</sup></td><td>3.277 × 10<sup>8</sup></td><td>12.368</td><td>2.391</td></tr><tr><td>lbg10</td><td>0.272</td><td>7.967 × 10<sup>10</sup></td><td>4.876 × 10<sup>10</sup></td><td>1.733 × 10<sup>10</sup></td><td>8.665 × 10<sup>8</sup></td><td>2.929 × 10<sup>7</sup></td><td>1.793 × 10<sup>7</sup></td><td>6.409 × 10<sup>8</sup></td><td>3.425 × 10<sup>8</sup></td><td>13.590</td><td>10.544</td></tr><tr><td>HZ4</td><td>0.545</td><td>1.965 × 10<sup>11</sup></td><td>6.329 × 10<sup>10</sup></td><td>2.171 × 10<sup>11</sup></td><td>1.085 × 10<sup>10</sup></td><td>7.224 × 10<sup>7</sup></td><td>2.326 × 10<sup>7</sup></td><td>1.169 × 10<sup>10</sup></td><td>2.468 × 10<sup>9</sup></td><td>36.933</td><td>7.462</td></tr><tr><td>HZ9</td><td>3.260</td><td>5.255 × 10<sup>11</sup></td><td>6.259 × 10<sup>10</sup></td><td>1.065 × 10<sup>11</sup></td><td>5.323 × 10<sup>9</sup></td><td>1.932 × 10<sup>8</sup></td><td>2.301 × 10<sup>7</sup></td><td>6.645 × 10<sup>9</sup></td><td>1.473 × 10<sup>9</sup></td><td>52.271</td><td>6.617</td></tr><tr><td>lbg53</td><td>1.134</td><td>4.525 × 10<sup>10</sup></td><td>3.716 × 10<sup>10</sup></td><td>3.521 × 10<sup>10</sup></td><td>1.760 × 10<sup>9</sup></td><td>1.663 × 10<sup>7</sup></td><td>1.366 × 10<sup>7</sup></td><td>1.144 × 10<sup>9</sup></td><td>4.967 × 10<sup>8</sup></td><td>8.347</td><td>6.023</td></tr><tr><td>MACSJD1149</td><td>0.470</td><td>2.963 × 10<sup>11</sup></td><td>9.010 × 10<sup>10</sup></td><td>9.870 × 10<sup>10</sup></td><td>4.935 × 10<sup>9</sup></td><td>1.089 × 10<sup>8</sup></td><td>3.312 × 10<sup>7</sup></td><td>1.038 × 10<sup>10</sup></td><td>2.595 × 10<sup>9</sup></td><td>30.810</td><td>7.835</td></tr><tr><td>566428</td><td>0.309</td><td>2.557 × 10<sup>11</sup></td><td>1.856 × 10<sup>11</sup></td><td>1.657 × 10<sup>11</sup></td><td>1.283 × 10<sup>10</sup></td><td>9.401 × 10<sup>7</sup></td><td>6.822 × 10<sup>7</sup></td><td>1.404 × 10<sup>10</sup></td><td>3.351 × 10<sup>9</sup></td><td>33.097</td><td>21.002</td></tr><tr><td>HZ2</td><td>2.282</td><td>9.420 × 10<sup>7</sup></td><td>1.189 × 10<sup>10</sup></td><td>1.603 × 10<sup>11</sup></td><td>8.016 × 10<sup>9</sup></td><td>3.463 × 10<sup>9</sup></td><td>4.371 × 10<sup>9</sup></td><td>1.900 × 10<sup>10</sup></td><td>3.587 × 10<sup>9</sup></td><td>15.596</td><td>2.741</td></tr><tr><td>HZ10</td><td>3.652</td><td>1.275 × 10<sup>12</sup></td><td>8.743 × 10<sup>10</sup></td><td>1.227 × 10<sup>11</sup></td><td>1.019 × 10<sup>10</sup></td><td>4.688 × 10<sup>8</sup></td><td>3.214 × 10<sup>7</sup></td><td>3.018 × 10<sup>7</sup></td><td>4.225 × 10<sup>9</sup></td><td>69.682</td><td>6.826</td></tr><tr><td>HZ8</td><td>1.735</td><td>9.815 × 10<sup>9</sup></td><td>1.134 × 10<sup>10</sup></td><td>1.747 × 10<sup>11</sup></td><td>8.734 × 10<sup>9</sup></td><td>3.608 × 10<sup>9</sup></td><td>4.168 × 10<sup>7</sup></td><td>1.322 × 10<sup>10</sup></td><td>2.055 × 10<sup>9</sup></td><td>12.006</td><td>2.347</td></tr><tr><td>HZ7</td><td>0.593</td><td>2.831 × 10<sup>10</sup></td><td>3.034 × 10<sup>10</sup></td><td>1.519 × 10<sup>11</sup></td><td>7.594 × 10<sup>9</sup></td><td>1.041 × 10<sup>7</sup></td><td>1.115 × 10<sup>7</sup></td><td>1.342 × 10<sup>10</sup></td><td>2.410 × 10<sup>9</sup></td><td>10.312</td><td>4.663</td></tr><tr><td>HZ1</td><td>0.484</td><td>1.298 × 10<sup>10</sup></td><td>1.544 × 10<sup>10</sup></td><td>1.739 × 10<sup>11</sup></td><td>8.694 × 10<sup>9</sup></td><td>4.773 × 10<sup>9</sup></td><td>5.676 × 10<sup>6</sup></td><td>2.423 × 10<sup>10</sup></td><td>4.157 × 10<sup>9</sup></td><td>11.801</td><td>2.932</td></tr><tr><td>HZ6</td><td>2.196</td><td>3.475 × 10<sup>11</sup></td><td>5.015 × 10<sup>10</sup></td><td>3.346 × 10<sup>10</sup></td><td>1.573 × 10<sup>10</sup></td><td>1.278 × 10<sup>8</sup></td><td>6.388 × 10<sup>6</sup></td><td>2.922 × 10<sup>10</sup></td><td>3.131 × 10<sup>9</sup></td><td>29.039</td><td>7.748</td></tr><tr><td>lbg48</td><td>0.216</td><td>3.614 × 10<sup>10</sup></td><td>5.828 × 10<sup>10</sup></td><td>2.056 × 10<sup>11</sup></td><td>1.051 × 10<sup>10</sup></td><td>1.575 × 10<sup>8</sup></td><td>1.843 × 10<sup>7</sup></td><td>5.096 × 10<sup>9</sup></td><td>1.021 × 10<sup>9</sup></td><td>4.106</td><td>5.002</td></tr><tr><td>WMH5</td><td>1.951</td><td>4.284 × 10<sup>11</sup></td><td>5.828 × 10<sup>10</sup></td><td>2.056 × 10<sup>11</sup></td><td>6.540 × 10<sup>9</sup></td><td>1.217 × 10<sup>7</sup></td><td>1.465 × 10<sup>7</sup></td><td>2.142 × 10<sup>7</sup></td><td>2.677 × 10<sup>10</sup></td><td>3.744 × 10<sup>9</sup></td><td>22.695</td></tr><tr><td>HZ3</td><td>0.837</td><td>3.310 × 10<sup>10</sup></td><td>3.986 × 10<sup>10</sup></td><td>1.308 × 10<sup>11</sup></td><td>2.109 × 10<sup>10</sup></td><td>2.859 × 10<sup>9</sup></td><td>1.807 × 10<sup>6</sup></td><td>2.451 × 10<sup>10</sup></td><td>6.059 × 10<sup>9</sup></td><td>7.356</td><td>4.944</td></tr><tr><td>CLM1</td><td>0.990</td><td>7.776 × 10<sup>9</sup></td><td>4.916 × 10<sup>9</sup></td><td>2.693 × 10<sup>11</sup></td><td></td><td></td><td></td><td></td><td></td><td></td><td>1.500</td></tr></tbody></table>

Notes. As in Table 3, the parameters for which the uncertainty (e.g. param_err) is ≥ parameter (i.e. param) should be considered as upper limits.
<figure><img src="image_7.png" /></figure>

**Fig. 5.** Hi-z LBGs (filled dots, top) and Low-zZ galaxies (filled boxes, panel) DFRD, $sM_{\text{dust}}$ vs. sSFR, colour-coded by age. We see a sequence from top right to bottom left. Using a colour-coding system based on the age of the dominant stellar population, we see that this sequence corresponds to galaxies getting older and older from top right to bottom left.

low IRX, low $sM_{\text{dust}}$, and low ages are found at high sSFR. We also find Hi-z LBGs with low IRX and even lower $sM_{\text{dust}}$, but large ages are found at low sSFR. In between, we see an apparent maximum observed in IRX and $sM_{\text{dust}}$ for the present Hi-z LBG sample at sSFR $\sim 10^{-8} \text{ yr}^{-1}$ that needs to be confirmed by future observational data. Both IRX and $sM_{\text{dust}}$ decline, and age increases to a locus at sSFR $\sim 10^{-10}-10^{-9} \text{ yr}^{-1}$, where we find the LBGs with upper limits like those discussed in Capak et al. (2015), Ferrara et al. (2017) and Faiset et al. (2017). This seems to define an evolutionary sequence from the right of the DFRD, where the first dust grains form to the left and where they are removed/destroyed. In Fig. 7, we confirm the trends observed in Figs. 5 and 6 by plotting $sM_{\text{dust}}$ vs. age$_{\text{main}}$. The decline of $sM_{\text{dust}}$ is interpreted as an efficient dust removal (destruction or outflow) with a timescale of $\sim 500$ Myrs. A word of caution, though: if the dust grains are carried away into the wider ISM, they could cool down or warm up from the injection of kinetic energy, and could therefore be more difficult to detect given our wavelength coverage (Fig. 1). We can try to quantify the removal trends: the distribution of Bayesian values $t/\tau$ found by the CIGALE SED fitting are in the range $0.5 \le \text{age}/\tau_{\text{main}} \le 7.5$. The galaxies with the maximum $sM_{\text{dust}}$ are in the range 0.4–2.9, corresponding to an integrated SFR = 12% and 80%, respectively. The maximum SFR occurs at age/$\tau_{\text{main}}$ = 1.0, where about 26% of the stars have formed. This means (predictably) that $sM_{\text{dust}}$ correlates with SFR and therefore with the dust production by stars and not the ISM.

If we now focus on the Low-zZ sample, our interpretation is that separate parallel $sM_{\text{dust}}$ vs. sSFR sequences could be stratified from bottom to top in addition to from right to left. Accounting for the horizontal age sequence for these low-zZ galaxies from the lower panel of Fig. 5, it could be possible that these local galaxies undergo an increase of $sM_{\text{dust}}$ as they get older. However, this interpretation is not valid for the Hi-z sample, where we do not see low-IRX objects at inter-

<figure><img src="image_8.png" /></figure>

**Fig. 6.** Hi-z LBGs (filled dots, top) and Low-zZ galaxies (filled boxes, middle) dust formation rate diagram, $sM_{\text{dust}}$ vs. sSFR, colour-coded in IRX. Like in Fig. 5, colour-coding the data points helps towards interpreting the sequence as an evolutionary sequence with low-IRX objects at high sSFR and again at low sSFR. In the middle range, IRX is higher. Low-zZ galaxies (filled dots, bottom) dust formation rate diagram, $sM_{\text{dust}}$ vs. sSFR, colour-coded in metallicity.

mediate sSFR, and there are no high-IRX objects on the left-hand side of the plot, at low sSFR, while we still observe the evolutionary sequence from right to left. Under this hypothesis, we can suggest that the Low-zZ and Hi-z LBG samples are not perfect analogues. It is possible that we observe the effect of another parameter, the metallicity, for the Low-zZ galaxies (Fig. 6, lower panel), because these galaxies formed at much later cosmic time than the Hi-z LBGs. We already knew that the metallicity and IRX of Low-zZ galaxies are correlated (e.g. Heckman et al. 1998; Boissier et al. 2004). In Fig. 8, the linear fit we estimate for the Low-zZ sample is close to, but different, from the starburst galaxy relation found by Heckman et al. (1998). We also added two data points corresponding to galaxies at z = 3.3 by Pannella et al. (2015). Using the updated derived relation between IRX and $12+\log_{10}(\text{O/H})$ for the Low-zZ sample ($\log_{10}(\text{IRX}) = 1.856 \times (12 + \log_{10}(\text{O/H}))$, we can estimate the metallicity of the Hi-z LBG sample. The maximum of the normalised distribution of the latter shown in the lower panel of Fig. 8 is offset to lower metallicities by about 0.3 compared to the former. So, the origin of the difference between the two samples in Figs. 5 and 6 might be due to the higher metallicity of the Low-zZ galaxies sample, where dust growth in the ISM could contribute to the dust budget, especially in the
<figure><img src="image_9.png" /><figcaption>**Fig. 7.** *sM<sub>dust</sub> vs. Age<sub>main</sub> diagrams colour-coded in sSFR for Hi-z LBGs. Low-sSFR objects are the oldest objects and low-*sM<sub>dust</sub>* objects are the youngest ones. The decline of *sM<sub>dust</sub>* is interpreted as an efficient dust removal (destruction or outflow) with a time-scale of ~500 Myrs. This diagram makes it possible to obtain an estimate of the dust removal grain in these Hi-z galaxies. Symbols coded as in Fig. 1. From left to right, the vertical dotted lines show 10%, 20%, 30%, 40%, 50%, 60%, 70%, 80%, 90%, and 99% of the integrated SFR, assuming a delayed SFH.</figcaption></figure>

<figure><img src="image_10.png" /><figcaption>**Fig. 9.** *Top:* the evolution of Hi-z LBGs in the IRX vs.β diagram is an evolutionary sequence, with younger galaxies being bluer and older galaxies redder. The continuous line is the Calzetti dust attenuation law, and the dotted line is the SMC one. The trend observed here explains the location of galaxies not detected by Capak et al. (2015). *Bottom:* the same objects are colour-coded with *sM<sub>dust</sub>*. Hi-z LBGs with low dust masses are found below Calzetti and SMC dust attenuation laws.</figcaption></figure>

<figure><img src="image_11.png" /><figcaption>**Fig. 8.** *Top:* we estimate the relation (black line) between the metallicity and IRX of our Low-zZ galaxies (black filled dots) $12+\log_{10}(\text{O/H}) = 0.2829 \times \log_{10}(\text{IRX}) + 7.8155$. The dashed red line is the relation for starburst galaxies from Heckman et al. (1998) and the two blue dots are two objects from Pannella et al. (2015) at $z \sim 3.3$. *Bottom:* from the derived relation between IRX and $12+\log_{10}(\text{O/H})$ for the Low-zZ sample (blue dashed histogramme), we estimate the metallicity of the Hi-z LBG sample (red filled histogramme).</figcaption></figure>

highest metallicity of these Low-zZ galaxies, as is suggested by several papers for the DGS sample (e.g. Asano et al. 2013). Asano et al. (2013) define a metallicity as the critical metallicity $Z_{cr}$, the switching metallicity point at which the increasing rate of dust mass by grain growth exceeds the dust production rate by stars. This $Z_{cr} \sim 0.1-0.5$, which approximately corresponds

to $12+\log_{10}(\text{O/H}) = 8.0-8.6$, in agreement with the Fig. 5 in Ma et al. (2016). In other words, the Low-zZ sample might not be as clean and homogeneous as the Hi-z LBG one, and they are not complete analogues to galaxies in the early universe.

In the upper panel of Fig. 9, we see that Hi-z LBGs with young stellar ages are found to the left of, and above, Meurer's law (the local starburst law linking IRX and β, Meurer et al. 1999) while older and lower *sM<sub>dust</sub>* LBGs are found below this law, in the region identified from ALMA-undetected LBGs by Capak et al. (2015). Therefore, the origin of the position of these Hi-z LBGs is consistent with having both low-IRX values and red UV slopes β.

## 5. Modelling the dust evolutionary sequence

It is beyond the scope of this paper to describe the models in detail. They will be fully explained in an associated paper (Nanni et al. 2020). In a nutshell, we performed the calculations for the metal evolution using the OMEGA code (One-zone Model for the Evolution of Galaxies, Côté et al. 2017). We assumed the metal yields for Type II SNe from Kobayashi et al. (2006) computed up to $40 M_{\odot}$, and from the FRUITY database for low-mass stars with $1.3 < M/M_{\odot} < 6$ evolving through the thermally pulsing AGB phase and developing stellar winds (Cristallo et al. 2011, 2015; Piersanti et al. 2013). The yields for population III stars are from Heger & Woosley (2010) and are limited to the mass range $10 < M/M_{\odot} < 30$ in OMEGA.

Dust removal from the galaxy through galactic outflow follows a rate proportional to the SFR through the “mass-loading factor”: $\text{ML} \times \text{SFR}$. This assumes that the galactic outflow is generated by the feedback of stars on the gas in the ISM (e.g. Murray et al. 2005). The SFH is a delayed one with $\tau$ assumed to be 83 Myrs (in agreement with the average value found from the SED fitting). Two kinds of IMFs are tested: top-heavy IMFs:
$m^{-\alpha}$ with $\alpha = 1.0, 1.35, 1.5$ defined above in this paper and a Chabrier IMF.

We assumed that a fraction of silicates (olivine and pyroxene), iron and carbon grains ejected in the ISM are condensed. In different works (Ventura et al. 2012; Nanni et al. 2013, 2014), it has been shown that for low-mass stars, the condensation fraction ($f_i$) of silicate and carbon dust can be up to 50–60% or slightly more during the super wind phase, when most of the mass is lost by the stars through stellar winds.

In a recent work in which dust dust condensation in SNe remnants is modelled (Marassi et al. 2019), a fair estimate of the mass of dust over the mass of metals is found to be $0.3 < f_i < 0.6$. In this calculation, however, the effect of the reverse shock on dust destruction is not taken into account. There is a tension on the amount of dust produced by SNe in models and what is actually estimated from observations. Bocchio et al. (2016) argued that only a small fraction (a few per cent) of the dust produced initially condensed in SNe remnants would survive the reverse shock, while Gall et al. (2014) found evidence that large dust grains (>1 µm) in SNe remnants might survive the passage of the reverse shock. Matsuura et al. (2019) also found that dust might be reformed after the passage of a shock.

In this work, we define the condensation fraction at each time-step as the ratio between the number of atoms locked into dust grains over the maximum that can condense according to its stoichiometric formula and to the chemical composition of the gas. On the basis of theoretical calculations, we assume, for TP-AGB stars: $f_{\text{sil}} = 0.6$ for silicates, $f_{\text{iron}} = 0.01$ for solid iron, and $f_{\text{car}} = 0.5$ for carbon dust grains. For SNe, we consider two cases: a very high value of the condensation fraction: $f_{\text{sil}} = 1.0$ for silicates, $f_{\text{iron}} = 1.0$ and $f_{\text{car}} = 0.5$ for carbon dust grains. Note that we only considered pyroxene in silicates (MgSiO3). We also add an additional test case with half of the previously assumed condensation fraction and IMF with $\alpha = 1$. These assumptions correspond to two very favourable scenarios for dust condensation in SNe remnants.

For each dust species, the destruction rate by SNe is computed as: $\text{SNe}_{\text{destr}} = M_{\text{dust},i}/\tau_{\text{destr}}$, where $M_{\text{dust},i}$ is the dust mass of each dust species, $\tau_{\text{destr}}$ is the destruction timescale evaluated as: $\tau_{\text{destr}} = M_{\text{gas}}/(\tau_{\text{destr}}R_{\text{SNe}}M_{\text{swept}})$, where $M_{\text{gas}}$ is the gas mass in the galaxy that changes as a function of the outflow and star formation, and $\epsilon_{\text{destr}}$ is the efficiency of the destruction. Here, $\epsilon_{\text{destr}} = 0.1$ (Hirashita & Aoyama 2019) $R_{\text{SNe}}$ is the SNe rate and $M_{\text{swept}}$ is the mass of gas swept by the SNe blast wave. We explore values for $M_{\text{swept}} = 1000$ and $6800 M_{\odot}$. The latter case represents the typical value assumed in the literature (e.g. Dwek et al. 2007), while the former assumption considers that much lower values for $M_{\text{swept}}$ are also possible (Nozawa et al. 2006). The average $\tau_{\text{SFH}} = 83$ Myrs found by fitting individual LBGs is used for the SFH in our modelling. The calculations were performed by normalising the total stellar mass formed in the galaxy after 13 Gyr to $1 M_{\odot}$.

The best chemical evolution models selected by computing the $\chi^2$ with respect to $sM_{\text{dust}}$, sSFR, and age are presented in Fig. 10. The blue models show that we can explain the structure of the DFRD diagram. However, even though the models are consistent with the data if we account for the uncertainties, we cannot reach the top data points. This may or may not mean that we are still missing some important physics about the dust cycle in these early galaxies. The fact that the models are low in the DFRD could be relevant, given the uncertainties in the theoretical metal yields presently available in the literature. Further observational constraints and deeper theoretical studies need to be devoted to the understanding of the dust cycle in the first galaxies.

<figure><img src="image_12.png" /><figcaption>**Fig. 10.** Left: dust formation rate diagram: $sM_{\text{dust}}$ vs. sSFR in the Hi-z LBGs (red data points). Only models with a delayed star formation history are plotted. Consistently, a delayed star formation history only was used in the SED fitting with CIGALE to estimate $sM_{\text{dust}}$ and sSFR. The models shown do not have any grain growth in the ISM. We select the best models (shown in blue) by computing the $\chi^2$ of the models with respect to $sM_{\text{dust}}$, sSFR, and age, and finally use the combined total $\chi^2$ as quality of the fit to select the best models. We also plot the other models in light grey to plot which parameter space the models cover. Right: specific dust mass vs. age of the dominant stellar population. Models with a larger condensation fraction provide better fits, on average. The Hi-z LBGs are young galaxies. The steep decline of $sM_{\text{dust}}$ for Hi-z LBGs with respect to Low-zZ galaxies suggest a fast removal (destruction + outflow) of grains in about 0.5 Gyr for Hi-z LBGs that can explain the various detectability of these Hi-z LBGs. Colours have the same meaning as in the left panel.</figcaption></figure>

We adopted a different initial total baryonic mass of the galaxy between $M_{\text{gas}} = 20 M_{\odot}$ and $100 M_{\odot}$. The larger the mass of the gas, the lower the dust destruction efficiency of SNe. Different values of the mass loading factor (ML) have also been tested, ranging from 0.5 to $0.98 \times M_{\text{gas}}$. For the best models shown in blue in Fig. 10, we computed the mean and standard deviation values: $M_{\text{gas}} = 57.9 \pm 26.7$, $\alpha = -1.14 \pm 0.18$, $ML = 38.3 \pm 18.2$, and $M_{\text{swept}} = 3951 \pm 2900$. In the dust removal budget, outflows are dominant ($\ge 70-80\%$) over SN destruction ($\le 20-30\%$).

## 6. Conclusions

The major result from this work is a consistent model that explains the (sub-)mm detections (or non-detections) of Hi-z LBGs in the EoR. The non-detections are explained by dust destruction from shocks produced by SNe plus removal in the circum- and intergalactic media by outflows (Fig. 10). Such large-scale outflows of interstellar material are a generic feature of LBGs (Shapley et al. 2003; Pettini et al. 2002). The models presented here do not include grain growth in the ISM of these galaxies since the efficiency of such a process can be highly inhibited at the very low temperature characterising the ISM (Ferrara et al. 2016; Ceccarelli et al. 2018). Our modelling suggests that a Chabrier IMF cannot produce enough dust mass to reach the upper part of the sequence. For these LBGs, a top-heavy IMF is preferred.

The possible rise, suggested from data, and decline of $sM_{\text{dust}}$ with sSFRs would be the formation and removal (dust destruction and outflow) of the first dust grains in the universe. If a coeval and massive burst of star formation occurred in the early universe (population III), we could predict an earlier bump at higher redshift in the evolution of the cosmic average dust attenuation (Burgarella et al. 2013). However, this bump is likely to be blurred by a non-coeval galaxy formation.
There are still large uncertainties in our modelling due to a lack of knowledge regarding the physical conditions in the early universe. It is crucial to further constrain the models by securing new observational data, and to invest in work on dust modelling to fully explain the dust cycle in Hi-z LBGs (and in other objects).

*Acknowledgements.* Part of this work was supported by the Centre National d'Etudes Spatiales (CNES) through a post-doctoral fellowship for AN.

## References

Álvarez-Márquez, J., Burgarella, D., Buat, V., Ilbert, O., & Pérez-González, P. G. 2019, A&A, 630, A153

Aravena, M., Decarli, R., Walter, F., et al. 2016, ApJ, 833, 71

Asano, R. S., Takeuchi, T. T., Hirashita, H., & Inoue, A. K. 2013, EPS, 65, 213

Beelen, A., Cox, P., Benford, D. J., et al. 2006, ApJ, 642, 694

Bertoldi, F., Carilli, C. L., Cox, P., et al. 2003, A&A, 406, L55

Bocchio, M., Marassi, S., Schneider, R., et al. 2016, A&A, 587, A157

Boissier, S., Boselli, A., Buat, V., Donas, J., & Milliard, B. 2004, A&A, 424, 465

Bolmer, J., Greiner, J., Krühler, T., et al. 2018, A&A, 609, A62

Boquien, M., Burgarella, D., Roehlly, Y., et al. 2019, A&A, 622, A103

Bouwens, R. J., Illingworth, G. D., Oesch, P. A., et al. 2011, ApJ, 737, 90

Bouwens, R. J., Aravena, M., Decarli, R., et al. 2016, ApJ, 833, 72

Bruzual, G., & Charlot, S. 2003, MNRAS, 344, 1000

Burgarella, D., Buat, V., & Iglesias-Páramo, J. 2005, MNRAS, 360, 1413

Burgarella, D., Le Floc'h, E., Takeuchi, T. T., et al. 2007, MNRAS, 380, 986

Burgarella, D., Buat, V., Gruppioni, C., et al. 2013, A&A, 554, A70

Calzetti, D., Kinney, A. L., & Storchi-Bergmann, T. 1994, ApJ, 429, 582

Capak, P. L., Carilli, C., Jones, G., et al. 2015, Nature, 522, 455

Casey, C. M. 2012, MNRAS, 425, 3094

Castellano, M., Sommariva, V., Fontana, A., et al. 2014, A&A, 566, A19

Ceccarelli, C., Viti, S., Balucani, N., & Taquet, V. 2018, MNRAS, 476, 1371

Chabrier, G. 2003, PASP, 115, 763

Ciesla, L., Boquien, M., Boselli, A., et al. 2014, A&A, 565, A128

Cormier, D., Abel, N. P., Hony, S., et al. 2019, A&A, 626, A23

Côté, B., O'Shea, W., Ritter, C., Herwig, F., & Venn, K. A. 2017, ApJ, 835, 128

Cristallo, S., Piersanti, L., Straniero, O., et al. 2011, ApJS, 197, 17

Cristallo, S., Straniero, O., Piersanti, L., & Gobrecht, D. 2015, ApJS, 219, 40

Cullen, F., McLure, R. J., Khochfar, S., Dunlop, J. S., & Dalla Vecchia, C. 2017, MNRAS, 470, 3006

Dabringhausen, J., Fellhauer, M., & Kroupa, P. 2010, MNRAS, 403, 1054

Dell'Agli, F., Valiante, R., Kamath, D., Ventura, P., & García-Hernández, D. A. 2019, MNRAS, 486, 4738

Draine, B. T. 2009, Space Sci. Rev., 143, 333

Draine, B. T., & Li, A. 2007, ApJ, 657, 810

Draine, B. T., Aniano, G., Krause, O., et al. 2014, ApJ, 780, 172

Dwek, E. 1998, ApJ, 501, 643

Dwek, E., & Cherchneff, I. 2011, ApJ, 727, 63

Dwek, E., Galliano, F., & Jones, A. P. 2007, ApJ, 662, 927

Dwek, E., Staguhn, J. G., Arendt, R. G., et al. 2011, ApJ, 738, 36

Engelbracht, C. W., Kundurthy, P., Gordon, K. D., et al. 2006, ApJ, 642, L127

Faisst, A. L., Capak, P. L., Yan, L., et al. 2017, ApJ, 847, 21

Ferland, G. J., Korista, K. T., Verner, D. A., et al. 1998, PASP, 110, 761

Ferrara, A., Viti, S., & Ceccarelli, C. 2016, MNRAS, 463, L112

Ferrara, A., Hirashita, H., Ouchi, M., & Fujimoto, S. 2017, MNRAS, 471, 5018

Finkelstein, S. L., Papovich, C., Giavalisco, M., et al. 2010, ApJ, 719, 1250

Gall, C., Andersen, A. C., & Hjorth, J. 2011, A&A, 528, A13

Gall, C., Hjorth, J., Watson, D., et al. 2014, Nature, 511, 326

Hashimoto, T., Laporte, N., Mawatari, K., et al. 2018, Nature, 557, 392

Heckman, T. M., Robert, C., Leitherer, C., Garnett, D. R., & van der Rydt, F. 1998, ApJ, 503, 646

Heger, A., & Woosley, S. E. 2010, ApJ, 724, 341

Hirashita, H., & Aoyama, S. 2019, MNRAS, 482, 2555

Hirashita, H., & Omukai, K. 2009, MNRAS, 399, 1795

Hirashita, H., Burgarella, D., & Bouwens, R. J. 2017, MNRAS, 472, 4587

Hjorth, J., Vreeswijk, P. M., Gall, C., & Watson, D. 2013, ApJ, 768, 173

Hou, K.-C., Aoyama, S., Hirashita, H., Nagamine, K., & Shimizu, I. 2019, MNRAS, 485, 1727

Inoue, A. K. 2011, MNRAS, 415, 2920

Jones, A. P., & Nuth, J. A. 2011, A&A, 530, A44

Jones, T., Stark, D. P., & Ellis, R. S. 2018, ApJ, 863, 191

Kobayashi, C., Umeda, H., Nomoto, K., Tominaga, N., & Ohkubo, T. 2006, ApJ, 653, 1145

Komatsu, E., Smith, K. M., Dunkley, J., et al. 2011, ApJS, 192, 18

Leśniewska, A., & Michalowski, M. J. 2019, A&A, 624, L13

Ma, X., Hopkins, P. F., Faucher-Giguère, C.-A., et al. 2016, MNRAS, 456, 2140

Madden, S. C., Rémy-Ruyer, A., Galametz, M., et al. 2013, PASP, 125, 600

Maiolino, R., Nagao, T., Grazian, A., et al. 2008, A&A, 488, 463

Marassi, S., Schneider, R., Limongi, M., et al. 2019, MNRAS, 484, 2587

Matsuura, M., Dwek, E., Barlow, M. J., et al. 2015, ApJ, 800, 50

Matsuura, M., De Buizer, J. M., Arendt, R. G., et al. 2019, MNRAS, 482, 1715

Meurer, G. R., Heckman, T. M., & Calzetti, D. 1999, ApJ, 521, 64

Michalowski, M. J. 2015, A&A, 577, A80

Michalowski, M. J., Watson, D., & Hjorth, J. 2010, ApJ, 712, 942

Mortlock, D. J., Warren, S. J., Venemans, B. P., et al. 2011, Nature, 474, 616

Murray, N., Quataert, E., & Thompson, T. A. 2005, ApJ, 618, 569

Nanni, A., Bressan, A., Marigo, P., & Girardi, L. 2013, MNRAS, 434, 2390

Nanni, A., Bressan, A., Marigo, P., & Girardi, L. 2014, MNRAS, 438, 2328

Nanni, A., Burgarella, D., Theulé, P., Côté, B., & Hirashita, H. 2020, A&A, submitted

Noll, S., Burgarella, D., Giovannoli, E., et al. 2009, A&A, 507, 1793

Nozawa, T., Kozasa, T., Umeda, H., Maeda, K., & Nomoto, K. 2003, ApJ, 598, 785

Nozawa, T., Kozasa, T., & Habe, A. 2006, ApJ, 648, 435

Ohyama, Y., Sakamoto, K., Aalto, S., & Gallagher, J. S. I. 2019, ApJ, 871, 191

Pannella, M., Elbaz, D., Daddi, E., et al. 2015, ApJ, 807, 141

Pavesi, R., Riechers, D. A., Capak, P. L., et al. 2016, ApJ, 832, 151

Pettini, M., Rix, S. A., Steidel, C. C., et al. 2002, Ap&SS, 281, 461

Piersanti, L., Cristallo, S., & Straniero, O. 2013, ApJ, 774, 98

Priddey, R. S., Isaak, K. G., McMahon, R. G., Robson, E. I., & Pearson, C. P. 2003, MNRAS, 344, L74

Rémy-Ruyer, A., Madden, S. C., Galliano, F., et al. 2013, A&A, 557, A95

Rémy-Ruyer, A., Madden, S. C., Galliano, F., et al. 2015, A&A, 582, A121

Riechers, D. A., Bradford, C. M., Clements, D. L., et al. 2013, Nature, 496, 329

Robson, I., Priddey, R. S., Isaak, K. G., & McMahon, R. G. 2004, MNRAS, 351, L29

Rupke, D. 2018, Galaxies, 6, 138

Scoville, N., Sheth, K., Aussel, H., et al. 2016, ApJ, 820, 83

Shapley, A. E., Steidel, C. C., Pettini, M., & Adelberger, K. L. 2003, ApJ, 588, 65

Slavin, J. D., Dwek, E., & Jones, A. P. 2015, ApJ, 803, 7

Stanway, E. R., & Eldridge, J. J. 2018, MNRAS, 479, 75

Stratta, G., Gallerani, S., & Maiolino, R. 2011, A&A, 532, A45

Todini, P., & Ferrara, A. 2001, MNRAS, 325, 726

Valiante, R., Schneider, R., Bianchi, S., & Andersen, A. C. 2009, MNRAS, 397, 1661

Venemans, B. P., McMahon, R. G., Walter, F., et al. 2012, ApJ, 751, L25

Ventura, P., di Criscienzo, M., Schneider, R., et al. 2012, MNRAS, 424, 2345

Watson, D., Christensen, L., Knudsen, K. K., et al. 2015, Nature, 519, 327

Willott, C. J., Carilli, C. L., Wagg, J., & Wang, R. 2015, ApJ, 807, 180

Yoshida, M., Kawabata, K. S., Ohyama, Y., Itoh, R., & Hattori, T. 2019, PASP, 71, 87

Yuan, F.-T., Burgarella, D., Corre, D., et al. 2019, A&A, 631, A123
Appendix A: Individual fits

<figure><img src="image_13.png" /><figcaption>Fig. A.1. Final fits for all Hi-z SEDs. Only fits with reduced $\chi^2 \le 5.0$ are kept.</figcaption></figure>
<figure><img src="image_14.png" /></figure>

Best model for HZ9 at z = 5.548. Reduced χ²=3.0

Best model for 566428 at z = 5.89. Reduced χ²=0.0

<figure><img src="image_15.png" /></figure>

Best model for CLM1 at z = 6.166. Reduced χ²=1.0

<figure><img src="image_16.png" /></figure>

Best model for WMH5 at z = 6.07. Reduced χ²=2.0

<figure><img src="image_17.png" /></figure>

Best model for MACSJD1149at z = 9.11. Reduced χ²=0.0

<figure><img src="image_18.png" /></figure>

<figure><img src="image_19.png" /></figure>

Fig. A.1. continued.
<figure><img src="image_20.png" /><figcaption>Fig. A.1. continued.</figcaption></figure>
Appendix B: Quality of the fits

<figure><img src="image_21.png" /><figcaption>Fig. B.1. Probability Distribution function (PDF) derived from the fit for each of the important parameters for this paper: $L_{\text{FUV}}$, $L_{\text{dust}}$, $M_{\text{dust}}$, $\text{Age}_{\text{main}}$, $M_{\text{star}}$ and SFR for each representative Hi-z LBGs.</figcaption></figure>
<figure><img src="image_22.png" /><figcaption>**Fig. B.2.** Mock analysis derived from the fit for each of the important parameters for this paper: $L_{\text{FUV}}$, $L_{\text{dust}}$, $M_{\text{dust}}$, $\text{Age}_{\text{main}}$, $M_{\text{star}}$, and SFR for each representative Hi-z LBG. To build these plots, we use the best-fit modelled SED for each of the Hi-z LBGs. Then, we add the observational noise to the observed data points and we refit the noised best models. Finally, we compare the input and the output parameters. These plots in this figure shows that CIGALE is able to re-derive the input parameters quite nicely and also provide an estimate of the degeneracy between these parameters.</figcaption></figure>