Components of the Diffuse Ultraviolet Radiation at High Latitudes

TL;DR

This study analyzes the diffuse ultraviolet background at high Galactic latitudes using GALEX data, identifying components and dust properties, and highlighting an unexplained zero-offset in UV radiation.

Contribution

It quantifies the components of UV background radiation and determines dust optical constants at high latitudes, revealing an unexplained zero-offset in UV measurements.

Findings

Zero-offset in UV background at high latitudes (~240-394 photon units)

Approximately half of the UV radiation is from extragalactic sources and dust-scattered starlight

Dust grain optical constants determined as a=0.4 and g=0.8 (FUV), a=0.4 and g=0.5 (NUV)

Abstract

We have used data from the Galaxy Evolution Explorer to study the different components of the diffuse ultraviolet background in the region between the Galactic latitudes 70-80 degree. We find an offset at zero dust column density (E(B - V) = 0) of $240 \pm 18$ photon units in the FUV (1539A) and $394 \pm 37$ photon units in the NUV (2316A). This is approximately half of the total observed radiation with the remainder divided between an extragalactic component of $114 \pm 18$ photon units in the FUV and $194 \pm 37$ photon units in the NUV and starlight scattered by Galactic dust at high latitudes. The optical constants of the dust grains were found to be a=0.4$\pm$0.1 and g=0.8$\pm$0.1 (FUV) and a=0.4$\pm$0.1 and g=0.5$\pm$0.1 (NUV). We cannot differentiate between a Galactic or extragalactic origin for the zero-offset but can affirm that it is not from any known source.

Tables

Table 1. Table 1: Correlation coefficients

Wavelength	r	a	b	${\chi }^2$
FUV–100$\mu$m (1.13 MJy sr-1)	0.03	7.97	353.89	0.84
FUV–100$\mu$m ($\ge$1.13 MJy sr-1)	0.46	110.35	238.20	0.70
FUV–E(B-V) ( 0.027 mag)	0.06	942.14	339.97	0.82
FUV–E(B-V) ($\ge$ 0.027 mag)	0.47	4871.51	233.88	0.69
NUV–100$\mu$m	0.25	41.63	587.47	0.94
NUV–E(B-V)	0.27	1677.25	589.33	0.52
NGP∗
FUV–100$\mu$m (1.08 MJy sr-1)	0.27	57.43	288.27	1.15
FUV–100$\mu$m ($>$1.08 MJy sr-1)	0.57	156.33	182.10	1.37
FUV–E(B-V)	0.52	4245.40	250.11	1.24
NUV–100$\mu$m	0.42	68	530.89	1.18
NUV–E(B-V)	0.40	2655.67	531.06	1.21
r - Spearman’s correlation coefficient (P 0.05 for all cases)
a - Scale factor
b - Offset (photon units)
∗ From Akshaya et al. (2018)

Table 2. Table 2: Albedo and phase function asymmetry factor from literature

References	Wavelength (Å)	$a$	$g$
Onaka & Kodaira (1991)	1450–1800	$>$0.32	$>$0.5
Henry & Murthy (1993)	1500	$>0.5$	$>0.7$
Witt et al. (1993)	1000–1600	0.42$\pm$0.04	0.75
Murthy et al. (1993)	912–1150	$>0.3$	0.8
Gordon et al. (1994)	1362	0.47–0.7	0.8
	1769	0.55–0.73	0.8
Witt & Petersohn (1994)	1500	0.5	0.9
Sasseen & Deharveng (1996)	1400–1800	0.3	0.8
Witt et al. (1997)	1400–1800	0.45$\pm$0.05	0.68$\pm$0.1
Schiminovich et al. (2001)	1740	0.45$\pm$0.05	0.77$\pm$0.1
Burgh et al. (2002)	900–1400	0.2–0.4	0.85
Mathis et al. (2002)	1300	$\ge$0.5	0.6–0.85
Shalima & Murthy (2004)	1100	0.4$\pm$0.2	-
Sujatha et al. (2005)	1100	0.4$\pm$0.1	0.55$\pm$0.25
Shalima et al. (2006)	900–1200	0.3–0.7	0.55–0.85
Lee et al. (2008)	1370–1670	0.36$\pm$0.2	0.52$\pm$0.22
Sujatha et al. (2009)	1350–1750	0.4	0.7
Puthiyaveettil et al. (2010)	1400–1900	0.6	0.8
Sujatha et al. (2010)	1350–1750	0.32$\pm$0.09	0.51$\pm$0.19
	1750–2850	0.45$\pm$0.08	0.56$\pm$0.1
