Open Access
Issue
A&A
Volume 689, September 2024
Article Number A102
Number of page(s) 18
Section Extragalactic astronomy
DOI https://doi.org/10.1051/0004-6361/202450373
Published online 11 September 2024

© The Authors 2024

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

The ultraviolet (UV) to near-infrared (NIR) spectral energy distribution enables the characterisation of stellar age distributions and the assessment of star formation histories (SFHs) of galaxies (e.g. Salim et al. 2005, 2007). The UV continuum is a tracer of young stars, while the infrared (IR) unveils the older stellar populations and dust-obscured regions within galaxies. Dust grains preferentially absorb shorter-wavelength light, leading to dimming effects that are prominent in the UV and optical regions. In contrast, the NIR is much less affected by dust attenuation making it essential to derive unbiased star formation rates (SFRs) and stellar masses (e.g. Calzetti et al. 2000; Leja et al. 2017). Emission lines, such as those arising from ionised hydrogen, provide crucial insights into gas properties and very young (≤10 Myrs) stellar objects (Kennicutt 1998; Boselli et al. 2009). A wide spectral coverage of both continuum and emission line features thus enables robust constraints on SFRs, gas metallicity, and ionisation conditions (Shivaei et al. 2020; Reddy et al. 2023).

Optical and NIR emission lines, driven primarily by stars with lifetimes of 3 to 10 Myrs serve as robust tracers of SFRs (Kennicutt et al. 2009). Notably, recombination lines of hydrogen, such as the Balmer series, are very valuable tools for tracing star formation activities due to their independence on the exact physical conditions of the gas (Osterbrock 1989). However, these diagnostic lines are affected by dust attenuation. A correction can be made by using the Balmer decrement (Hα/Hβ) but this method is hindered by the spatial resolution of the emission and attenuation (Kennicutt & Evans 2012; Keel et al. 2023; Robertson et al. 2024) as the latter may not be uniform over the emission region and high optical depths to Balmer emission.

Paschen lines, which are found in the NIR and are thus suffering less dust attenuation than the optical Balmer lines, can also be used to measure SFRs in star-forming regions obscured by dust clouds (Alonso-Herrero et al. 2006; Calzetti et al. 2007; Kessler et al. 2020; Cleri et al. 2022). When combined with Balmer lines, they can also be used to constrain nebular reddening and also provide insights into the complex interplay between stars, gas, and dust within galaxies (Calzetti et al. 2000; Prescott et al. 2022; Reddy et al. 2023).

For the last two decades, astronomers gathered data from ground-based multi-object NIR spectrographs, such as Keck/MOSFIRE and VLT/KMOS, and space-based facilities like the HST/WFC3 grism, which provided extensive measurements of hydrogen recombination lines and nebular reddening for thousands of galaxies up to a redshift of z ∼ 2.6 (e.g. Kashino et al. 2013; Reddy et al. 2015; Shivaei et al. 2020; Rezaee et al. 2021; Battisti et al. 2022). While using these instruments, only a handful of Paschen lines have been robustly detected in z > 2 lensed galaxies (Papovich et al. 2009; Finkelstein et al. 2011). With the advent of the James Webb Space Telescope (JWST, Gardner et al. 2006, 2023) equipped with its NIRSpec instrument, Paschen lines are now becoming accessible for intermediate redshift sources (Paα is detectable until z ∼ 1.7 and Paβ until z ∼ 3).

The comparison of the combination of stellar, nebular and dust models predicting the full spectral energy distribution (SED), with observed broadband fluxes from both continuum and line emissions has proven to be highly effective in deducing physical parameters within populations of star-forming galaxies (Boselli et al. 2016; Fossati et al. 2018; Buat et al. 2018; Villa-Vélez et al. 2021; Tacchella et al. 2023; Larson et al. 2023; Arrabal Haro et al. 2023a). Even so, the addition of spectroscopic data implies a careful consideration in the choice of models and their impact on the output parameters.

Combining nebular and stellar emissions also requires to account for a different attenuation recipe for both the continuum and the line emission. This differential attenuation can be effectively modelled with the theoretical framework of Charlot & Fall (2000) as it allows for two independent attenuation recipes, one for the interstellar medium and one for the current birth clouds.

Emission line fluxes often suffer from slit losses which needs to be addressed by either using the photometric data on the same spatial aperture or estimating the emission line fluxes beyond the slit through other methods. In this work, we propose to fit simultaneously photometric broadband fluxes and emission line equivalent widths (EWs) to get around the issue posed by slit losses. Boselli et al. (2016) has shown the feasibility and robustness of this spectro-photometric combination. We develop a new method and apply it on a sample of 13 galaxies observed by JWST in the Cosmic Evolution Early Release Science (CEERS) field. These 1 < z< 3 sources are selected to have at least one Balmer line and one Paschen line detected as well as ancillary data from the Hubble Space Telescope (HST) to cover the rest-frame UV range. We assess the feasibility and reliability of combining spectro-photometric data to derive parameters such as the SFR, stellar mass and dust attenuation parameters. We quantify the impact of including the binary stellar populations from the BPASS models (Eldridge et al. 2017; Stanway & Eldridge 2018) into our models on the aforementioned output parameters. We compare the quality of the same fit done on spectro-photometric and photometric data only to quantify the uncertainty on the SFR and stellar mass in the absence of spectroscopic input data. Finally, we measure the impact of using a flexible attenuation recipe inspired by the original Charlot & Fall (2000) (CF00) model on the SED fitting output parameters.

Our paper is organised as follows. The selection of the sample used for this work is presented in section 2. The theoretical framework linking line ratios and nebular dust reddening is outlined in Section 3. Section 4 is dedicated to our SED fitting process. We present our results on the impact of the choice of data and models in Section 5. The results pertaining to attenuation are shown in Section 6. Finally, we discuss the impact of the choice of models of stellar populations, input data and attenuation recipes in Section 7. We adopt a cosmology with H0 = 70 km/s/Mpc, ΩΛ = 0.7, and Ωm = 0.3.

2. Data and sample

The data for this work are taken from the Cosmic Evolution Early Release Science Survey (CEERS; ERS 1345, PI: S. Finkelstein) in the CANDELS survey (Grogin et al. 2011; Koekemoer et al. 2011) of the Extended Groth (EGS, Davis et al. 2007; Noeske et al. 2007) field. The main features of the CEERS program are presented in Arrabal Haro et al. (2023b) and Finkelstein et al. (2023), and will be described in more detail in Finkelstein et al. (in prep.). For each object in this study, we need spectroscopic data to trace the hydrogen line emission of the gas ionised by the very young stars and as much photometric data as we can gather to trace the stellar continuum flux.

2.1. Spectroscopy

The galaxies studied in this work have been observed with JWST/NIRSpec (Jakobsen et al. 2022) multi-object shutter configurations. These were taken with the Micro Shutter Array (MSA; Ferruit et al. 2022) of size 0.2″× 0.46″during the CEERS epoch 2 observations in December 2022.

These NIRSpec observations are split into six different MSA pointings, each of them observed with the G140M/F100LP, G235M/F170LP and G395M/F290LP medium resolution (R = 1000) gratings. The total wavelength coverage of the NIRSpec instrument is from ∼1 to ∼5.2 μm. The MSA is configured to use three shutter slitlets thus enabling a three point nodding pattern, shifting the pointing by a shutter length plus the size of the bar between the shutters in each direction for background subtraction. The total exposure time per grating is 3107s for each MSA pointing.

The spectral reduction of NIRSpec observations is described in Arrabal Haro et al. (2023a) and will be further explained in Arrabal Haro et al. (in prep.). Here, we choose to focus on the main processing steps based on the STScI Calibration Pipeline1 and on the Calibration Reference Data System (CRDS) mapping 1061.

The first step is to correct the raw images for the detector 1/f noise, dark current, bias and “snowball” contaminations (trails produced by cosmic rays). Count-rate maps are then used to create two-dimensional (2D) spectra for each source.

The second step is to subtract the background by using the three point nodding pattern, apply the flat-field correction, and calibrate the photometry and wavelength of these spectra. They are then resampled in case of spectral trace distortions. The multiple background subtracted 2D spectra (one for each nod position) are then combined to create the final 2D spectrum of the source.

A post-combination inspection check is done on the reduced 2D spectra in which hot pixels, detector gaps and reduction artefacts are masked. The noise spectrum is automatically calculated by the JWST pipeline and then corrected as described in Arrabal Haro et al. (2023a) to take into account the spectral resampling. In the pipeline a slit loss correction is performed by default in the pathloss step.

The line fluxes and EWs are measured using the beta release of LIME (Fernández et al. 2023, 2024). The fluxes are measured assuming a Gaussian distribution and taking into consideration the pixel uncertainty produced by the NIRSpec reduction described in Arrabal Haro et al. (2023a) for the CEERS 0.7 data release. The reported flux uncertainty comes from the standard error definition of the least-squares trust-region method. This uncertainty is propagated to the equivalent width calculation. In this computation the linear continuum flux and its standard deviation are derived from two adjacent continua bands (see Fig. 1 for an example of line fitting).

thumbnail Fig. 1.

Example of the LIME fitting process on the Hα + [NII] doublet of ID 6563 (z ∼ 1.7). The top figure shows the Gaussian profile fittings while the lower figure shows the residual between the observation and the theoretical fits. The red shade in the bottom plot shows the calibration pixel error while the yellow shade shows the standard deviation from the continua bands. The gap between the red and yellow regions is intended to avoid potential contamination of uncertain line wings and artefacts.

