We investigate the effects of dust on Ly{\alpha} photons emergent from an optically thick medium by solving the integro-differential equation of the radiative transfer of resonant photons. To solve the differential equations numerically we use the Weighted Essentially Non-oscillatory method (WENO). Although the effects of dust on radia-tive transfer is well known, the resonant scattering of Ly{\alpha} photons makes the problem non-trivial. For instance, if the medium has the optical depth of dust absorption and scattering to be {\tau}a>>1, {\tau}>>1, and {\tau}>>{\tau}a, the effective absorption optical depth in a random walk scenario would be equal to \surd {\tau}a({\tau}a+{\tau}). We show, however, that for a resonant scattering at frequency {\nu}0, the effective absorption optical depth would be even larger than {\tau}({\nu}0). If the cross section of dust scattering and absorption is frequency-independent, the double-peaked structure of the frequency profile given by the resonant scattering is basically dust-independent. That is, dust causes neither narrowing nor widening of the width of the double peaked profile. One more result is that the time scales of the Ly{\alpha} photon transfer in the optically thick halo are also basically independent of the dust scattering, even when the scattering is anisotropic. This is because those time scales are mainly determined by the transfer in the frequency space, while dust scattering, either isotropic or anisotropic, does not affect the behavior of the transfer in the frequency space when the cross section of scattering is wavelength-independent. This result does not support the speculation that dust will lead to the smoothing of the brightness distribution of Ly{\alpha} photon source with optical thick halo.

With a state-of-the-art numerical method for solving the integral-differential equation of radiative transfer, we investigate the flux of the Ly$\alpha$ photon $\nu_0$ emergent from an optically thick halo containing a central light source. Our focus is on the time-dependent effects of the resonant scattering. We first show that the frequency distribution of photons in the halo are quickly approaching to a locally thermalized state around the resonant frequency, even when the mean intensity of the radiation is highly time-dependent. Since initial conditions are forgotten during the thermalization, some features of the flux, such as the two peak structure of its profile, actually are independent of the intrinsic width and time behavior of the central source, if the emergent photons are mainly from photons in the thermalized state. In this case, the difference $|\nu_{\pm}-\nu_0|$, where $\nu_{\pm}$ are the frequencies of the two peaks of the flux, cannot be less than $2$ times of Doppler broadening. We then study the radiative transfer in the case where the light emitted from the central source is a flash. We calculate the light curves of the flux from the halo. It shows that the flux is still a flash. The time duration of the flash for the flux, however, is independent of the original time duration of the light source but depends on the optical depth of the halo. Therefore, the spatial transfer of resonant photons is a diffusion process, even though it is not a purely Brownian diffusion. This property enables an optically thick halo to trap and store thermalized photons around $\nu_0$ for a long time after the cease of the central source emission. The photons trapped in the halo can yield delayed emission, of which the profile also shows typical two peak structure as that from locally thermalized photons. Possible applications of these results are addressed.

The 21 cm emission and absorption from gaseous halos around the first generation of star depend on the Wouthuysen-Field (W-F) coupling, which relates the spin temperature with the kinetic temperature of hydrogen gas via the resonant scattering between Lyman alpha photons and neutral hydrogen. Although the center object generally is a strong source of these photons, the transfer of these photons in the 21 cm region is inefficient, as the optical depth of the photons is large. Consequently, these photons from the source may not be able to transfer to the entire 21 cm region timely to provide the W-F coupling. This problem is important because the lifetime of first stars generally is short. The problem is investigated with numerical solution of the integro-differential equation, which describes the kinetics of these resonant photons in both physical and frequency spaces. We show that the photon transfer process in the physical space is actually coupled to that in the frequency space. Firstly diffusion in the frequency space provides a shortcut for the diffusion in the physical space. It makes the mean time for the escape of the resonant photon in optical depth \tau media roughly proportional to the optical depth \tau, not \tau^2. Secondly the resonant scattering is effective in bouncing photons with a frequency which is not equal to initial frequency back to the initial frequency. This process can restore initial frequency photons and establish the local Boltzmann distribution of the photon spectrum around the initial frequency. Therefore, the mechanism of 'escape via shortcut' plus 'bounce back' enables W-F coupling to be properly realized in the 21 cm region around first stars. This mechanism also works for photons injected into the 21 cm region by redshift.

We develop a numerical solver for the integral-differential equations, which describes the radiative transfer of photon distribution in the frequency space with resonant scattering of Lyalpha photons by hydrogen gas in the early universe. The time-dependent solutions of this equation is crucial to the estimation of the effect of the Wouthuysen-Field (WF) coupling in relation to the 21 cm emission and absorption at the epoch of reionization. The resonant scattering leads to the photon distribution in the frequency space to be piecewise smooth containing sharp changes. The weighted essentially nonoscillatory (WENO) scheme is suitable to handle this problem, as this algorithm has been found to be highly stable and robust for solving Boltzmann equation. We test this numerical solver against analytic solutions of the evolution of the photon distribution in rest background, analytic solution in expanding background without resonant scattering and formation of local Boltzmann distribution around the resonant frequency with the temperature same as that of atom for recoil. We find that evolution of photon distribution undergoes three phases; profile is similar to the initial one, a flat plateau (without recoil) or local Boltzmann distribution (with recoil) forms around the resonant frequency, and finally the distribution around the resonant frequency is saturated when the photons from the source is balanced by the redshift of the expansion. This result indicates that the onset of the W-F coupling should not be determined by the third phase, but by the time scale of the second phase. We found that the time scale of the W-F coupling is equal to about a few hundreds of the mean free flight time of photons with resonant frequency, and is independent of the Sobolev parameter if this parameter is much less than 1.

We study the Wouthuysen-Field coupling at early universe with numerical solutions of the integrodifferential equation describing the kinetics of photons undergoing resonant scattering. The numerical solver is developed based on the weighted essentially non-oscillatory (WENO) scheme for the Boltzmann-like integrodifferential equation. We focus on the time evolution of the Wouthuysen-Field (W-F) coupling in relation to the 21 cm emission and absorption at the epoch of reionization. We show that a local Boltzmann distribution will be formed if photons with frequency \sim \nu_0 have undergone a ten thousand or more times of scattering, which corresponds to the order of 10^3 yrs for neutral hydrogen density of the concordance \Lambda CDM model. The time evolution of the shape and width of the local Boltzmann distribution actually doesn't dependent on the details of atomic recoil, photon sources, or initial conditions very much. However, the intensity of photon flux at the local Boltzmann distribution is substantially time-dependent. The time scale of approaching the saturated intensity can be as long as 10^5-10^6 yrs for typical parameters of the \Lambda CDM model. The intensity of the local Boltzmann distribution at time less than 10^5 yrs is significantly lower than that of the saturation state. Therefore, it may not be always reasonable to assume that the deviation of the spin temperature of 21 cm energy states from cosmic background temperature is mainly due to the W-F coupling if first stars or their emission/absorption regions evolved with a time scale equal to or less than Myrs.