Jo et al. (2012)	1350–1750	0.42$\pm$0.03	0.45$\pm$0.02
Choi et al. (2013)	1330–1780	0.38$\pm$0.06	0.46$\pm$0.06
Hamden et al. (2013)	1344–1786	0.62$\pm$0.04	0.78$\pm$0.05
Lim et al. (2013)	1360–1680	0.42$\pm$0.05	0.2–0.58
Murthy (2016)	1500	0.4	0.6
	2300	0.4	0.6
Mattila et al. (2018)	3500	0.58$\pm$0.05	0.6$\pm$0.1

Full text

††pubyear: 2019††pagerange: Components of the Diffuse Ultraviolet Radiation at High Latitudes–Components of the Diffuse Ultraviolet Radiation at High Latitudes

1 Introduction

The diffuse ultraviolet radiation field at the Galactic poles has been assumed to be a combination of diffuse Galactic light from the scattering of starlight by high latitude dust and extragalactic radiation but with the relative contributions unknown (Bowyer, 1991; Henry, 1991). Later observations and a better understanding of the known Galactic and extragalactic components found that they were not sufficient to account for the observed light and that there was a leftover component of 200 – 300 photon units111photons s*-1* cm*-2* sr*-1* Å*-1* (Onaka & Kodaira, 1991; Hamden et al., 2013; Henry et al., 2015; Boissier et al., 2015).

Akshaya et al. (2018) used observations from the Galaxy Evolution Explorer (GALEX) in the far ultraviolet (FUV: 1344 – 1786 Å) and the near ultraviolet (NUV: 1771 – 2831 Å) to explore the contributions of the different components at both Galactic poles ($|b|\geq 80\degree$) in conjunction with a model for the dust scattered light (Murthy, 2016). They found that there was an offset at zero E(B-V) of about 150 photon units in the FUV and 350 photon units in the NUV over and beyond the known extragalactic contributors. The dust scattered light was never more than about 30% of the total signal, even at high column densities. Molecular hydrogen fluorescence kicked in for column densities of greater than log N(H) = 20.2 (where N(H) is the column density in cm*-2*), significantly lower than the canonical values of log N(H)$>$20.6 (Savage et al., 1977; Franco & Cox, 1986; Reach et al., 1994).

We have chosen to study the region between 70° – 80° in Galactic latitude (Fig. 1). The original motivation was to relook at a cirrus cloud discovered by Haikala et al. (1995) (Fig. 2) using Far Ultraviolet Space Telescope (FAUST: Bowyer et al. (1993)) data. These points occupy a limited range with 100µm values between 1 – 4 MJy sr*-1*(0.2 < E(B - V) < 0.12 mag) and, based on a suggestion from an anonymous referee, we expanded our area of interest to the entire range between 70$\degree$ and 80$\degree$. The Haikala cloud may be seen to be part of much larger complex which includes Markkanen’s Cloud (Markkanen, 1979) near the North Galactic Pole (NGP). Cirrus clouds at high latitudes stand out against the general infrared background at 60 and 100µm (Low et al., 1984). They scatter light from Galactic plane stars (Jura, 1979) and block extragalactic radiation (Mattila et al., 2017; Mattila & Väisänen, 2019). Denser regions of the clouds may contain molecular hydrogen (Weiland et al., 1986; Deul & Burton, 1990), which emits in the FUV (Martin et al., 1990; Hurwitz, 1998). We have decomposed the different components of the diffuse UV radiation through a multi-wavelength analysis of the observed radiation at these latitudes. Although there is still not much interstellar dust at these latitudes, the peak column densities are higher than in the NGP and the effects of Galactic absorption and emission are more readily observed.

2 Data

We have collected archival data on the diffuse emission in this field from the UV to the IR. The UV data are from Murthy (2014) who extracted the diffuse emission in the NUV and FUV bands from GALEX (Martin et al., 2005); the IR data at 100µm from Schlegel et al. (1998) who used the data from the Infrared Astronomy Satellite (IRAS: Neugebauer et al. (1984)) and the Cosmic Background Explorer (COBE: Boggess et al. (1992)) to produce a reprocessed 100µm all sky map; and the dust extinction from the $Planck$ mission (Planck Collaboration et al., 2014). We have placed all the data on a common reference frame with 6*′* bins and displayed them in Fig. 1. The cirrus features in this field are distinctly seen in the GALEX FUV, 100µm and E(B-V) maps but not in the NUV, suggesting that dust scattering from the clouds is a small part of the total emission.