These continua regions are manually adjusted for every line to guarantee a representative linear continuum measurement and to avoid artefacts. We take special care to confirm that the pixel uncertainty was of the same order than the flux standard deviation from the adjacent continua. The LIME fittings for the 0.7 release spectra used in this work can be found at the interactive repository2 using the corresponding NIRSpec MSA IDs. This site will be updated as the JWST pipeline and flux calibration keep improving. Additionally, this manuscript line measurements can be found on the supplementary online material as a FITS file, where each page contains the measurements from one galaxy. The complete measurements description can be found in the online documentation3.

2.2. Photometry

Our sample galaxies were also observed with JWST/Near Infrared Camera (NIRCam, Rieke et al. 2023) in December 2022. The ten pointings and the reduction steps are described in Bagley et al. (2023). We adopt the multi-band SExtractor catalog described in Finkelstein et al. (2023). This catalog includes all available HST/ACS and WFC3 data (v1.9) and all seven JWST/NIRCam bands (v0.51).

Data from the Canada–France–Hawaii Telescope (CFHT)/Megacam and from the Spitzer/InfraRed Array Camera (IRAC) channels are taken from the multi-wavelength photometric catalog of the EGS field of Stefanon et al. (2017).

Fluxes from the Multiband Imaging Photometer for Spitzer/MIPS, 24μm band, the 450 and 850 μm James Clerk Maxwell Telescope (JCMT)/SCUBA-2 bands are retrieved from the Super-deblended Catalog of the EGS field of Le Bail et al. (2024).

In cases where two data points are overlapping due to their wavelength proximity (CFHT/Megacam r with HST F606W, CFHT/Megacam i and z with HST F814W, NIRCam F115W with HST F125W, NIRCam F150W with HST F140W and F160W, NIRCam F356W with Spitzer/IRAC 1 and NIRCam F444W with Spitzer/IRAC 2), only the flux with the best signal to noise ratio (S/N), (i.e. HST over CFHT, NIRCam over HST and IRAC) is kept in the SED fitting process.

The final photometry catalog includes measurements over the full CEERS NIRCam wavelength range in the F115W, F150W, F200W, F277W, F356W, F410M, and F444W filters, HST/CANDELS ACS_WFC3 F606W and F814W bands, CFHT/Megacam u* and g′, Spitzer/IRAC channels 3 and 4, Spitzer/MIPS 24 μm band, JCMT/SCUBA-2 450 and 850 μm bands. In the catalog, we only keep values with a S/N ≥2 for CFHT, HST, IRAC and MIPS and with a S/N ≥1 for PACS and SCUBA-2.

Table 1.

Available photometric bands.

2.3. Sample selection

Galaxies included in our study are selected based on several criteria. The first criterion is to have data from NIRCam and NIRSpec gratings since the spatial coverage of both instruments does not completely overlap (this gives us 126 sources). As we also combine NIRCam data with photometry obtained by other facilities, we need all of them to be in agreement. The overlaps between HST/NIRCam are checked to only keep galaxies for which fluxes measured in the same aperture of nearby bands (HST/F160W and NIRCam/F150W) were within 1σ (122 objects remain at this stage). We present the photometric data used in the study in Table 1.

Balmer and Paschen lines are necessary as they probe stars located in vastly different environments so we require our objects to have at least one Balmer and one Paschen line with an equivalent width S/N ≥5 (other Balmer and Paschen lines with a S/N ≥3 are also added, see Table 2). Getting the lines to compute the ratio of a Paschen line to a Balmer line (Pa/H hereafter) described in Section 3.1 restricts our sample to the redshift range 1 < z< 3. All CEERS spectra are inspected to identify and secure a spectroscopic redshift between 1 and 3 which amounts to 40 unique sources. The requirement of a S/N ≥5 for at least one Balmer line and one Paschen line further restricts our sample down to 17 objects.

Table 2.

Characteristics of the objects.

Finally, we want to avoid galaxies with an active galactic nucleus (AGN) as their presence would require much more complex models beyond the scope of this paper. To check if our sources could host AGNs, we compute the ratio of the [NII]-6585 over Hα if the [NII]-6585 line is detected (a non-detection of the [NII]-6585 line is indicated by a cross in Table 2). [NII]-6585/Hα is a linear function of the nebular metallicity and it presents a saturation point so a further increase in the [NII]-6585/Hα value is only due to AGN contribution (Stasińska et al. 2006). Following the authors’ prescription, we only keep galaxies with log([NII]-6585/Hα) < −0.2 (four objects were dropped at this stage). Some sources (IDs 8515, 8588 and 8710) are considered “hidden AGNs” by Calabrò et al. (2023) following their diagnostics based on IR emission lines. These AGNs diagnostics are BPT-like diagrams (Baldwin et al. 1981) of ratios of [CI]-9850, [PII] 1.188 μm or [FeII] 1.257 μm to Paschen lines (either Paβ or Paγ) versus a ratio of [SIII]-9530 to Paβ or Paγ. As we do not detect [PII] or [CI] in these three galaxies and two of them have a [FeII] detection with a S/N < 3 σ ([FeII] is not detected for 8588), we decide to still include them as their exact nature as AGNs still needs to be confirmed. Using all of these criteria, we end up with 13 galaxies in our sample as requiring an observation with NIRCam and the NIRSpec gratings drastically reduces the number of available sources.

3. Line ratios as an indicator of differential attenuation

3.1. Theoretical line ratios

Within our sample, we are able to get multiple Balmer and Paschen lines. In this section, we compute theoretical ratios of lines, i.e. ratios expected for a case B recombination and a temperature of 10 000 K (Osterbrock 1989), as well as model the reddening of the stellar emission following a CF00 attenuation recipe. The key feature of this model is the computation of two complementary recipes, one for the birth clouds (BCs) and one for the interstellar medium (ISM) In this way, we are able to properly account for an age-dependent attenuation where not only the total amount of dust attenuation changes as a function of the stellar age but also the way the stellar light is obscured, as originally introduced by CF00.

In this model, young stars up to an age of 10 Myrs are surrounded by BCs which progressively dissipate over time. Dust attenuation works differentially, with an efficient attenuation in the BCs combined with a lower attenuation in the surrounding ISM. The double component attenuation CF00 recipe is critical as it means that the emission lines and UV light emitted by young stars are attenuated by both the BCs and the ISM. However, once stars have migrated outside the BCs, the emitting radiation is subject only to the diffuse ISM attenuation. Throughout this study, we refer to the results or models of dust attenuation as “recipes” and not “curves”. Indeed, the latter would suggest that a unique or universal result exists while the real effective result depends on the SFH and attenuation parameters of each individual galaxy.

We utilise the CF00 model in which two different power-law attenuation recipes are used to compute the total attenuation, Aλ:

A λ BC = A V BC ( λ / 0.55 ) n BC , $$ \begin{aligned} A_\lambda ^{\mathrm{BC} }&= A_{\rm V}^{\mathrm{BC} } (\lambda /0.55)^{n^{\mathrm{BC} }}, \end{aligned} $$(1)

A λ ISM = A V ISM ( λ / 0.55 ) n ISM , $$ \begin{aligned} A_\lambda ^{\mathrm{ISM} }&= A_{\rm V}^{\mathrm{ISM} } (\lambda /0.55)^{n^{\mathrm{ISM} }},\end{aligned} $$(2)

A λ = A λ BC + A λ ISM $$ \begin{aligned} A_\lambda&= A_\lambda ^{\mathrm{BC} } + A_\lambda ^{\mathrm{ISM} } \end{aligned} $$(3)

where nISM is the slope of the attenuation curve for the ISM and nBC the slope of the attenuation curve for the birth clouds. Both power-law exponents are fixed to −0.7 in the original recipe of CF00. However, these parameters have been shown to vary among galaxies (e.g. Buat et al. 2012; Chevallard et al. 2013; Kriek & Conroy 2013; Battisti et al. 2017; Lo Faro et al. 2017; Salim et al. 2018; Trayford et al. 2020; Pantoni et al. 2021; Boquien et al. 2022) so we decide to keep these exponents as free parameters. The attenuation for the birth clouds, A V BC $ A_{\mathrm{V}}^{\mathrm{BC}} $, is computed using the μ parameter:

μ = A V ISM A V ISM + A V BC $$ \begin{aligned} \mu = \frac{A_{\rm V}^{\mathrm{ISM} }}{A_{\rm V}^{\mathrm{ISM} }+A_{\rm V}^{\mathrm{BC} }} \end{aligned} $$(4)

In the absence of dust, the flux ratios of the different recombination lines are fixed, for a given temperature and electron density. In this case, comparing the observed line ratios with the intrinsic ones allows us to infer the reddening coming from the nebular regions.

We use the Balmer decrement (Hα/Hβ) as it is available for all our sources with fluxes measured with a high S/N. For the Pa/H lines ratio, we have many different combinations available thanks to the large wavelength coverage of the three gratings. We check which combination gives us the best chance at discriminating between two different exponents for a power law which follows Eq. 1 of the CF00 recipe by plotting the Pa/H ratios at our disposal in Fig. 2. The two ratios with the widest gap are Paα/Hγ and Paβ/Hδ, they are however, only available for a few sources in our sample (see Table 2). Because of this issue, we settle for the two next best ratios: Paα/Hβ and Paβ/Hγ (in cases where both ratios are available, both will be included in the analysis). Finally we note that the Paβ/Hβ ratio needs very high levels of attenuation to start being discriminant between both exponents, so we are not using this ratio in this section.

thumbnail Fig. 2.

Theoretical line ratios for a case B recombination computed following the BCs component of the Charlot & Fall (2000) attenuation recipe and varying the exponents presented in Eq. 1. nBC = −0.7 corresponds to a shallower curve (solid lines) while nBC = −1.3 corresponds to a steeper curve (dashed lines).

3.2. Analysis of observed line ratios

