A Novel Solution for Resonant Scattering Using Self-Consistent Boundary Conditions

Abstract

We present two novel additions to the semi-analytic solution of Lyman $\alpha$ (Ly$\alpha$) radiative transfer in spherical geometry: (1) implementation of the correct boundary condition for a steady source, and (2) solution of the time-dependent problem for an impulsive source. For the steady-state problem, the solution can be represented as a sum of two terms: a previously-known analytic solution of the equation with mean intensity $J=0$ at the surface, and a novel, semi-analytic solution which enforces the correct boundary condition of zero-ingoing intensity at the surface. This solution is compared to that of the Monte Carlo method, which is valid at arbitrary optical depth. It is shown that the size of the correction is of order unity when the spectral peaks approach the Doppler core and decreases slowly with line center optical depth, specifically as $(a{\tau }_0{)}^{-1/3}$, which may explain discrepancies seen in previous studies. For the impulsive problem, the time, spatial, and frequency dependence of the solution are expressed using an eigenfunction expansion in order to characterize the escape time distribution and emergent spectra of photons. It is shown that the lowest-order eigenfrequency agrees well with the decay rate found in the Monte Carlo escape time distribution at sufficiently large line-center optical depths. The characterization of the escape-time distribution highlights the potential for a Monte Carlo acceleration method, which would sample photon escape properties from distributions rather than calculating every photon scattering, thereby reducing computational demand.

23 pages, 10 figures, 2 appendices, resubmitted to ApJ after review

References in corpus (16)

• Array Programming with NumPy
• A giant comet-like cloud of hydrogen escaping the warm Neptune-mass exoplanet GJ 436b
• Radiative transfer in very optically thick circumstellar disks
• Hubble PanCET: An extended upper atmosphere of neutral hydrogen around the warm Neptune GJ 3470 b
• Reconnaissance of the TRAPPIST-1 exoplanet system in the Lyman-$\alpha$ line
• The Lyman-alpha signature of the first galaxies
• A Model of the H$\alpha$ and Na Transmission Spectrum of HD 189733b
• No hydrogen exosphere detected around the super-Earth HD97658 b
• Lyman-alpha radiation hydrodynamics of galactic winds before cosmic reionization
• Is pi Men c's atmosphere hydrogen-dominated? Insights from a non-detecton of Hy Ly-alpha absorption
• Ly-alpha Radiative Transfer: Monte-Carlo Simulation of the Wouthuysen-Field Effect
• Lyman-alpha in the GJ 1132 System: Stellar Emission and Planetary Atmospheric Evolution
• Resonant-line radiative transfer within power-law density profiles
• Lyman-alpha radiation pressure: an analytical exploration
• Atmosphere escape inferred from modelling the H$\alpha$ transmission spectrum of WASP-121b
• The Hubble PanCET program: Long-term chromospheric evolution and flaring activity of the M dwarf host GJ 3470