Draft version June 25, 2021 Typeset using LATEX twocolumn style in AASTeX63
An Ultra-long Wavelength Sky Model with Absorption Effect
Yanping Cong,1, 2 Bin Yue∗,1 Yidong Xu∗,1 Qizhi Huang,1 Shifan Zuo,3, 1 and Xuelei Chen∗1, 2, 4
1National Astronomical Observatories, Chinese Academy of Sciences, 20A, Datun Road, Chaoyang District, Beijing 100101, China 2School of Astronomy and Space Science, University of Chinese Academy of Sciences, Beijing 100049, China 3Department of Astronomy, Tsinghua University, Beijing 100084, China 4Center of High Energy Physics, Peking University, Beijing 100871, China
ABSTRACT The radio sky at frequencies below 10 MHz is still largely unknown, this remains the last unex- ∼ plored part of the electromagnetic spectrum in astronomy. The upcoming space experiments aiming at such low frequencies (ultra-long wavelength or ultra-low frequency) would benefit from reasonable expectations of the sky brightness distribution at relevant frequencies. In this work, we develop a radio sky model that is valid down to 1 MHz. In addition to the discrete HII objects, we take into account ∼ the free-free absorption by thermal electrons in the Milky Way’s warm ionized medium (WIM). This absorption effect becomes obvious at . 10 MHz, and could make the global radio spectrum turn over at 3 MHz. Our sky map shows unique features at the ultra-long wavelengths, including a darker ∼ Galactic plane in contrast to the sky at higher frequencies, and the huge shadows of the spiral arms on the sky map. It would be a useful guidance for designing the future ultra-long wavelength observa- tions. Our Ultralong-wavelength Sky Model with Absorption (ULSA) model could be downloaded at DOI:10.5281/zenodo.4454153.
Keywords: Radio continuum emission; Interstellar emissions; Galaxy structure; HII regions; Extra- galactic radio sources
1. INTRODUCTION rate. These includes the Global Sky Model (GSM) (de The radio sky has been surveyed from 10 MHz up Oliveira-Costa et al. 2008), its improvements (Danny ∼ to THz. Based on data of surveys at the relevant C. Price 2016; Zheng et al. 2017; Sathyanarayana Rao ∼ frequencies, full-sky maps were produced at 408 MHz et al. 2017; Kim et al. 2018), and the Self-consistent (Haslam et al. 1982; Remazeilles et al. 2015), 1.42 GHz Sky Model (SSM) (Huang et al. 2019). They are very (Reich 1982; Reich & Reich 1986; Reich et al. 2001), useful in the design of new instruments, study of obser- and a number of higher frequencies from 22.8 GHz to vation strategies, and testing foreground removal/miti- 857 GHz by WMAP and Planck satellites, see Hinshaw gation methods (e.g., Shaw et al. 2014; Zuo et al. 2019). et al.(2009); Planck Collaboration et al.(2016a). In ad- However, the extrapolation becomes trickier at fre- dition, there are also surveys which cover a large portion quencies . 10 MHz, where the absorption of the inter- but not the full sky. In order to produce sky maps at stellar medium (ISM) becomes significant, and at the frequencies in which the full sky have not be surveyed, a same time there is a dearth of observation data. Ob- arXiv:2104.03170v2 [astro-ph.GA] 24 Jun 2021 number of sky models have been developed, based essen- servations below 10 MHz is severely hampered by the tially on interpolation or extrapolation of the available absorption and distortions of the ionosphere (Jester & data, and statistical modeling of the sky intensity dis- Falcke 2009). Until now there are only a few ground- based observations performed during 1950s 1960s tribution. Over much of the radio wave band, the sky ∼ − intensity has a nearly power law spectrum, making the (George et al. 2015; Cane & Whitham 1977), and a few interpolation/extrapolation relatively simple and accu- space-based low-resolution observations performed dur- ing 1960s 1970s (Alexander et al. 1969; Brown 1973; ∼ − Alexander & Novaco 1974; Alexander et al. 1975; Cane Corresponding author: Bin Yue; Yidong Xu; Xuelei Chen 1979). Below we shall call this wave band the ultra-long [email protected]; [email protected]; wavelength or ultra-low frequency band, though these [email protected] are not the name used by the radio engineers. The 2 ultra-long wavelength sky is still a largely unexplored absorption is probably the reason of this downturn in regime in the astronomical observations. To open this spectrum. However, it is also possible that the intrin- new window in the electromagnetic spectrum, a num- sic energy spectrum of the cosmic ray particles becomes ber of space missions, e.g. the lunar-orbit array for softer at low energy (Strong & Moskalenko 1998a). More Discovering the Sky at the Longest wavelength (DSL) references about early low frequency observations could (Chen et al. 2019), and the lunar-surface based FAR- be found in Cane(1979). SIDE (Burns et al. 2019) have been proposed, with the For dense electron clumps, such as those in discrete aim of obtaining high resolution1 full-sky maps at ultra- HII regions, the free-free absorption could make these long wavelengths. objects opaque even for radio signal at higher frequen- However, to design such low frequency experiments cies (e.g. Odegard 1986; Kassim 1988). This effect can and observe the sky at this largely unexplored part of be used to separate the synchrotron emissivity behind electromagnetic spectrum, one needs a reasonable esti- and in front of such regions. If there are sufficient num- mate of the sky brightness at these frequencies, and a ber of such HII regions with known distance, the 3D fair expectation of the structures, so that the required distribution of the synchrotron emissivity could be re- system gain and dynamic range can be derived. More constructed. This has been successfully put into prac- importantly, at present or in the near future, a lunar- tice, see e.g. Nord et al.(2006); Hindson et al.(2016); Su orbit or lunar-surface based array necessarily has a lim- et al.(2017, 2018); Polderman et al.(2019) and refer- ited number of antenna elements (though the orbital ences therein. Since the synchrotron radiation is pro- precession could improve the uvw coverage at the price duced by the cosmic ray electrons in magnetic field, of longer observation time), and the antennas are all-sky the cosmic ray electron distribution can be inferred if sensitive, these make the imaging process quite challeng- the magnetic field is known by other means (Polderman ing, especially when the sky brightness itself varies dra- et al. 2020). The degeneracy between the cosmic ray matically with frequency and direction, as we will show and the magnetic field could also be broken by combin- below. It is therefore crucial to have a reasonable full- ing with gamma-ray observations (Nord et al. 2006), as sky model for these unexplored frequencies, which could the Galactic gamma-ray is mainly from the collisions be- serve as a starting point for the deconvolution process, tween the cosmic ray particles and ISM particles (Nava for both end-to-end simulations in designing such exper- et al. 2017), and the ISM density profile could be con- iments and the upcoming real observations. structed from 21cm, CO, and dust thermal emission ob- At . 10 MHz, the existing sky models based on ex- servations (Ackermann et al. 2012). trapolation are grossly inadequate, because the free-free A sky model for the ultra-long wavelengths has to absorption and synchrotron self-absorption (Orlando & take into account the free-free absorption so as to Strong 2013; Ghisellini 2013) become very significant. give a reasonable prediction for the upcoming obser- In the Galactic diffuse ISM the former is much more vations (Reynolds 1990). In this work, we develop an important (Scheuer & Ryle 1953; Shain 1959; Hurley- observation-based sky model that is still valid below Walker et al. 2017). Free-free absorption is the inverse 10 MHz, by including the free-free absorption effect. ∼ process of bremsstrahlung (free-free) radiation, so it is We start with a description of our methods in Sec.2, proportional to the square of free electron density. The which incorporate a fitting result for the Galactic syn- absorption level depends not only on the smoothed dis- chrotron spectral index, a fitted Galactic synchrotron tribution, but also on the small-scale turbulence. emissivity, and an adopted electron density distribution Since as early as 1960s, it has been discovered that the that is responsible for the free-free absorption. The global radio background (the sum of Galactic and extra- main results are given in Sec.3, including the model galactic radiation) spectrum has a downturn at 3 5 with constant spectral index, the model with frequency- ∼ − MHz, by ground-based telescope (e.g. Ellis & Hamil- dependent spectral index and the one with direction- ton 1966), by radio-astronomy rocket (e.g. Alexander dependent spectral index. At last, we summarize our & Stone 1965), and by space satellites (e.g. Ariel II main conclusions in Sec.4. A brief description of our Smith 1965; RAE-1 Alexander et al. 1969; IMP-6 Brown software is presented in AppendixA. 1973; and RAE-2 Novaco & Brown 1978), and free-free
1 At such low frequencies, the resolution is limited by the scatter- 2. METHODS ing in the interstellar medium and interplanetary medium. This The free-free absorption becomes significant at low fre- angular scattering limit is roughly arc-minute level at 1 MHz and roughly scales with ∝ ν−2, see Jester & Falcke(2009). quencies. The optical depth can be written as the inte- gration of the absorption coefficient, κν , along the line 3 of sight, Z s 0.05 0 1 MHz τ(ν, s) = κν (ne,Te) ds , (1) 3 MHz 0 0.04 10 MHz where ne and Te are the electron density and tempera- 0 ture at s respectively. For constant electron tempera- 0.03 ) < z
ture, the following approximation holds (Condon & Ran- ( IGM som 2016), τ 0.02 −1.35 −2.1 −7 Te ν EM τν 3.28 10 , 0.01 ≈ × 104 K GHz pc cm−6 (2) 0.00 where the emission measure (EM) is the integration of 0 1 2 3 4 5 6 7 z electron density squared at all scales, Figure 1. The IGM mean optical depth for ν = 1, 3, 10 MHz EM Z n 2 ds = e . (3) to redshift z. pc cm−6 cm−3 pc The ISM inside the Milky Way, the circum-galactic We use the tanh model to describe the reionization medium (CGM), and the intergalactic medium (IGM) history for both HI and HeI, and HeII, i.e. can all absorb low frequency radio waves. We can make " # some simple estimates to assess the absorption on the HI,HeI 1 + nHe/nH yre y fe = 1 + tanh − , (6) various astrophysical scales. 2 ∆y Regarding the CGM, van de Voort et al.(2019) has in- vestigated the distribution of the hydrogen surrounding " # Milky Way-size galaxies from high resolution zoom-in yre y f HeII = 0.5 1 + tanh − , (7) cosmological hydrodynamic simulations with star for- e ∆y mation. They presented the column number density of the hydrogen as a function of r in their Figure 2. We and adopt a similar radial density profile, and assume that HI,HeI HeII the free electron density follows the same distribution, fe = fe + (nHe/nH)fe . (8) with a clumping factor of 3 and CGM electron tem- 5 ∼ 3/2 perature of 10 K. We find τCGM = 0.9, 0.09 and 0.007 Here y(z) = (1 + z) , yre = y(zre), ∆y = 3/2(1 + ∼ 1/2 for 1, 3 and 10 MHz respectively. If, however, there are zre) ∆z, and nHe/nH 0.083. For HI and HeI, we ≈ some cooler dense clumps in the CGM, the absorption assume zre = 7.68 and ∆z = 0.5, and for HeII, we adopt would be stronger. Thus, the CGM is nearly transparent zre = 3.5 and ∆z = 0.5 (Planck Collaboration et al. above 1 MHz, but becomes opaque below 1 MHz. 2020). We use the clumping factor model that fits reion- The EM of the IGM up to a redshift z is ization simulations presented in Iliev et al.(2007), but Z z renormalize it to be 3.0 at z = 5, i.e. κν0 0 2 0 0 EM = CIGM(z )¯ne (z )drp, (4) 0 κν 2 CIGM(z) = 6.8345 exp( 0.1822z + 0.003505z ). (9) 0 0 0 − where CIGM(z ) is the clumping factor at z , and drp = cdz0/[(1 + z0)H(z0)] is the proper distance element at z0. The IGM free-free absorption optical depths at frequen- 3 Alson ¯e(z) = fe(z)¯nH(1 + z) is the mean free electrons cies 1, 3 and 10 MHz are shown in Fig.1 as a function −7 −3 density at z, withn ¯H 1.89 10 cm the mean of z, it is not large even at 1 MHz. ≈ × hydrogen number density at present. Note that to use Next we consider the ISM, which is much denser than the EM as in Eq. (2), we add the factor κν0 /κν to correct the IGM and CGM. The observed sky brightness tem- for the redshift-dependence of the absorption for a given perature in the direction given by Galactic coordinates −2.1 observing frequency at ν. Since κν ν (Condon & (l, b) is the sum of the Galactic radiation TG(ν, l, b) and 0 0 ∝ Ransom 2016) and ν = ν(1 + z ), we have κν0 /κν = the extragalactic background TE(ν), (1 + z0)−2.1, hence T (ν, l, b) = TG(ν, l, b) + TE(ν) Z z cdz0 EM = (1 + z0)−2.1C (z0)¯n2(z0) . Z sG IGM e 0 0 −τ(ν,s) iso −τ(ν,sG) 0 (1 + z )H(z ) = (ν, R, Z, φ)e ds + TE (ν)e , 0 (5) (10) 4 where sG is the maximum distance to the edge of the 2011), minus the cosmic microwave background (CMB). Milky Way along any line of sight, is the three- This is in good agreement with the fit obtained from dimensional radiation emissivity in Galactic-centric the Long Wavelength Array (LWA) (Dowell & Taylor cylindrical coordinates (R, Z, φ), and τ(ν, s) is the free- 2018) below 100 MHz. Strictly speaking, “background” free absorption optical depth integrated up to a distance is the cumulative radiation from unresolved sources, so iso of s. Here we set sG = 50 kpc. TE is the extragalactic it depends on the sensitivity and angular resolution of radio background for which we assume to be isotropic in the telescope. For simplicity we directly use the back- the first approximation. With models of an extragalac- ground fit given by the ARCADE-2 experiment, since in iso tic radio background TE at ultra-long wavelengths, a our model we do not specify any particular instrument. three-dimensional emissivity (ν, R, Z, φ), and a three- In addition to the absorption by our Milky Way, the dimensional distribution of free electrons in the Milky extragalactic radio background would also be absorbed Way ne, we can derive the sky map at the required fre- by the ISM in the host galaxies. Moreover, the spectrum quencies from Eq. (10). In the following subsections, we of relativistic electrons could become softer at lower en- iso describe the modeling of TE , , and ne respectively. ergy, resulting in softer extragalactic radio background (Protheroe & Biermann 1996). In a recent estimate for 2.1. The isotropic extragalactic background the extragalactic radio background, Nit¸u et al.(2021) The extragalactic background is assumed to be largely found the extragalactic radio background deviates from isotropic, and a summary of some past results is given power-law below 1 MHz. In Sec. 3.3 we investigate ∼ in Guzm´anet al. 2011. Two methods are commonly the case that the extragalactic radio background is not used to extract it from the observed sky maps (Kogut a power-law form. et al. 2011; Seiffert et al. 2011). In the first approach the Galactic radiation is described by a plane-parallel 2.2. The Galactic free-free emission model, so the total radiation is The free-free emission also contributes to the low fre- iso quency sky brightness and should be subtracted from T (ν) = TE (ν) + TG(ν) csc b , (11) × | | the observed maps for the purpose of assessing Galactic iso where both TE and TG depend only on frequency. In synchrotron. A free-free template could be constructed the second approach, the Galactic radiation is derived by using the Hα emission line as a tracer. After cor- from a tracer, usually the Galactic [CII] emission, say recting the absorption and scattering, the intrinsic Hα 2 flux could be converted into free-free flux (Dickinson X iso p T (ν) = T (ν) + ai(ν) ICII, (12) et al. 2003; Lian et al. 2020). Moreover, it can also E,i × i=1 be constructed from the multi-frequency microwave ob- servations. In this work, we use the EM and Te maps where ICII is from the CORE/FIRAS measurements. constructed from Planck multi-frequency observations In practice, all coefficients are fitted from observed ra- (Planck Collaboration et al. 2016b)2 to derive the free- dio maps and these two methods give consistent re- free map, say sults (Kogut et al. 2011; Dowell & Taylor 2018). How- ever, whether the extracted isotropic background really " # √ originates from extragalactic sources is still being de- 3 −1.5 gff (ν) = ln exp 5.960 ln ν9 T4 + exp(1) bated. Seiffert et al.(2011) extracted a background − π which is larger than the integrated radio emission of ex- −1.5 −2 τff (ν) = 0.05468 Te ν9 EM gff ternal galaxies (Gervasi et al. 2008). The excess could 6 Tff (ν) = 10 Te[1 exp( τff )], (14) be from unknown extragalactic source populations, or − − could be unaccounted Galactic foreground (see also Sub- where ν is the frequency in units 109 Hz and T is the rahmanyan & Cowsik 2013). 9 4 electron temperature in units 104 K(Draine 2011). In Here we assume that the isotropic background is from the top panel of Fig.2, we show the 408 MHz Galac- extragalactic sources. We adopt the following simple tic free-free emission map constructed from the Planck model of isotropic extragalactic background, observations. The bottom panel shows the Galactic free- ν −2.58 free spectrum. Clearly, most of the free-free emission is T iso(ν) = 1.2 [K]. (13) E 1 GHz from the Galactic plane. This is the fit obtained by the ARCADE-2 balloon ex- periment, which observed the absolute brightness tem- 2 http://pla.esac.esa.int/pla/aio/product-action?MAP.MAP ID= perature of the sky between 3 and 90 GHz (Seiffert et al. COM CompMap freefree-commander 0256 R2.00.fits 5
T (408MHz) ff 2012; McKinley et al. 2018),
ν βG TG(ν) = TG(ν∗) , (15) ν∗
where ν∗ = 408 MHz. We determine βG by minimizing
2 X X 1 χ = 2 T (νj) νj i i 2 " βG # νj TE(νj) + TG,i(ν∗) + Tff,i(ν) Ti(νj) , × ν∗ − -5.0 -4.0 -3.0 -2.0 -1.0 0.0 1.0 2.0 log(T/K) (16)
where Ti(νj) is the observed brightness temperature of the i-th pixel in the sky-map at frequency νj, TG,i(ν∗) 103 is the Galactic radiation for the i-th pixel at ν∗. The summation is performed for all available pixels at all observed frequencies of our selected data. ] 102 K [
f Below several hundred MHz, there are observed sky f T maps, each of which covers a part of the sky. In deter- 101 mining the spectral index, in addition to the Haslam 408 MHz map, we use a selection of observations that have large sky coverages: the 45 MHz sky map in combina- 100 tion of two observations (Guzm´anet al. 2011), and the 10 1 100 101 102 3 [MHz] 35 MHz , 38 MHz, 40 MHz, 50 MHz, 60 MHz, 70 MHz, 74 MHz, and 80 MHz maps observed by the LWA (Dow- Figure 2. Top: the Galactic free-free emission map at 408 ell et al. 2017). The Guzman 45 MHz map covers 96% of MHz. Bottom: the spectrum of the mean Galactic free-free the sky, with resolutions of 4.6◦ 2.4◦ and 3.6◦ 3.6◦ for × × emission. the southern and northern sky respectively. It combines two surveys performed by devices of similar character- Comparing the above free-free map with observations istics. All of the LWA maps cover a total of about 80% at 408 MHz (see next subsection), we find that the free- of the sky, with resolutions ranging from 4.8◦ 4.5◦ to × free emission contributes 4% of the mean observed sky 2.1◦ 2.0◦. As the different maps have different angular ◦ × brightness. At b > 10 , the free-free fraction decreases resolutions, we smooth all the maps with a beam with | | ◦ to 1%. If we smooth both maps with a resolution of 10 , FWHM = 5◦ when deriving the power-law index. For then the r.m.s. contributed by the free-free emission is all the maps we subtract the free-free emission given in 15% of the total, while at b > 10◦ the r.m.s. fraction | | Sec. 2.2. Our fiducial model assumes that the spectrum decreases to 6%. This is consistent with Dickinson et al. index is a constant, but we will also investigate a model (2003). with a frequency-dependent spectral index in Sec. 3.3 and a model with direction dependence in the spectral index in Sec. 3.4. 2.3. The absorption-free sky map Using the above method, we derive a spectral index βG = 2.51 for the Galactic synchrotron radiation. This − We first derive the synchrotron spectral index from is actually quite similar to the spectral index of the ex- the observed sky maps at frequencies where absorption tragalactic background, βE = 2.58. In Fig.3 we plot − is assumed to be negligible. At low frequencies, there the 1 MHz sky map that is extrapolated from free-free are full-sky maps at 408 MHz (Haslam et al. 1974, 1981, 1982; Remazeilles et al. 2015) and 1.42 GHz (Reich 1982; Reich & Reich 1986; Reich et al. 2001). Here we use the 3 Using Eq. (2) and electron model in Sec. 2.5, we check that at free-free subtracted Haslam 408 MHz map as the base, this frequency, for most of the sky regions the free-free absorption is negligible, except at a very thin plane with b ∼ 0◦ and near assuming that in the absence of absorption, the bright- the Galactic center where the optical depth could be up to ∼ 1 ness temperature of the Galactic synchrotron radiation , so it is safe to include this low frequency observations when follows a power-law form (Platania et al. 1998; Kogut constructing the absorption-free map. 6
Table 1. Fitted Galactic emissivity model parameters
60° A 43.10 K kpc−1
30° R0 3.41 kpc α 0.46
150° 120° 90° 60° 30° 0° 330° 300° 270° 240° 210° 180° Z0 1.12 kpc γ 1.23 30°
60°