Impact of galactic and intergalactic dust on the stellar EBL

Current theories assume that the low intensity of the stellar extragalactic background light (stellar EBL) is caused by finite age of the Universe because the finite-age factor limits the number of photons that have been pumped into the space by galaxies and thus the sky is dark in the night. We oppose this opinion and show that two main factors are responsible for the extremely low intensity of the observed stellar EBL. The first factor is a low mean surface brightness of galaxies, which causes a low luminosity density in the local Universe. The second factor is light extinction due to absorption by galactic and intergalactic dust. Dust produces a partial opacity of galaxies and of the Universe. The galactic opacity reduces the intensity of light from more distant background galaxies obscured by foreground galaxies. The inclination-averaged values of the effective extinction AV\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$A_{V}$\end{document} for light passing through a galaxy is about 0.2 mag. This causes that distant background galaxies become apparently faint and do not contribute to the EBL significantly. In addition, light of distant galaxies is dimmed due to absorption by intergalactic dust. Even a minute intergalactic opacity of 1×10−2mag\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$1 \times 10^{-2}~\mbox{mag}$\end{document} per Gpc is high enough to produce significant effects on the EBL. As a consequence, the EBL is comparable with or lower than the mean surface brightness of galaxies. Comparing both extinction effects, the impact of the intergalactic opacity on the EBL is more significant than the obscuration of distant galaxies by partially opaque foreground galaxies by factor of 10 or more. The absorbed starlight heats up the galactic and intergalactic dust and is further re-radiated at IR, FIR and micro-wave spectrum. Assuming static infinite universe with no galactic or intergalactic dust, the stellar EBL should be as high as the surface brightness of stars. However, if dust is considered, the predicted stellar EBL is about 290nWm−2sr−1\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$290~\mbox{nW}\,\mbox{m}^{-2}\,\mbox{sr}^{-1}$\end{document}, which is only 5 times higher than the observed value. Hence, the presence of dust has higher impact on the EBL than currently assumed. In the expanding universe, the calculated value of the EBL is further decreased, because the obscuration effect and intergalactic absorption become more pronounced at high redshifts when the matter was concentrated at smaller volume than at present.


