Issue |
A&A
Volume 663, July 2022
|
|
---|---|---|
Article Number | A120 | |
Number of page(s) | 16 | |
Section | Cosmology (including clusters of galaxies) | |
DOI | https://doi.org/10.1051/0004-6361/202243238 | |
Published online | 19 July 2022 |
The first molecules in the intergalactic medium and halos of the Dark Ages and Cosmic Dawn
1
Ivan Franko National University of Lviv, Kyryla i Methodia str., 8 Lviv, 79005, Ukraine
e-mail: bohdan.novosyadlyj@lnu.edu.ua
2
Jilin University, Qianjin Street 2699, Changchun 130012, PR China
3
Institute of Radio Astronomy of NASU, 4 Mystetstv str., 61002 Kharkiv, Ukraine
e-mail: shulga@rian.kharkov.ua
Received:
1
February
2022
Accepted:
11
April
2022
We study the formation and destruction of the first molecules at the epochs of the Dark Ages and Cosmic Dawn to evaluate the luminosity of the protogalaxy clumps (halos) in the molecular lines. The cosmological recombination is described using the RecFast model of an effective three-level atom, while the chemistry of the molecules is examined using the relevant basic kinetic equations. We then studied the effect of collisional and radiative excitation of molecules on the intensity of molecular emission in both warm and hot halos. Using the Planck data on the reionization of the intergalactic medium at z ∼ 6−8, we evaluated the upper limits of the light energy density for four models of thermal light from the first sources that appeared in the Cosmic Dawn epoch. Assuming that in the halos, the light energy density may essentially be even higher, we estimated the impact of the light from the first sources (the first light) on the formation and destruction of the first molecules in them as well as between them. We show that the molecules H2 and HD are destroyed by photodissociation processes shortly before the full reionization in the inter-halo medium, in the medium of both types of halos and for all models of the first light. At the same time, the number density of helium hydride ions, HeH+, shows essentially more complicated dependences on the kinetic temperature of halos and the models of the first light. These features characterizing the abundance of molecules also determine the intensity of the halos luminescence during their evolution. Furthermore, we calculated the evolution of the brightness temperature of the individual halo in the rotational lines of H2, HD and HeH+ molecules relative to the temperature of the cosmic microwave background at redshifts corresponding to the Dark Ages and Cosmic Dawn epochs. It does not exceed the microkelvin, but its detection may be an important source of information about the physical processes taking place at the beginning of the formation of the first stars and galaxies at the epochs of the Dark Ages and Cosmic Dawn.
Key words: cosmology: dark ages / reionization / first stars
© B. Novosyadlyj et al. 2022
Open 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 upcoming decade promises to be a time of significant breakthroughs in the observational investigation exploring the epochs of the Dark Ages and Cosmic Dawn. These studies are expected to fill the gap in our understanding of how the first stars and galaxies formed. A few teams of cosmologists and observers are competing to detect the signal in the hyperfine 21 cm line of atomic hydrogen that is redshifted according to these epochs. The first detection of the absorption in this line at z = 15−20 was announced by the EDGES1 team a few years ago (Bowman et al. 2018; Hills et al. 2018). Other sources of important information about the formation of the first stars and galaxies are the first molecules, which absorb the cosmic microwave background radiation (CMBR) in their ro-vibrational transitions in cold baryonic matter clouds, while emitting additional radiation in hot ones or resonantly scattering it in moving clouds (Dubrovich 1977; Maoli et al. 1994, 1996; Kamaya & Silk 2002, 2003; Ripamonti et al. 2002; Omukai & Kitayama 2003; Mizusawa et al. 2004, 2005; Dubrovich et al. 2008; Liu et al. 2019; Novosyadlyj et al. 2020; Maio et al. 2022). Unfortunately, predictions of the amplitudes of such distortions of CMBR are at the level of microkelvins and too low for their detection by current telescopes. The upper observational limits obtained over the last few decades for such signals are at the level of a few microkelvins (de Bernardis et al. 1993; Gosachinskij et al. 2002; Persson et al. 2010). The detectability of primordial pre-galactic clouds in the molecular lines has been discussed by some authors over the two past decades (Kamaya & Silk 2002, 2003; Ripamonti et al. 2002; Omukai & Kitayama 2003; Mizusawa et al. 2004, 2005; Liu et al. 2019; Novosyadlyj et al. 2020; Maio et al. 2022) and a marginal possibility has been predicted for modern and projected telescopes such as ALMA2, JWST3, SAFIR Observary4, SPICA5, and OST6.
In ours previous papers (Novosyadlyj et al. 2017, 2018, 2020; Kulinich et al. 2020) we studied the formation and destruction of the first molecules, H2, HD, and HeH+, in the halos of the Dark Ages and their luminescence in the rotational lines of these molecules. Here, we elucidate the matter of how the radiation of the first sources, arising at Cosmic Dawn, influences the abundances of the first molecules and their luminescence in the same rotational lines. For simplicity, we suppose the thermal energy distribution of the first sources and consider the different evolution of its spectral density and radiation temperature at the moment of cosmological reionization, zrei. The latter we assume to be between 6 ≤ z ≤ 8, according to the observational data coming from the measurements of the polarization of CMBR at low spherical harmonics (Planck Collaboration I 2020; Planck Collaboration VI 2020), Gunn-Peterson trough in the spectra of high-redshift quasars (Bouwens et al. 2015; Banados et al. 2018; Davies et al. 2018), and the Lyman Break galaxies (Mason et al. 2018).
All computations in the paper were performed for consistent values of the main parameters of the cosmological model following the final data release of the Planck Space Observatory (Planck Collaboration VI 2020): the Hubble constant H0 = 67.4 km s−1 Mpc−1, the mean density of baryonic matter in the units of critical one Ωb = 0.0493, the mean density of dark matter Ωdm = 0.2657, the mean density of dark energy Ωde = 0.685, its equation of state parameter wde = −1.03, the spectral index of the scalar mode of cosmological perturbations ns = 0.965 and the matter fluctuation amplitude σ8 = 0.811, and the current temperature of cosmic microwave radiation T0 = 2.7255 K. We use also the primordial helium abundance Yp = 0.2446 (Peimbert et al. 2016) and deuterium fraction yDp = 2.527 × 10−5 (Cooke et al. 2018), which are in good agreement with the posterior means from Planck Collaboration VI (2020).
The outline of this paper is as follows. In Sect. 2, we present the network of chemical reactions leading to the formation and destruction of molecules and the recombination and photoionization of atoms as well as the corresponding kinetics equations and energy balance. In Sect. 3, we describe the models of the energy distribution for the first sources of light (the first light) and estimate the value of the dilution coefficient for each from them. In Sect. 4, we compute the kinetics of chemical reactions in the inter-halo medium and analyze the influence of the first light on the number density of the first molecules H2, HD and HeH+. The intra-halo chemistry and luminosity of halos in the ro-vibrational molecular lines are analyzed in Sect. 5. To do this, together with the equations of the kinetics of chemical reactions and the kinetics of excitations and de-excitations of the ro-vibrational energy levels of the first molecules, we solve the equations for the evolution of the density and temperature of matter. The top-hat halos with hot and warm baryonic gas inside and with masses that are typical for ΛCDM/ΛWDM models of the Universe are considered. As a result, we obtain the evolution of the number density of molecules H2, HD, and HeH+ at different stages of halo formation up to the reionization epoch for different models of the first light. Our conclusions are given in Sect. 6. Useful approximations for the rates of photorecombination and photoionization for HI and DI are presented in Appendix A, while the rates of photodissociation for molecules H2 and HD are presented in Appendix B.
2. Chemical network
We analyze the formation and destruction of the most abundant molecules in the inter-proto-galaxy medium of the Dark Ages and Cosmic Dawn epochs: H2, HD, and HeH+. The minimal set of reactions important for epochs starting from cosmological recombination, Dark Ages, Cosmic Dawn, and reionization is presented in Table 1. The rates and fitting formulas for them are taken from the sources cited in the last column. Most of fitting formulas are collected in Tables 1–3 of Abel et al. (1997), Tables 1–3 of Galli & Palla (1998), and Table B.1 of Schleicher et al. (2008). We also keep here the notations of reactions (in parentheses) introduced there. In comparison with a minimal set of reactions important for the Dark Ages epoch (Galli & Palla 1998; Novosyadlyj et al. 2018), we have added the important reactions for the Cosmic Dawn epoch, when the first light would ionize atoms and molecules and then photodissociate the molecules. These are: H6, H11, H12, H18, H23, D25, He10, He12, and He15. Since the ionization of hydrogen is important for the formation and destruction of molecules, we propose the more accurate analytical approximation of effective cross-sections of hydrogen photoionization from the base level (case A) and the first metastable 2s one (case B) in Appendix A. The direct reaction of photodissociation H23, as well as that of D25, becomes crucial for the destruction of hydrogen and hydrogen deuteride molecules by first light before reionization (Glover & Jappsen 2007; Schleicher et al. 2008). We propose the new approximation of the rate of photodissociation of H2 (HD) using recent and accurate computations of effective cross-section for this purpose from Heays et al. (2017), presented in Appendix B.
Main chemical reactions and references for their rates.
We computed the evolution of fractions, xi, of neutral atoms, molecules, and ions from the cosmological recombination up to hydrogen reionization at zrei. The index “i” notes any atom, molecule, or ion presented in Table 1. The kinetic equations for chemical reactions in the general form are as follows (Puy et al. 1993; Galli & Palla 1998; Vonlanthen et al. 2009; Novosyadlyj et al. 2018):
where z is redshift; is the rate of collision reaction with reactants, m and n, which leads to the formation of the atom, molecule, ion, i; is the rate of radiative reaction for reactant, m, which leads to the formation of an atom, molecule, ion, i; Cin is the rate of destruction of i-reactant by collision with an n-reactant; and αi is the photorecombination rate of the reactant, i. The total fraction of the main chemical elements in the unit number density of hydrogen, denoted as fm, namely, fHe = nHe/nH for a reactant, m, containing helium; fD = nD/nH for a reactant, m, containing deuterium; and fH ≡ 1 for reactant, m, containing hydrogen only. For chemical species containing only hydrogen, the fraction, m, is xm = nm/nH, where nm is the number density of species, m, and nH is the total number density of hydrogen; for species containing deuterium and helium, xm = nm/nD and xm = nm/nHe, respectively, where nD and nHe are the total number densities of deuterium and helium.
Equation (1) for the reactions presented in Table 1 together with equations of cosmological recombinations for hydrogen, deuterium, and helium (Seager et al. 1999, 2000; Novosyadlyj et al. 2017) compose a system of kinetic equations which we solve numerically together with the equations of expansion of the Universe and energy balance for the baryonic component:
where σT is the Thomson scattering cross-section, ar is the radiative constant, and me is the mass of an electron. We set the initial conditions for these reactions at the early epoch before the cosmological recombination when all ingredients were completely ionized and a Saha approximation was applicable (see details in Appendix A in Novosyadlyj et al. 2017). The publicly available codes RecFast7 and DDRIV18 have been used in the general code CDhalo.f, which was designed for integrating the system of equations (1) in the expanding Universe over cosmological recombination, Dark Ages, and Cosmic Dawn, when the first light becomes important for the ionization and dissociation of atoms and molecules. The last Eq. (3) is used at z ≤ 850, at higher zTb = Tr = T0(1 + z). It takes into account the main mechanisms of cooling or heating of diluted gas in the Dark Ages and Cosmic Dawn: adiabatic cooling and heating by CMBR at the early epoch and by the first light in the late ones.
3. Models of the first light
In the epoch of Cosmic Dawn, the radiation of the first sources (warm or hot halos, the first dwarf galaxies with Pop III stars, and shock waves), referred to here as the first light, starts to impact the kinetics of the chemical reactions, exciting the ro-vibrational levels of molecules; and at the end of Cosmic Dawn, it starts to ionize the hydrogen, marking the beginning of the epoch of reionization. To indicate the moment of reionization, we assume that this happens when the fraction of ionized hydrogen reaches a value of 0.5. According to the Planck 2018 results on the CMBR anisotropy, this takes place at 6 ≤ z ≤ 8 (see Fig. 45 in Planck Collaboration VI 2020). Other observations, such as the Gunn-Peterson trough in the spectra of high-redshift quasars (Bouwens et al. 2015; Banados et al. 2018; Davies et al. 2018) and Lyman Break galaxies (Mason et al. 2018), support this value. Thus, here we analyze models of the first light that satisfy the following condition: xHI = xHII = 0.5 at 6 ≤ z ≤ 8.
To analyze the allowable levels of illumination in the inter-proto-galaxy medium of the Cosmic Dawn epoch, we assume that sources of the first light are thermal. We consider here the thermal models of the first light described by the Planck function with temperature, Tfl, which is some smooth function of redshift and dilution coefficient αfl: Ifl(ν) = αflBν(Tfl).
So, in each moment of the Dark Ages and Cosmic Dawn, the total spectral energy density of radiation is the sum of the CMBR energy density and the first light one:
Since the rates of ionization or dissociation of atoms and molecules are computed by integration over the energy of product of cross-section and intensity per photon energy (see Eq. (A1)) in Appendix), they can be presented as the following sum:
First, we estimated which of the cases of photoionization-recombination, A or B, is realized in the Dark Ages and Cosmic Dawn. For this purpose, we compared the characteristic lifetime of the hydrogen atom in the excited states, 2s and 2p (, ), with the time intervals between following photo recombinations to the ground state and excited levels (, ) and photoionization from them by CMBR, and the first light photons (, ). The collision excitation, ionization, and recombination of atomic hydrogen are insignificant in these epochs (Peebles 1968; Matsuda et al. 1971; Grachev & Dubrovich 1991; Abel et al. 1997) and we omitted them in this estimation. Since the Lyα line is narrow, the medium is rarefied and expanding, while the Lyα photons escape from the processes of scattering due to the Doppler shift of the frequency (ibid). Thus, the comparison of characteristic times is reasonable. The results are shown in Fig. 1. The characteristic times of photoionization by CMBR are shown as black lines at z > 200, and those by the first light are shown by color lines at z < 30.
Fig. 1. Characteristic lifetime of the hydrogen atom in the excited states 2s and 2p (, ), time intervals between following recombinations to the base and excited levels (, ), and the resulting photoionization (, ) by CMBR (at z > 200) as well as the first light photons (z < 30). The computations are made for the cosmological model with Planck 2018 parameters and thermal models of the first light. |
Case B is realized when the characteristic time of photoionization is lower or close to the lifetime of hydrogen atoms on the metastable level of 2s. We can see that this condition is satisfied at z ≳ 1000 when the temperature of the CMBR was ≳3000 K. The energy density of the first light that appeared at Cosmic Dawn, increases but does not achieve the value when the population of the metastable level becomes important. Therefore, all computations here and below are made for case A of the recombination-photoionization of atomic hydrogen. We supported this conclusion by computations of the populations for hydrogen levels using the publicly available code Cloudy (Ferland et al. 2017): the values of hydrogen ionization fractions obtained by Cloudy and our code agree well.
We consider here the simple evolution model of the radiation temperature of the first light:
with such values of parameters: (a) T* = 5000 K, n = 2 (Ia), (b) T* = 10 000 K, n = 2, (Ib), (c) T* = 5000 K, n = 1 (Ic), and (d) T* = 10 000 K, n = 1 (Id). So, at the moment of reionization, zrei, the radiation temperature of the first light is 5000 K or 10 000 K. For the dilution factor, αfl, we set the constant and estimate its value using the condition of reionization. The condition xHI = xHII = 0.5 at 6 ≤ zrei ≤ 8, which we suppose here, gives 3.3 × 10−10 ≤ αfl ≤ 4.85 × 10−10 for model Ia, 6.0 × 10−17 ≤ αfl ≤ 8.6 × 10−17 for model Ib, 1.9 × 10−10 ≤ αfl ≤ 3.0 × 10−10 for Ic, and 4.5 × 10−17 ≤ αfl ≤ 6.3 × 10−17 for Id. So, for the same T* the range of values of αfl is narrow for such a range of redshift of reionization, but the value strongly depends on T*.
The spectral energy densities for all these models are shown in Fig. 2. The potential of ionization of hydrogen atoms is in the Wien range of the Planck distribution. Thus, the higher the radiation temperature at the moment when xHII ≈ 0.5, the lower the dilution coefficient should be. For comparison, we present in this figure the spectral energy density of CMBR at the same redshifts as well as at the current epoch (z = 0). The observable cosmic infrared background (CIB; Kogut et al. 2019) is shown as well. Such a comparison shows that at frequencies below 2.5 × 104 GHz (λ > 0.1 μm), the intensity of such a first light is essentially lower than CMBR and CIB. However, it can excite the ro-vibrational levels of the most abundant molecules to produce the spectral distortions of the CMBR in the shortwave range.
Fig. 2. Spectral energy density of the first light and CMBR at z = 8 (solid lines) and at z = 6 (dashed lines) when the ionization of hydrogen is ≈0.5 for models Ia and Ib. The vertical dotted lines show the thresholds of dissociation of molecules H−, HeH+, H, H2, and HD (overlapped with H2) and the ionization of HI, correspondingly. The dotted lines show these spectral energy densities at the current epoch. The spectral energy density of cosmic infrared background (dotted green line) is shown for comparison. In the models Ic and Id, the corresponding lines are lower by the ratios of dilution coefficients. |
In the left panel of Fig. 3, we present the evolution of the hydrogen ionized fraction xHII for the models of the first light Ia-Id. The helium atoms are represented by the neutral fraction. The values of dilution coefficients, αr, are adjusted so that xHII ≈ 0.5 on z = 8 and 6. We can see that the data in this figure well match the data in Fig. 45 of Planck Collaboration VI (2020). It means that these models of the first light are in agreement with the observational data on the optical depth of reionization.
Fig. 3. Reionization of hydrogen by the first light at the end of the Cosmic Dawn epoch. Left: evolution of the ionized hydrogen fraction xHII for the first light models Ia-Id, which provide xHII ≈ 0.5 at 6 ≤ z ≤ 8. Right: ranges of dilution coefficient against the temperature of thermal radiation for two evolution models of the first light. |
In the right panel of Fig. 3, we present the ranges of dilution coefficient against the temperature of thermal radiation for two evolution models of the first light, which provide xHII ≈ 0.5 at zrei = 8 (solid line) and zrei = 6 (dashed line).
4. Molecules in the inter-proto-galactic medium
The evolution of the number density H, D, and He and the temperature of CMBR and baryonic matter through the Dark Ages and Cosmic Dawn up to reionization are shown in Fig. 4 (top panel and bottom one correspondingly). In the bottom panel, we also show the evolution of radiation temperature of the first light. The correspondence of the lines of the radiation temperature to the models of the first light is the same as in Fig. 2. The heating of the baryonic gas by the first light is noticeable. We can see that the dilution factor is important in the heating of gas by radiation: the thermal radiation with lower Tfl but higher αfl (model Ia and Ic) heats more than the thermal radiation with higher Tfl but lower αfl (models Ic and Id). The temperature of baryonic matter is approximately the same (∼2000 K) in the cases of models Ia and Ic at the moment of reionization z = 8 and 6, as it is shown in Fig. 4 (solid red lines of Tb). At the same time, the temperature of the first light is 5000 K. For the model of the first light with radiation temperature Tfl = 10 000 K at the moment of reionization, the dilution factor is ∼10−17 and the baryonic matter is heated by radiation only by a few degrees (solid blue lines of Tb).
Fig. 4. Evolution of the number density of H, D and He (top) and temperature of CMBR and baryonic matter through the Dark Ages to Cosmic Dawn (bottom). |
We computed the evolution of the number densities of atoms and molecules, ni, starting from the early epoch before the cosmological recombination (z = 10 000) up to the beginning of the epoch of reionization (z = 6) for different models of the first light. In Figs. 5 and 6, we present the results for the Ia – Id models of the first light. The black lines show the evolution of the number densities of atoms, ions, neutral and ionized molecules without the first light, , which well agree with results by Galli & Palla (1998), Schleicher et al. (2008); and Vonlanthen et al. (2009). Recent studies by Maio et al. (2022) of H2 gas evolution over cosmic time (2 ≤ z ≤ 17) in the cosmological simulation ColdSIM do not contradict the semi-analytical ones presented here. The fact that the H2 molecules at z ≤ 8−9 do not disappear due to photodissociation can be explained by the effect of self-shielding of dense or extended regions of molecular hydrogen, which is taken into account in those studies. The lines of the same atoms and molecules in the epochs preceding the Cosmic Dawn are superimposed and black. The color lines show the influence of the first light on the number density of the first molecules.
We can see that negative hydrogen ions H− are most sensitive to the first light. This is expected since its threshold of dissociation is low, namely, 0.75 eV. The photodissociation reaction H− + γ → H + e− effectively decreases the number density of H− up to the moment when the ionization of hydrogen starts to increase. Increasing the number density of free electrons makes raises the importance of the reaction H + e− → H− + γ, which leads to the growth in the number density of H−. This is not for long, however, as the rapid increase of the number density of photons subsequently makes the process of photodissociation dominant.
The evolution of the number densities of the molecular ions and HeH+ is similar for the high energy density of the first light because their photodissociation thresholds (2.77 eV and 1.85 eV accordingly) are essentially lower than the photoionization potential of atomic hydrogen (13.6 eV). In both cases, an increase in the value of nHII causes a brief increase in the number density of these molecular ions, as well as a subsequent decrease due to the growth of the photodissociation. For the low energy density of the first light (bottom panels in the right columns of Figs. 5 and 6), the number density of hydride helium ions HeH+ increases when the ionization of atomic hydrogen increases because the reaction He + H+ → HeH+ + γ prevails over the reaction HeH+ + γ → He + H+.
Fig. 5. Evolution of the number densities of atoms, ions, neutral, and ionized molecules ni(z) from cosmological recombination through the Dark Ages and Cosmic Dawn up to reionization for the models of the first light Ia and Ib: the red lines for zrei = 6, the blue lines for zrei = 8. The black lines are the number density of species for the case without the first light, . |
Fig. 6. Evolution of the number densities of atoms, ions, and neutral and ionized molecules ni(z) from cosmological recombination through the Dark Ages and Cosmic Dawn up to reionization for the models of the first light Ic and Id: the red lines for zrei = 6, the blue lines for zrei = 8. The black lines are the number density of species for the case without the first light, . |
The number density of the most abundant molecules H2 (threshold of dissociation is 4.53 eV) drastically decreases when the number density of photons with energy above the threshold becomes larger than a specific value when the rate of direct dissociation of H2 becomes ∼10−18 − 10−19 cm3 s−1. Thus, the spectral energy distribution of the first light and its time dependence crucially impact the formation and dissociation of molecules during the Cosmic Dawn. To demonstrate this impact more clearly, we present in Fig. 7 the ratios of the number densities of molecules in the models with the first light, ni(z) to the number density of molecules without it, .
Fig. 7. Impact of the first light on the number densities of molecules at Cosmic Dawn for the different evolution of the energy density of thermal radiation of the first sources (models Ia-Id) (lines for H2 and HD superimposed). |
In this study, we supposed the primordial chemical content of hydrogen, helium, and deuterium. Taking into account the PopII stellar nucleosynthesis, metal enrichment will change the relationship between them, which will somewhat affect the number densities of the first molecules in the halos. Obviously, these changes will not be decisive at z ≥ 10.
5. Molecules in the intra-proto-galactic medium and their luminescence
We went on to repeat the same computations for the intra-proto-galactic medium. We supposed that the proto-galaxies are spherical top-hat halos in the multicomponent Universe, which are filled with dynamical dark energy, dark matter, baryon matter, and two “sorts” of thermal radiation, namely, the CMBR and the first light, which appeared in the Cosmic Dawn epoch. Since the time dependence of the density and temperature of the baryonic matter in the halos at stages before, during, and after virialization differ from the background ones, the number densities of the first molecules are different too. Moreover, the warm and hot halos can slightly shine at the cosmic microwave background (Novosyadlyj et al. 2020; Kulinich et al. 2020). To analyze the impact of the first light on the abundance of the first molecules and their luminescence in the halos, we use the theory of cosmological perturbations and the halo model to produce accurate computations of physical conditions and integrations of the corresponding kinetics equations.
5.1. Dynamical model of halo
We computed the physical conditions and chemistry of the halo by modeling the evolution of a single spherical perturbation in the four-component Universe (cold dark matter, baryon matter, dark energy, and thermal relict radiation) starting from the linear stage at the early epoch, long before the cosmological recombination, through the quasi-linear stage, turnaround point and collapse up to virialized state (Novosyadlyj et al. 2016, 2018). We assume that mass of a halo Mh ∼ 106 M⊙, which is the typical one in the cosmological model with cold dark matter and dark energy. The initial amplitudes of adiabatic perturbation in the conformal Newtonian gauge for this scale are taken as follows: for thermal radiation δtr = 1.5 × 10−4, for the matter components δcdm = δb = 1.125 × 10−4, and for phantom dark energy δde = −3.376 × 10−6 at zinit = 108, that correspond to the root mean square values at this scale in the cosmological model with parameters presented in the Introduction. At the redshift zinit, such a perturbation is the superhorizon one, which enters the particle horizon at zeh ≈ 2.5 × 106. The integration of the equation system that follows from those of Einstein and energy conservation (Novosyadlyj et al. 2016, 2018) gives the moment of virialization zv = 23.8 when the overdensity reaches the value Δv = 178 (top panel of Fig. 8).
Fig. 8. Formation of a spherical halo from a specific cosmological perturbation. Top panel: evolution of number densities of H, D, and He in the top-hat halo with k = 80 Mpc (Mh = 1.26 × 10−6 M⊙) and initial amplitude of adiabatic scalar perturbation δr = 1.5 × 10−4 at z = 108). Bottom panel: evolution of the temperature of CMBR (green line) and kinetic temperature of the top-hat halo with adiabatic heating (dark line) and virial temperature (8; red line). |
The halo mass Mh [M⊙], the comoving wave number of linear (seed) perturbation k [Mpc−1] and the comoving radius rh [kpc] after virialization are related as:
where and h ≡ H0/100 km s−1 Mpc−1. In the cosmological model with Ωm = 0.315 and h = 0.674 for k = 80 Mpc−1, we obtain the value Mh ≈ 1.26 × 106 M⊙ and rh ≈ 0.15 kpc. The angular size of such halo at zv ÷ zrei is in the range of 0.07 ÷ 0.03 arcsec. We suppose that a halo is spherical and homogeneous (top-hat) with values of matter density and kinetic temperature unchanged after virialization (Fig. 8). The total number density of baryon particles is unchanged too and equals 0.68 cm−3.
The kinetic temperature of baryon gas in the halo during the epochs of the Dark Ages and Cosmic Dawn is determined by its interaction with thermal relic radiation and the first light and the adiabatic heating over compression. The bremsstrahlung, molecular cooling and other cooling mechanisms are not effective due to the low gas density (top panels of Figs. 4 and 8). So, the evolution of baryon gas temperature Tb at 6 ≤ z ≤ 850 can be obtained by the integration of the energy balance equation for the intra-halo baryon matter:
The kinetic temperature of the gas in such halo after virialization is about 280 K and does not depend on halo mass. The halo with such temperature we call the “warm” halo.
The gas in the halo can also be heated by shocks in the processes of violent relaxation to the virial temperature (Barkana & Loeb 2001; Bromm & Yoshida 2011):
where μH is the mass per H atom. In our case, the temperature of such halo is 2280 K and we call it the “hot” halo.
In the bottom panel of Fig. 8, the temperature of the gas in the virialized halo is shown for two cases: when it is heated by adiabatic compression only (dark solid line) and when it is heated in the processes of violent relaxation to virial temperature (8; red solid line). In the first case, the final temperature is a result of the integration of Eq. (7), in the second one, we set it equal to by hand at z ≥ zv.
5.2. Number density of molecules
In the halo model described above, the number density of baryon (and lepton) particles is unchanged after virialization, while the number densities of molecules, protons (np) and electrons (ne) are changed via recombination and photoionization processes of atoms and molecules as well as the formation and destruction reactions of molecules.
In Fig. 9, we show how the number densities of molecules H2, HD, and HeH+ are changed in the warm and hot halos in comparison with the background ones (inter-proto-galaxy medium) without the first light. We can see that increasing of density and temperature of baryon gas leads to the increasing of the number density of molecules H2 and HD even after virialization essentially more than Δv times. Moreover, in the hot halos, the number density of these molecules is higher than in the warm ones. The number density of HeH+ molecular ions behaves differently. In the warm halo, the nHeH+ decreases so that the ratio is almost constant. In the hot halo, the number density of these molecules sharply decreases since destruction reactions He11 and He12 (Table 1) dominate.
Fig. 9. Ratio of the number densities of molecules in the warm and hot halo () to the number density of molecules at the cosmological background () without the first light. |
To estimate the impact of the first light on the chemistry of the intra-halo medium, it is necessary to note that the spectral density of radiation in the halo is larger than the background one and this excess increases over time quickly during the Cosmic Dawn epoch. We analyzed the impact of the first light with different values of dilution coefficients on the halo chemistry: (i) αfl, estimated in Sect. 3; (ii) 103αfl; and (iii) 106αfl. We computed the ratios of the number densities of molecules in the halo with the first light to the number densities of molecules in the halo without it. The results for the first light models Ia and Ib are presented in Fig. 10. In the top panels the models of the first light in the halo are the same, which reionize the inter-halo medium at z = 6 (dotted lines) and z = 8 (solid lines). In the left panel, the radiation temperature of the first light at the moment of reionization is 5000 K, with dilution coefficients αIa6 = 3.3 × 10−10 and αIa8 = 4.85 × 10−10 for zrei = 6 and 8, correspondingly. In the right panel, the radiation temperature of the first light at the moment of reionization is 10 000 K, with dilution coefficients αIb6 = 6.0 × 10−17 and αIb8 = 8.6 × 10−17 for zrei = 6 and 8, correspondingly. In the bottom panels the energy densities of the first light in the halos are 103 (dotted line) and 106 (solid lines) times higher than the energy densities in the corresponding models in the top panels shown by solid lines. In all the panels, the lines for H2 and HD are superimposed. We can see that these molecules are quite destroyed before reionization. The larger is the spectral density of the first light, the earlier they disappear because of photodissociation.
Fig. 10. Impact of the first light on the number densities of molecules in the halos of Cosmic Dawn. In the top row, the models of the first light are the same as in background ones (left – Ia, right – Ib), which reionize the medium at z = 6 (dotted lines) and z = 8 (solid lines). In the bottom row, the energy density of the first light in the halos are 103 (dotted line) and 106 (solid lines) times higher than energy density in the corresponding models in the top panels shown by solid lines. In all the panels, the lines for H2 and HD are superimposed. |
The number density of HeH+ molecular ions in the halo evolves quite differently. It starts to change early over the low threshold of dissociation and changes nonmonotonically over the competition of formation and destruction reactions (He8, He10, He11, He12, He14, and He15 in Table 1) in the medium with the growth of nHII, ne, and nγ. Moreover, in the models of the first light with α ≤ 10−14 and T* ∼ 10 000 K, the number density of these molecules increases by 102 − 103 times on zrei, when other molecules are destroyed in their entirety. It gives us the chance to observe the epoch of reionization in the lines of helium hydride ions.
In order not to overload the paper we do not present the figures for models of the first light Ic and Id, in which the spectral energy distributions on zrei are the same as in the models Ia and Ib but with slower growth of temperature and energy density rate. As in the case of the inter-halo medium, the qualitative behavior is similar to the cases with Ia and Ib models, but the effect of the first light in the models Ic and Id is already noticeable at larger redshifts because its energy density is higher there.
5.3. Luminescence of halo in the molecular lines
We went on to analyze the impact of the first light on the luminescence of a halo in the lines of the first molecules. Since the first light with temperature, Tfl ∼ 103 − 104, can excite the rotational-vibrational energy levels of the first molecules the luminescence of halos in the Cosmic Dawn would complicatedly depend on the spectral energy density of background radiation and the number density of the molecules.
To evaluate the luminescence of the halo in the molecular lines, we must compute the population of rotational levels. We integrate for that the kinetic equations of population or depopulation of rotational levels:
where the indices, i, and j, mark the rotational levels, ni is the number density of molecules at the rotational level i, and
are the rates of transitions between them. Here, i > j, Aij is the rate of spontaneous transition from i-level to j-level, Bij and Bji are the rates of induced transition and absorption at frequency νij, accordingly, Cji and Cij are the rates of collisional excitation of i-level from j-level and de-excitation from the i-level to the j-level, is the total energy density of radiation at any redshift from cosmological recombination to reionization. We used the same approximations for excitation and de-excitation rates of rotational levels of molecules H2, HD, and HeH+ as in our previous papers, Novosyadlyj et al. (2020) and Kulinich et al. (2020), and references therein.
Solution of the system of equations for populations of rotational levels with J = 0−5 gives us the possibility to compute the excitation temperature from the i-level to the j-level:
where hP is the Planck constant, kB is the Boltzmann constant, and gi, gj are statistical weights of levels. This, in turn, makes it possible to calculate the optical depths in the rotational lines, νij, along the diameter of the spherical homogeneous top-hat halo,
where mA is the atomic number of a molecule, rh is the radius of the halo in units of kpc, and the rest values are in the CGS system. The results show that the maximum values of , , and after virialization are in the range of 10−7 ÷ 10−9, which means that such halos are optically thin. Our estimation shows that they are also optically thin in the frequencies of photo-dissociations of molecules (τpd ≤ 10−6). Thus, the use of the optically thin halo approximation in the molecular lines and in the continuum is justified. Since the column densities of the molecules N(H2), N(HD), N(HeH+)≪1020 cm2, the self-shielding effect is not important for such halos (Draine & Bertoldi 1996).
Traditionally, the intensity of microwave radiation is expressed in units of brightness temperature:
Since the luminescence of the Dark Ages and Cosmic Dawn halos in the lines of the first molecules occurs on the cosmic microwave background at the same frequencies, it is reasonable to calculate the differential brightness temperature:
In the case of an optically thin top-hat halo the radiative transfer equation gives the next expression for the differential brightness temperature; for details, see Novosyadlyj et al. (2020):
Using it we compute the evolution of luminescence of warm and hot halos in the rotational lines of molecules H2, HD, and HeH+ during the Dark Ages and Cosmic Dawn with the different energy densities of the first light. The results are presented in Figs. 11 and 12 for the same models of the first light as in Fig. 10.
Fig. 11. Differential brightness temperature of warm halos (Tb ≈ 280 K) in the base rotational lines of molecules H2, HD, and HeH+ (from top to down) with and without of the first light: solid dark line – without first light, solid red and blue lines with background first light which reionize the medium at z = 6 and z = 8 accordingly, dotted and dashed blue lines with the first light with energy density 103 and 106 times higher than the background one (zrei = 8). In the left column, the first light is Ia, in the right one it is Ib. |
Fig. 12. Differential brightness temperature of hot halos (Tb ≈ 2280 K) in the base rotational lines H2, HD and HeH+ (from top to down) with and without of the first light: solid dark line is without the first light; solid red and blue lines with background are the first light that reionize the medium at z = 6 and z = 8, accordingly; dotted and dashed blue lines with the first light with an energy density that is 103 and 106 times higher than the background one (zrei = 8). In the left column, the first light is Ia, in the right one – Ib. |
We can see that the luminescence of a single halo in the rotational lines of H2 and HD molecules appears after the turn-around point of seed perturbation, when adiabatic heating of baryon matter becomes noticeable and disappears before hydrogen reionization of the intra-halo medium caused by photodissociation of molecules. It has formed a plateau with sharp uphill and downhill. Their detection could be very informative about the first halos, first molecules, and first light. Unfortunately, the amplitude of the differential brightness temperature is too low to be detected by current telescopes. However, in the case of the “hotter” first light with lower energy density, the differential brightness temperature of hot halos can reach a few hundredths or even tenths of a microkelvin (Fig. 12, right panels in the top and middle rows), at the sensitivity level of existing and planned advanced telescopes for the microwave range.
The time dependence of the luminescence of a single halo in the rotational lines of HeH+ molecular ions is as interesting with regard to both warm halos as well as hot ones. It has three phases of sharp changes: increasing after the turn-around point of seed perturbation, decreasing during the Cosmic Dawn caused by photodissociation of molecules in the case of the warm halo (Fig. 11, bottom row), or by collision in the case of the hot halo and a noticeable increase during the reionization. Thus, the halo could exhibit luminescence at the reionization epoch in the rotational lines of HeH+ molecular ions only (from the studied ones, of course). Again, however, the amplitude of differential brightness temperature in the basic rotational lines is too low to be detected by existing telescopes.
The magnitude of fluxes in the rotational lines of the first molecules obtained in this work and in Novosyadlyj et al. (2020) without reionizing by the first light practically coincides with that predicted in the early works (Dubrovich 1977; Maoli et al. 1994, 1996; Kamaya & Silk 2002, 2003; Ripamonti et al. 2002; Omukai & Kitayama 2003; Mizusawa et al. 2004, 2005; Dubrovich et al. 2008) despite completely different source models. The main conclusion about their marginal detectability by modern and planned telescopes coincides with the conclusions of these and other authors. The novel result presented here is the luminescence of H2 and HD molecules in the lines long before the reionization of atomic hydrogen and their disappearance due to photodissociation of these molecules shortly prior. The halo luminescence with M ∝ 1010 − 1013 M⊙ at 7 ≤ z ≤ 24 has been estimated using GIZMO cosmological simulations (see, e.g., Liu et al. 2019). The interpolation of their results to smaller halo masses shows a good agreement (up to orders of magnitude of fluxes) with our results in the models without the first light. The fact that there is no disappearance of H2 emission at z ≤ 8−9 due to photodissociation can be explained by the underestimation of UV radiation, since the model of reionization is not self-consistent in this case.
The results presented here were obtained within the framework of the simplest models of halo and first light, which can be far from real, as well as a limited set of chemical reactions and energy levels of molecules. Complicating the models by expanding the number of reactions and energy transitions in the molecules will be a task for further research.
6. Discussion and conclusions
We analyzed the ionization of the hydrogen atoms and the formation and destruction of the first molecules in the inter-proto-galactic medium of Cosmic Dawn, assuming the thermal energy distribution of the background radiation of the first sources of light. The spectral energy distribution at those epochs was presented by the Planck function with a certain temperature dependent on redshift, along with a dilution factor constrained by data about the redshift of reionization. The latter was assumed to be xHII = 0.5 at 6 ≤ z ≤ 8. We used the simplest model of the first light with a minimal number of parameters (temperature and dilution coefficient). Such radiation is not in thermodynamic equilibrium but features Planck energy distribution, which essentially simplifies the computations of all activation and deactivation and photodissociation of atoms and molecules, according to the generalized formula (5). It corresponds to the following scenario: the sources of the first light, proto-, and PopIII stars, evolve, while the average energy density and temperature of their radiation increase according to Eqs. (6) and (4). The parameters in Eq. (6) could be constrained using the observational data on the luminescence of halos in the lines of molecules, which we have evaluated here. However, in this paper, we set them arbitrarily in order to clarify the effect of the star formation rate on our estimates. Mathematically, the problem has consisted of the numerical integration of the system of kinetic equations of the ionization and recombination of atoms and molecules, as well as the formation and destruction of molecules and their ions together with the equations of expansion of the Universe and the thermal energy balance for the baryonic component. We only took into account the most important reactions (Table 1) that are connected with the most abundant atoms and molecules: atomic hydrogen H, helium He, deuterium D, molecular hydrogen H2, hydrogen deuteride molecule HD, and helium hydride ion HeH+. The obtained results lead to conclusions that provide a clearer understanding of the possible physical conditions at Cosmic Dawn and the epoch of reionization. Our conclusions are summarized as follows:
1. The number densities of hydrogen atoms, kinetic temperature, and the radiation field in the inter-proto-galaxy medium of Cosmic Dawn are such that the atoms are in the base state mainly and the case of A of photoionization is realized (Fig. 1).
2. The thermal radiation with a temperature of Tfl = 5000 K and dilution factor of αfl ∼ 2 ÷ 5 × 10−10 ionizes the atomic hydrogen up to xHII = 0.5 at 6 ≤ z ≤ 8. At the current epoch, the spectral energy density of such radiation is a few orders lower than the spectral energy density of CIB at frequencies of < 30 THz (λ > 10 μm). At higher frequencies, it would be detected by the current radio telescopes. For higher radiation temperatures, the dilution factor is lower and the current energy density of such radiation becomes too low for its detection. For the first light with comoving radiation temperature, Tfl = 10 000 K, which ionizes the inter-proto-galaxy medium at 6 ≤ z ≤ 8, the dilution factor is ∼10−17. Its spectral flux is in the range of the frequency maximum at the current epoch, which is ∼Jy, and masks in the Wien range of CIB. The range of values for the dilution factor αfl is narrow for any fixed radiation temperature of the first light, Tfl, at the moment of reionization, but the dependence αfl(Tfl) is steep since the potential of ionization of the atomic hydrogen is in the Wien range of spectral energy distribution (Figs. 2 and 3).
3. The first light heats the baryonic gas in the inter-proto-galaxy medium up to dozens or hundreds of kelvins at zrei, which depends on its radiation temperature and dilution factor. In any case, the change of baryonic temperature Tb becomes noticeable not long before the ionized hydrogen fraction reaches 1/2. The different history of increasing of radiation temperature, Tfl, forms the temperature history of the baryonic component at that time (Fig. 4).
4. The number densities of negative hydrogen ion H− and molecular ions H start to decrease long before reionization, since the thresholds energy of photodissociation for them are essentially lower than the potential of ionization of the atomic hydrogen (Figs. 5 and 6).
5. The most abundant molecules in the Dark Ages, namely, H2 and HD, are destroyed over the course of the Cosmic Dawn not long before hydrogen reionization since the threshold of photodissociation for them is lower than the potential of ionization of atomic hydrogen (Figs. 5–7).
6. The dependence of the number density of the helium hydride ion HeH+ on redshift in the Cosmic Dawn is intriguing via the competition of photodissociation by the first light and its formation in the reaction He + H+ → HeH+ + γ when the number density of H+ increase. The number density of HeH+ decreases in the models of the first light with a large number density of photons, and increases for diluted first light with higher radiation temperature (Figs. 5–7). We consider whether these molecules can survive in the remote corners of the intergalactic medium? This is an interesting task for further research with more realistic models of the spectral energy distribution of radiation in the different periods of the evolution of the Universe.
7. The time dependence of the luminescence of a single top-hat halo in the rotational lines of H2 and HD molecules has a form of plateau with sharp uphill and downhill and maximal amplitude of a differential brightness temperature of ∼10−10 − 10−8 K for warm halos and ∼10−8 − 10−6 K for hot ones (Figs. 11 and 12). The beginning of the plateau is connected with the formation and adiabatic heating of the forming halo, while the ending of the plateau is caused by photodissociation of the molecules before reionization.
8. The time dependence of the luminescence of a single halo in the rotational lines of HeH+ molecular ions has three phases of sharp changes: increasing during virialization halo (heating), decreasing during Cosmic Dawn (photodissociation or collision destruction of molecules), and increasing during the reionization epoch. The maximal amplitude of differential brightness temperature does not exceed 10−7 K.
9. The detection of the luminescence of the Dark Ages and Cosmic Dawn halos in the rotational lines of the first molecules could provide very important information about formation of the first halos, first molecules, and the first sources of light.
Acknowledgments
This work was performed within the project of Ministry of Education and Science of Ukraine “Modeling the luminosity of elements of the large-scale structure of the early universe and the remnants of galactic supernovae and the observation of variable stars” (state registration number 0122U001834) and was supported also by the International Center of Future Science of Jilin University. We would also like to thank the Armed Forces of Ukraine for their resilience and courage which saved us and provided security to finalize this work during the unprovoked invasion russian army into Ukraine.
References
- Abel, T., Anninos, P., Zhang, Yu, & Norman, M. L. 1997, New Astron., 2, 181 [NASA ADS] [CrossRef] [Google Scholar]
- Banados, E., Venemans, B. P., Mazzucchelli, C., et al. 2018, Nature, 553, 473 [NASA ADS] [CrossRef] [Google Scholar]
- Barkana, R., & Loeb, A. 2001, Phys. Rep., 349, 125 [NASA ADS] [CrossRef] [Google Scholar]
- Bates, D. R., & Opik, U. 1968, J. Phys. B, 1, 543 [NASA ADS] [CrossRef] [Google Scholar]
- Bauman, R. P., Ferland, G. J., & MacAdam, K. B. 2001, BAAS, 33, 1331 [NASA ADS] [Google Scholar]
- de Bernardis, P., Dubrovich, V., Encrenaz, P., et al. 1993, A&A, 269, 1 [NASA ADS] [Google Scholar]
- Bowman, J. D., Rogers, A. E. E., Monsalve, R. A., et al. 2018, Nature, 555, 67 [NASA ADS] [CrossRef] [Google Scholar]
- Bouwens, R. J., Illingworth, G. D., Oesch, P. A., et al. 2015, ApJ, 811, 140 [NASA ADS] [CrossRef] [Google Scholar]
- Bromm, V., & Yoshida, N. 2011, ARA&A, 49, 373 [CrossRef] [Google Scholar]
- Cooke, R. J., Pettini, M., & Steidel, C. C. 2018, ApJ, 855, 102 [Google Scholar]
- Davies, F. B., Hennawi, J. F., Bañados, E., et al. 2018, ApJ, 864, 142 [NASA ADS] [CrossRef] [Google Scholar]
- Draine, B. T., & Bertoldi, F. 1996, ApJ, 468, 269 [Google Scholar]
- Dubrovich, V. K. 1977, Sov. Astron. Lett., 3, 128 [NASA ADS] [Google Scholar]
- Dubrovich, V., Bajkova, A., & Khaikin, V. B. 2008, New Astron., 13, 28 [NASA ADS] [CrossRef] [Google Scholar]
- Ferland, G. J., Peterson, B. M., Horne, K., et al. 1992, ApJ, 387, 95 [NASA ADS] [CrossRef] [Google Scholar]
- Ferland, G. J., Chatzikos, M., Guzmán, F., et al. 2017, Rev. Mex. Astron. Astrof., 53, 385 [Google Scholar]
- de Jong, T. 1972, A&A, 20, 263 [NASA ADS] [Google Scholar]
- Galli, D., & Palla, F. 1998, A&A, 335, 403 [NASA ADS] [Google Scholar]
- Glover, S. C., & Jappsen, A.-K. 2007, ApJ, 666, 1 [NASA ADS] [CrossRef] [Google Scholar]
- Gosachinskij, I. V., Dubrovich, V. K., Zhelenkov, S. R., Il’in, G. N., & Prozorov, V. A. 2002, Astron. Rep., 46, 543 [NASA ADS] [CrossRef] [Google Scholar]
- Grachev, S. I., & Dubrovich, V. K. 1991, Astrophysics, 34, 124 [Google Scholar]
- Heays, A. N., Bosman, A. D., & van Dishoeck, E. F. 2017, A&A, 602, A102 [Google Scholar]
- Hills, R., Kulkarni, G., Meerburg, P. D., & Puchwein, E. 2018, Nature, 564, E32 [Google Scholar]
- Holliday, M. G., Muckerman, J. T., & Friedman, L. 1971, J. Chem. Phys., 54, 1058 [NASA ADS] [CrossRef] [Google Scholar]
- Kamaya, H., & Silk, J. 2002, MNRAS, 332, 251 [NASA ADS] [CrossRef] [Google Scholar]
- Kamaya, H., & Silk, J. 2003, MNRAS, 339, 1256 [NASA ADS] [CrossRef] [Google Scholar]
- Karpas, Z., Anicich, V., & Huntress, W. T. 1979, J. Chem. Phys., 70, 2877 [NASA ADS] [CrossRef] [Google Scholar]
- Kogut, A., Abitbol, M. H., Chluba, J., et al. 2019, BAAS, 51, 113 [Google Scholar]
- Kulinich, Yu, Novosyadlyj, B., Shulga, V., & Han, W. 2020, Phys. Rev. D, 101, 083519 [NASA ADS] [CrossRef] [Google Scholar]
- Launay, J. M., Le Dourneuf, M., & Zeippen, C. J. 1991, A&A, 252, 842 [NASA ADS] [Google Scholar]
- Liu, B., Jaacks, J., Finkelstein, S. L., & Bromm, V. 2019, MNRAS, 486, 3617 [NASA ADS] [CrossRef] [Google Scholar]
- Maio, U., Péroux, C., & Ciardi, B. 2022, A&A, 657, A47 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Maoli, R., Melchiorri, F., & Tosti, D. 1994, ApJ, 425, 372 [NASA ADS] [CrossRef] [Google Scholar]
- Maoli, R., Ferrucci, V., Melchiorri, F., & Tosti, D. 1996, ApJ, 457, 1 [NASA ADS] [CrossRef] [Google Scholar]
- Mason, C. A., Treu, T., Dijkstra, M., et al. 2018, ApJ, 856, 2 [Google Scholar]
- Matsuda, T., Sato, H., & Takeda, H. 1971, Prog. Theor. Phys., 46, 416 [CrossRef] [Google Scholar]
- Mizusawa, H., Nishi, R., & Omukai, K. 2004, PASJ, 56, 487 [NASA ADS] [Google Scholar]
- Mizusawa, H., Omukai, K., & Nishi, R. 2005, PASJ, 57, 951 [NASA ADS] [Google Scholar]
- Moseley, J., Aberth, W., & Peterson, J. A. 1970, Phys. Rev. Lett., 24, 435 [NASA ADS] [CrossRef] [Google Scholar]
- Novosyadlyj, B., Tsizh, M., & Kulinich, Yu 2016, Gen. Rel. Grav., 48, 30 [NASA ADS] [CrossRef] [Google Scholar]
- Novosyadlyj, B., Sergijenko, O., & Shulga, V. M. 2017, Kinematics Phys. Celestial Bodies, 33, 255 [NASA ADS] [CrossRef] [Google Scholar]
- Novosyadlyj, B., Shulga, V., Han, W., Kulinich, Yu, & Tsizh, M. 2018, ApJ, 865, 38 [NASA ADS] [CrossRef] [Google Scholar]
- Novosyadlyj, B., Shulga, V., Kulinich, Yu, & Han, W. 2020, ApJ, 888, 27 [Google Scholar]
- Omukai, K., & Kitayama, T. 2003, ApJ, 599, 738 [NASA ADS] [CrossRef] [Google Scholar]
- ONeil, S. V., & Reinhardt, W. P. 1978, J. Chem. Phys., 69, 2126 [NASA ADS] [CrossRef] [Google Scholar]
- Osterbrock, D. E. 1989, Astrophysics of Gaseous Nebulae and Active Galactic Nuclei (University Science Books), 408 [Google Scholar]
- Peebles, P. J. E. 1968, ApJ, 153, 1 [NASA ADS] [CrossRef] [Google Scholar]
- Peimbert, A., Peimbert, M., & Luridiana, V. 2016, Rev. Mex. Astron. Astrofis., 52, 419 [NASA ADS] [Google Scholar]
- Persson, C. M., Maoli, R., Encrenaz, P., et al. 2010, A&A, 515, A72 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Planck Collaboration I. 2020, A&A, 641, A1 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Planck Collaboration VI. 2020, A&A, 641, A6 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Poulaert, G., Brouillard, F., Claeys, W., McGowan, J. W., & van Wassenhove, G. 1978, J. Phys. B, 11, L671 [NASA ADS] [CrossRef] [Google Scholar]
- Puy, D., Alecian, G., Le Bourlot, J., et al. 1993, A&A, 267, 337 [NASA ADS] [Google Scholar]
- Ripamonti, E., Haardt, F., Ferrara, A., & Colpi, M. 2002, MNRAS, 334, 401 [NASA ADS] [CrossRef] [Google Scholar]
- Roberge, W., & Dalgarno, A. 1982, ApJ, 255, 489 [NASA ADS] [CrossRef] [Google Scholar]
- Schleicher, D. R. G., Galli, D., Palla, F., et al. 2008, A&A, 490, 521 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Schneider, I. F., Dulieu, O., Giusti-Suzor, A., & Roueff, E. 1994, ApJ, 424, 983 [NASA ADS] [CrossRef] [Google Scholar]
- Seager, S., Sasselov, D. D., & Scott, D. 1999, ApJ, 523, L1 [NASA ADS] [CrossRef] [Google Scholar]
- Seager, S., Sasselov, D. D., & Scott, D. 2000, ApJS, 128, 407 [NASA ADS] [CrossRef] [Google Scholar]
- Smith, D., Adams, N. G., & Alge, E. 1982, ApJ, 263, 123 [NASA ADS] [CrossRef] [Google Scholar]
- Stancil, P. C. 1994, ApJ, 430, 360 [NASA ADS] [CrossRef] [Google Scholar]
- Stancil, P. C., Babb, J. F., & Dalgarno, A. 1993, ApJ, 414, 672 [NASA ADS] [CrossRef] [Google Scholar]
- Verner, D. A., & Ferland, G. J. 1996, ApJS, 103, 467 [NASA ADS] [CrossRef] [Google Scholar]
- Verner, D. A., Ferland, G. J., Korista, K. T., & Yakovlev, D. G. 1996, ApJ, 465, 487 [Google Scholar]
- Vonlanthen, P., Rauscher, T., Winteler, C., et al. 2009, A&A, 503, 47 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Watson, W. D., Christensen, R. B., & Deissler, R. J. 1978, A&A, 69, 159 [NASA ADS] [Google Scholar]
- Yousif, F. B., & Mitchell, J. B. A. 1989, Phys. Rev. A, 40, 4318 [NASA ADS] [CrossRef] [Google Scholar]
- Zygelman, B., & Dalgarno, A. 1990, ApJ, 365, 239 [NASA ADS] [CrossRef] [Google Scholar]
Appendix A: Rates of photorecombination and effective cross-sections of photoionization for HI and DI
Before the epoch of reionization, during epochs of the Dark Ages and Cosmic Dawn, hydrogen atoms were mainly in the ground electron state. To obtain the total rate of electron recombination on the second level, we use the analytical approximation of the total rate of recombination on all levels by Verner & Ferland (1996):
with the coefficients a = 7.982 ⋅ 10−11, b = 0.748, T1 = 3.148, T2 = 7.036 ⋅ 105, and the approximation of the rate of recombination on the base level by Ferland et al. (1992):
with refined values of approximation coefficients from the publicly available code Cloudy (Ferland et al. 2017) a1 = −10.781359, b1 = −0.38790005, c1 = 4.6748385, d1 = 0.063898347, e1 = −0.87297469, f1 = −5.0923340 ⋅ 10−3, g1 = 0.081445619, h1 = 2.4935490 ⋅ 10−4, and i1 = −3.9043790 ⋅ 10−3. So, the expected value of the recombination rate coefficient for case B is as follows:
The testing shows that this approximation of the case B of recombination coefficient αrB(T) matches the values in the range of electron temperature 3 ≤ T ≤ 105K computed by Cloudy with an accuracy of about 0.01%.
The rate of photoionization of atomic hydrogen from the level n by radiation with intensity, I(ν), we can compute, for the known effective cross-section of photoionization from this level σn(ν), according to Abel et al. (1997), as follows:
where νnc is the threshold frequency of photoionization from level n. For the intergalactic cold rarefied medium of the Dark Ages and Cosmic Dawn, photoionization goes mainly from base (n = 1, l = 0) or second metastable levels (n = 2, l = 0). We use the analytical approximation proposed by Verner et al. (1996) with a small correction:
where , σ0, E = hν, E0, P, ya and α are the best-fit parameters, which we re-determined using the open source codes by Bauman et al. (2001) for computations the cross-sections of photoionization of hydrogen from the base level (n = 1, l = 0) and the second one (n = 2, l = 0) and (n = 2, l = 1) shown in Fig. A.1 as solid lines. The best-fit coefficients and accuracy of approximation are presented in Table A.1.
The effective cross-sections of photoionization of HI from the base and second levels are shown in Fig. A.1 by the dotted red line. We assume that the rates of photorecombination and effective cross-sections of photoionization for deuterium atoms are the same as for atomic hydrogen.
Fig. A.1. Effective cross-sections of photoionization of hydrogen atoms H from the base level (n = 1, l = 0) and second one (n = 2, l = 0), and (n = 2, l = 1). Solid lines show the numerical computations, dotted ones the analytical approximations. |
Coefficients of approximation (A2) for cross-sections of photoionization of HI from the first and second energy levels.
Appendix B: Rates of photodissociation of molecules H2 and HD
At the Cosmic Dawn, energetic photons appear and they are capable to ionize of atoms and molecules and dissociate the last. The threshold of direct photodissociation of H2 (H2 + γ → H + H) is 4.53 eV, so, it can be important long before the reionization. The physics and computations of the cross-section of this reaction are well described in Glover & Jappsen (2007), Heays et al. (2017) (see also citations therein). Here, we use the photodissociation cross-section of molecule H2 computed by Heays et al. (2017) and presented on the site of Leiden University9 in the table form with high wavelength resolution (4 ⋅ 10−5μm). We illustrate these data in Figure B.1 using low wavelength resolution table data (0.1μm). We computed the rate of dissociation according to (A1) with I(ν) = Bν(T) for range of temperatures 103 < T < 2 ⋅ 105 K. The result is shown in the right panel of Fig. B.1 by open circles. We also propose the analytical approximation of this dependence by the formula:
with the best-fit parameters (a, b, c, d, f) = (2394.9, 0.62681, 2.4635, 0.56957, 140050) for T ≤ 15000 K and (a, b, c, d, f) = (6.4471 ⋅ 107, 0.322, 6.5341 ⋅ 10−8, 1.4752, 153570) for T > 15000 K. Its deviation from the numerical result is not larger than 1% in the interested temperature range (bottom small panel of the left figure). We present also in the main panel of the right figure the analytical approximation by Glover & Jappsen (2007), Schleicher et al. (2008) for comparison. We can see that in the low temperature range (T < 5000 K), they are close. We assume here that the rate of direct photodissociation of molecule HD (HD + γ → D + H) is the same as for molecule H2.
Fig. B.1. Effective cross-sections of hydrogen molecule H2 photodissociation from Heays et al. (2017) shown in the left panel and the rate of dissociation computed here in the right panel. The solid red line is the analytical approximation (B1), the dotted blue line is from Glover & Jappsen (2007), Schleicher et al. (2008). The deviation of analytical approximation (B1) from the numerical results is shown in the low panel. |
All Tables
Coefficients of approximation (A2) for cross-sections of photoionization of HI from the first and second energy levels.
All Figures
Fig. 1. Characteristic lifetime of the hydrogen atom in the excited states 2s and 2p (, ), time intervals between following recombinations to the base and excited levels (, ), and the resulting photoionization (, ) by CMBR (at z > 200) as well as the first light photons (z < 30). The computations are made for the cosmological model with Planck 2018 parameters and thermal models of the first light. |
|
In the text |
Fig. 2. Spectral energy density of the first light and CMBR at z = 8 (solid lines) and at z = 6 (dashed lines) when the ionization of hydrogen is ≈0.5 for models Ia and Ib. The vertical dotted lines show the thresholds of dissociation of molecules H−, HeH+, H, H2, and HD (overlapped with H2) and the ionization of HI, correspondingly. The dotted lines show these spectral energy densities at the current epoch. The spectral energy density of cosmic infrared background (dotted green line) is shown for comparison. In the models Ic and Id, the corresponding lines are lower by the ratios of dilution coefficients. |
|
In the text |
Fig. 3. Reionization of hydrogen by the first light at the end of the Cosmic Dawn epoch. Left: evolution of the ionized hydrogen fraction xHII for the first light models Ia-Id, which provide xHII ≈ 0.5 at 6 ≤ z ≤ 8. Right: ranges of dilution coefficient against the temperature of thermal radiation for two evolution models of the first light. |
|
In the text |
Fig. 4. Evolution of the number density of H, D and He (top) and temperature of CMBR and baryonic matter through the Dark Ages to Cosmic Dawn (bottom). |
|
In the text |
Fig. 5. Evolution of the number densities of atoms, ions, neutral, and ionized molecules ni(z) from cosmological recombination through the Dark Ages and Cosmic Dawn up to reionization for the models of the first light Ia and Ib: the red lines for zrei = 6, the blue lines for zrei = 8. The black lines are the number density of species for the case without the first light, . |
|
In the text |
Fig. 6. Evolution of the number densities of atoms, ions, and neutral and ionized molecules ni(z) from cosmological recombination through the Dark Ages and Cosmic Dawn up to reionization for the models of the first light Ic and Id: the red lines for zrei = 6, the blue lines for zrei = 8. The black lines are the number density of species for the case without the first light, . |
|
In the text |
Fig. 7. Impact of the first light on the number densities of molecules at Cosmic Dawn for the different evolution of the energy density of thermal radiation of the first sources (models Ia-Id) (lines for H2 and HD superimposed). |
|
In the text |
Fig. 8. Formation of a spherical halo from a specific cosmological perturbation. Top panel: evolution of number densities of H, D, and He in the top-hat halo with k = 80 Mpc (Mh = 1.26 × 10−6 M⊙) and initial amplitude of adiabatic scalar perturbation δr = 1.5 × 10−4 at z = 108). Bottom panel: evolution of the temperature of CMBR (green line) and kinetic temperature of the top-hat halo with adiabatic heating (dark line) and virial temperature (8; red line). |
|
In the text |
Fig. 9. Ratio of the number densities of molecules in the warm and hot halo () to the number density of molecules at the cosmological background () without the first light. |
|
In the text |
Fig. 10. Impact of the first light on the number densities of molecules in the halos of Cosmic Dawn. In the top row, the models of the first light are the same as in background ones (left – Ia, right – Ib), which reionize the medium at z = 6 (dotted lines) and z = 8 (solid lines). In the bottom row, the energy density of the first light in the halos are 103 (dotted line) and 106 (solid lines) times higher than energy density in the corresponding models in the top panels shown by solid lines. In all the panels, the lines for H2 and HD are superimposed. |
|
In the text |
Fig. 11. Differential brightness temperature of warm halos (Tb ≈ 280 K) in the base rotational lines of molecules H2, HD, and HeH+ (from top to down) with and without of the first light: solid dark line – without first light, solid red and blue lines with background first light which reionize the medium at z = 6 and z = 8 accordingly, dotted and dashed blue lines with the first light with energy density 103 and 106 times higher than the background one (zrei = 8). In the left column, the first light is Ia, in the right one it is Ib. |
|
In the text |
Fig. 12. Differential brightness temperature of hot halos (Tb ≈ 2280 K) in the base rotational lines H2, HD and HeH+ (from top to down) with and without of the first light: solid dark line is without the first light; solid red and blue lines with background are the first light that reionize the medium at z = 6 and z = 8, accordingly; dotted and dashed blue lines with the first light with an energy density that is 103 and 106 times higher than the background one (zrei = 8). In the left column, the first light is Ia, in the right one – Ib. |
|
In the text |
Fig. A.1. Effective cross-sections of photoionization of hydrogen atoms H from the base level (n = 1, l = 0) and second one (n = 2, l = 0), and (n = 2, l = 1). Solid lines show the numerical computations, dotted ones the analytical approximations. |
|
In the text |
Fig. B.1. Effective cross-sections of hydrogen molecule H2 photodissociation from Heays et al. (2017) shown in the left panel and the rate of dissociation computed here in the right panel. The solid red line is the analytical approximation (B1), the dotted blue line is from Glover & Jappsen (2007), Schleicher et al. (2008). The deviation of analytical approximation (B1) from the numerical results is shown in the low panel. |
|
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.