Context. Metal-poor massive stars are supposed to be progenitors of certain supernovae, gamma-ray bursts and compact object mergers, potentially contributing to the early epochs of the Universe with their strong ionizing radiation. However, they remain mainly theoretical as individual spectroscopic observations of such objects have rarely been carried out below the metallicity of the SMC. Aims. This work aims at exploring what our state-of-the-art theories of stellar evolution combined with those of stellar atmospheres predict about a certain type of metal-poor (0.02 Z$_{\odot}$) hot massive stars, the chemically homogeneously evolving ones, called TWUIN stars. Methods. Synthetic spectra corresponding to a broad range in masses (20-130 M$_{\odot}$) and covering several evolutionary phases from the zero-age main-sequence up to the core helium-burning stage were computed. Results. We find that TWUIN stars show almost no emission lines during most of their {core hydrogen-burning} lifetimes. Most metal lines are completely absent, including nitrogen. During their core helium-burning stage, lines switch to emission and even some metal lines (oxygen and carbon, but still almost no nitrogen) show up. Mass loss and clumping play a significant role in line-formation in later evolutionary phases, particularly during core helium-burning. Most of our spectra are classified as an early O type giant or supergiant, and we find Wolf-Rayet stars of type WO in the core helium-burning phase. Conclusions. An extremely hot, early O type star observed in a low-metallicity galaxy could be the outcome of chemically homogeneous evolution $-$ and therefore the progenitor of a long-duration gamma-ray burst or a type Ic supernova. TWUIN stars may play an important role in reionizing the Universe due to their being hot without showing prominent emission lines during the majority of their lifetimes.

{We aim to demonstrate the effect of atmospheric inhomogeneities on the emergent specific intensity and radiation flux of a spectral line radiation.} {We self-consistently solve the NLTE problem for a two-level atom in a 3D atmosphere using the Cartesian grid. For that purpose, we use the 3D radiative transfer code PORTA. By examining simple examples, we study cases where the temperature inhomogeneities in the atmosphere models lead to the modification of the emergent radiation.} {We show that specific temperature inhomogeneities in the model atmospheres influence the emerging radiation, and that interpretation of the stellar spectra based on a plane-parallel atmosphere models can lead to erroneous conclusions about the atmospheric structure.}

Eclipsing binaries with a $\delta$ Sct component are powerful tools to derive the fundamental parameters and probe the internal structure of stars. In this study, spectral analysis of 6 primary $\delta$ Sct components in eclipsing binaries has been performed. Values of $T_{\rm eff}$, $v \sin i$, and metallicity for the stars have been derived from medium-resolution spectroscopy. Additionally, a revised list of $\delta$ Sct stars in eclipsing binaries is presented. In this list, we have only given the $\delta$ Sct stars in eclipsing binaries to show the effects of the secondary components and tidal-locking on the pulsations of primary $\delta$ Sct components. The stellar pulsation, atmospheric and fundamental parameters (e.g., mass, radius) of 92 $\delta$ Sct stars in eclipsing binaries have been gathered. Comparison of the properties of single and eclipsing binary member $\delta$ Sct stars has been made. We find that single $\delta$ Sct stars pulsate in longer periods and with higher amplitudes than the primary $\delta$ Sct components in eclipsing binaries. The $v \sin i$ of $\delta$ Sct components is found to be significantly lower than that of single $\delta$ Sct stars. Relationships between the pulsation periods, amplitudes, and stellar parameters in our list have been examined. Significant correlations between the pulsation periods and the orbital periods, $T_{\rm eff}$, $\log g$, radius, mass ratio, $v \sin i$, and the filling factor have been found.