Introduction
The extragalactic background light (EBL) covers the nearultraviolet, visible and infrared wavelengths from 0.1 to 1000 µm. Measurements of the EBL are provided by data from the Cosmic Background Explorer (COBE) mission, by the Infrared Space Observatory (ISO) instruments and by the Submillimeter Common User Bolometer Array (SCUBA) instrument (for reviews, see Hauser and Dwek 2001;Lagache et al. 2005;Cooray 2016). The direct measurements are supplemented by analysing integrated light from extragalactic source counts which provide a lower limit on the EBL (Madau and Pozzetti 2000;Hauser and Dwek 2001). The upper limits are provided by attenuation of gamma rays from distant blazars due to scattering on the EBL (Kneiske et al. 2004;Dwek and Krennrich 2005;Primack et al. 2011;Gilmore et al. 2012;Biteau and Williams 2015). The spectral energy distribution of the EBL has two distinct maxima: at visible-to-near-infrared wavelengths in the range from 0.7 to 2 µm associated with the radiation of stars,  Domínguez et al. 2011) and at far-infrared wavelengths from 100 to 200 µm associated with the thermal radiation of cold and warm dust in galaxies (Schlegel et al. 1998;Calzetti et al. 2000). Despite the extensive number of measurements of the EBL, the uncertainties in their peak values are still large (see Fig. 1). The total EBL from 0.1 to 1000 µm lies roughly in the range from 40 to 200 nW m −2 sr −1 and half of this value comes from the visible-to-near-infrared part of the spectrum (Hauser and Dwek 2001;Bernstein et al. 2002aBernstein et al. , 2002bBernstein et al. , 2002cMatsumoto et al. 2005;Bernstein 2007;Dwek and Krennrich 2013;Driver et al. 2016). This value is quite low and reflects the fact that sky is dark in the night.
Current theories assume that the low intensity of the EBL is caused primarily by finite age of the Universe and by its expansion (Bondi 1961;Wesson et al. 1987;Wesson 1989;Knutsen 1997). It is argued according to the Olbers' paradox that the infinite static universe predicts a bright sky with intensity of the EBL comparable with surface brightness of stars being thus higher by more than 10 orders than the observed value. Consequently, the dark sky is taken as an important evidence for expanding universe of finite age. The finite age of the Universe implies that galaxies have not had time to populate the intergalactic space with enough photons to make it bright (Wesson 1989). This argument is not, however, fully correct because it neglects an impact of light absorption by interstellar and intergalactic dust on the intensity of the EBL.
The light extinction due to presence of absorbing interstellar dust has been observed, measured and numerically modelled by many authors (Mathis 1990;Charlot and Fall 2000;Draine 2003Draine , 2011Tuffs et al. 2004;Draine and Li 2007;da Cunha et al. 2008;Somerville et al. 2012;Popescu et al. 2011). The rate of light extinction is roughly 0.7-1.0 mag/kpc in the Milky Way (Milne and Aller 1980;Koppen and Vergely 1998), but this value can vary being traced, for example, by dust temperature mapping (Bernard et al. 2010). Obviously, light extinction is observed also in other galaxies and its rate depends on the type of the galaxy, its dust content and the galaxy inclination (Goudfrooij et al. 1994;Calzetti 2001;Holwerda et al. 2005bHolwerda et al. , 2007Lisenfeld et al. 2008;Finkelman et al. 2008Finkelman et al. , 2010. The light extinction due to interstellar dust affects the EBL in two ways. First, it reduces light radiated by galaxies and subsequently their surface brightness. Second, presence of dust in nearby foreground galaxies causes obscuration of distant background galaxies (González et al. 1998;Alton et al. 2001;Holwerda et al. 2005aHolwerda et al. , 2005b. Consequently, the distant obscured galaxies become faint and do not contribute to the EBL significantly. In addition, the brightness of all distant galaxies is decreased because of light absorption by intergalactic dust. Obviously, accurate calculations of the EBL should take into account these effects. In this paper, we calculate the impact of light absorption by the galactic and intergalactic dust on the intensity of the stellar EBL. We show that the key factor responsible for the observed low stellar EBL is not the finite age or the expansion of the Universe (Harrison 1984(Harrison , 1990Wesson 1989;Knutsen 1997) but a low surface brightness of galaxies and a partial opacity of galaxies and of the Universe due to absorbing dust.

Light extinction
Energy emitted by light sources and received at the Earth's surface is controlled by three basic factors: distance of light sources, extinction of light along a light ray by scattering and dust absorption, and obscuration of distant light sources by those at closer distance. The obscuration depends on the number of light sources in unit volume, their size and transparency. For fully opaque sources like stars, the obscuration is most effective, for partially opaque sources like galaxies, the obscuration is suppressed. The energy received at the Earth's surface is calculated by summing contributions of all light sources in a specified universe model.
Energy flux I received per unit area and time from light sources in the static universe is expressed as follows where r is the distance, n is the mean number density of light sources (i.e., the mean number of light sources per unit volume), L is the mean energy radiated by a light source per time (in W), κ is the mean opacity of a light source, λ is the mean absorption coefficient along a ray path, dV is the volume element and coefficient γ is the mean free path of a light ray between light sources (i.e., the mean travelling distance of a photon emitted by one light source and reaching another light source) where a is the mean radius of light sources. The opacity κ in Eq. (1) quantifies how much energy is absorbed by a light source when external light goes through the source. Hence, κ is 1 for a fully opaque galaxy and 0 for a fully transparent galaxy. Terms e −λr and e −(κ/γ )r in Eq. (1) describe the light extinction and the obscuration (Harrison 1990;Knutsen 1997) weighted by the opacity. Integrating Eq. (1) we get and is the luminosity density (in W m −3 ). Alternatively, Eq. (4) can be expressed as is the mean surface energy density (in W m −2 ) radiated by a light source. For a universe with uniformly distributed stars (opacity κ of stars is 1) and for zero interstellar absorption, λ = 0, Eq. (7) yields the equality between the energy received at the unit area of the Earth l E and the mean surface energy l S radiated by a star (in W m −2 ) which is the mathematical formulation of the well-known Olbers' paradox (Harrison 1990;Knutsen 1997). For a universe with uniformly distributed galaxies with opacity κ and for zero intergalactic absorption, λ = 0, Eq. (7) yields where l G is the mean surface energy radiated by a galaxy. Note that factor 4 in Eq. (7) is missing in Eqs. (9) and (10) because l E means the flux coming from the upper hemisphere and received at a unit area with a fixed (vertical) normal. Hence the integration in Eq. (1) is slightly different than in Eqs. (9) and (10) and yields a value which is four times lower (Knutsen 1997).
In the case of fully transparent galaxies (κ = 0) in the fully transparent static universe (λ = 0), the received energy in Eq. (10) diverges. Obviously, this is not realistic, because stars, interstellar dust and intergalactic dust are fully or partially opaque, so κ and λ are apparently non-zero and cause the intensity of light observed at the Earth being low. For example, if 1 % of light energy is absorbed when the ray passes through a partially opaque galaxy and no energy is lost in the intergalactic space, Eq. (10) predicts the observed intensity of the EBL to be 100 times higher than the mean surface energy radiated by galaxies. Since the mean surface brightness of galaxies is extremely low, such intensity of the EBL corresponds effectively to the dark sky in the night.
In order to calculate an accurate value of the intensity of the stellar EBL using Eqs. (4)-(6) we need values of number density n, radius of galaxies a, luminosity density j , galactic opacity κ and intergalactic absorption λ (also called the intergalactic or universe opacity). We shortly review estimates of these parameters based on observations in the next sections.

Number density, galaxy size and luminosity density
The number density is fairly variable because of galaxy clustering and presence of voids in the universe (Peebles 2001;Jones et al. 2004;von Benda-Beckmann and Müller 2008). The number density might be ten times higher or more in clusters, at distances up to 15-20 Mpc than the density averaged over larger distances. The mean value of the number density over hundreds of Mpc is, however, stable. It is derived from the Schechter luminosity function (Schechter 1976) being in the range of 0.010-0.025 h 3 Mpc −3 (Peebles 1993;Peacock 1999;Blanton et al. 2001Blanton et al. , 2003. The size distribution of galaxies is dependent on their luminosity, stellar mass and the morphological type. Observed galaxies cover luminosities between ∼ 10 8 L and ∼ 10 12 L with effective radii between ∼ 0.1 h −1 kpc and ∼ 10 h −1 kpc. For late-type galaxies, the characteristic luminosity in the R-band is −20.5 (Shen et al. 2003). The corresponding Petrosian half-light radius is ∼ 2.5-3.0 h −1 kpc and the R90 radius is about 3 times larger, R90 = 7.5-9 h −1 kpc (Graham and Driver 2005). This is close to a commonly assumed value R = 10 h −1 kpc (Peebles 1993;Peacock 1999).
The luminosity density is a fundamental quantity in observational cosmology standing in the Schechter luminosity function (Schechter 1976). The most recent determination of the optical luminosity function come from large flux-limited redshift surveys such as the Two Degree Field Galaxy Redshift Survey (2dFGRS; Cross et al. 2001), the Sloan Digital Sky Survey (SDSS; Blanton et al. 2001Blanton et al. , 2003 or Century Survey (CS; Geller et al. 1997;Brown et al. 2001). Independent estimates of the luminosity function in the R-band are well consistent being (1.84 ± 0.04) × 10 8 h L Mpc −3 for the SDSS data  and (1.9 ± 0.6) × 10 8 h L Mpc −3 for the CS data (Brown et al. 2001).

Galactic opacity
The galactic opacity can be measured by a variety of methods usually applied to a large set of samples. The most widely used methods test dependence of the surface brightness on inclination, multi-wavelength comparisons, and statistical analysis of the colour and number count variations induced by a foreground galaxy onto background sources (Calzetti 2001).
The most transparent galaxies are elliptical. Goudfrooij et al. (1994) and Goudfrooij and de Jong (1995) found that the observed infrared luminosities are compatible with central optical depths of the diffuse component τ (0) ≤ 0.7, with a typical value of τ (0) ∼ 0.2. The corresponding effective extinction A V is 0.04-0.08 mag. The giant elliptical galaxies found at the centres of cooling flow clusters are often surrounded by extended (≈ 10-100 kpc) Note: A V is the effective inclination-averaged visual extinction, κ is the mean visual galactic opacity of the individual galaxy types, and w is the relative frequency of the galaxy types taken from Table 4 (typical galactic content of regular clusters) of Bahcall (1999). and dusty emission-line nebulae (Voit and Donahue 1997;Donahue et al. 2000). Optical and UV emission-line studies give extinction values in the range A V ≈ 0.3-2.0 mag for the dust associated with the nebula. However, such galaxies are not statistically significant in the population of elliptical galaxies (Calzetti 2001). The dust extinction in spiral and irregular galaxies is higher than in elliptical galaxies. For estimating the extinction by dust in spiral galaxies, Holwerda et al. (2005a) used the so-called Synthetic field method (SFM) which counts the number of background galaxies seen through a foreground galaxy (González et al. 1998;Holwerda et al. 2005a). The advantage of the SFM is that it yields the average opacity for the area of a galaxy disk without making assumptions about either the distribution of absorbers or of the disk starlight. Holwerda et al. (2005b) found that the dust opacity of the disk in the face-on view apparently arises from two distinct components: an optically thicker component (A I = 0.5-4 mag) associated with the spiral arms and a relatively constant optically thinner disk (A I = 0.5 mag). The earlytype spiral disks show less extinction than the later types. As regards the inclination-averaged extinction, the typical values are according to Calzetti (2001): 0.5-0.75 mag for Sa-Sab galaxies, 0.65-0.95 mag for the Sb-Scd galaxies and 0.3-0.4 mag for the irregular galaxies at the B-band.
In summary, galaxies in the local universe are moderately opaque, and extreme values of the opacity are found only in the statistically insignificant more active systems (Calzetti 2001). Adopting estimates of the relative frequency of specific galaxy types in the Universe and their average extinctions (see Table 1), we can calculate their mean visual opacities and finally the overall mean galactic opacity using the weighted average where subscripts 1, 2, and 3 stay for quantities of the elliptical, spiral and lenticular galaxies. According to Eq. (12) Note: n is the number density of galaxies, γ is the mean free path of light between galaxies defined in Eq. (3), κ is the mean opacity of galaxies, λ is the intergalactic absorption, and j R is the R-band luminosity density , I theor is the predicted intensity of the stellar EBL, and I obs is the observed intensity of the stellar EBL. The mean effective radius of galaxies a is considered to be 10 kpc.
and Table 1, the average value of visual opacity κ is about 0.22 ± 0.08. A more accurate approach should take into account the statistical distributions of the galaxy size and of the mean galaxy surface brightness for individual types of galaxies.

Intergalactic opacity
Observations indicate that the absorption of light is not limited to the interstellar medium within galaxies but is present also in the intergalactic space (Nickerson and Partridge 1971;Margolis and Schramm 1977;Chelouche et al. 2007). It is lower by several orders than in galaxies and depends on distance from galaxies. High values of attenuation are detected in galaxy halos (Ménard et al. 2010) and in cluster centres. The attenuation due to dust in the galaxy clusters has been measured by reddening of background objects behind the clusters (Chelouche et al. 2007;Bovy et al. 2008;Muller et al. 2008;Ménard et al. 2010). The attenuation can also be investigated by correlations between the positions of high-redshift QSOs and low-redshift galaxies using catalogues of UVX objects. The excess of high-redshift QSOs around low-redshift galaxies is explained by a model in which dust situated in foreground clusters of galaxies obscures the QSOs lying behind them (Boyle et al. 1988;Romani and Maoz 1992). Based on constructing extinction curves around galaxies, Ménard et al. (2010) found visual attenuation of A V = (1.3 ± 0.1) × 10 −2 mag at distance from a galaxy up to 170 kpc and A V = (1.3 ± 0.3) × 10 −3 mag on large scale at distance up to 1.7 Mpc. Values of the same order are reported for an average visual extinction by intracluster dust also by Muller et al. (2008) and Chelouche et al. (2007). In addition, consistent opacity was recently reported by Xie et al. (2015) who studied the luminosity and redshifts of the quasar continuum at the data sample of ∼ 90,000 objects and estimated the effective dust density nσ V ≈ 0.02 h Gpc −1 at z < 1.5. However, the intergalactic absorption is redshift dependent. According to Davies et al. (1997) the intergalactic extinction increases with redshift and transparent uni-verse becomes significantly opaque (optically thick) at redshifts of z = 1-3. The increase of intergalactic extinction with redshift is confirmed by results of Ménard et al. (2010) who estimated A V to about 0.03 mag at z = 0.5 but to about 0.05-0.09 mag at z = 1.

Predicted and observed stellar EBL
Taking into account estimates of the galactic and intergalactic opacity and other cosmological parameters (see Table 2), the intensity of the stellar EBL calculated by Eqs. (4) and (5) lies for wavelengths between 300 and 3500 nm in the range of 190-560 nW m −2 sr −1 with the optimum value of I theor ≈ 290 nW m −2 sr −1 .
(13) Figure 2 indicates that the EBL is rather insensitive to mean free path between galaxies γ but quite sensitive to intergalactic opacity λ and luminosity density j . Obviously, high values of the EBL are produced for high luminosity density and low intergalactic opacity. The observed intensity of the stellar EBL (Hauser and Dwek 2001;Bernstein et al. 2002aBernstein et al. , 2002bBernstein et al. , 2002cBernstein 2007) lies in the range of 20-140 nW m −2 sr −1 (see Fig. 1) with the optimum value of Hence the predicted stellar EBL is about 5 times higher than the observed EBL. This result is surprising because it is commonly assumed that the EBL calculated for the infinite static universe must be higher than the observed EBL by more than 10 orders. The rather low value of the EBL in Eq. (13) evidences that the key factor for a successful prediction of the EBL is including the effects of the galactic and intergalactic absorption of light by dust. In fact, considering an expanding universe in the EBL calculations is not essential. Substituting the infinite static universe by expanding universe of finite age in the EBL calculations further reduces the predicted EBL but by factor of 5 only. Fig. 2 (a) The stellar EBL as a function of intergalactic opacity λ and mean free path between galaxies γ . (b) The stellar EBL as a function of the intergalactic opacity λ and relative luminosity density j . The relative luminosity density is normalized to its optimum value (see Table 2). The black open circles mark the positions of optimally chosen parameters. The colour-coded EBL is in W m −2 sr −1

Obscuration by galaxies versus intergalactic opacity
As shown above, the low value of the stellar EBL is caused by two effects: (1) obscuration of background galaxies by partially opaque foreground galaxies, and (2) the intergalactic opacity. Considering observations of the number density of galaxies and of the galactic and intergalactic opacity, we can determine which effect has a more significant impact on the EBL. We calculate ratio k that is higher (lower) than 1 for the intergalactic opacity (obscuration of galaxies) being predominant. Figure 3 shows ratio k as a function of intergalactic opacity λ and mean free path γ . The ratio is calculated for three values of the mean galactic opacity, κ = 0.30, 0.22 and 0.14. The corresponding values of k for optimum values of λ and κ are 10, 13 and 21. This evidences that the EBL is affected predominantly by the intergalactic opacity. The impact of the obscuration effect on the EBL is almost negligible for the majority of combinations of realistic values of λ and κ.

Discussion
The calculations prove that the impact of light absorption by galactic and intergalactic dust on the stellar EBL is significant. It is even more significant than the expansion, redshift or the finite age of the Universe. The absorbed starlight heats up the dust and is further re-radiated at the IR, FIR and micro-wave spectrum. Although the presented calculations are rough approximations, the estimate of the stellar EBL is robust and reliable. It is based on: (1) the luminosity density measurements (Blanton et al. 2001Brown et al. 2001), (2) the estimate of the effective inclination-averaged Fig. 3 Ratio k evaluating the impact of the intergalactic opacity and the obscuration of galaxies on the stellar EBL. The ratio is shown as a function of intergalactic opacity λ and mean free path between galaxies γ . The galactic opacity is: (a) κ = 0.30, (b) κ = 0.22, and (c) κ = 0.14. Considering optimum values of λ and γ , ratio k is: (a) k = 10, (b) k = 13, and (c) k = 21 galactic opacity (Calzetti 2001), and (3) the estimate of the intergalactic opacity (Xie et al. 2015). Even a minute inter-galactic opacity of 1 × 10 −2 mag Gpc −1 is high enough to produce significant effects on the EBL. In addition, galaxies are partially opaque due to galactic dust. Consequently, the opacity of foreground galaxies reduces the intensity of light radiated by distant background galaxies. The background galaxies become apparently faint and do not contribute to the stellar EBL significantly. As a result, the stellar EBL is comparable with or even lower than the mean surface brightness of galaxies. Comparing the impact of the intergalactic opacity and obscuration of galaxies on the EBL, the intergalactic opacity is more significant by factor of 10 or more. If the intergalactic opacity and the obscuration effects are neglected, the EBL predicted for the static universe should be according to the Olbers' paradox as high as the surface brightness of stars. This value is higher by almost 14 orders than the observed one. Hence, the galactic and intergalactic absorption is the most important factor for the observed low stellar EBL. The corrections for the expansion of the Universe and for the redshift alter the predictions by less than one order only.