We want to ascertain if some objects in our sample have hidden attenuation before going further into our analysis. Indeed, a galaxy with a high Pa/H ratio compared to a Balmer decrement close to 2.86 would indicate a hidden and highly attenuated component. To do this, we compare the observed line fluxes ratios to the theoretical ones. However, this comparison comes with a few caveats, the first one being slit losses. Following (Napolitano et al. 2024) we assume that the slit losses affect all lines equally. The second issue comes from the underlying absorption lines of the stellar continuum, especially for Hβ, Hγ and Hδ which have an absorption line EW ∼ 5Å(see Table 3). We correct the fluxes for these absorptions by computing the absorbed flux from the best model of CIGALE (Sect. 3.3).

Table 3.

Absorption line corrections in Angstroms.

We plot Hα/Hβ against Paα/Hβ or Paβ/Hγ, depending on the redshift and availability of the lines along with the theoretical predictions for each model along with the evolution of the attenuation in Fig. 3. For reference, we also add a Milky Way (Cardelli et al. 1989) and a Small Magellanic Cloud (Gordon et al. 2003) attenuation recipe on top of several power laws with varying exponents.

thumbnail Fig. 3.

Observed fluxes ratios corrected for underlying stellar absorption. The thin coloured curves are power laws following the CF00 framework with varying exponents, from top to bottom: −2.0, −1.3, −1.0, −0.7, −0.4.

We find that our objects encounter at most a very moderate hidden attenuation, a result comparable to what Reddy et al. (2023) found for a sample of similar redshift CEERS sources. Taking into account the uncertainties on the flux measurements, our galaxies have either a very low attenuation or are compatible with the attenuation model presented in Sect. 3.1. In the left panel of Fig. 3, four objects are located to the right of the model lines (IDs 8736, 8515, 10293 and 5409). The Hα/Hβ of IDs 8736 and 10293 shows that both objects are consistent with the models and their Pa/H ratio is compatible with this result within less than 1σ (0.8 and 0.5 respectively). ID 5409 has a large uncertainty on its Pa/H ratio which also makes this galaxy compatible with the models within 1 σ. Finally, the Paα/Hβ ratio of ID 8515 is compatible with the models within 2 σ while its other Pa/H ratio, Paβ/Hγ, is compatible with the models and a very low amount of attenuation within 0.1 σ (see the right panel of Fig. 3). From Fig. 3, we conclude that the modified CF00 recipe is able to reproduce the attenuation processes at work in obscured galaxies.

3.3. Absorption line EWs

One common difficulty of measuring hydrogen emission line EWs is the presence of an underlying absorption line. In our fitting method, presented in Sect. 4.3, we reproduce the EW measurements carried out on observed spectra (Sect. 2.1), avoiding any correction for an underlying absorption either on observed or model spectra.

Absorption lines are part of the stellar emission templates used to fit the galaxies and in this section we make use of the best fit model to compute seven hydrogen absorption line EWs for each galaxy of our sample. We use the binary models of BPASS as even the resolution of the high resolution models from Bruzual & Charlot (2003) (BC03) is not enough to accurately model these absorption lines. We show an example of the best model of ID 5430 in Fig. 4 and the values of the EWs are reported in Table 3.

thumbnail Fig. 4.

Snapshots of the stellar continuum of the best model from the fit of spectro-photometric data made with BPASS for ID 5430. All absorption line EWs for this object can be found in Table 3.

Table 4.

CIGALE modules and input parameters.

In this Table, we notice a trend of decreasing EWs with the increase of the wavelength of the line. It is in agreement with previous studies (e.g. Groves et al. 2012; Domínguez et al. 2013; Valentino et al. 2015; Calabrò et al. 2019). We find that the EWs are highly dependant on the galaxy and its specific SFH and care should be taken when correcting with averaged EWs.

Domínguez et al. (2013) and Valentino et al. (2015) did not find any absorption EWs for Paschen lines while Calabrò et al. (2019) finds slightly higher values (no uncertainties are given) than our average measurements for Paγ (2.50 Å compared to 1.10 ± 0.19 Å) and Paβ (2.00 Å compared to 1.30 ± 0.23 Å) and no absorption for Paα (against 0.20 ± 0.08 Å). Part of these differences may be explained by the sample selection as their study is focusing solely on massive dusty starbursts at z < 1.

Our average value for Hα (2.30 ± 0.27 Å) agrees with the value reported by Domínguez et al. (2013) and Calabrò et al. (2019)(2.70 ± 0.52 Å, and 2.50 Å respectively). Valentino et al. (2015) finds a value of 3.50 Å through SED fitting but do not provide an associated uncertainty. EWs of Hγ and Hδ have not been computed in these other studies so we cannot compare them.

The main driver of the strength of absorption line EWs is likely to be the SFH of the galaxy. We find a continuous decrease of the strength of the absorption lines with increasing wavelengths from Paα to Hδ. We explain this result by a larger contribution of old stars to the stellar emission at longer wavelength, reducing the impact of absorption lines from late B and A stars on the stellar continuum.

4. Spectral energy distribution fitting

To fit the SEDs from the combination of broadband fluxes and EWs, we use the modelling software CIGALE. For a complete description of CIGALE and its functionalities see Boquien et al. (2019). The models are built by successively calling modules, each corresponding to a single physical component or process. CIGALE combines a stellar SED built with an SFH, stellar population models and nebular emission with dust absorption and emission components. The energy balance between stellar and nebular dust absorption and dust re-emission is conserved. Energetic photons produced by massive stars ionise the surrounding gas which re–emits the energy in the form of a series of emission lines and a nebular continuum.

The quality of the fit is assessed by the value of the χ2 and the values of the physical parameters and their corresponding uncertainties are estimated as the likelihood-weighted means and standard deviations. Below we briefly present the modules we use and our input parameter values which are summarised in Table 4.

Our SED fitting combines photometric fluxes and EWs as spectroscopic data. When combining integrated photometric data and EWs measured on the NIRSpec spectra we assume that these EWs are representative of the whole galaxies.

Along with the SED fitting process with CIGALE it is possible to build and analyse a mock catalogue related to the observed dataset (Buat et al. 2014; Ciesla et al. 2015). CIGALE uses the best-fit model of each of the objects obtained with the SED-fitting procedure. The flux densities and EWs of the mock SEDs are then computed by randomly picking a flux or EW value from the normal distribution generated using the best model value and the error of the input data as standard deviation. CIGALE is then run on this artificial catalogue with the same configuration as for the first run in order to compare the exact values of the physical parameters corresponding to the artificial SEDs to the parameters estimated by the code with the probability density function of each parameter. Below we use this comparison to check the robustness of the output SED fitting parameters of interest for our study and our ability to constrain them.

4.1. Stellar population models and star formation history

We first define an SFH to compute the stellar spectrum and we choose a delayed exponential SFH in the form of t × exp(−t/τ) where τ is the e-folding time of the main stellar population. We add an on-going burst of constant star formation, defined with the mass fraction created during the burst (fburst from 0 to 20%) and with its age (ageburst from 10 to 100 Myr). We refer to the age of the onset of the delayed star formation as agemain, this parameter ranges from 2000 Myrs to 5000 Myrs and is always lower than the age of the universe at the redshift of the galaxy.

We consider two models of stellar populations: the widely used BC03 models and the BPASS models (version 2.2) which account for binary stars (Eldridge et al. 2017; Stanway & Eldridge 2018). The latter has been shown to be very efficient when fitting UV and optical photometry with spectroscopic data (e.g. Larson et al. 2023; Bolamperti et al. 2023; Rezaee et al. 2023; Tacchella et al. 2023; Lecroq et al. 2024) as they were built with a focus on best reproducing the emission from young and massive stars (e.g. Eldridge et al. 2017; Stanway & Eldridge 2018; Reddy et al. 2022; Marchant & Bodensteiner 2023). The BPASS binary models have known issues when fitting stellar populations older than a few gigayears (Eldridge et al. 2017; Stanway & Eldridge 2018; Han & Han 2019), the first one being that models remain bluer compared to single star models. Two reasons are given by Eldridge et al. (2017) to explain this effect: the fixed binary fraction regardless of the mass of the star and the width of the time bins at late ages combined with a relatively sparse sampling of stellar masses. Binary models also take a longer time to reach the quiescence (Eldridge et al. 2017). This leads to an underestimation of the number of old and small binary stars.

In both cases, we choose the Initial Mass Function of Chabrier (2003) and we assume that the metallicity is fixed to the solar value: Z = 0.02.

4.2. Dust attenuation and emission

To account for the dust attenuation, we apply the model presented in Sect. 3.1 which corresponds to the modified_CF00 module in CIGALE. In this module, A V ISM $ A_{\mathrm{V}}^{\mathrm{ISM}} $, μ, nBC and nISM (Eq. 1, 2, and 4) are free parameters. We sample the parameter space of A V ISM $ A_{\mathrm{V}}^{\mathrm{ISM}} $, μ and nISM very thoroughly (with 0.1 step). nISM ranges from −0.4 to −1.5 to represent the broad distributions of values found from both observations (e.g. Lo Faro et al. 2017; Buat et al. 2019; Pantoni et al. 2021; Seillé et al. 2022) and simulations (e.g. Chevallard et al. 2013; Roebuck et al. 2019; Trayford et al. 2020). The tests we perform on mock catalogues show that nBC is not constrained. We thus introduce only three values for nBC: the widely used value of −1.3 (Charlot & Fall 2000; da Cunha et al. 2008; Chevallard et al. 2013), −0.7 as proposed in the original CF00 recipe (Charlot & Fall 2000) and an in-between value of −1.0.

