A Model of the Hα and Na Transmission Spectrum of HD 189733b

This paper presents a detailed hydrostatic model of the upper atmosphere of HD 189733b, with the goal of constraining its temperature, particle densities, and radiation field over the pressure range 10(-4) to 10/mu bar, where the observed Ha transmission spectrum is produced. The atomic hydrogen level population is computed including both collisional and radiative transition rates. The Ly alpha resonant scattering is computed using a Monte Carlo simulation. The model transmission spectra are in broad agreement with the data. Excitation of the H(2l) population is mainly by Ly alpha radiative excitation due to the large Ly alpha intensity. The density of H(2f) is nearly flat over two decades in pressure and is optically thick to H alpha. Additional models computed for a range of the stellar Lyman continuum (LyC) flux suggest that the variability in H alpha transit depth may be due to the variability in the stellar LyC. Since metal lines provide the dominant cooling of this part of the atmosphere, the atmosphere structure is sensitive to the density of species such as Mg and Na, which may themselves be constrained by observations. Since the H alpha and Na D lines have comparable absorption depths, we argue that the centers of the Na D lines are also formed in the atomic layer where the H alpha line is formed.