Emission features formed through Raman scattering with atomic hydrogen provide unique and crucial information to probe the distribution and kinematics of a thick neutral region illuminated by a strong far-ultraviolet radiation source. We introduce a new 3-dimensional Monte-Carlo code in order to describe the radiative transfer of line photons that are subject to Raman and Rayleigh scattering with atomic hydrogen. In our Sejong Radiative Transfer through Raman and Rayleigh Scattering (STaRS) code, the position, direction, wavelength, and polarization of each photon is traced until escape. The thick neutral scattering region is divided into multiple cells with each cell being characterized by its velocity and density, which ensures exibility of the code in analyzing Raman-scattered features formed in a neutral region with complicated kinematics and density distribution. To test the code, we revisit the formation of Balmer wings through Raman scattering of the far-UV continuum near Lyβ and Lyγ in a static neutral region. An additional check is made to investigate Raman scattering of Ovi in an expanding neutral medium. We find a good agreement of our results with previous works, demonstrating the capability of dealing with radiative transfer modeling that can be applied to spectropolarimetric imaging observations of various objects including symbiotic stars, young planetary nebulae, and active galactic nuclei.

We investigate the escape of Lyβ from emission nebulae with a significant population of excited hydrogen atoms in the level n = 2, rendering them optically thick in Hα. The transfer of Lyβ line photons in these optically thick regions is complicated by the presence of another scattering channel leading to re-emission of Hα, alternating their identities between Lyβ and Hα. In this work, we develop a Monte Carlo code to simulate the transfer of Lyβ line photons incorporating the scattering channel into Hα. Both Hα and Lyβ lines are formed through diffusion in frequency space, where a line photon enters the wing regime after a fairly large number of resonance scatterings with hydrogen atoms. Various line profiles of Hα and Lyβ emergent from our model nebulae are presented. It is argued that the electron temperature is a critical parameter which controls the flux ratio of emergent Lyβ and Hα. Specifically for T = 3 × 104 K and Hα line center optical depth α = 10, the number flux ratio of emergent Lyβ and Hα is ∼ 49 percent, which is quite significant. We propose that the leaking Lyβ can be an interesting source for the formation of Hα wings observed in many symbiotic stars and active galactic nuclei. Similar broad Hα wings are also expected in Lyα emitting halos found in the early universe, which can be potentially probed by the James Webb Telescope in the future.