Rayleigh scattering is a result of an interaction of photons with bound electrons. Rayleigh scattering is mostly neglected in calculations of hot star model atmospheres because most of the hydrogen atoms are ionized and the heavier elements have a lower abundance than hydrogen. In atmospheres of some chemically peculiar stars, helium overabundant regions containing singly ionized helium are present and Rayleigh scattering can be a significant opacity source. We evaluate the contribution of Rayleigh scattering by neutral hydrogen and singly ionized helium in the atmospheres of hot stars with solar composition and in the atmospheres of helium overabundant stars. We computed several series of model atmospheres using the TLUSTY code and emergent fluxes using the SYNSPEC code. These models describe atmospheres of main sequence B-type stars with different helium abundance. We used an existing grid of models for atmospheres with solar chemical composition and we calculated an additional grid for helium-rich stars with $N$(He)/$N$(H)=10. Rayleigh scattering by neutral hydrogen can be neglected in atmospheres of hot stars, while Rayleigh scattering by singly ionized helium can be a non-negligible opacity source in some hot stars, especially in helium-rich stars.

During the summer of 2013, a 4-month spectroscopic campaign took place to observe the variabilities in three Wolf-Rayet stars. The spectroscopic data have been analyzed for WR 134 (WN6b), to better understand its behaviour and long-term periodicity, which we interpret as arising from corotating interaction regions (CIRs) in the wind. By analyzing the variability of the He II $\lambda$5411 emission line, the previously identified period was refined to P = 2.255 $\pm$ 0.008 (s.d.) days. The coherency time of the variability, which we associate with the lifetime of the CIRs in the wind, was deduced to be 40 $\pm$ 6 days, or $\sim$ 18 cycles, by cross-correlating the variability patterns as a function of time. When comparing the phased observational grayscale difference images with theoretical grayscales previously calculated from models including CIRs in an optically thin stellar wind, we find that two CIRs were likely present. A separation in longitude of $\Delta \phi \simeq$ 90$^{\circ}$ was determined between the two CIRs and we suggest that the different maximum velocities that they reach indicate that they emerge from different latitudes. We have also been able to detect observational signatures of the CIRs in other spectral lines (C IV $\lambda\lambda$5802,5812 and He I $\lambda$5876). Furthermore, a DAC was found to be present simultaneously with the CIR signatures detected in the He I $\lambda$5876 emission line which is consistent with the proposed geometry of the large-scale structures in the wind. Small-scale structures also show a presence in the wind, simultaneously with the larger scale structures, showing that they do in fact co-exist.

Stable shell stars are ideal objects for studying basic physical principles of the formation of disks in Be stars. If these stars have a close unresolved visual companion, its contribution toward the modelling of the disk cannot be overlooked, as is sometimes done. The study aims to spectroscopically resolve close visual binary Be (shell) star 1 Del, which up to now was only resolved by speckle or micrometric measurements. The integral field spectroscopy obtained by the SINFONI spectrograph at the VLT telescope in the European Southern Observatory (ESO) in the infrared region was used; we supplemented these observations with visual spectroscopy with the Perek Telescope at the Ondrejov Observatory. Spectra of 1 Del were successfully resolved, and, for the first time, spectra of 1 Del B were obtained. We found that 1 Del A is a Be/shell star, while 1 Del B is not an emission-line object.

Various available codes for NLTE modeling and analysis of hot star spectra are reviewed. Generalizations of standard equations of kinetic equilibrium and their consequences are discussed.

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.

Basic assumptions of the NLTE approximation in stellar atmospheres are summarized. The assumptions of thermodynamic equilibrium, local thermodynamic equilibrium (LTE), and non-LTE (NLTE) are compared. It is emphasized that LTE is a poor approximation if radiative transitions dominate in stellar atmospheres. The equations of kinetic equilibrium and methods of their solution are discussed.

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.