We have plotted the UV observations as a function of the IR and the E(B - V) in Fig. 3 and tabulated the correlations in Table 1. The NUV surface brightness is linearly correlated with both the 100µm and the E(B -V). The FUV surface brightness is also correlated with the IR and the E(B - V) but with an inflection point at $1.13\pm 0.13$ MJy sr*-1* (0.027 mag in E(B - V)) where the slope increases. This inflection point is also seen in the FUV/NUV ratio (Fig. 4). The UV radiation and the IR are both due to the amount of dust along the line of sight and one might expect them all to be tightly correlated, especially at low column densities where the optical depth is low both in the UV and the IR. The excess radiation in the FUV is most likely due to Lyman band fluorescent emission from molecular hydrogen (Martin et al., 1990; Hurwitz, 1998). The inflection point represents the column density above which self-shielding protects the hydrogen molecules from photodissociation by interstellar UV photons and corresponds to a column density of log N(H) = 20.2. Canonically, this critical column density has been taken as log N(H)$>$20.6 (Savage et al., 1977; Franco & Cox, 1986; Reach et al., 1994) but has been found to be lower in high latitudes cirrus clouds (Gillmon & Shull, 2006; Planck Collaboration et al., 2011), where Gillmon et al. (2006) suggests that H2 formation is more efficient due to compression from dynamical processes at the disk-halo interface.

The intercepts observed represents the surface brightness at zero column density where we would expect no contribution from dust scattering. We will discuss its possible components below but note that they are in agreement with the spectroscopic data of Anderson et al. (1979); Tennyson et al. (1988); Seon et al. (2011) who found similar offsets after correcting for airglow lines in the spectrum.

3 Modelling and Results

We have assumed four independent sources for the UV radiation in this field:

Extragalactic background light (EBL). 2. 2.

Starlight scattered from interstellar dust. 3. 3.

Molecular hydrogen fluorescence. 4. 4.

Unexplained offset.

We will discuss each of these in the following sections.

3.1 EBL

Akshaya et al. (2018) tabulated the different components of the known EBL:

• •

60 – 81 photon units in FUV and 121 – 181 photon units in NUV from the integrated light from galaxies (Gardner et al., 2000; Xu et al., 2005; Voyer et al., 2011; Driver et al., 2016);

• •

16 – 30 photon units from QSOs (Madau, 1992);

• •

$<$20 photon units from the intergalactic medium (IGM) (Martin et al., 1991);

for a total of $114\pm 18$ photon units in the FUV and $194\pm 37$ photon units in the NUV. We have extincted the EBL using the known E(B - V) and accounted for the scattered EBL photons for different values of albedo and phase function asymmetry factor using Mattila (1976). This is shown in Fig. 5 (green line).

3.2 Dust-scattered Starlight

We have used the same Monte Carlo model for the dust-scattered starlight as did Akshaya et al. (2018). This model is fully described in Murthy (2016) and may be downloaded from Murthy (2015). The model uses the known positions of hot, UV stars from the Hipparcos star catalog (Perryman et al., 1997) with model spectra from Castelli & Kurucz (2003) to predict the number of photons at any position in the Galaxy. The amount and location of the interstellar dust was derived from Schlegel et al. (1998) reddening maps using a scale height of 125pc (Marshall et al., 2006). The stellar photons were scattered from the dust using the Henyey-Greenstein scattering function (Henyey & Greenstein, 1941) and assuming a range of values for albedo ($a$) and phase function asymmetry factor ($g$).

We derived the dust scattered light by subtracting the zero-offset (Table 1) and EBL with extinction from the total observed radiation. We then fit this remainder with the model predictions for different values of $a$ and $g$ for the entire data in NUV and for E(B-V)$<$0.027 in case of FUV. As there is lot of scatter in the model they were compared with the data at a resolution of 12*′*. We have only considered those bins with more than five data points to compare with the model. We find 1$\sigma$ confidence interval of $a=0.4\pm 0.1$ and $g=0.8\pm 0.1$ in the FUV and $a=0.4\pm 0.1$ and $g=0.5\pm 0.1$ in the NUV using the formulation given by Lampton et al. (1976). These are in agreement with the theoretical predictions of Draine (2003) and various observational estimates (Table 2). The models fit the data well in the NUV but fall short in the FUV for E(B - V) $>$ 0.027, where we believe that the fluorescent emission from molecular hydrogen begins to be important.

3.3 Molecular Hydrogen Fluorescence

