Methods of calculation of NLTE model atmosphere are discussed. The NLTE trace element procedure is compared with the full NLTE model atmosphere calculation. Differences between LTE and NLTE atmosphere modeling are evaluated. The ways of model atom construction are discussed. Finally, modelling of expanding atmospheres of hot stars with winds is briefly reviewed.

Recent studies of O-type stars demonstrated that discrepant mass-loss rates are obtained when different diagnostic methods are employed - fitting the unsaturated UV resonance lines (e.g. P v) gives drastically lower values than obtained from the H{\alpha} emission. Wind clumping may be the main cause for this discrepancy. In a previous paper, we have presented 3-D Monte-Carlo calculations for the formation of scattering lines in a clumped stellar wind. In the present paper we select five O-type supergiants (from O4 to O7) and test whether the reported discrepancies can be resolved this way. In the first step, the analyses start with simulating the observed spectra with Potsdam Wolf-Rayet (PoWR) non-LTE model atmospheres. The mass-loss rates are adjusted to fit best to the observed H{\alpha} emission lines. For the unsaturated UV resonance lines (i.e. P v) we then apply our 3-D Monte-Carlo code, which can account for wind clumps of any optical depths, a non-void inter-clump medium, and a velocity dispersion inside the clumps. The ionization stratifications and underlying photospheric spectra are adopted from the PoWR models. From fitting the observed resonance line profiles, the properties of the wind clumps are constrained. Our results show that with the mass-loss rates that fit H{\alpha} (and other Balmer and He II lines), the UV resonance lines (especially the unsaturated doublet of P v) can also be reproduced without problem when macroclumping is taken into account. There is no need to artificially reduce the mass-loss rates, nor to assume a sub-solar phosphorus abundance or an extremely high clumping factor, contrary to what was claimed by other authors. These consistent mass-loss rates are lower by a factor of 1.3 to 2.6, compared to the mass-loss rate recipe from Vink et al. Macroclumping resolves the previously reported discrepancy between H{\alpha} and P v mass-loss diagnostics.

Clumping in stellar winds of hot stars is a possible consequence of radiative-acoustic instability appearing in solutions of radiative-hydrodynamical equations. However, clumping is usually included to stellar atmosphere modeling and radiative transfer calculations in a highly approximate way via a global free parameter called the clumping factor. Using different values of clumping factors many researchers succeeded to fit the observed spectra better and to correct empirical mass loss rates. This usually leads to a conclusion that the stellar wind is clumped. To understand how clumping may influence theoretical predictions of mass-loss rates, different clumping properties have to be taken into account. If clumping appears already below the critical point, the mass-loss rates is changed.

The true mass-loss rates from massive stars are important for many branches of astrophysics. For the correct modeling of the resonance lines, which are among the key diagnostics of stellar mass-loss, the stellar wind clumping turned out to be very important. In order to incorporate clumping into radiative transfer calculation, 3-D models are required. Various properties of the clumps may have strong impact on the resonance line formation and, therefore, on the determination of empirical mass-loss rates. We incorporate the 3-D nature of the stellar wind clumping into radiative transfer calculations and investigate how different model parameters influence the resonance line formation. We develop a full 3-D Monte Carlo radiative transfer code for inhomogeneous expanding stellar winds. The number density of clumps follows the mass conservation. For the first time, realistic 3-D models that describe the dense as well as the tenuous wind components are used to model the formation of resonance lines in a clumped stellar wind. At the same time, non-monotonic velocity fields are accounted for. The 3-D density and velocity wind inhomogeneities show very strong impact on the resonance line formation. The different parameters describing the clumping and the velocity field results in different line strengths and profiles. We present a set of representative models for various sets of model parameters and investigate how the resonance lines are affected. Our 3-D models show that the line opacity is reduced for larger clump separation and for more shallow velocity gradients within the clumps. Our new model demonstrates that to obtain empirically correct mass-loss rates from the UV resonance lines, the wind clumping and its 3-D nature must be taken into account.

We study the effects of optically thick clumps, non-void inter-clump medium, variation of the onset of clumping, and velocity dispersion inside clumps on the formation of resonance lines. For this purpose we developed a full 3-D Monte Carlo Radiative Transfer (MCRT) code that is able to handle 3-D shapes of clumps and arbitrary 3-D velocity fields. The method we developed allows us to take into account contributions from density and velocity wind inhomogeneities to the total opacity very precisely. The first comparison with observation shows that 3-D density and velocity wind inhomogeneities have a very strong influence on the resonance line formation, and that they have to be accounted for in order to obtain reliable mass-loss rate determinations.

We have studied the Residual Blaze Functions (RBF) resulting from division of individual echelle orders by extracted flat-field in spectra obtained by slit-fed OES spectrograph of 2m telescope of Ond\v{r}ejov observatory, Czech Republic. We have eliminated the dependence on target and observation conditions by semiautomatic fitting of global response function, thus getting the instrument-only dependent part, which may be easily incorporated into data reduction pipeline. The improvement of reliability of estimation of continuum on spectra of targets with wide and shallow lines is noticeable and the merging of all orders into the one long spectrum gives much more reliable results.