Niebla: an open-source code for modeling the extragalactic background light
Pith reviewed 2026-05-18 04:24 UTC · model grok-4.3
The pith
Niebla is an open-source code that models the extragalactic background light using stellar evolution and dust re-emission to show its effect on very-high-energy gamma-ray spectra.
A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.
Core claim
Niebla computes the EBL by evolving a single stellar population spectrum with redshift, mean metallicity and star-formation-rate density, then absorbs the optical emissivity with dust and re-emits it in the infrared using multiple prescriptions. Three resulting EBL models fitted to data are presented, and a simulated LHAASO observation of Markarian 501 in a high-flux state shows that the VHE spectrum is highly sensitive to infrared EBL photon density, enabling distinction between dust reemission models and constraints on EBL parameters with future observations.
What carries the argument
The Niebla code, which generates EBL photon density by evolving single stellar population spectra and applying selectable dust re-emission prescriptions to determine gamma-ray optical depth.
If this is right
- Future VHE gamma-ray observations can distinguish between different dust re-emission models.
- EBL parameters can be constrained using data from instruments such as LHAASO.
- The code enables detailed modeling of EBL optical depth effects on gamma-ray spectra.
- Users can add extra EBL sources such as axion decay or intra-halo light to test their influence.
Where Pith is reading between the lines
- The infrared sensitivity may allow gamma-ray data to indirectly test dust properties at high redshift.
- Open availability of the code supports community refinement of EBL models as more observations accumulate.
- The framework could extend to other high-energy processes attenuated by background photon fields.
Load-bearing premise
The single stellar population evolution with mean metallicity and star-formation-rate density, together with the chosen dust re-emission prescriptions, sufficiently captures the true EBL without major unaccounted contributions or incorrect evolution assumptions.
What would settle it
A measured very-high-energy spectrum from a blazar such as Markarian 501 whose attenuation does not match predictions from any of the provided dust re-emission models at infrared wavelengths.
Figures
read the original abstract
The flux of extragalactic gamma rays is attenuated through interactions with optical and infrared photons of the extragalactic background light (EBL). The EBL is an isotropic, diffuse photon field that is difficult to measure directly at these wavelengths due to strong foreground emission. We present niebla, the first open-source code to compute the EBL from optical to far-infrared wavelengths using a phenomenological approach that accepts fully customizable inputs. This software enables a detailed modeling of the influence of EBL optical depth on gamma-ray observations and facilitates the distinction between different dust reemission models. The code models the optical background primarily from stellar emission, by evolving the spectrum of a single stellar population as a function of redshift, considering mean metallicity evolution and star formation rate density. Additional sources to the EBL can be provided by the user. The code already includes optional contributions from, e.g., stripped stars, intra-halo light, or the decay of axion dark matter. The optical emissivity is then absorbed by interstellar dust and reemitted in the infrared regime. We provide multiple prescriptions to model this process, using spectral dust templates or a combination of blackbodies. We provide three EBL models calculated with different dust reemission prescriptions, which have been fitted to various observational data sets. In addition, we showcase the versatility of our model through a simulated observation of the blazar Markarian 501 in a high-flux state with the LHAASO array. We find that the simulated VHE spectrum is highly sensitive to the photon density of the EBL at infrared wavelengths. Our model will therefore allow the community to distinguish between different dust reemission models and constrain EBL parameters with future observations.
Editorial analysis
A structured set of objections, weighed in public.
Referee Report
Summary. The manuscript introduces niebla, the first open-source code for computing the extragalactic background light (EBL) from optical to far-infrared wavelengths via a phenomenological model. Optical emissivity is generated by evolving a single stellar population spectrum with mean metallicity evolution and a star-formation rate density SFRD(z); this emissivity is absorbed by dust and re-emitted in the IR using user-selectable prescriptions (spectral templates or blackbody combinations). Three specific EBL models are supplied after fitting to observational data sets, with optional user-added components such as stripped stars, intra-halo light, or axion decay. The code is demonstrated through a simulated LHAASO observation of the blazar Markarian 501 in a high state, where the VHE spectrum is shown to be sensitive to the IR EBL photon density; the authors conclude that the framework will enable the community to distinguish dust re-emission models and constrain EBL parameters with future observations.
Significance. If the modeling assumptions are shown to be robust, niebla would provide a valuable, customizable open-source tool for gamma-ray astrophysics and cosmology, allowing reproducible exploration of EBL attenuation effects on VHE spectra. The provision of multiple dust prescriptions, fully customizable inputs, and optional non-stellar contributions (e.g., axion decay) are clear strengths that facilitate community use and falsifiable tests. The LHAASO simulation illustrates practical utility for future observations. Significance is reduced, however, until the relative size of systematics from the single-population stellar modeling versus dust-model variations is quantified.
major comments (1)
- [§2] §2 (Phenomenological modeling of stellar emission): The adoption of a single stellar population evolved with mean metallicity and a single SFRD(z) to generate optical emissivity omits the dispersion in star-formation histories, metallicities, and dust geometries across the galaxy population. Because the VHE attenuation kernel is most sensitive in the 10–100 µm range, any resulting systematic offset in the absorbed energy budget and IR photon density could be comparable to the template-to-blackbody variations highlighted in the three EBL models. This directly affects the central claim (abstract and §4 simulation) that future LHAASO-type spectra will cleanly separate the dust re-emission models; a quantitative comparison to multi-population or semi-analytic EBL calculations is needed to establish that the dust-model differences exceed these systematics.
minor comments (3)
- [Abstract and §3] Abstract and §3 (fitting procedure): No quantitative validation metrics (e.g., χ² values, residual plots, or posterior uncertainties) are reported for the fits of the three EBL models to the observational data sets; this makes it difficult to assess how post-hoc choices in dust parameters propagate into the final photon densities.
- [§4] §4 (LHAASO simulation): The simulated VHE spectrum sensitivity is presented without error bands arising from stellar-model uncertainties or from the fitting procedure; adding these would clarify whether the dust-model separation remains distinguishable.
- Code and data availability: While the manuscript states that niebla is open-source, the repository link, version number, and exact input files used for the three published EBL models should be provided to ensure full reproducibility.
Simulated Author's Rebuttal
We thank the referee for their constructive review and for highlighting an important limitation in our phenomenological modeling approach. We address the major comment in detail below, have revised the manuscript to incorporate appropriate caveats and clarifications, and note one aspect that remains a standing limitation.
read point-by-point responses
-
Referee: [§2] §2 (Phenomenological modeling of stellar emission): The adoption of a single stellar population evolved with mean metallicity and a single SFRD(z) to generate optical emissivity omits the dispersion in star-formation histories, metallicities, and dust geometries across the galaxy population. Because the VHE attenuation kernel is most sensitive in the 10–100 µm range, any resulting systematic offset in the absorbed energy budget and IR photon density could be comparable to the template-to-blackbody variations highlighted in the three EBL models. This directly affects the central claim (abstract and §4 simulation) that future LHAASO-type spectra will cleanly separate the dust re-emission models; a quantitative comparison to multi-population or semi-analytic EBL calculations is needed to establish that the dust-model differences exceed these systematics.
Authors: We agree that the single stellar population with mean metallicity and a single SFRD(z) is a simplification that does not capture the full dispersion in star-formation histories, metallicities, and dust geometries present in the real galaxy population. This choice was made to keep the model computationally lightweight, fully customizable, and focused on the dust re-emission component, which is the primary novelty of the code. The manuscript already emphasizes that additional non-stellar contributions can be added by the user; the same modular structure allows replacement or extension of the stellar emissivity module. In the revised version we have added an explicit limitations paragraph in §2 and a corresponding caveat in §4 and the abstract, stating that the reported differences between the three dust prescriptions are shown under the adopted single-population stellar assumptions and that systematic uncertainties arising from stellar-population dispersion could affect the absolute IR photon density at a level comparable to or larger than the template variations. We have also revised the language of the central claim to indicate that the framework enables the community to explore and distinguish dust re-emission models once the stellar modeling uncertainties are separately quantified or marginalized. A full quantitative error-budget comparison against multi-population or semi-analytic EBL calculations (e.g., those based on GALFORM or other SAMs) would require new integrations and is not performed in the present work. revision: partial
- Quantitative comparison of the single-population EBL models against multi-population or semi-analytic calculations, including a direct assessment of whether dust-model differences exceed stellar-population systematics in the 10–100 µm range relevant for VHE attenuation.
Circularity Check
No significant circularity in the phenomenological EBL modeling chain
full rationale
The niebla framework is explicitly constructed as a customizable, open-source code that accepts arbitrary user inputs for stellar population evolution, metallicity, SFRD, dust re-emission prescriptions, and optional additional components. The three provided EBL models are stated to have been fitted to observational datasets, but this is presented as an application of the tool rather than a derivation that reduces to its own inputs by construction. The simulated VHE spectrum for Markarian 501 is a forward-modeling exercise demonstrating sensitivity to IR photon density; it does not constitute a tautological prediction forced by prior fits. No self-definitional equations, load-bearing self-citations, uniqueness theorems, or ansatz smuggling appear in the described methodology. The derivation from single-population emissivity through dust absorption/re-emission remains independent and falsifiable against external data.
Axiom & Free-Parameter Ledger
free parameters (3)
- star formation rate density evolution
- mean metallicity evolution parameters
- dust re-emission template or blackbody parameters
axioms (2)
- domain assumption Optical background is modeled primarily from evolving stellar emission
- domain assumption Optical emissivity is absorbed by interstellar dust and re-emitted in the infrared
Lean theorems connected to this paper
-
IndisputableMonolith/Cost/FunctionalEquation.leanwashburn_uniqueness_aczel unclear?
unclearRelation between the paper passage and the cited Recognition theorem.
The code models the optical background primarily from stellar emission, by evolving the spectrum of a single stellar population as a function of redshift, considering mean metallicity evolution and star formation rate density.
-
IndisputableMonolith/Foundation/AlphaCoordinateFixation.leanalpha_pin_under_high_calibration unclear?
unclearRelation between the paper passage and the cited Recognition theorem.
We provide multiple prescriptions to model this process, using spectral dust templates or a combination of blackbodies.
What do these tags mean?
- matches
- The paper's claim is directly supported by a theorem in the formal canon.
- supports
- The theorem supports part of the paper's argument, but the paper may add assumptions or extra steps.
- extends
- The paper goes beyond the formal theorem; the theorem is a base layer rather than the whole result.
- uses
- The paper appears to rely on the theorem as machinery.
- contradicts
- The paper's claim conflicts with a theorem or certificate in the canon.
- unclear
- Pith found a possible connection, but the passage is too broad, indirect, or ambiguous to say the theorem truly supports the claim.
Reference graph
Works this paper leans on
-
[1]
Stellar contribution A simple stellar population (SSP) is a group of stars born from a single gas cloud. The time scale for this creation is of the order of 10 7 years, which can be ap- proximated as instantaneous in cosmic terms. Following this first burst no more star formation occurs. The emis- sivity of a stellar population can be calculated as ενstel...
-
[2]
Intra-halo light The IHL is the light caused by stars that have been ejected from their galaxies by dynamical events. These stars would create a faint halo surrounding their host galaxies, but would not have been used to calculate the ρ⋆. We follow the parametrization described by Bernal et al. [8] to calculate the emissivity ελ,IHL(λ, z) = Z Mmax Mmin dn...
-
[3]
Cosmic ALP decay Axions are hypothetical pseudo-Nambu-Goldstone bosons, which are a consequence of the Peccei-Quinn so- lution to the strong CP problem in QCD. In contrast to axions, ALPs do not solve the strong CP problem and their massm a is independent of the coupling constant to photonsg aγ. Both axions and ALPs are candidates to make up cold dark mat...
-
[4]
Other contributions The list of sources introduced above is not exhaustive. Additional contributions to the COB could come, e.g., from AGN as their spectra span the entire electromag- netic spectrum. However their contribution to the EBL intensity is as low as∼10% of the stellar one [5, 34]. We do not include a parametrization for this source in niebla. H...
-
[5]
TheCharyspectral models are normalized to the total luminosity of ULIRGs
Spectral templates The first library we use is a set of synthetic galactic spectra derived from observations of ultra luminous in- frared galaxies (ULIRGs) [13], which we refer to asChary from here on. TheCharyspectral models are normalized to the total luminosity of ULIRGs. We need to rescale them to the luminosity of a single SSP, which we do by introdu...
-
[6]
Black body spectra The second method to calculate the infrared emission assumes different dust populations, where the reemitted light follows a black-body spectrum for each population peaking at different temperatures. For example, poly- cyclic aromatic hydrocarbon and very small grains are heated in galaxies toT∼ O(100 K), and therefore emit in the mid i...
-
[7]
The metallicity evolution is taken from Tanikawa et al
to parametrize the SFRD. The metallicity evolution is taken from Tanikawa et al. [55]. We assume the same functional form for the dust absorption as Finke et al. [20]. The mathematical formulas for these models are compiled in Appendix A. Each of the dust reemission models has a different number of free parameters. We introduce two free pa- rameters forCh...
work page 2004
-
[8]
to characterize the detector. A flare in the energy range between a few hundred GeV and 20 TeV would be observed with the water Cherenkov detector array (WCDA), which is sensitive to energies 100 GeV-30 TeV. The number of counts that would be observed is given by µ(∆Eobs) = ∆T Z ∆Eobs dEobs × Z ∞ 0 R(E obs, Etrue)A eff (Etrue)ϕ(E true)dE true, (18) where ...
work page 1997
-
[9]
Abdalla, H., Abe, H., Acero, F., et al. 2021, J. Cosmology Astropart. Phys., 2021, 048
work page 2021
-
[10]
Ackermann, M., Ajello, M., Allafort, A., et al. 2012, Sci- ence, 338, 1190
work page 2012
-
[11]
Aharonian, F. A., Akhperjanian, A. G., Barrio, J. A., et al. 1999, A&A, 349, 11
work page 1999
-
[12]
Aharonian, F. A., Akhperjanian, A. G., Barrio, J. A., et al. 2001, A&A, 366, 62
work page 2001
-
[13]
Andrews, S. K., Driver, S. P., Davies, L. J. M., Lagos, C. d. P., & Robotham, A. S. G. 2017, Monthly Notices of the Royal Astronomical Society, 474, 898
work page 2017
-
[14]
Asplund, M., Grevesse, N., Sauval, A. J., & Scott, P. 2009, ARA&A, 47, 481
work page 2009
-
[15]
2022, Chinese Physics C, 46, 035001
Bai, X., Bi, B., Bi, X., et al. 2022, Chinese Physics C, 46, 035001
work page 2022
-
[16]
L., Caputo, A., Sato-Polito, G., Mirocha, J., & Kamionkowski, M
Bernal, J. L., Caputo, A., Sato-Polito, G., Mirocha, J., & Kamionkowski, M. 2023, Phys. Rev. D, 107, 103046
work page 2023
-
[17]
Biasuzzi, B., Hervet, O., Williams, D. A., & Biteau, J. 2019, A&A, 627, A110
work page 2019
-
[18]
2023, The Multi-Messenger Extragalactic Spectrum
Biteau, J. 2023, The Multi-Messenger Extragalactic Spectrum
work page 2023
- [19]
- [20]
- [21]
-
[22]
Cherenkov Telescope Array Consortium, Acharya, B. S., Agudo, I., et al. 2019, Science with the Cherenkov Tele- scope Array
work page 2019
-
[23]
2011, European Physical Journal C, 71, 1554
Cowan, G., Cranmer, K., Gross, E., & Vitells, O. 2011, European Physical Journal C, 71, 1554
work page 2011
-
[24]
de Matos Pimentel, D. R. & Moura-Santos, E. 2019, J. Cosmology Astropart. Phys., 2019, 043
work page 2019
-
[25]
Dembinski, H. et al. 2020
work page 2020
-
[26]
Driver, S. P., Andrews, S. K., Davies, L. J., et al. 2016, The Astrophysical Journal, 827, 108, aDS Bibcode: 2016ApJ...827..108D
work page 2016
-
[27]
Fermi-LAT Collaboration, Abdollahi, S., Ackermann, M., et al. 2018, Science, 362, 1031
work page 2018
-
[28]
D., Ajello, M., Dom´ ınguez, A., et al
Finke, J. D., Ajello, M., Dom´ ınguez, A., et al. 2022, ApJ, 941, 33
work page 2022
-
[29]
Finke, J. D., Razzaque, S., & Dermer, C. D. 2010, ApJ, 712, 238
work page 2010
-
[30]
Fioc, M. & Rocca-Volmerange, B. 2019, Astronomy & Astrophysics, 623, A143
work page 2019
-
[31]
Forbes, D. A. & Ponman, T. J. 1999, Monthly Notices of the Royal Astronomical Society, 309, 623
work page 1999
- [32]
-
[33]
Franceschini, A., Rodighiero, G., & Vaccari, M. 2008, A&A, 487, 837
work page 2008
-
[34]
Gould, R. J. & Schr´ eder, G. P. 1967, Physical Review, 155, 1408
work page 1967
-
[35]
Gould, R. J. & Schr´ eder, G. P. 1967, Physical Review, 155, 1404
work page 1967
-
[36]
2024, The Astrophysical Journal Letters, 975, L18
Gr´ eaux, L., Biteau, J., & Nievas Rosillo, M. 2024, The Astrophysical Journal Letters, 975, L18
work page 2024
-
[37]
G¨ otberg, Y., de Mink, S. E., Groh, J. H., Leitherer, C., & Norman, C. 2019, A&A, 629, A134
work page 2019
-
[38]
H. E. S. S. Collaboration, Abramowski, A., Acero, F., et al. 2013, A&A, 550, A4
work page 2013
-
[39]
Haardt, F. & Madau, P. 2012, The Astrophysical Journal, 746, 125, aDS Bibcode: 2012ApJ...746..125H
work page 2012
- [40]
-
[41]
Hill, R., Masui, K. W., & Scott, D. 2018, Applied Spec- troscopy, 72, 663, arXiv:1802.03694 [astro-ph]
work page internal anchor Pith review Pith/arXiv arXiv 2018
-
[42]
Khaire, V. & Srianand, R. 2019, Monthly Notices of the Royal Astronomical Society, 484, 4174
work page 2019
-
[43]
Kneiske, T. M. & Dole, H. 2010, Astronomy and Astro- physics, 515, A19, arXiv:1001.2132 [astro-ph]
work page internal anchor Pith review Pith/arXiv arXiv 2010
-
[44]
M., Mannheim, K., & Hartmann, D
Kneiske, T. M., Mannheim, K., & Hartmann, D. H. 2002, Astronomy & Astrophysics, 386, 1, arXiv:astro- ph/0202104
- [45]
-
[46]
2014, The Astrophysical Journal Supplement Series, 212, 14, aDS Bibcode: 2014ApJS..212...14L
Leitherer, C., Ekstr¨ om, S., Meynet, G., et al. 2014, The Astrophysical Journal Supplement Series, 212, 14, aDS Bibcode: 2014ApJS..212...14L
work page 2014
- [47]
-
[48]
Madau, P. & Dickinson, M. 2014, Annual Review of Astronomy and Astrophysics, 52, 415, arXiv:1403.0007 [astro-ph]
work page internal anchor Pith review Pith/arXiv arXiv 2014
-
[49]
Matsumoto, T., Kim, M. G., Pyo, J., & Tsumura, K. 2015, ApJ, 807, 57
work page 2015
- [50]
- [51]
- [52]
-
[53]
Moll´ a, M., Garc´ ıa-Vargas, M. L., & Bressan, A. 2009, MNRAS, 398, 451
work page 2009
-
[54]
Murray, S. G., Power, C., & Robotham, A. S. G. 2013, Astronomy and Computing, 3, 23
work page 2013
- [55]
-
[56]
Overduin, J. M. & Wesson, P. S. 2004, Physics Reports, 402, 267, arXiv:astro-ph/0407207
work page internal anchor Pith review Pith/arXiv arXiv 2004
- [57]
- [58]
-
[59]
Postman, M., Lauer, T. R., Parker, J. W., et al. 2024, arXiv e-prints, arXiv:2407.06273
-
[60]
Probing the peak of the star formation rate density with the extragalactic background light
Raue, M. & Meyer, M. 2012, Monthly Notices of the Royal Astronomical Society, 426, 1097, arXiv:1203.0310 [astro-ph]
work page internal anchor Pith review Pith/arXiv arXiv 2012
-
[61]
Razzaque, S., Dermer, C. D., & Finke, J. D. 2009, The Astrophysical Journal, 697, 483
work page 2009
- [62]
-
[63]
Tanikawa, A., Yoshida, T., Kinugawa, T., et al. 2022, ApJ, 926, 83
work page 2022
- [64]
-
[65]
Windhorst, R. A., Cohen, S. H., Jansen, R. A., et al. 2023, The Astronomical Journal, 165, 13 Appendix A: Parametrizations already listed in the EBL code We present a list of all the parametrizations present in the EBL code in Table V. The possibility to input custom parameters and functions is also available for all the expressions that we list here. App...
work page 2023
-
[66]
use. We also share the dust absorption methodol- ogy, so the spectral shapes are very similar in the optical regime. Moreover, our 2BB model shares their methodol- ogy of dust reemission in the MIR, so these two are sim- ilar in this regime. Saldana-Lopez et al. [54] derive the EBL from observed spectral energy densities of galaxies without fits to the ac...
work page 1997
discussion (0)
Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.