Two methods are developed for solving the steady-state spherically symmetric radiative transfer equation for resonance line radiation emitted by a point source in the Intergalactic Medium. One method is based on solving the ray and moment equations using finite differences. The second uses a Monte Carlo approach incorporating methods that greatly improve the accuracy compared with previous approaches in this context. Several applications are presented serving as test problems for both a static medium and an expanding medium, including inhomogeneities in the density and velocity fields. Solutions are obtained in the coherent scattering limit and for Doppler RII redistribution with and without recoils. We find generally that the radiation intensity is linear in the cosine of the azimuthal angle with respect to radius to high accuracy over a broad frequency region across the line centre for both linear and perturbed velocity fields, yielding the Eddington factors f(nu) = 1/3 and g(nu) = 3/5. We show the radiation field produced by a point source divides into three spatial regimes for a uniformly expanding homogeneous medium: at radii r small compared with a characteristic radius r*, the mean intensity near line centre varies as 1/ r^(7/3), while at r > r* it approaches 1/ r^2; for r << r* it is modified by frequency redistribution. Before the reionization epoch, r* takes on the universal value 1.1 Mpc, independent of redshift. The mean intensity and scattering rate are found to be very sensitive to the gradient of the velocity field, growing exponentially with the amplitude of the perturbation as the limit of a vanishing velocity gradient is approached near the source. We expect the 21cm signal from the Epoch of Reionization to thus be a sensitive probe of both the density and the peculiar velocity fields.

With the coming data deluge from synoptic surveys, there is a growing need for frameworks that can quickly and automatically produce calibrated classification probabilities for newly-observed variables based on a small number of time-series measurements. In this paper, we introduce a methodology for variable-star classification, drawing from modern machine-learning techniques. We describe how to homogenize the information gleaned from light curves by selection and computation of real-numbered metrics ("feature"), detail methods to robustly estimate periodic light-curve features, introduce tree-ensemble methods for accurate variable star classification, and show how to rigorously evaluate the classification results using cross validation. On a 25-class data set of 1542 well-studied variable stars, we achieve a 22.8% overall classification error using the random forest classifier; this represents a 24% improvement over the best previous classifier on these data. This methodology is effective for identifying samples of specific science classes: for pulsational variables used in Milky Way tomography we obtain a discovery efficiency of 98.2% and for eclipsing systems we find an efficiency of 99.1%, both at 95% purity. We show that the random forest (RF) classifier is superior to other machine-learned methods in terms of accuracy, speed, and relative immunity to features with no useful class information; the RF classifier can also be used to estimate the importance of each feature in classification. Additionally, we present the first astronomical use of hierarchical classification methods to incorporate a known class taxonomy in the classifier, which further reduces the catastrophic error rate to 7.8%. Excluding low-amplitude sources, our overall error rate improves to 14%, with a catastrophic error rate of 3.5%.

We show that, due to the high optical depth of the intergalactic medium to Lyman-alpha photons before the Epoch of Reionization, the Lyman-alpha scattering rate responsible for the Wouthuysen-Field effect from an isolated source will be negligible unless (1) there is sufficient time for the scattering photons to establish a steady state, or (2) the scattering gas is undergoing internal expansion or has a peculiar motion of tens to hundreds of km/s away from the source. We present steady-state solutions in the radiative diffusion approximation for the radiation field trapped in a clump of gas and show that this may result in an enhancement, by a factor of up to 10^6, of the strength of the Wouthuysen-Field effect over that obtained from the free-streaming limit. Solutions to the time-dependent diffusion equation, however, suggest that the timescales required to reach such a steady state will generally exceed the source lifetimes. In the presence of internal expansion, a steady state may be established as photons are redshifted into the red wing, and significant enhancement in the scattering rate may again be produced. Alternatively, a substantial scattering rate may arise in systems with a peculiar motion away from the source that redshifts the received radiation into the resonance line centre. As a consequence, at epochs z<30, when collisional decoupling is small except in dense regions, and prior to the establishment of any large-scale diffuse radiation field of resonance line photons, the 21cm signature from the Intergalactic Medium produced by the Wouthuysen-Field effect will in general trace the peculiar velocity field of the gas in addition to its density structure.