2008A&A...478..543Z

We calculate the hydrogen and helium ionization in B[e] envelopes and explore their dependence on mass-loss and effective temperature. We also present simulated observations of the Hα emission line and the CIV λλ1550 doublet, and study their behavior. This paper reports our first results in an ongoing study of B[e] supergiants, and provides a glimpse on the ionization of the most important elements in self-consistent numerical simulations. Our newly developed 2D stellar atmosphere code, ASTAROTH, was used for the numerical simulations. The code self-consistently solves for the continuum radiation, non-LTE level populations, and electron temperature in axi-symmetric stellar envelopes. Observed profiles were calculated by an auxiliary program developed separately from ASTAROTH. In all but one of our models, H remained fully ionized - only for {dot}(M)>10^–5M_☉/yr and T_eff≤18000K did we obtain a neutral H disk, and then only for radii beyond 3R_*. Due to ionizations from excited states it is much more difficult to get a H neutral disk than indicated by previous analytical calculations. Near the poles, the ionization is high in all models, while helium recombined in the equatorial regions for all but our lowest mass-loss rate (10^–6M_☉/yr). Although the model parameters were not adjusted to provide fits to any particular star, the theoretical profiles show some features seen in the profiles of R126. These include the partially resolved double peaked profile of Hα, and the weak emission associated with the UV CIV resonance line.