We have plotted the contributions to the observed data from the EBL, the dust-scattered starlight and the offset in Fig. 5. The NUV data are fit well by the model as are the FUV data for E(B - V) $<$ 0.3, corresponding to log N(H) = 20.2, assuming N(H)/E(B-V) = 5.8$\times$1021 atoms cm*-2* mag*-1* (Bohlin et al., 1978). We have attributed this excess emission to fluorescence in the Lyman band (1400 – 1700 Å) of molecular hydrogen (Duley & Williams, 1980; Martin et al., 1990; Hurwitz, 1998). The level of this emission averages out to about 100 photon units, consistent with the surface brightness found by Jo et al. (2017). Using the formulation of Martin et al. (1990), we can derive the local gas density in the cloud from our observations assuming a constant value of R = $1\times 10^{-17}$ cm3 s*-1* for the formation rate of H2 (Fig. 6). Adding this emission to the other components gives the solid line in Fig. 5.

3.4 Offsets

It has long been noted (Henry, 1991; Hamden et al., 2013; Henry et al., 2015) that the diffuse background, particularly at high latitudes, cannot be fully explained by known sources (dust scattered radiation and EBL) without an additional offset at zero column density. Akshaya et al. (2018) found this offset to be $174\pm 18$ (FUV) and $337\pm 37$ (NUV) photon units in the North Galactic Pole and $127\pm 18$ (FUV) and $385\pm 37$ (NUV) photon units in the South Galactic Pole. In this work, we find that the offsets are $240\pm 18$ and $394\pm 37$ photon units in the FUV and NUV, respectively, after subtracting the known EBL. These offsets are slightly higher than those in the polar regions which may be due to uncertainties in the data.

Unfortunately, we cannot distinguish residual airglow from other diffuse sources in our study but other, spectroscopic, studies have found similar offsets at the NGP (Anderson et al., 1979; Tennyson et al., 1988; Seon et al., 2011). Brune et al. (1978) has measured the contribution from the 1356 Å OI line, the only significant contributor in the GALEX FUV band, to be less than 50 photon units. We are therefore confident that the offset is predominantly not due to airglow.

4 Conclusions

We have studied the latitude range between 70° and 80° as a continuation of our characterisation of the diffuse UV sky. This region has a greater amount of nebulosity and hence higher column densities than our earlier study at the Galactic Poles (Akshaya et al., 2018). We confirm the presence of an offset at zero column density (E(B -V) = 0) with a level of $240\pm 18$ photon units in the FUV and $394\pm 37$ photon units in the NUV, slightly higher than in the polar regions (Akshaya et al., 2018). Henry et al. (2015) has suggested that this is due to an unknown Galactic component of the DGL but it may be due to an unidentified EBL component. This may be resolved as we extend our study to lower Galactic latitudes where the dust component is enhanced and the extragalactic component is further extincted.

We find a component due to scattering from interstellar dust rising to a level of about 100 photon units in the FUV and 200 photon units in the NUV, about the same level as the EBL. We were able to set 1$\sigma$ limits on $a$ to be 0.4$\pm$0.1 in both FUV and NUV and the limits on $g$ were found to be $g=0.8\pm 0.1$ in FUV and $g=0.5\pm 0.1$ in NUV. We find an increase over the dust scattered emission in the FUV for column densities greater than E(B - V) = 0.027 mag (log N(H) = 20.2) which we interpret as molecular hydrogen fluorescent emission in the Lyman band.

Molecular hydrogen fluorescence emission was observed in the NGP by Akshaya et al. (2018) for column densities greater than log N(H)=20.2. We also observe this emission at the same column density for the region 70°$\leq$GLAT$\leq 80$°. Other studies have found a higher value (log N(H) $>$ 20.6) for the critical density at which molecular hydrogen can form through self-shielding from the interstellar UV photons (Savage et al., 1977; Franco & Cox, 1986; Reach et al., 1994) but Gillmon et al. (2006) suggest that the critical density is less at high Galactic latitudes.

Acknowledgements

Support for JM has come from proposal EMR/2016/001450 from DST/SERB. We have used the Gnu Data Language (http://gnudatalanguage.sourceforge.net/index.php) for the analysis of these data. The data presented in this paper were obtained from the Mikulski Archive for Space Telescopes (MAST). STScI is operated by the Association of Universities for research in Astronomy, Inc., under the NASA contract NAS5-26555. Support for MAST for non-HST data is provided by the NASA Office of Space Science via grant NNX09AF08G and by other grants and contracts.