2014A&A...565A.105B

Transit observations in the MgI line of HD209458b revealed signatures of neutral magnesium escaping the upper atmosphere of the planet, while no atmospheric absorption was found in the MgII doublet. Here we present a 3D particle model of the dynamics of neutral and ionized magnesium populations, coupled with an analytical modeling of the atmosphere below the exobase. Theoretical MgI absorption line profiles are directly compared with the absorption observed in the blue wing of the line during the planet transit. Observations are well-fitted with an escape rate of neutral magnesium {dot}(M)_Mg0=2.9⁺0.5-0.9_x10⁷g/s, an exobase close to the Roche lobe R_exo=3^+1.3_–0.9R_p, where R_p is the planet radius) and a planetary wind velocity at the exobase v_pl–wind=25km/s. The observed velocities of the planet-escaping magnesium up to -60km/s are well explained by radiation pressure acceleration, provided that UV-photoionization be compensated for by electron recombination up to ∼13R_p. If the exobase properties are constrained to values given by theoretical models of the deeper atmosphere (R_exo=2R_p and v_pl–wind=10km/s), the best fit to the observations is found at a similar electron density and escape rate within 2σ. In all cases, the mean temperature of the atmosphere below the exobase must be higher than ∼6100K. Simulations predict a redward expansion of the absorption profile from the beginning to the end of the transit. The spatial and spectral structure of the extended atmosphere is the result of complex interactions between radiation pressure, planetary gravity, and self-shielding, and can be probed through the analysis of transit absorption profiles in the Mgi line.