A tutorial associated with the distribution of 2D non-LTE multilevel atom radiative transfer codes.

HD 5980 is a unique system containing one massive star (star A) that is apparently entering the luminous blue variable phase, and an eclipsing companion (star B) that may have already evolved beyond this phase to become a Wolf-Rayet star. In this paper we present the results from FUSE observations obtained in 1999, 2000, and 2002 and one far-UV observation obtained by ORFEUS/BEFS in 1993 shortly before the first eruption of HD 5980. The eight phase-resolved spectra obtained by FUSE in 2002 are analyzed in the context of a wind-eclipse model. This analysis shows that the wind of the eruptor obeyed a very fast velocity law in 2002, which is consistent with the line-driving mechanism. Large amplitude line-profile variations on the orbital period are shown to be due to the eclipse of star B by the wind of star A, although the eclipse due to gas flowing in the direction of star B is absent. This can only be explained if the wind of star A is not spherically symmetric, or if the eclipsed line radiation is "filled-in" by emission originating from somewhere else in the system, e.g., in the wind-wind collision region. Except for a slightly lower wind speed, the ORFEUS/BEFS spectrum is very similar to the spectrum obtained by FUSE at the same orbital phase: there is no indication of the impending eruption. However, the trend for decreasing wind velocity suggests the occurrence of the "bi-stability" mechanism, which in turn implies that the restructuring of the circumbinary environment caused by the transition from "fast, rarefied wind" to "slow, dense wind" was observed as the eruptive event. The underlying mechanism responsible for the long-term decrease in wind velocity that precipitated this change remains an open issue.

We attempt to fit observations with 5" resolution of the J=2-1 transition of CS in the directions of H II regions A, B, and G of W49A North as well as observations with 20" resolution of the J=2-1, 3-2, 5-4, and 7-6 transitions in the directions of H II regions A and G by using radiative transfer calculations. These calculations predict the intensity profiles resulting from several spherical clouds along the line of sight. We consider three models: global collapse of a very large (5 pc radius) cloud, localized collapse from smaller (1 pc) clouds around individual H II regions, and multiple, static clouds. For all three models we can find combinations of parameters that reproduce the CS profiles reasonably well provided that the component clouds have a core-envelope structure with a temperature gradient. Cores with high temperature and high molecular hydrogen density are needed to match the higher transitions (e.g. J=7-6) observed towards A and G. The lower temperature, low density gas needed to create the inverse P-Cygni profile seen in the CS J=2-1 line (with 5" beam) towards H II region G arises from different components in the 3 models. The infalling envelope of cloud G plus cloud B creates the absorption in global collapse, cloud B is responsible in local collapse, and a separate cloud, G', is needed in the case of many static clouds. The exact nature of the velocity field in the envelopes for the case of local collapse is not important as long as it is in the range of 1 to 5 km/s for a turbulent velocity of about 6 km/s. High resolution observations of the J=1-0 and 5-4 transitions of CS and C34S may distinguish between these three models. Modeling existing observations of HCO+ and C18O does not allow one to distinguish between the three models but does indicate the existence of a bipolar outflow.