A&A 475, 109-114 (2007)
DOI: 10.1051/0004-6361:20077333
J. Chluba1 - R. A. Sunyaev1,2
1 - Max-Planck-Institut für Astrophysik, Karl-Schwarzschild-Str. 1,
85741 Garching bei München, Germany
2 -
Space Research Institute, Russian Academy of Sciences, Profsoyuznaya 84/32,
117997 Moscow, Russia
Received 20 February 2007 / Accepted 16 September 2007
Abstract
We compute the corrections to the cosmological hydrogen
recombination history due to delayed feedback of Lyman-series photons
and the escape in the Lyman-continuum. The former
process is expected to slightly delay recombination, while the latter
should allow the medium to recombine a bit faster.
It is shown that the subsequent feedback of released Lyman-n photons on the lower lying Lyman-(n-1) transitions yields a maximal
correction of
at
.
Including only Lyman-
feedback onto the Lyman-
transition, accounts for most of the effect.
We find corrections to the cosmic microwave background TT and EE power spectra with typical peak to peak amplitude
and
at
.
The escape in the Lyman-continuum and feedback of Lyman-
photons on
the photoionization rate of the second shell lead to modifications of
the ionization history which are very small (less than
).
Key words: cosmic microwave background - cosmology: miscellaneous
The tremendous advances in observations of the cosmic microwave background
(CMB) temperature and polarization angular fluctuations
(e.g. Hinshaw et al. 2006; Page et al. 2006) and the prospects with the PLANCK
Surveyor have motivated several groups to
re-examine the problem of cosmological hydrogen recombination,
with the aim of identifying previously neglected physical processes
which could affect the ionization history of the Universe at the level of
,
and may lead to modifications of the cosmological hydrogen
recombination spectrum
(Chluba et al. 2007; Beigman & Sunyaev 1978; Burgin 2003; Kholupenko et al. 2005; Rybicki & dell'Antonio 1993; Dubrovich 1975; Rubiño-Martín et al. 2006; Wong et al. 2006; Chluba & Sunyaev 2006a; Dubrovich & Stolyarov 1995; Dubrovich & Shakhvorostova 2004), which could become observable in the future.
For example, effects connected with the two-photon transitions from
high s and d-states to the ground state (Dubrovich & Grachev 2005; Wong & Scott 2007; Chluba & Sunyaev 2007), the induced 2s-two-photon decay (Chluba & Sunyaev 2006b), the
increase in the 2s-two-photon absorption rate due to the large
Lyman-
spectral distortion (Kholupenko & Ivanchik 2006), and details in the
evolution of the populations of the angular momentum sub-states
(Chluba et al. 2007; Rubiño-Martín et al. 2006) were discussed.
All these studies show that sub-percent-level corrections to the ionization
history do exist, which in principle could bias the values of the key
cosmological parameters (Lewis et al. 2006).
It is clear, that when reaching percent-level accuracy in the determination of
the key cosmological parameters or when considering signatures from inflation,
e.g. the possibility of a running spectral index, accurate theoretical
predictions of the ionization history are required.
Here we examine the effects due to feedback of hydrogen Lyman-series photons
and the escape of photons in the Lyman-continuum on the ionization history
during the epoch of hydrogen recombination, and the impact of this
process on the Lyman-
distortion of the CMB blackbody spectrum.
The strongest distortions of the CMB spectrum arising during the epoch of
hydrogen recombination are due to the Lyman-
transition and the 2s
two-photon decay (Peebles 1968; Zeldovich et al. 1968).
The feedback of these excess photons on the photoionization
rates of the second shell has been considered by Seager et al. (1999,2000),
with no significant changes for the number density of free electrons,
.
We confirm this statement for the feedback from Lyman-
photons and
find a maximal correction of
at
,
the redshift which roughly corresponds to the time, when the
maximum of the CMB spectral distortion due to the Lyman-
transition is
reaching the Balmer-continuum frequency.
However, it has been shown by Kholupenko & Ivanchik (2006) that the huge excess of
photons in the Wien-tail of the CMB due to the Lyman-
distortion leads
to an increase of the
two-photon absorption
rate, which delays recombination and introduces corrections to the ionization
history at the percent-level.
Similarly, one expects some feedback of escaping Lyman-series photons on
lower lying Lyman-transitions and eventually the Lyman-
transition.
Recently, this process was also considered by Switzer & Hirata (2007) for both
helium and hydrogen recombination, using an iterative approach.
Here we directly compute the escaping Lyman-n radiation within the Sobolev
approximation (e.g. see Chluba et al. 2007; Rubiño-Martín et al. 2006) and include the
feedback on the Lyman-(n-1) transition at all times by evaluating the
distorted spectrum at the corresponding frequency.
In addition, we assume that due to the huge optical depth in all the
Lyman-lines, at the end only the Lyman-
line remains and all the
Lyman-n lines will be completely re-processed by the closest lower lying
Lyman-(n-1) transition. Therefore, feedback always only works as
.
It is also important to ensure that the line is not producing any
feedback on itself.
In addition, we estimate the escape probability for the Lyman-continuum and include the approximate net rate for direct recombinations into our multi-level code. As argued earlier (Peebles 1968; Zeldovich et al. 1968), we find that the modification of the ionization history due to this process is indeed completely negligible.
Table 1:
Appearance of the first few Lyman lines for the computations not including any
feedback. In particular we give the approximate redshift,
,
at
which
is maximal (see Fig. 1), and the
redshift at which the peak of the line reaches the next lower Lyman resonance,
,
where
is the resonance frequency of the Lyman-k transition. Also
we give the approximate total number ofescaping photons per hydrogen nucleus for each transition.
![]() |
Figure 1:
Spectral distortions due to the first few Lyman-transitions at z = 0.
No feedback
![]() ![]() |
Open with DEXTER |
![]() |
Figure 2:
Changes in the free electron fraction due to inclusion of Lyman-series
feedback relative to the reference model without feedback.
The computations were performed for a 15-shell atom where for the Lyman-series
the escape of photons was modeled using the Sobolev approximation. The curves
are labeled according to the sequence of feedback that was included. For
example,
![]() ![]() ![]() ![]() ![]() |
Open with DEXTER |
In Fig. 2 we present the results for detailed computations of
feedback-induced changes in the free electron fraction. Including
only Lyman-
feedback on the Lyman-
transition accounts for most
of the effect, leading to a maximal difference of
at
.
When including in addition the feedback of Lyman-
on Lyman-
,
the
maximal relative difference decreases to
,
now peaking at
.
This behavior is expected, since the Lyman-
photon reduces the
transition rate in the Lyman-
channel such that its feedback on the
Lyman-
line should become slightly smaller. In addition, the number of
electrons reaching the ground state via the Lyman-
transition will be
reduced, which very likely leads to the small increase of the correction at
.
Subsequent inclusion of more feedback results in an alternating behavior of
.
We found that for 30-shells, convergence is reached when including feedback
from np-levels with
.
In Fig. 2 we also give the
full result including all feedback within a 30-shell hydrogen atom.
In addition we found that the correction barely depends on the total number of
included shells.
Therefore we expect that even for computations with up to 100-shells the total
feedback-induced relative difference will not change significantly.
Examining the final differences in the Lyman-
distortion and
the feedback-induced modifications of the Lyman-
and
lines
before their absorption within the corresponding lower lying resonance
(see Fig. 3) shows that the situation is a bit more
complicated.
Including only Lyman-
feedback the amount of emission in the
Lyman-
line reduces by
at
.
Note that this
redshift is very close to the expected value,
,
for the
maximal feedback (see Table 1).
On the other hand, due to the feedback-induced small changes in the
populations of the levels, the net 2s-1s two-photon transition rate also
increases with a maximum of roughly
at
(see Fig. 4), which partially cancels the delaying feedback effect on the
Lyman-
transition.
When including Lyman-
feedback, as expected, the emission in the
Lyman-
line reduces, with a maximal difference at the expected value
.
Adding Lyman-
feedback reduces the strength of the
feedback of Lyman-
on the Lyman-
transition by a significant
fraction. Comparing the relative strength of the Lyman-
at its maximum
with the Lyman-
line at
(Fig. 1) yields
,
which indicates that Lyman-
should be able to affect Lyman-
strongly.
Looking at the relative difference in the Lyman-
when including
Lyman-
feedback shows that at
the line is
reduced by
.
![]() |
Figure 3:
Feedback-induced changes of
![]() ![]() ![]() ![]() |
Open with DEXTER |
![]() |
Figure 4:
Changes in the net 2s-1s two-photon decay rate due to the inclusion of
![]() |
Open with DEXTER |
With this in mind, in the vicinity of the Lyman-continuum the photon number
only changes by recombination and ionization to the ground state (we neglect
electron scattering):
The optical depth for absorption of a photon, which has been emitted at
redshift
and is observed at frequency
at redshift
,
within the Lyman-continuum is given by
It is possible to simplify Eq. (5) even further. First,
since the emission profile is very narrow one may use
,
with the
corresponding normalization constant
.
Assuming
one then obtains
![]() |
Figure 5:
Escape probability for the Lyman-continuum as given by
Eq. (5) and Eq. (6). For
comparison we also show the Sobolev escape probability for the
Lyman-![]() |
Open with DEXTER |
![]() |
Figure 6:
Changes in the free electron fraction due to inclusion of Lyman-continuum
escape relative to the reference model without any direct recombinations to
1s. The computations were performed for a 5-shell atom where for the
Lyman-series escape the Sobolev approximation was used. Note that the absolute
value of
![]() ![]() |
Open with DEXTER |
Using Eq. (6) it is easy to include the possibility of
direct recombinations to the ground state connected with the escape of
continuum photons. For this,one should add
Comparing the escape probability for the Lyman-continuum with the Sobolev
escape probability for the Lyman-
photons (see Fig. 5) one
can see that
is roughly 10-100 times
larger at most times.
However, including the possibility of direct recombinations to the ground
state in our recombination code yields a tiny correction to
(see
Fig. 6), which shows that
as given by
Eq. (7) is still many orders of magnitude smaller than the net
Lyman-
transition rate. One can safely neglect direct recombinations
to the ground state of hydrogen for computations of the recombination history.
![]() |
Figure 7: Changes in the CMB temperature (TT) and polarization (EE) power spectra. The differences were computed using our modified versions of CMBEASY (Doran 2005), which allows loading of pre-calculated recombination histories, and where the corresponding RECFAST-routine was improved to achieve higher numerical accuracy with solvers from the NAG-library. |
Open with DEXTER |
As shown by Chluba & Sunyaev (2007) due to two-photon processes
one can expect some asymmetries in the emission profiles of the higher
Lyman-series lines. This could modify the amount and time-dependence for the
feedback process and may lead to
additional differences in the results presented above. However, because the
bulk of photons that escape in the red wing of a particular Lyman-transition
is expected to come from close to the line center, this may be of minor
importance. Here the most interesting aspect may be the changes in the
time-dependence of the feedback, but more detailed computations would be
necessary to understand this problem.
In addition one should look at the feedback-induced corrections in the
hydrogen recombination history due to the
spectral
distortions, in particular due to the photons appearing in the
resonance transition, the
intercombination line and the
two-photon continuum.
It is clear that all high frequency
lines will be re-processed
during
recombination, and correspondingly can only
affect the ionization history during that epoch.
However, for the feedback of the
lines on hydrogen recombination
it will be important to compute the re-processing of photons in the
Lyman-continuum and all the subsequent Lyman-series transitions. Since all of these transitions are very optically thick a huge part of the
photons released during
and
recombination will never reach the observer today, but due to
re-absorption by neutral hydrogen, they are fully converted to hydrogen
Lyman-
and 2s-1s continuum photons (see
also Kholupenko et al. 2007; Rubiño-Martín et al. in preparation).
Since most of the helium
photons are
released at
,
they should be re-absorbed by neutral hydrogen
atoms at the early stages
.
One can also conclude this from the paper of Kholupenko et al. (2007), where the
re-processed
photons appear on the red
side of the hydrogen Lyman-
distortion (most hydrogen Lyman-
photons are released at
). Therefore the helium photons affect the
hydrogen recombination history well before the maximum of the
Thomson-visibility function (
), and hence should have a rather
small impact on the CMB power spectra. However, a more careful computation is
required and will be described in a forthcoming paper (Rubiño-Martín et al. in preparation).
Acknowledgements
The authors thank the anonymous referee for his useful comments.