We present results of our calculations of NLTE model stellar atmospheres for hot Population III stars composed of hydrogen and helium. We use our own computer code for calculation of spherically symmetric NLTE model atmospheres in hydrostatic and radiative equilibrium. The model atmospheres are then used for calculation of emergent fluxes. These fluxes serve for the evaluation of the flow of high-energy photons for energies higher than ionization energies of hydrogen and helium, the so-called ionizing photon fluxes. We also present the time evolution of the ionizing photon fluxes.

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.

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.

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 analyse new spectra of the multiple system SZ Cam because previous studies found different values of the primary radial velocity amplitude. The older solutions of light curves also have different ratios of secondary to primary luminosity as inferred from the observed equivalent widths of spectral lines. We therefore reanalyse the light curves of the eclipsing pair. Only the light curve derived by Wesselink has a solution that agrees with the observed equivalent width ratio. The resulting parameters of the binary are discussed. Masses of $M_1=16.6$ and $M_2=11.9$ M$_{\odot}$, and radii $R_1=9.4$ and $R_2=5.4$ R$_{\odot}$ are derived. We point out that radial velocities measured with the CCF method can be misleading when the method is applied to multiple systems with complex line blends. New radial velocities are also obtained for the visual component ADS 2984 A (HD 25639).

We studied both components of a slightly overlooked visual binary HR 1847 spectroscopically to determine its basic physical and orbital parameters. Basic stellar parameters were determined by comparing synthetic spectra to the observed echelle spectra, which cover both the optical and near-IR regions. New observations of this system used the Ond\v{r}ejov and Rozhen 2-m telescopes and their coud\'e spectrographs. Radial velocities from individual spectra were measured and then analysed with the code {\FOTEL} to determine orbital parameters. The spectroscopic orbit of HR 1847A is presented for the first time. It is a single-lined spectroscopic binary with a B-type primary, a period of 719.79 days, and a highly eccentric orbit with e=0.7. We confirmed that HR 1847B is a Be star. Its H\alpha emission significantly decreased from 2003 to 2008. Both components have a spectral type B7-8 and luminosity class IV-V.

Radial-velocity variations of the H-alpha emission measured on the steep wings of the H-alpha line, prewhitened for the long-time changes, vary periodically with a period of (218.025 +/- 0.022)d, confirming the suspected binary nature of the bright Be star Pleione, a member of the Pleiades cluster. The orbit seems to have a high eccentricity over 0.7, but we also briefly discuss the possibility that the true orbit is circular and that the eccentricity is spurious owing to the phase-dependent effects of the circumstellar matter. The projected angular separation of the spectroscopic orbit is large enough to allow the detection of the binary with large optical interferometers, provided the magnitude difference primary - secondary is not too large. Since our data cover the onset of a new shell phase up to development of a metallic shell spectrum, we also briefly discuss the recent long-term changes. We confirm the formation of a new envelope, coexisting with the previous one, at the onset of the new shell phase. We find that the full width at half maximum of the H-alpha profile has been decreasing with time for both envelopes. In this connection, we briefly discuss Hirata's hypothesis of precessing gaseous disk and possible alternative scenarios of the observed long-term changes.

We have developed a time-dependent two-component hydrodynamics code to simulate radiatively-driven stellar winds from hot stars. We use a time-explicit van Leer scheme to solve the hydrodynamic equations of a two-component stellar wind. Dynamical friction due to Coulomb collisions between the passive bulk plasma and the line-scattering ions is treated by a time-implicit, semi-analytic method using a polynomial fit to the Chandrasekhar function. This gives stable results despite the stiffness of the problem. This method was applied to model stars with winds that are both poorly and well-coupled. While for the former case we reproduce the mCAK solution, for the latter case our solution leads to wind decoupling.

A new method for the formal solution of the 2D radiative transfer equation in axial symmetry in the presence of arbitrary velocity fields is presented. The combination of long and short characteristics methods is used to solve the radiative transfer equation. We include the velocity field in detail using the Local Lorentz Transformation. This allows us to obtain a significantly better description of the photospheric region, where the gradient of the global velocity is too small for the Sobolev approximation to be valid. Sample test calculations for the case of a stellar wind and a rotating atmosphere are presented.