To model the dust re-emission of the five sources that have far-IR data, we use the models from Draine & Li (2007), updated in Draine et al. (2014) which can be found in the dl2014 module in CIGALE. These models reproduce the diffuse dust emission heated by the general stellar population which means that the dust is illuminated with a single radiation field Umin. The models also include a component to recreate the emission of the dust tightly linked to star-forming regions. In that case, the dust is illuminated with a variable radiation field intensity superior to Umin The dust mass fraction of dust linked to the star–forming regions (respectively diffuse emission) is referred to as γ (respectively 1 – γ). The polycyclic aromatic hydrocarbon (PAH) abundance, qPAH is also a free parameter in the models.

4.3. Nebular emission and EWs

CIGALE models the emission of the ionised gas with CLOUDY 17.01 (Ferland et al. 1998, 2017) with the configuration described in Villa-Vélez et al. (2021) and Theulé et al. (2024). 129 nebular emission lines intensities are calculated and re-scaled with the number of Lyman continuum photons from the stellar emission of the galaxy. We are only interested in hydrogen lines which makes the impact of the input parameters on the emission line intensities negligible. Therefore we decide to fix all parameters pertaining to the nebular emission to the values mentioned below. The radiation field intensity is given by the ionisation parameter:

log U = log ( n γ n H ) $$ \begin{aligned} \log U = \log (\frac{n_{\gamma }}{n_H}) \end{aligned} $$(5)

where nγ is the number density of photons capable of ionising hydrogen and nH the number density of hydrogen. The ionisation parameter is fixed to −2.0, a value typical of extra galactic HII regions (Osterbrock & Ferland 2006). These configurations are also parameterised according to the gas-phase metallicity, Zgas, which is fixed to the value of the stellar metallicity (Z = 0.02). The electron density is set to 100 cm−3, a typical density for the ionised parts of a diffuse nebula (Osterbrock & Ferland 2006). The ionising spectrum is modelled as a constant burst of star formation of 10 Myrs (Boquien et al. 2019; Villa-Vélez et al. 2021) computed with either BC03 or BPASS models.

We model the Paschen and Balmer line EWs4 so as to reproduce the measurements described in Sect. 2.1. For each emission line of each galaxy, the input wavelengths (a six values “array”) defined with LIMe are used to delimit the line region and two nearby featureless continua on each side of the line. The code then performs a linear fit of the overall continuum which will be subtracted from the line. LIMe deblends the Hα and [NII] doublet of our galaxy sample (see Fig. 1) and these lines are also considered separately in CIGALE. Finally as pointed out in Sect 3.2 these EWs are not corrected for underlying absorption as they are computed and thus affected in the same way as the EWs measured with LIMe.

5. Impact of the choice of data and models on the fits

5.1. Fit of the spectro-photometric data of BC03 vs BPASS

In Sect. 4.1, we outlined the widely acknowledged performance of the BPASS models in reproducing the emission of young stars but also some caveats which could impact the results of the SED fitting process. In this subsection, we compare the quality of the fits of the spectro-photometric data made with either BPASS or the more all-purpose BC03 models.

In Fig. 5 we show a representative example of fits performed with BPASS and BC03. The photometric data appear to be better fitted with the BC03 models. We check if we can statistically (from the quality of the fit) select the best model of stellar populations to fit each object from our sample. For this, we use the Bayesian Information Criterium (BIC, adapted from Salmon et al. 2016), the best model corresponding to the lowest BIC. The two fits being performed with the same number of free parameters, the difference between the two BIC, ΔBIC, is equal to the difference between the two best χ2 obtained for the fits with BC03 and BPASS. The χ2 is computed as:

χ 2 = i ( f i α × m i σ i ) 2 + j ( f j m j σ j ) 2 $$ \begin{aligned} \chi ^{2} = \sum _{i} \left( \frac{f_i-\alpha \times m_i}{\sigma _i}\right)^{2} + \sum _{j} \left(\frac{f_j-m_j}{\sigma _j}\right)^{2} \end{aligned} $$(6)

thumbnail Fig. 5.

Comparison of the SEDs of the galaxy ID 5430 fitted with BC03 (left panel) and with BPASS (right panel).

with α a rescaling parameter (see Boquien et al. 2019, for a detailed explanation), fi and mi being the observed and model fluxes, fj and mj being the observed and model EWs and σ the corresponding observational uncertainties.

We adopt the limit of |ΔBIC|> 6 to put a strong preference on a model, a positive preference being defined as 2 < |ΔBIC|< 6 (Salmon et al. 2016). The Δχ2 distribution is shown in Fig. 6. Eleven objects present a strong preference (and one, ID 10293, a positive preference) for the BC03 models while only one galaxy (ID 23542) is found with a strong preference for the BPASS models. For this object the photometric data is better fitted with the BC03 models (see the SEDs in Appendix A) but the EWs are better fitted with the BPASS models. The uncertainties on the measured EWs are very small for this galaxy (∼10 times lower than for the other sources) which leads to a higher impact of spectroscopic data on the calculation of the χ2 value.

thumbnail Fig. 6.

Difference of the χ2 between the BC03 and BPASS stellar population models for all galaxies in our sample. A negative value favors BC03 while a positive one favors BPASS.

We present a more detailed comparison of the quality of the fit of the photometric data with BPASS and BC03 models in Fig. 7. For each filter and both models we compute the ratio between the bayesian estimate of the flux given by CIGALE over the measured flux. While the CFHT and HST fluxes do not show a significant difference of quality between both models, there is a clear trend for the four bluest NIRCam filters to underestimate the model fluxes with respect to the observed values when the binary BPASS models are used. The CFHT and HST filters cover the rest-frame UV to optical spectrum while the JWST filters cover the rest-frame optical to NIR range. As pointed out in Sect. 4.1, the BPASS models are not optimised to reproduce older stars which we suspect to cause the underestimation of the model fluxes in these specific filters, a result supported by Osborne & Salim (2024). According to this test on our full sample we conclude that the BC03 models better reproduce the photometric measurements.

thumbnail Fig. 7.

Quality assessment of the photometry fitting for the fits of the spectro-photometric data with BC03 (red) and BPASS (blue). For each model, the median of each filter is represented as a square.

We also compare the quality of the fit of the spectroscopic data for both BC03 and BPASS models by comparing the EWs measured on the observed spectra to their bayesian estimation with the CIGALE fit. As can be seen in Fig. 8 both models are equally good at reproducing the observed EW larger than ∼10 Åwhile fainter emission lines (EW < 10 Å) tend to be slightly underestimated.

thumbnail Fig. 8.

Quality assessment of the EWs fitting for the fits of the spectro-photometric data. Filled dots are used for BC03 fits and empty dots for BPASS fits. The linear fit of the Hα points is plotted as a solid magenta line for BC03, dashed magenta for BPASS and the black solid line is the 1:1 line.

To conclude both stellar models are able to fit well the EWs but photometric data are better fitted with BC03 models.

5.2. Impact of fitting only photometric data and line emission contributions

When a substantial amount of nebular emission contaminates the broadband data (Anders & Fritze-Alvensleben 2003; Zackrisson et al. 2008; Pacifici et al. 2015), it can lead to large uncertainties on the results of the SED fitting as the code lacks the information on the relative contribution of the line to the broadband emission. Notorious examples of erroneous interpretations due to intense nebular emission are broadband measurements made with Spitzer/IRAC (Fazio et al. 2004) interpreted as a Balmer break (e.g. Chary et al. 2005; Zackrisson et al. 2008; van der Wel et al. 2011) and more recently measurements with JWST which mimic a Lyman break leading to incorrect estimates of redshifts (e.g. Donnan et al. 2023; Arrabal Haro et al. 2023a; Zavala et al. 2023; Davis et al. 2023). Below we focus on the consequences of fitting only photometric data which accounts for the vast majority of studies making use of SED fitting. We also measure the contribution of line emission in some of the NIRCam broadband filters obtained with the fits of both photometric fluxes and EWs. For these studies we fit our 13 galaxies with both BC03 and BPASS models despite the results of the BIC test. We choose to use both as more and more studies include the BPASS models.

Table 5.

Hα contribution to the NIRCam filter.

5.2.1. Fitting only photometric data

We start by looking at the quality of the fit of photometric data only. We obtain the same trend as in Fig. 7 i.e. an overestimation of the fluxes in the four bluest NIRCam filters with the BPASS models. Next we compare the estimations of the EWs of hydrogen recombination lines with the measured EWs and we present the results in Fig. 9. As expected these estimations are less accurate than the ones obtained when fitting simultaneously photometric fluxes and EWs (Fig. 8). Both stellar population models tend to overestimate the EWs of the lines as CIGALE is favouring stronger bursts of star formation compared to the ones obtained when fitting the spectro-photometric data: the stellar mass fraction created during the burst increases from an average value of 0.05 to 0.13. This overestimation is much more prevalent with the BPASS models. This larger contribution of the lines we find for the fits with BPASS (average of ∼45%) compared to those made with BC03 (average of ∼15%) is explained by the addition of a more intense and younger burst of star formation with median values of fburst and ageburst equal to 0.2 and 50 Myrs for BPASS against 0.1 and 100 Myrs for BC03.

thumbnail Fig. 9.

Quality assessment of the EW estimations when fitting photometric data only. Filled dots are used for BC03 fits and empty dots for BPASS fits. The linear fit of the Hα points is plotted as a solid magenta line for BC03, dashed magenta for BPASS and the black solid line is the 1:1 line.

5.2.2. Contribution of line emission in broad-band filters

