Lyman alpha radiative transfer at the epoch of cosmic reionization

We present a numerical code for the random scattering histories of Lyman alpha photons in the intergalactic medium. The numerical code calculates the radiative transfer under generic three dimensional density, ionization fraction, and peculiar velocity fields based on N-body + radiation transfer simulations of the epoch of reionization. The code is tested with models having analytical solutions, which have idealized geometry and simplified velocity fields. The emergent line profiles can give constraints to the ionization structure around Lyman alpha sources in the early universe.