To estimate the contribution of the emissions lines to the flux measured in each of our broadband filters, we go back to the fit of photometric data and EWs and we consider the best model obtained with the BC03 and BPASS models. We compute the ratio of the flux coming from the nebular emission and the flux of the continuum within the broadband filter. At the redshift of our galaxies, strong emission lines (Hα, Hβ, Hγ, [OIII]-4959, [OIII]-5007, [NII]-6548 and [NII]-6583) are lying within the bandpasses of the NIRCam F150W and F200W filters. Out of all these lines, Hα is the most intense which makes its contribution easier to track. In Table 5 we provide the contribution of the nebular emission5 in the NIRCam broadband filter containing the Hα line for the BC03 and BPASS models. The filter to which Hα is contributing changes between F150W and F200W depending on the redshift of each source. We plot the percentage of nebular emission contribution in the filter containing Hα in Fig. 10 against the specific SFRs (sSFRs) of the galaxies. The sSFRs are obtained by dividing the bayesian estimation of the SFRs averaged over the last 10 Myrs with the bayesian estimation of the stellar mass both obtained from the output parameters of the best model from the SED fitting. Four galaxies (IDs 8515, 8736, 16991 and 18294) are not included since the Hα line is redshifted in-between NIRCam filters.

thumbnail Fig. 10.

Correlation between the sSFR and the level of contamination in the NIRCam filter containing the Hα line. The error on the y-axis comes from the uncertainty on the attenuation. The dashed lines represent the best fit of each sample of data points.

The contribution of the Hα line is found to increase with the sSFR for both models. We find a stronger average contribution from the nebular emission with BPASS (42%) than with BC03 (36%) as well as a larger average sSFR of 7.4 × 10−9yr−1 for BPASS and 5.4 × 10−9yr−1 for BC03. McKinney et al. (2023) find an average contribution of the Hα+[NII] doublet of 60% for their sample composed of hot, dust-obscured galaxies and dusty, star-forming galaxies at 1 < z < 4 and include the [NII] doublet which may explain the higher percentage they found. We observe a large scatter between the values obtained with BC03 and BPASS.

This dispersion reflects the inherent differences between both models already discussed in Sect. 4.1 such as the more intense emission lines generated by BPASS or the “missing” low mass stars. We further discuss the impact of the contribution of nebular emission on the stellar mass and SFR with both BPASS and BC03 models in Sect. 7.

6. Dust attenuation

In this section, we study the output dust attenuation parameters and for a few cases where the amount of attenuation is high enough, we also reconstruct the total effective dust attenuation curve of each object. We only consider the run performed with the BC03 models since the fit of photometric fluxes is crucial to measure the dust attenuation of the stellar continuum.

6.1. Dust attenuation parameters

We start by discussing the output parameters computed by CIGALE: AV, nISM and μ defined in the CF00 recipe (Section 3.1). The two quantities nISM and μ cannot be accurately estimated if the amount of dust in the ISM is too low. Following Corre et al. (2018) we select a subsample of fitted sources with a required minimum AV> 0.3 mag. Six galaxies fulfill this condition: ID 5409, ID 5430, ID 6563, ID 8588, ID 8710 and ID 23542. Only two objects have AV> 1 mag, the four other ones have AV < 0.6 mag.

Table 6.

Attenuation parameters of our subsample.

In Table 6, we present the results of the SED fitting for the attenuation parameters of our subsample. The average values of μ and nISM (nISM = −0.73 ± 0.15 and μ = 0.37 ± 0.08) are in good agreement with the original values of CF00 (μ = 0.3 and nISM = −0.7). However, similarly to Battisti et al. (2017, 2020), we find that μ can vary a lot (from 0.1 to 0.7) between objects and fixing it to the standard CF00 value of 0.3 can change the shape of the SED, the SFRs and stellar masses (see Sect. 7.3).

Several studies have shown a link between the redshift and a higher colour excess (as defined with the attenuation recipe of Calzetti et al. 2000) needed for emission lines than for the stellar continuum (e.g. Yoshikawa et al. 2010; Whitaker et al. 2014). For our subsample of galaxies, we conclude that there is no major trend of μ increasing with redshift between 1 and 3.

6.2. Derivation of the effective attenuation curves

In the CF00 recipe, the interplay between A v ISM $ A_v^{\mathrm{ISM}} $ and A v BC $ A_v^{\mathrm{BC}} $ (through μ) depends on the star formation history (Charlot & Fall 2000). As a consequence, no universal attenuation curve can be defined.

With CIGALE, we can derive the effective attenuation curve Aλ/AV for each galaxy of our subsample, AV being the total attenuation in the V band. We select fifteen wavelength bands (from the rest-frame far-UV to the near-IR) as the code is able to compute the attenuation in filters different from the ones used to fit the input data. With this feature in mind we add fifteen custom-made filters to the code. We choose to create very narrow filters (10 nm each) to carefully select which part of the spectrum is going to be sampled. The reason for this choice is two-fold: we can avoid contamination from emission and absorption lines but also pinpoint a specific part of the spectrum. We use the bayesian estimations of CIGALE for fifteen photometric bands to reconstruct the total effective attenuation curve of our subsample of galaxies.

We show the attenuation curves of our subsample. All curves except that of ID 8588 agree well on the UV to optical part of the spectrum but all six of them are much more spread out once we reach the near-IR.

We fit each attenuation curve with a single powerlaw:

A λ = A V ( λ / 0.55 ) m $$ \begin{aligned} A_\lambda = A_{\rm V} ({\lambda }/0.55)^{m} \end{aligned} $$(7)

and we report the resulting exponents in Table 6. We also fit our whole subsample with a single powerlaw and find an exponent of m = −0.96± 0.10. By colour-coding each curve with its associated value of μ we see a trend between the decrease in μ and the steepening of the curve (i.e. the exponent m is closer to −1.3). A high μ equals a negligible increase of attenuation for stars in birth clouds compared to the ISM meaning a single powerlaw (for ID 8588). The differential attenuation found for the galaxies with μ < 0.35 mag leads to a steeper attenuation law. Using radiative transfer modelling, it has been shown that the effective attenuation curve depends on the amount of dust attenuation (e.g. Pierini et al. 2004; Saftly et al. 2015; Feldmann et al. 2017; Roebuck et al. 2019; Trayford et al. 2020). As predicted by the models the IR part of the curves do not follow the model of Calzetti et al. (2000) and ID 8588 is the only galaxy to follow the model in the UV-optical which shows the need for a two-components dust attenuation model. Finally we compared the values of μ obtained through a fit made with μ, A v ISM $ A_v^{\mathrm{ISM}} $, nISM and nBC being free and a fit made with μ and A v ISM $ A_v^{\mathrm{ISM}} $ being free and both exponents fixed to −0.96. Both sets of values are consistent within their respective errorbars.

7. Impact of input data type and models on SFRs and stellar mass determinations

The precision and confidence in constraining the physical parameters for our sample of galaxies through statistical fits of observations rely on two key factors: the nature of the observations (whether photometry alone or in conjunction with EWs) and the spectral libraries used to model them. Changing one or both of these factors can lead to major modifications in the SFH which may result in different output parameters and especially SFRs and stellar masses.

As getting robust SFRs and stellar masses is crucial in numerous contexts (e.g. Pacucci et al. 2023; Juodžbalis et al. 2024), we explore the effects of the choice of stellar populations to fit the spectro-photometric data on both aforementioned output parameters. Then we quantify the variation of SFRs and stellar masses between a fit performed with broad-band photometry alone and with spectro-photometric data, the latter being used as our reference fit. Finally we assess the impacts of the change of the attenuation recipe on the SFRs and stellar masses when fitting spectro-photometric data with the BC03 models.

7.1. Fits with BPASS or BC03 models

Here we compare SFRs and stellar mass estimations obtained with either BPASS or BC03 models. We start with SFRs estimations. Fig. 12 shows a clear trend of increased SFRs when using BPASS compared to BC03 models. Excluding the three outliers (IDs 10293, 8736 and 5409 in red in Fig. 12) the SFRs derived through SED fitting with BPASS are found on average 0.13 dex larger than with BC03. The value of the current SFR is expected to be sensitive to the burst fraction (fburst) added to the delayed SFH. The average value of fburst with BPASS is found to be almost double the average fburst value with BC03 (0.15 and 0.08 respectively) which may explain the higher SFRs we get while fitting our sample with the BPASS models. For the three outliers the difference is even higher as one of the two fburst values is found to be null. The two galaxies with a much higher SFR obtained with BPASS are fitted with no burst with BC03 and with a fburst value equal to 0.20 and 0.12 for ID 10293 and ID 8736 respectively with BPASS. Conversely for the galaxy with a much higher SFR with BC03 (ID 5409): no burst is found with BPASS while a small but positive burst fraction (0.01) is obtained with BC03. All in all the difference in SFRs obtained with both stellar population models seems to be correlated with the burstiness of the model rather than with an intrinsic difference between single and binary stars. Indeed Eldridge et al. (2017) and Wilkins et al. (2019) found that the SFRs obtained with BPASS are ∼ 0.1 dex only lower than those obtained with BC03 when the SFH is fixed.

thumbnail Fig. 11.

Total attenuation curves of the subsample with a fit of the mean curve in a black solid line. On the rightmost measurement from top to bottom: IDs 8588, 8710, 5409, 23542, 6563 and 5430. The classical CF00 recipe as well as the Calzetti law have been added for comparison. Each total attenuation curve is color-coded with its value of μ presented in Table 6.

thumbnail Fig. 12.

Comparison of the SFR (left panel) and stellar mass (right panel) obtained through fitting spectro-photometric data in a run made with BC03 and a run made with BPASS. The solid line represents the 1:1 relation while the dashed line represents the mean deviation of our sample (without the three outliers in the case of the SFR).

We find higher stellar masses with BC03 models as shown in Fig. 12. The mass derived through SED fitting with the BPASS models is on average 0.18 dex lower than with the BC03 models.Osborne & Salim (2024) also report lower masses (by 0.17 dex) when using the BPASS models compared to the BC03 models for a sample of CANDELS sources with redshift between 0.7 and 2.3. We attribute these lower masses to the deficit in low mass stars of the BPASS models addressed in Sect 4.2.

7.2. Fits of photometric data or spectro-photometric data

Here we assess the impact of using only photometry with both stellar population models on the stellar mass and SFR estimations. In Fig. 13 we compare the estimates of the SFR obtained with a fit of photometric data to those from a fit of spectro-photometric data with both BC03 and BPASS. The SFR derived with the fit of photometric data only is on average 0.16 dex higher with BC03 and 0.03 dex higher with BPASS compared to values obtained with spectro-photometric measurements. Considering the scatter found in this comparison we also provide the maximal and minimal absolute value of the difference in log(SFR). For BC03 (BPASS) the largest discrepancy is of 0.42 (0.60) dex and the lowest is of 0.04 (0.09) dex. Errors on the SFRs also become much lower when fitting the spectro-photometric data. SFRs lower than 30 M/yr tend to be systematically overestimated when only photometric fluxes are fitted by 0.18 dex with BC03 and by 0.09 dex with BPASS. We attribute this increase to the stronger emission lines computed by the code in the absence of spectroscopic data as shown in Fig. 9 which are due to a more active star formation process. Our results are in agreement with Pacifici et al. (2015) who also find a higher SFR (by 0.13 dex) while using BC03 and photometric data only for a sample of 364 3D-HST galaxies at z < 3.

thumbnail Fig. 13.

Comparison of the SFR (left panel) and stellar mass (right panel) of BC03 and BPASS when fitting photometry only (y-axis) and spectro-photometry (x-axis).

The stellar masses are found larger when EWs are added to the fit by ∼0.15 dex regardless of the chosen stellar populations as shown in Fig. 13. The lower masses obtained with only photometric data are due to an overestimation of the line contribution (see Fig. 9): the contribution of young stars to the SED is higher leading to a lower contribution of older stars which make the bulk of the galaxy mass. Our results are in agreement with the findings of Schaerer & de Barros (2010) who fits z ∼ 7 sources while Pacifici et al. (2015) finds no noticeable difference in the stellar mass output with the addition of spectroscopic data to the fits for their sample of z < 3 galaxies.

These findings highlight the importance of considering both spectroscopic and photometric data in SED fitting to obtain accurate estimates of star formation rates and stellar masses, particularly for galaxies with lower star formation rates. They also show the impact of emission line contributions on parameter estimation emphasising that care needs to be taken when considering the results for sources with strong emission lines.

7.3. Fits with different attenuation recipes

Dust attenuation has a direct impact on the derived SFRs and to a lesser extent on stellar masses. Any change in the shape of the attenuation curve can affect (albeit in varying degrees) the estimations of these quantities as shown in (e.g. Schaerer & de Barros 2010; Pacifici et al. 2015; Lo Faro et al. 2017; Buat et al. 2018, 2019). In the following analysis we compare the effects of changing the attenuation curve parameters in the framework of the CF00 recipe.

We test three different attenuation curves with separate runs using the BC03 models only and the full spctro-photometric dataset. We do not include the BPASS models in these tests as the attenuation curves presented in Sect. 6.2 cannot be reconstructed with the binary stellar population models. The first run called Runref is performed with the four parameters of the modified CF00 recipe taken free as presented in Table 4; it will be considered the reference run. The second run is performed with the original CF00 recipe (RunCF00) in which μ is fixed to 0.3 and both exponents are equal to −0.7 while the third one makes use of the single average power-law introduced in Sect. 6.2; μ is taken free and we call this run Runfix.

As already noted in Sect. 6.2 we do not find noticeable changes in μ between Runref and Runfix. The SFRs are found consistent with a slightly larger scatter with CF00 than with the single power-law as shown in Fig. 14. The agreement between the three stellar mass estimations is very good with a scatter of 0.03 and 0.06 dex for Runfix and RunCF00 respectively when they are compared to Runref (see Fig. 14).

thumbnail Fig. 14.

Comparison of the SFR (left panel) and stellar mass (right panel) obtained with a fit in which both exponents are fixed to −0.96 (dots) and a fit following the original CF00 recipe (squares) with respect to the one we get while using the parameters presented in Sect 4.

From the results shown in this section we conclude that fixing the slope of the attenuation curve (and even μ) is much less impactful on the SFRs and stellar masses than including binary stellar populations or fitting only photometric data. However this lesser effect could be a result of our sample being mostly composed of objects with low dust content. Indeed the most discrepant objects in Fig. 14 have an amount of dust Av > 1 mag.

8. Conclusions

In this work we assessed the robustness of the estimates of the SFRs, stellar masses and dust attenuation parameters of 13 sources in the CEERS field observed with HST, JWST/NIRCam and NIRSpec between redshifts 1 and 3. These objects are selected based on the availability of at least one Balmer and one Paschen line with S/N ≥ 5 to allow for a solid determination of the amount of differential attenuation at play in our sample. The aforementioned physical parameters are computed using a newly developed version of CIGALE which combines photometric fluxes and emission line EWs as input data. We summarise our main results as follows:

  • The absorption line EWs of seven hydrogen lines from Paα (up to 2 Å) to Hδ (up to 7 Å) are measured. This is the first time an absorption EW for the Paα line is measured and a steady trend of decreasing EWs with longer wavelengths is found.

  • The photometric data is better fitted with BC03 compared to BPASS while the quality is similar for EWs. Meanwhile the SED fitting code overestimates the EWs when fitting photometry only (by an average of ∼15%) an effect we find to be even more pronounced when using the BPASS models (average of ∼45%).

  • We find higher SFRs (0.13 dex) and lower stellar masses (0.17 dex) with the fit of spectro-photometric data using BPASS compared to BC03. There is an overestimation of the SFRs (0.16 dex with BC03 and 0.03 with BPASS) and an underestimation of stellar masses (0.15 dex) when fitting only photometry compared to fitting spectro-photometric data.

  • From the best model computed by CIGALE for each source we assess the contribution of the Hα line emission to the measurements made with the NIRCam filters (F150W or F200W depending on the redshift). We find Hα to be a major contributor to the broadband photometric fluxes (average of 36% and up to 68%) even more so while using the BPASS models (average of 42% and up to 76%).

  • We find that using a double power law as introduced by the CF00 recipe to model the dust attenuation allows for solid constraints on the differential attenuation. The attenuation parameters derived for a subsample of sources with Av > 0.3 mag are on average consistent with the CF00 recipe, especially the differential attenuation (μ = 0.37 ± 0.08). The change of flexibility in the dust attenuation recipe has no effect on SFRs and stellar masses but this could be different for more attenuated objects. We reconstruct the effective attenuation curves of our subsample and derive a unique power law having an exponent of −0.96.

Data availability

FITS files are available at the CDS via anonymous ftp to cdsarc.cds.unistra.fr (130.79.128.5) or via https://cdsarc.cds.unistra.fr/viz-bin/cat/J/A+A/689/A102


4

The EWs computation used in this work is not included in the current released version of CIGALE, 2022.1. The newest version of the EWs computation is available on request.

5

Other lines than these strong emission lines, as well as the nebular continuum, are included by CIGALE in the nebular emission, but their contribution is found negligible.

Acknowledgments

Lise-Marie Seillé gratefully acknowledges support from the PNCG; the Programme Investissements d’Avenir; and the Initiative d’excellence A*Midex (doctoral grant AMX-21-HAN-01).

References

  1. Alonso-Herrero, A., Rieke, G. H., Rieke, M. J., et al. 2006, ApJ, 650, 835 [Google Scholar]
  2. Anders, P., & Fritze-Alvensleben, U. 2003, A&A, 401, 1063 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  3. Arrabal Haro, P., Dickinson, M., Finkelstein, S. L., et al. 2023a, ApJ, 951, L22 [NASA ADS] [CrossRef] [Google Scholar]
  4. Arrabal Haro, P., Dickinson, M., Kartaltepe, J., et al. 2023b, Am. Astron. Soc. Meeting Abstr., 55, 149.06 [Google Scholar]
  5. Bagley, M. B., Finkelstein, S. L., Koekemoer, A. M., et al. 2023, ApJ, 946, L12 [NASA ADS] [CrossRef] [Google Scholar]
  6. Baldwin, J. A., Phillips, M. M., & Terlevich, R. 1981, PASP, 93, 5 [Google Scholar]
  7. Battisti, A. J., Calzetti, D., & Chary, R. R. 2017, ApJ, 840, 109 [NASA ADS] [CrossRef] [Google Scholar]
  8. Battisti, A. J., Cunha, E. D., Shivaei, I., & Calzetti, D. 2020, ApJ, 888, 108 [NASA ADS] [CrossRef] [Google Scholar]
  9. Battisti, A. J., Bagley, M. B., Baronchelli, I., et al. 2022, MNRAS, 513, 4431 [NASA ADS] [CrossRef] [Google Scholar]
  10. Bolamperti, A., Zanella, A., Meštrić, U., et al. 2023, MNRAS, 526, 5263 [NASA ADS] [CrossRef] [Google Scholar]
  11. Boquien, M., Burgarella, D., Roehlly, Y., et al. 2019, A&A, 622, A103 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  12. Boquien, M., Buat, V., Burgarella, D., et al. 2022, A&A, 663, A50 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  13. Boselli, A., Boissier, S., Cortese, L., et al. 2009, ApJ, 706, 1527 [Google Scholar]
  14. Boselli, A., Roehlly, Y., Fossati, M., et al. 2016, A&A, 596, A11 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  15. Bruzual, G., & Charlot, S. 2003, MNRAS, 344, 1000 [NASA ADS] [CrossRef] [Google Scholar]
  16. Buat, V., Noll, S., Burgarella, D., et al. 2012, A&A, 545, A141 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  17. Buat, V., Heinis, S., Boquien, M., et al. 2014, A&A, 561, A39 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  18. Buat, V., Boquien, M., Małek, K., et al. 2018, A&A, 619, A135 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  19. Buat, V., Ciesla, L., Boquien, M., Małek, K., & Burgarella, D. 2019, A&A, 632, A79 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  20. Calabrò, A., Daddi, E., Puglisi, A., et al. 2019, A&A, 623, A64 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  21. Calabrò, A., Pentericci, L., Feltre, A., et al. 2023, A&A, 679, A80 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  22. Calzetti, D., Armus, L., Bohlin, R. C., et al. 2000, ApJ, 533, 682 [NASA ADS] [CrossRef] [Google Scholar]
  23. Calzetti, D., Kennicutt, R. C., Engelbracht, C. W., et al. 2007, ApJ, 666, 870 [Google Scholar]
  24. Cardelli, J. A., Clayton, G. C., & Mathis, J. S. 1989, ApJ, 345, 245 [Google Scholar]
  25. Chabrier, G. 2003, PASP, 115, 763 [Google Scholar]
  26. Charlot, S., & Fall, S. M. 2000, ApJ, 539, 718 [Google Scholar]
  27. Chary, R.-R., Stern, D., & Eisenhardt, P. 2005, ApJ, 635, L5 [CrossRef] [Google Scholar]
  28. Chevallard, J., Charlot, S., Wandelt, B., & Wild, V. 2013, MNRAS, 432, 2061 [CrossRef] [Google Scholar]
  29. Ciesla, L., Charmandaris, V., Georgakakis, A., et al. 2015, A&A, 576, A10 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  30. Cleri, N. J., Trump, J. R., Backhaus, B. E., et al. 2022, ApJ, 929, 3 [NASA ADS] [CrossRef] [Google Scholar]
  31. Corre, D., Buat, V., Basa, S., et al. 2018, A&A, 617, A141 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  32. da Cunha, E., Charlot, S., & Elbaz, D. 2008, MNRAS, 388, 1595 [Google Scholar]
  33. Davis, M., Guhathakurta, P., Konidaris, N. P., et al. 2007, ApJ, 660, L1 [NASA ADS] [CrossRef] [Google Scholar]
  34. Davis, K., Trump, J. R., Simons, R. C., et al. 2023, ApJ, submitted [arXiv:2312.07799] [Google Scholar]
  35. Domínguez, A., Siana, B., Henry, A. L., et al. 2013, ApJ, 763, 145 [CrossRef] [Google Scholar]
  36. Donnan, C. T., McLeod, D. J., Dunlop, J. S., et al. 2023, MNRAS, 518, 6011 [Google Scholar]
  37. Draine, B. T., & Li, A. 2007, ApJ, 657, 810 [CrossRef] [Google Scholar]
  38. Draine, B. T., Aniano, G., Krause, O., et al. 2014, ApJ, 780, 172 [Google Scholar]
  39. Eldridge, J. J., Stanway, E. R., Xiao, L., et al. 2017, PASA, 34, e058 [Google Scholar]
  40. Fazio, G. G., Hora, J. L., Allen, L. E., et al. 2004, ApJS, 154, 10 [Google Scholar]
  41. Feldmann, R., Quataert, E., Hopkins, P. F., Faucher-Giguère, C.-A., & Kereš, D. 2017, MNRAS, 470, 1050 [NASA ADS] [CrossRef] [Google Scholar]
  42. Ferland, G. J., Korista, K. T., Verner, D. A., et al. 1998, PASP, 110, 761 [Google Scholar]
  43. Ferland, G. J., Chatzikos, M., Guzmán, F., et al. 2017, Rev. Mex. Astron. Astrofis., 53, 385 [NASA ADS] [Google Scholar]
  44. Fernández, V., Amorín, R., Sanchez-Janssen, R., del Valle-Espinosa, M. G., & Papaderos, P. 2023, MNRAS, 520, 3576 [CrossRef] [Google Scholar]
  45. Fernández, V., Amorín, R., Firpo, V., & Morisset, C. 2024, A&A, 688, A69 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  46. Ferruit, P., Jakobsen, P., Giardino, G., et al. 2022, A&A, 661, A81 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  47. Finkelstein, K. D., Papovich, C., Finkelstein, S. L., et al. 2011, ApJ, 742, 108 [CrossRef] [Google Scholar]
  48. Finkelstein, S. L., Bagley, M. B., Ferguson, H. C., Wilkins, S. M., & Kartaltepe, J. S. 2023, ApJ, 946, L13 [NASA ADS] [CrossRef] [Google Scholar]
  49. Fossati, M., Mendel, J. T., Boselli, A., et al. 2018, A&A, 614, A57 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  50. Gardner, J. P., Mather, J. C., Clampin, M., et al. 2006, Space Sci. Rev., 123, 485 [Google Scholar]
  51. Gardner, J. P., Mather, J. C., Abbott, R., et al. 2023, PASP, 135, 068001 [NASA ADS] [CrossRef] [Google Scholar]
  52. Gordon, K. D., Clayton, G. C., Misselt, K. A., Landolt, A. U., & Wolff, M. J. 2003, ApJ, 594, 279 [NASA ADS] [CrossRef] [Google Scholar]
  53. Grogin, N. A., Kocevski, D. D., Faber, S. M., et al. 2011, ApJS, 197, 35 [NASA ADS] [CrossRef] [Google Scholar]
  54. Groves, B., Brinchmann, J., & Walcher, C. J. 2012, MNRAS, 419, 1402 [NASA ADS] [CrossRef] [Google Scholar]
  55. Han, Y., & Han, Z. 2019, ApJS, 240, 3 [NASA ADS] [CrossRef] [Google Scholar]
  56. Jakobsen, P., Ferruit, P., Alves de Oliveira, C., et al. 2022, A&A, 661, A80 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  57. Juodžbalis, I., Maiolino, R., Baker, W. M., et al. 2024, ArXiv e-prints [arXiv:2403.03872] [Google Scholar]
  58. Kashino, D., Silverman, J. D., Rodighiero, G., et al. 2013, ApJ, 777, L8 [Google Scholar]
  59. Keel, W. C., Windhorst, R. A., Jansen, R. A., et al. 2023, AJ, 165, 166 [NASA ADS] [CrossRef] [Google Scholar]
  60. Kennicutt, R. C., Jr 1998, ARA&A, 36, 189 [NASA ADS] [CrossRef] [Google Scholar]
  61. Kennicutt, R. C., & Evans, N. J. 2012, ARA&A, 50, 531 [NASA ADS] [CrossRef] [Google Scholar]
  62. Kennicutt, R. C., Jr, Hao, C.-N., Calzetti, D., et al. 2009, ApJ, 703, 1672 [NASA ADS] [CrossRef] [Google Scholar]
  63. Kessler, S., Leroy, A., Querejeta, M., et al. 2020, ApJ, 892, 23 [NASA ADS] [CrossRef] [Google Scholar]
  64. Koekemoer, A. M., Faber, S. M., Ferguson, H. C., et al. 2011, ApJS, 197, 36 [NASA ADS] [CrossRef] [Google Scholar]
  65. Kriek, M., & Conroy, C. 2013, ApJ, 775, L16 [NASA ADS] [CrossRef] [Google Scholar]
  66. Larson, R. L., Finkelstein, S. L., Kocevski, D. D., et al. 2023, ApJ, 953, L29 [NASA ADS] [CrossRef] [Google Scholar]
  67. Le Bail, A., Daddi, E., Elbaz, D., et al. 2024, A&A, 688, A53 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  68. Lecroq, M., Charlot, S., Bressan, A., et al. 2024, MNRAS, 527, 9480 [Google Scholar]
  69. Leja, J., Johnson, B. D., Conroy, C., van Dokkum, P. G., & Byler, N. 2017, ApJ, 837, 170 [NASA ADS] [CrossRef] [Google Scholar]
  70. Lo Faro, B., Buat, V., Roehlly, Y., et al. 2017, MNRAS, 472, 1372 [NASA ADS] [CrossRef] [Google Scholar]
  71. Marchant, P., & Bodensteiner, J. 2023, ARA&A, submitted [arXiv:2311.01865] [Google Scholar]
  72. McKinney, J., Finnerty, L., Casey, C. M., et al. 2023, ApJ, 946, L39 [NASA ADS] [CrossRef] [Google Scholar]
  73. Napolitano, L., Pentericci, L., Santini, P., et al. 2024, A&A, 688, A106 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  74. Noeske, K. G., Weiner, B. J., Faber, S. M., et al. 2007, ApJ, 660, L43 [CrossRef] [Google Scholar]
  75. Osborne, C., & Salim, S. 2024, ApJ, 962, 59 [NASA ADS] [CrossRef] [Google Scholar]
  76. Osterbrock, D. E. 1989, Astrophysics of Gaseous Nebulae and Active Galactic Nuclei (University Science Books) [Google Scholar]
  77. Osterbrock, D. E., & Ferland, G. J. 2006, Astrophysics of Gaseous Nebulae and Active Galactic Nuclei (Sausalito, CA: University Science Books) [Google Scholar]
  78. Pacifici, C., da Cunha, E., Charlot, S., et al. 2015, MNRAS, 447, 786 [NASA ADS] [CrossRef] [Google Scholar]
  79. Pacucci, F., Nguyen, B., Carniani, S., Maiolino, R., & Fan, X. 2023, ApJ, 957, L3 [CrossRef] [Google Scholar]
  80. Pantoni, L., Lapi, A., Massardi, M., et al. 2021, MNRAS, 504, 928 [NASA ADS] [CrossRef] [Google Scholar]
  81. Papovich, C., Rudnick, G., Rigby, J. R., et al. 2009, ApJ, 704, 1506 [CrossRef] [Google Scholar]
  82. Pierini, D., Gordon, K. D., Witt, A. N., & Madsen, G. J. 2004, ApJ, 617, 1022 [NASA ADS] [CrossRef] [Google Scholar]
  83. Prescott, M. K. M., Finlator, K. M., Cleri, N. J., Trump, J. R., & Papovich, C. 2022, ApJ, 928, 71 [NASA ADS] [CrossRef] [Google Scholar]
  84. Reddy, N. A., Kriek, M., Shapley, A. E., et al. 2015, ApJ, 806, 259 [NASA ADS] [CrossRef] [Google Scholar]
  85. Reddy, N. A., Topping, M. W., Shapley, A. E., et al. 2022, ApJ, 926, 31 [NASA ADS] [CrossRef] [Google Scholar]
  86. Reddy, N. A., Topping, M. W., Sanders, R. L., Shapley, A. E., & Brammer, G. 2023, ApJ, 952, 167 [CrossRef] [Google Scholar]
  87. Rezaee, S., Reddy, N., Shivaei, I., et al. 2021, MNRAS, 506, 3588 [CrossRef] [Google Scholar]
  88. Rezaee, S., Reddy, N. A., Topping, M. W., et al. 2023, MNRAS, 526, 1512 [NASA ADS] [CrossRef] [Google Scholar]
  89. Rieke, M. J., Kelly, D. M., Misselt, K., et al. 2023, PASP, 135, 028001 [CrossRef] [Google Scholar]
  90. Robertson, C., Holwerda, B. W., Young, J., et al. 2024, AJ, 167, 263 [CrossRef] [Google Scholar]
  91. Roebuck, E., Sajina, A., Hayward, C. C., et al. 2019, ApJ, 881, 18 [NASA ADS] [CrossRef] [Google Scholar]
  92. Saftly, W., Baes, M., De Geyter, G., et al. 2015, A&A, 576, A31 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  93. Salim, S., Charlot, S., Rich, R. M., et al. 2005, Am. Astron. Soc. Meeting Abstr., 207, 52.06 [NASA ADS] [Google Scholar]
  94. Salim, S., Rich, R. M., Charlot, S., et al. 2007, ApJS, 173, 267 [NASA ADS] [CrossRef] [Google Scholar]
  95. Salim, S., Boquien, M., & Lee, J. C. 2018, ApJ, 859, 11 [Google Scholar]
  96. Salmon, B., Papovich, C., Long, J., et al. 2016, ApJ, 827, 20 [NASA ADS] [CrossRef] [Google Scholar]
  97. Schaerer, D., & de Barros, S. 2010, A&A, 515, A73 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  98. Seillé, L. M., Buat, V., Haddad, W., et al. 2022, A&A, 665, A137 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  99. Shivaei, I., Reddy, N., Rieke, G., et al. 2020, ApJ, 899, 117 [NASA ADS] [CrossRef] [Google Scholar]
  100. Stanway, E. R., & Eldridge, J. J. 2018, MNRAS, 479, 75 [NASA ADS] [CrossRef] [Google Scholar]
  101. Stasińska, G., Cid Fernandes, R., Mateus, A., Sodré, L., & Asari, N. V. 2006, MNRAS, 371, 972 [Google Scholar]
  102. Stefanon, M., Yan, H., Mobasher, B., et al. 2017, ApJS, 229, 32 [NASA ADS] [CrossRef] [Google Scholar]
  103. Theulé, P., Burgarella, D., Buat, V., et al. 2024, A&A, 682, A119 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  104. Tacchella, S., Johnson, B. D., Robertson, B. E., et al. 2023, MNRAS, 522, 6236 [NASA ADS] [CrossRef] [Google Scholar]
  105. Trayford, J. W., Lagos, C. P., Robotham, A. S. G., & Obreschkow, D. 2020, MNRAS, 491, 3937 [NASA ADS] [CrossRef] [Google Scholar]
  106. Valentino, F., Daddi, E., Strazzullo, V., et al. 2015, ApJ, 801, 132 [NASA ADS] [CrossRef] [Google Scholar]
  107. van der Wel, A., Straughn, A. N., Rix, H. W., et al. 2011, ApJ, 742, 111 [NASA ADS] [CrossRef] [Google Scholar]
  108. Villa-Vélez, J. A., Buat, V., Theulé, P., Boquien, M., & Burgarella, D. 2021, A&A, 654, A153 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  109. Whitaker, K. E., Rigby, J. R., Brammer, G. B., et al. 2014, ApJ, 790, 143 [NASA ADS] [CrossRef] [Google Scholar]
  110. Wilkins, S. M., Lovell, C. C., & Stanway, E. R. 2019, MNRAS, 490, 5359 [NASA ADS] [CrossRef] [Google Scholar]
  111. Yoshikawa, T., Akiyama, M., Kajisawa, M., et al. 2010, ApJ, 718, 112 [NASA ADS] [CrossRef] [Google Scholar]
  112. Zackrisson, E., Bergvall, N., & Leitet, E. 2008, ApJ, 676, L9 [NASA ADS] [CrossRef] [Google Scholar]
  113. Zavala, J. A., Buat, V., Casey, C. M., et al. 2023, ApJ, 943, L9 [NASA ADS] [CrossRef] [Google Scholar]

Appendix A: Showcasing the sample SEDs

thumbnail Fig. A.1.

SEDs of the spectro-photometric data with the BC03 and BPASS models. The legend is identical to that of Fig 5.

thumbnail Fig. A.1.

continued.

thumbnail Fig. A.1.

continued.

All Tables

Table 1.

Available photometric bands.

Table 2.

Characteristics of the objects.

Table 3.

Absorption line corrections in Angstroms.

Table 4.

CIGALE modules and input parameters.

Table 5.

Hα contribution to the NIRCam filter.

Table 6.

Attenuation parameters of our subsample.

All Figures

thumbnail Fig. 1.

Example of the LIME fitting process on the Hα + [NII] doublet of ID 6563 (z ∼ 1.7). The top figure shows the Gaussian profile fittings while the lower figure shows the residual between the observation and the theoretical fits. The red shade in the bottom plot shows the calibration pixel error while the yellow shade shows the standard deviation from the continua bands. The gap between the red and yellow regions is intended to avoid potential contamination of uncertain line wings and artefacts.

In the text
thumbnail Fig. 2.

Theoretical line ratios for a case B recombination computed following the BCs component of the Charlot & Fall (2000) attenuation recipe and varying the exponents presented in Eq. 1. nBC = −0.7 corresponds to a shallower curve (solid lines) while nBC = −1.3 corresponds to a steeper curve (dashed lines).

In the text
thumbnail Fig. 3.

Observed fluxes ratios corrected for underlying stellar absorption. The thin coloured curves are power laws following the CF00 framework with varying exponents, from top to bottom: −2.0, −1.3, −1.0, −0.7, −0.4.

In the text
thumbnail Fig. 4.

Snapshots of the stellar continuum of the best model from the fit of spectro-photometric data made with BPASS for ID 5430. All absorption line EWs for this object can be found in Table 3.

In the text
thumbnail Fig. 5.

Comparison of the SEDs of the galaxy ID 5430 fitted with BC03 (left panel) and with BPASS (right panel).

In the text
thumbnail Fig. 6.

Difference of the χ2 between the BC03 and BPASS stellar population models for all galaxies in our sample. A negative value favors BC03 while a positive one favors BPASS.

In the text
thumbnail Fig. 7.

Quality assessment of the photometry fitting for the fits of the spectro-photometric data with BC03 (red) and BPASS (blue). For each model, the median of each filter is represented as a square.

In the text
thumbnail Fig. 8.

Quality assessment of the EWs fitting for the fits of the spectro-photometric data. Filled dots are used for BC03 fits and empty dots for BPASS fits. The linear fit of the Hα points is plotted as a solid magenta line for BC03, dashed magenta for BPASS and the black solid line is the 1:1 line.

In the text
thumbnail Fig. 9.

Quality assessment of the EW estimations when fitting photometric data only. Filled dots are used for BC03 fits and empty dots for BPASS fits. The linear fit of the Hα points is plotted as a solid magenta line for BC03, dashed magenta for BPASS and the black solid line is the 1:1 line.

In the text
thumbnail Fig. 10.

Correlation between the sSFR and the level of contamination in the NIRCam filter containing the Hα line. The error on the y-axis comes from the uncertainty on the attenuation. The dashed lines represent the best fit of each sample of data points.

In the text
thumbnail Fig. 11.

Total attenuation curves of the subsample with a fit of the mean curve in a black solid line. On the rightmost measurement from top to bottom: IDs 8588, 8710, 5409, 23542, 6563 and 5430. The classical CF00 recipe as well as the Calzetti law have been added for comparison. Each total attenuation curve is color-coded with its value of μ presented in Table 6.

In the text
thumbnail Fig. 12.

Comparison of the SFR (left panel) and stellar mass (right panel) obtained through fitting spectro-photometric data in a run made with BC03 and a run made with BPASS. The solid line represents the 1:1 relation while the dashed line represents the mean deviation of our sample (without the three outliers in the case of the SFR).

In the text
thumbnail Fig. 13.

Comparison of the SFR (left panel) and stellar mass (right panel) of BC03 and BPASS when fitting photometry only (y-axis) and spectro-photometry (x-axis).

In the text
thumbnail Fig. 14.

Comparison of the SFR (left panel) and stellar mass (right panel) obtained with a fit in which both exponents are fixed to −0.96 (dots) and a fit following the original CF00 recipe (squares) with respect to the one we get while using the parameters presented in Sect 4.

In the text
thumbnail Fig. A.1.

SEDs of the spectro-photometric data with the BC03 and BPASS models. The legend is identical to that of Fig 5.

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.