Ultraprecise space photometry enables us to reveal light variability even in stars that were previously deemed constant. A large group of such stars show variations that may be rotationally modulated. This type of light variability is of special interest because it provides precise estimates of rotational rates. We aim to understand the origin of the light variability of K2 targets that show signatures of rotational modulation. We used phase-resolved medium-resolution XSHOOTER spectroscopy to understand the light variability of the stars KIC~250152017 and KIC~249660366, which are possibly rotationally modulated. We determined the atmospheric parameters at individual phases and tested the presence of the rotational modulation in the spectra. KIC 250152017 is a HgMn star, whose light variability is caused by the inhomogeneous surface distribution of manganese and iron. It is only the second HgMn star whose light variability is well understood. KIC 249660366 is a He-weak, high-velocity horizontal branch star with overabundances of silicon and argon. The light variability of this star is likely caused by a reflection effect in this post-common envelope binary.

Fast line-driven stellar winds play an important role in the evolution of planetary nebulae. We provide global hot star wind models of central stars of planetary nebulae. The models predict wind structure including the mass-loss rates, terminal velocities, and emergent fluxes from basic stellar parameters. We applied our wind code for parameters corresponding to evolutionary stages between the asymptotic giant branch and white dwarf phases. We study the influence of metallicity and wind inhomogeneities (clumping) on the wind properties. Line-driven winds appear very early after the star leaves the asymptotic giant branch (at the latest for $T_\rm{eff}\approx10\,$kK) and fade away at the white dwarf cooling track (below $T_\rm{eff}=105\,$kK). Their mass-loss rate mostly scales with the stellar luminosity and, consequently, the mass-loss rate only varies slightly during the transition from the red to the blue part of the Hertzsprung-Russell diagram. There are the following two exceptions to the monotonic behavior: a bistability jump at around $20\,$kK, where the mass-loss rate decreases by a factor of a few (during evolution) due to a change in iron ionization, and an additional maximum at about $T_\rm{eff}=40-50\,$kK. On the other hand, the terminal velocity increases from about a few hundreds of $\rm{km}\,\rm{s}^{-1}$ to a few thousands of $\rm{km}\,\rm{s}^{-1}$ during the transition as a result of stellar radius decrease. The wind terminal velocity also significantly increases at the bistability jump. Derived wind parameters reasonably agree with observations. The effect of clumping is stronger at the hot side of the bistability jump than at the cool side. Derived fits to wind parameters can be used in evolutionary models and in studies of planetary nebula formation. A predicted bistability jump in mass-loss rates can cause the appearance of an additional shell of planetary nebula.

We present new ATCA multi-wavelength radio measurements (range 2.1-21.2 GHz) of the early-type magnetic star rho Oph A, performed in March 2019 during 3 different observing sessions. These new ATCA observations evidence a clear rotational modulation of the stellar radio emission and the detection of coherent auroral radio emission from rho Oph A at 2.1 GHz. We collected high-resolution optical spectra of rho Oph A acquired by several instruments over a time span of about ten years. We also report new magnetic field measurements of rho Oph A that, together with the radio light curves and the temporal variation of the equivalent width of the HeI line (lambda=5015 Angstrom), were used to constrain the rotation period and the stellar magnetic field geometry. The above results have been used to model the stellar radio emission, modelling that allowed us to constrain the physical condition of rho Oph A's magnetosphere. Past XMM measurements showed periodic X-ray pulses from rho Oph A. We correlate the X-ray light curve with the magnetic field geometry of rho Oph A. The already published XMM data have been re-analyzed showing that the X-ray spectra of rho Oph A are compatible with the presence of a non-thermal X-ray component. We discuss a scenario where the emission phenomena occurring at the extremes of the electromagnetic spectrum, radio and X-ray, are directly induced by the same plasma process. We interpret the observed X-ray and radio features of rho Oph A as having an auroral origin.

Helium-peculiar magnetic chemically peculiar stars show variations of helium abundance across their surfaces. As a result of associated atmospheric scale height variations, the stellar surface becomes distorted, with helium-rich regions dented inwards. Effectively, on top of flux variations due to opacity effects, the depressed helium-rich surface regions become less bright in the optical regions and brighter in the ultraviolet. We study the observational effects of the aspherical surface on the light curves of a Cen. We simulate the light curves of this star adopting surface distributions of He, N, O, Si, and Fe derived from Doppler mapping and introducing the effect of distortion proportional to helium abundance. We show that while most of the optical and UV variations of this star result from flux redistribution due to the non-uniform surface distributions of helium and iron, the reduction of light variations due to the helium-related surface distortion leads to a better agreement between simulated optical light curves and the light curves observed with the BRITE satellites.

Hot subdwarfs are helium burning objects in late stages of their evolution. These stars can develop winds driven by light absorption in the lines of heavier elements. The wind strength depends on chemical composition which can significantly vary from star to star. We aim to understand the influence of metallicity on the strength of the winds of the hot hydrogen-rich subdwarfs HD 49798 and BD+18$^\circ\,$2647. We used UV and optical spectra to derive stellar parameters and abundances. For derived stellar parameters, we predicted wind structure (including mass-loss rates and terminal velocities) with our METUJE code. We derived effective temperature $T_\text{eff}=45\,900\,$K and mass $M=1.46\,M_\odot$ for HD 49798 and $T_\text{eff}=73\,000\,$K and $M=0.38\,M_\odot$ for BD+18$^\circ\,$2647. The abundances can be interpreted as a result of interplay between stellar evolution and diffusion. HD 49798 has a strong wind that does not allow for chemical separation and consequently it shows solar chemical composition modified by hydrogen burning. On the other hand, we did not find any wind in BD+18$^\circ\,$2647 and its abundances are therefore most likely affected by radiative diffusion. Accurate abundances do not lead to a significant modification of wind mass-loss rate for HD 49798, because the increase of the contribution of Fe and Ni to the radiative force is compensated by the decrease of the force due to other elements. The resulting wind mass-loss rate $\dot M=2.1\times10^{-9}\,M_\odot\,\text{yr}^{-1}$ predicts an X-ray light curve during the eclipse which closely agrees with observations. On the other hand, the absence of the wind in BD+18$^\circ\,$2647 for accurate abundances is a result of its peculiar chemical composition. Wind models with accurate abundances provide more reliable wind parameters, but the influence of abundances on the wind parameters is limited in many cases.

CU Vir has been the first main sequence star that showed regular radio pulses that persist for decades, resembling the radio lighthouse of pulsars and interpreted as auroral radio emission similar to that found in planets. The star belongs to a rare group of magnetic chemically peculiar stars with variable rotational period. We study the ultraviolet (UV) spectrum of CU Vir obtained using STIS spectrograph onboard the Hubble Space Telescope (HST) to search for the source of radio emission and to test the model of the rotational period evolution. We used our own far-UV and visual photometric observations supplemented with the archival data to improve the parameters of the quasisinusoidal long-term variations of the rotational period. We predict the flux variations of CU Vir from surface abundance maps and compare these variations with UV flux distribution. We searched for wind, auroral, and interstellar lines in the spectra. The UV and visual light curves display the same long-term period variations supporting their common origin. New updated abundance maps provide better agreement with the observed flux distribution. The upper limit of the wind mass-loss rate is about $10^{-12}\,M_\odot\,\rm{yr}^{-1}$. We do not find any auroral lines. We find rotationally modulated variability of interstellar lines, which is most likely of instrumental origin. Our analysis supports the flux redistribution from far-UV to near-UV and visual domains originating in surface abundance spots as the main cause of the flux variability in chemically peculiar stars. Therefore, UV and optical variations are related and the structures leading to these variations are rigidly confined to the stellar surface. The radio emission of CU Vir is most likely powered by a very weak presumably purely metallic wind, which leaves no imprint in spectra.

In wind-powered X-ray binaries, the radiatively driven stellar wind from the primary may be inhibited by the X-ray irradiation. This creates the feedback that limits the X-ray luminosity of the compact secondary. Wind inhibition might be weakened by the effect of small-scale wind inhomogeneities (clumping) possibly affecting the limiting X-ray luminosity. We study the influence of X-ray irradiation on the stellar wind for different radial distributions of clumping. We calculate hot star wind models with external irradiation and clumping using our global wind code. The models are calculated for different parameters of the binary. We determine the parameters for which the X-ray wind ionization leads to a decrease of the radiative force. This causes a decrease of the wind velocity and even of the mass-loss rate in the case of extreme X-ray irradiation. Clumping weakens the effect of X-ray irradiation because it favours recombination and leads to an increase of the wind mass-loss rate. The best match between the models and observed properties of high-mass X-ray binaries (HMXB) is derived with radially variable clumping. We describe the influence of X-ray irradiation on the terminal velocity and on the mass-loss rate in a parametric way. The X-ray luminosities predicted within the Bondi theory agree nicely with observations when accounting for X-ray irradiation. The ionizing feedback regulates the accretion onto the compact companion resulting in a relatively stable X-ray source. The wind-powered accretion model can account for large luminosities in HMXBs only when introducing the ionizing feedback. There are two possible states following from the dependence of X-ray luminosity on the wind terminal velocity and mass-loss rate. One state has low X-ray luminosity and a nearly undisturbed wind, and the second state has high X-ray luminosity and exhibits a strong influence of X-rays on the flow.

Fraction of hot stars posses strong magnetic fields that channel their radiatively driven outflows. We study the influence of line splitting in the magnetic field (Zeeman effect) on the wind properties. We use our own global wind code with radiative transfer in the comoving frame to understand the influence of the Zeeman splitting on the line force. We show that the Zeeman splitting has a negligible influence on the line force for magnetic fields that are weaker than about 100~kG. This means that the wind mass-loss rates and terminal velocities are not affected by the magnetic line splitting for magnetic fields as are typically found on the surface of nondegenerate stars. Neither have we found any strong flux variability that would be due to the magnetically split line blanketing.

A small fraction of the radiative flux emitted by hot stars is absorbed by their winds and redistributed towards longer wavelengths. This effect, which leads also to the heating of the stellar photosphere, is termed wind blanketing. For stars with variable winds, the effect of wind blanketing may lead to the photometric variability. We have studied the consequences of line driven wind instability and wind blanketing for the light variability of O stars. We combined the results of wind hydrodynamic simulations and of global wind models to predict the light variability of hot stars due to the wind blanketing and instability. The wind instability causes stochastic light variability with amplitude of the order of tens of millimagnitudes and a typical timescale of the order of hours for spatially coherent wind structure. The amplitude is of the order of millimagnitudes when assuming that the wind consists of large number of independent concentric cones. The variability with such amplitude is observable using present space borne photometers. We show that the simulated light curve is similar to the light curves of O stars obtained using BRITE and CoRoT satellites.

Stars with B[e] phenomenon comprise a very diverse group of objects in a different evolutionary status. These objects show common spectral characteristics, including presence of Balmer lines in emission, forbidden lines, and strong infrared excess due to the dust. The observations of emission lines indicate the illumination by ultraviolet ionizing source, which is a key part to understand the elusive nature of these objects. We study the ultraviolet variability of many B[e] stars to specify the geometry of the circumstellar environment and its variability. We analyse massive hot B[e] stars from our Galaxy and from Magellanic Clouds. We study the ultraviolet broad-band variability derived from the flux-calibrated data. We determine variations of individual lines and its correlation with the total flux variability. We detected variability of the spectral energy distribution and of the line profiles. The variability has several sources of origin, including the light absorption by the disk, pulsations, LBV-type variations, and eclipses in the case of binaries. The stellar radiation of most of B[e] stars is heavily obscured by circumstellar material. This suggests that the circumstellar material is not present only in the disk but also above its plane. The flux and line variability is consistent with a two-component model of circumstellar environment composed of the dense disk and ionized envelope. The observations of B[e] supergiants show that many of these stars have nearly the same luminosity of about $1.9\times10^{5}\,L_\odot$ and similar effective temperature.

We provide mass-loss rate predictions for O stars from Large and Small Magellanic Clouds. We calculate global (unified, hydrodynamic) model atmospheres of main sequence, giant, and supergiant stars for chemical composition corresponding to Magellanic Clouds. The models solve radiative transfer equation in comoving frame, kinetic equilibrium equations (also known as NLTE equations), and hydrodynamical equations from (quasi-)hydrostatic atmosphere to expanding stellar wind. The models allow us to predict wind density, velocity, and temperature (consequently also the terminal wind velocity and the mass-loss rate) just from basic global stellar parameters. As a result of their lower metallicity, the line radiative driving is weaker leading to lower wind mass-loss rates with respect to the Galactic stars. We provide a formula that fits the mass-loss rate predicted by our models as a function of stellar luminosity and metallicity. On average, the mass-loss rate scales with metallicity as $ \dot M\sim Z^{0.59}$. The predicted mass-loss rates are lower than mass-loss rates derived from H$\alpha$ diagnostics and can be reconciled with observational results assuming clumping factor $C_\text{c}=9$. On the other hand, the predicted mass-loss rates either agree or are slightly higher than the mass-loss rates derived from ultraviolet wind line profiles. The calculated \ion{P}{v} ionization fractions also agree with values derived from observations for LMC stars with $T_\text{eff}\leq40\,000\,$K. Taken together, our theoretical predictions provide reasonable models with consistent mass-loss rate determination, which can be used for quantitative study of stars from Magellanic Clouds.

We calculate global (unified) wind models of main-sequence, giant, and supergiant O stars from our Galaxy. The models are calculated by solving hydrodynamic, kinetic equilibrium (also known as NLTE) and comoving-frame (CMF) radiative transfer equations from the (nearly) hydrostatic photosphere to the supersonic wind. For given stellar parameters, our models predict the photosphere and wind structure and in particular the wind mass-loss rates without any free parameters. Our predicted mass-loss rates are by a factor of 2--5 lower than the commonly used predictions. A possible cause of the difference is abandoning of the Sobolev approximation for the calculation of the radiative force, because our models agree with predictions of CMF NLTE radiative transfer codes. Our predicted mass-loss rates agree nicely with the mass-loss rates derived from observed near-infrared and X-ray line profiles and are slightly lower than mass-loss rates derived from combined UV and H$\alpha$ diagnostics. The empirical mass-loss rate estimates corrected for clumping may therefore be reconciled with theoretical predictions in such a way that the average ratio between individual mass-loss rate estimates is not higher than about $ 1.6 $. On the other hand, our predictions are by factor of $ 4.7 $ lower than pure H$\alpha$ mass-loss rate estimates and can be reconciled with these values only assuming a microclumping factor of at least eight.

The magnetic chemically peculiar (mCP) stars of the upper main sequence exhibit periodic light, magnetic, radio, and spectroscopic variations that can be adequately explained by a model of a rigidly rotating magnetized star with persistent surface structures. The majority of mCP stars rotate at strictly constant periods. However, there are a few mCP stars whose rotation periods vary on timescales of decades while the shape of their phase curves remains unchanged. In the case of CU Vir and V901 Ori, we have detected cyclic period variations. We demonstrate that the period oscillations of CU Vir may be a consequence of the interaction of the internal magnetic field and differential rotation.

Magnetic B-type stars exhibit photometric variability due to diverse causes, and consequently on a variety of timescales. In this paper we describe interpretation of BRITE photometry and related ground-based observations of 4 magnetic B-type systems: $\epsilon$ Lupi, $\tau$ Sco, a Cen and $\epsilon$ CMa.

Some chemically peculiar stars in the upper main sequence show rotational period variations of unknown origin. We propose these variations are a consequence of the propagation of internal waves in magnetic rotating stars that lead to the torsional oscillations of the star. We simulate the magnetohydrodynamic waves and calculate resonant frequencies for two stars that show rotational variations: CU Vir and HD 37776. We provide updated analyses of rotational period variations in these stars and compare our results with numerical models. For CU Vir, the length of the observed rotational-period cycle, $\mathit\Pi=67.6(5)$ yr, can be well reproduced by the models, which predict a cycle length of 51 yr. However, for HD 37776, the observed lower limit of the cycle length, $\mathit\Pi\geq100$ yr, is significantly longer than the numerical models predict. We conclude that torsional oscillations provide a reasonable explanation at least for the observed period variations in CU Vir.

A small fraction of hot OBA stars host global magnetic fields with field strengths of the order of 0.1-10 kG. This leads to the creation of persistent surface structures (spots) in stars with sufficiently weak winds as a result of the radiative diffusion. These spots become evident in spectroscopic and photometric variability. This type of variability is not expected in stars with strong winds, where the wind inhibits the radiative diffusion. Therefore, a weak photometric variability of the magnetic Of?p star HD 191612 is attributed to the light absorption in the circumstellar clouds. We study the nature of the photometric variability of HD 191612. We assume that the variability results from variable wind blanketing induced by surface variations of the magnetic field tilt and modulated by stellar rotation. We used our global kinetic equilibrium (NLTE) wind models with radiative force determined from the radiative transfer equation in the comoving frame (CMF) to predict the stellar emergent flux. Our models describe the stellar atmosphere in a unified manner and account for the influence of the wind on the atmosphere. The models are calculated for different wind mass-loss rates to mimic the effect of magnetic field tilt on the emergent fluxes. We integrate the emergent fluxes over the visible stellar surface for individual rotational phases, and calculate the rotationally modulated light curve of HD 191612. The wind blanketing that varies across surface is able to explain a part of the observed light variability in this star. The mechanism is able to operate even at relatively low mass-loss rates. The remaining variability is most likely caused by the flux absorption in circumstellar clouds. The variable wind blanketing is an additional source of the light variability in massive stars. The presence of the rotational light variability may serve as a proxy for the magnetic field.

Mass-loss rate is one of the most important stellar parameters. We aim to provide mass-loss rates as a function of subdwarf parameters and to apply the formula for individual subdwarfs, to predict the wind terminal velocities, to estimate the influence of the magnetic field and X-ray ionization on the stellar wind, and to study the interaction of subdwarf wind with mass loss from Be and cool companions. We used our kinetic equilibrium (NLTE) wind models with the radiative force determined from the radiative transfer equation in the comoving frame (CMF) to predict the wind structure of subluminous hot stars. Our models solve stationary hydrodynamical equations, that is the equation of continuity, equation of motion, and energy equation and predict basic wind parameters. We predicted the wind mass-loss rate as a function of stellar parameters, namely the stellar luminosity, effective temperature, and metallicity. The derived wind parameters (mass-loss rates and terminal velocities) agree with the values derived from the observations. The radiative force is not able to accelerate the homogeneous wind for stars with low effective temperatures and high surface gravities. We discussed the properties of winds of individual subdwarfs. The X-ray irradiation may inhibit the flow in binaries with compact components. In binaries with Be components, the winds interact with the disk of the Be star. Stellar winds exist in subluminous stars with low gravities or high effective temperatures. Despite their low mass-loss rates, they are detectable in the ultraviolet spectrum and cause X-ray emission. Subdwarf stars may lose a significant part of their mass during the evolution. The angular momentum loss in magnetic subdwarfs with wind may explain their low rotational velocities. Stellar winds are especially important in binaries, where they may be accreted on a compact or cool companion. (abridged)

Hot stars emit large amounts of X-rays, which are assumed to originate in the supersonic stellar wind. Part of the emitted X-rays is subsequently absorbed in the wind and influences its ionization state. Because hot star winds are driven radiatively, the modified ionization equilibrium affects the radiative force. We review the recent progress in modelling the influence of X-rays on the radiative equilibrium and on the radiative force. We focus particularly on single stars with X-rays produced in wind shocks and on binaries with massive components, which belong to the most luminous objects in X-rays.

The initial success of the Rigidly Rotating Magnetosphere (RRM) model application to the B2Vp star sigma OriE by Townsend, Owocki & Groote (2005) triggered a renewed era of observational monitoring of this archetypal object. We utilize high-resolution spectropolarimetry and the magnetic Doppler imaging (MDI) technique to simultaneously determine the magnetic configuration, which is predominately dipolar, with a polar strength Bd = 7.3-7.8 kG and a smaller non-axisymmetric quadrupolar contribution, as well as the surface distribution of abundance of He, Fe, C, and Si. We describe a revised RRM model that now accepts an arbitrary surface magnetic field configuration, with the field topology from the MDI models used as input. The resulting synthetic Ha emission and broadband photometric observations generally agree with observations, however, several features are poorly fit. To explore the possibility of a photospheric contribution to the observed photometric variability, the MDI abundance maps were used to compute a synthetic photospheric light curve to determine the effect of the surface inhomogeneities. Including the computed photospheric brightness modulation fails to improve the agreement between the observed and computed photometry. We conclude that the discrepancies cannot be explained as an effect of inhomogeneous surface abundance. Analysis of the UV light variability shows good agreement between observed variability and computed light curves, supporting the accuracy of the photospheric light variation calculation. We thus conclude that significant additional physics is necessary for the RRM model to acceptably reproduce observations of not only sigma Ori E, but also other similar stars with significant stellar wind-magnetic field interactions.

Binaries with hot massive components are strong X-ray sources. Besides the intrinsic X-ray emission of individual binary members originating in their winds, X-ray emission stems from the accretion on the compact companion or from wind collision. Since hot star winds are driven by the light absorption in the lines of heavier elements, wind acceleration is sensitive to the ionization state. Therefore, the over-ionization induced by external X-ray source strongly influences the winds of individual components. We studied the effect of external X-ray irradiation on hot star winds. We used our kinetic equilibrium (NLTE) wind models to estimate the influence of external X-ray ionization for different X-ray luminosities and source distances. The models are calculated for parameters typical of O stars. The influence of X-rays is given by the X-ray luminosity, by the optical depth between a given point and the X-ray source, and by a distance to the X-ray source. Therefore, the results can be interpreted in the diagrams of X-ray luminosity vs.~the optical depth parameter. X-rays are negligible in binaries with low X-ray luminosities or at large distances from the X-ray source. The influence of X-rays is stronger for higher X-ray luminosities and in closer proximity of the X-ray source. There is a forbidden area with high X-ray luminosities and low optical depth parameters, where the X-ray ionization leads to wind inhibition. There is excellent agreement between the positions of observed stars in these diagrams and our predictions. All wind-powered high-mass X-ray binary (HMXB) primaries lie outside the forbidden area. Many of them lie close to the border of the forbidden area, indicating that their X-ray luminosities are self-regulated. We discuss the implications of our work for other binary types.

BS Cir is a representative of moderately cool magnetic chemically peculiar stars which displays very strong light variations in Stroemgren index c1 indicating large changes in the height of the Balmer jump. We present two-spot model of light variations fitting successfully all of nine light curves obtained in the spectral region 335-750 nm. We also discuss the nature of the observed variations of intensities of Fe, Cr, Ti, Si, Mg and RE spectral lines and possible mechanisms matching the observed light variations. It was confirmed that the observed period of BS Cir 2.204 d is rising with the rate of dP/dt=5.4(4)x10^-9. The found minor secular changes in the shape of light curve should be compatible with the period changes caused by precessional motion due to magnetic distortion of the star.

Chemically peculiar stars of the upper part of the main sequence show periodical variability in line intensities and continua, modulated by the stellar rotation, which is attributed to the existence of chemical spots on the surface of these stars. The flux variability is caused by the changing redistribution rate of the radiative flux predominantly from the short-wavelength part of the spectra to the long-wavelength part, which is a result of abundance anomalies. We study the nature of the multi-spectral variability of one of the brightest chemically peculiar stars, $\theta$ Aur. We predict the flux variability of $\theta$ Aur from the emerging intensities calculated for individual surface elements of the star taking into account horizontal variation of chemical composition derived from Doppler abundance maps. The simulated optical variability in the Str\"omgren photometric system and the ultraviolet flux variability agree well with observations. The IUE flux distribution is reproduced in great detail by our models. The resonance lines of magnesium and possibly also some lines of silicon are relatively weak in the ultraviolet domain, which indicates non-negligible vertical abundance gradients in the atmosphere. We also derive a new period of the star, $P=3.618\,664(10)$ d, from all available photometric and magnetic measurements and show that the observed rotational period is constant over decades. The ultraviolet and visual variability of $\theta$ Aur is mostly caused by silicon bound-free absorption and chromium and iron line absorption. These elements redistribute the flux mainly from the far-ultraviolet region to the near-ultraviolet and optical regions in the surface abundance spots. The light variability is modulated by the stellar rotation. The ultraviolet domain is key for understanding the properties of chemically peculiar stars. (abridged)

Evolutionary models of fast-rotating stars show that the stellar rotational velocity may approach the critical speed. Critically rotating stars cannot spin up more, therefore they lose their excess angular momentum through an equatorial outflowing disk. The radial extension of such disks is unknown, partly because we lack information about the radial variations of the viscosity. We study the magnetorotational instability, which is considered to be the origin of anomalous viscosity in outflowing disks. We used analytic calculations to study the stability of outflowing disks submerged in the magnetic field. The magnetorotational instability develops close to the star if the plasma parameter is large enough. At large radii the instability disappears in the region where the disk orbital velocity is roughly equal to the sound speed. The magnetorotational instability is a plausible source of anomalous viscosity in outflowing disks. This is also true in the region where the disk radial velocity approaches the sound speed. The disk sonic radius can therefore be roughly considered as an effective outer disk radius, although disk material may escape from the star to the insterstellar medium. The radial profile of the angular momentum-loss rate already flattens there, consequently, the disk mass-loss rate can be calculated with the sonic radius as the effective disk outer radius. We discuss a possible observation determination of the outer disk radius by using Be and Be/X-ray binaries.

Hot star wind mass-loss rates depend on the abundance of individual elements. This dependence is usually accounted for assuming scaled solar chemical composition. However, this approach may not be justified in evolved rotating stars. The rotational mixing brings CNO-processed material to the stellar surface, increasing the abundance of nitrogen at the expense of carbon and oxygen, which potentially influences the mass-loss rates. We study the influence of the modified chemical composition resulting from the rotational mixing on the wind parameters, particularly the wind mass-loss rates. We use our NLTE wind code to predict the wind structure and compare the calculated wind mass-loss rate for the case of scaled solar chemical composition and the composition affected by the CNO cycle. We show that for a higher mass-fraction of heavier elements $Z/Z_\odot\gtrsim0.1$ the change of chemical composition from the scaled solar to the CNO-processed scaled solar composition does not significantly affect the wind mass-loss rates. The missing line force caused by carbon and oxygen is compensated for by nitrogen line force. However, for a very low-mass fraction of heavier elements $Z/Z_\odot\lesssim0.1$ the rotational mixing significantly affects the wind mass-loss rates. Moreover, the decrease of the mass-loss rate with metallicity is stronger at such low metallicities. We study the relevance of the wind momentum-luminosity relationship for different metallicities and show that for a metallicity $Z/Z_\odot\lesssim0.1$ the relationship displays a large scatter, which depreciates the use of this relationship at the lowest metallicities.

We calculate radiatively driven wind models of main-sequence B stars and provide the wind mass-loss rates and terminal velocities. The main-sequence mass-loss rate strongly depends on the stellar effective temperature. For the hottest B stars the mass-loss rate amounts to $10^{-9}\,\text{M}_\odot\,\text{year}^{-1}$, while for the cooler ones themass-loss rate is lower by more than three orders of magnitude. Main-sequence B stars with solar abundance and effective temperatures lower than about$15\,000\,\text{K}$ (later than spectral type B5) do not have any homogeneous line-driven wind. We predict the wind mass-loss rates for the solar chemical composition and for the modified abundance of heavier elements to study the winds of chemically peculiar stars. The mass-loss rate may either increase or decrease with increasing abundance, depending on the importance of the induced emergent flux redistribution. Stars with overabundant silicon may have homogeneous winds even below the solar abundance wind limit at $15\,000\,\text{K}$. The winds of main-sequence B stars lie below the static limit, that is, a static atmosphere solution is also possible. This points to an important problem regarding the initiation of these winds. We discuss the implications of our models for rotational braking, filling the magnetosphere of Bp stars, and for chemically peculiar stars.

The galactic magnetic chemically peculiar (mCP) stars of the upper main sequence are well known as periodic spectral and light variables. The observed variability is obviously caused by the uneven distribution of overabundant chemical elements on the surfaces of rigidly rotating stars. The mechanism causing the clustering of some chemical elements into disparate structures on mCP stars has not been fully understood up to now. The observations of light changes of mCP candidates recently revealed in the nearby Large Magellanic Cloud (LMC) should provide us with information about their rotational periods and about the distribution of optically active elements on mCP stars born in other galaxies. We queried for photometry at the Optical Gravitational Lensing Experiment (OGLE)-III survey of published mCP candidates selected because of the presence of the characteristic 5200A flux depression. In total, the intersection of both sources resulted in twelve stars. We performed our own and standard periodogram time series analyses of all available data. The final results are, amongst others, the frequency of the maximum peak and the bootstrap probability of its reality. We detected that only two mCP candidates, 190.1 1581 and 190.1 15527, may show some weak rotationally modulated light variations with periods of 1.23 and 0.49 days; however, the 49% and 32% probabilities of their reality are not very satisfying. The variability of the other 10 mCP candidates is too low to be detectable by their V and I OGLE photometry. The relatively low amplitude variability of the studied LMC mCP candidates sample can be explained by the absence of photometric spots of overabundant optically active chemical elements. The unexpected LMC mCPs behaviour is probably caused by different conditions during the star formation in the LMC and the Galaxy.

The light variability of hot magnetic chemically peculiar stars is typically caused by the flux redistribution in spots with peculiar abundance. This raises the question why some stars with surface abundance spots show significant rotational light variability, while others do not. We study the Bp star HD 64740 to investigate how its remarkable inhomogeneities in the surface distribution of helium and silicon, and the corresponding strong variability of many spectral lines, can result in one of the faintest photometric variabilities among the Bp stars. We used model atmospheres and synthetic spectra calculated for the silicon and helium abundances from surface abundance maps to predict the ultraviolet and visual light and line variability of HD 64740. The predicted fluxes and line profiles were compared with the observed ones derived with the IUE, HST, and Hipparcos satellites and with spectra acquired using the FEROS spectrograph at the 2.2m MPG/ESO telescope. We are able to reproduce the observed visual light curve of HD 64740 assuming an inhomogeneous distribution of iron correlated with silicon distribution. The light variations in the ultraviolet are hardly detectable. We detect the variability of many ultraviolet lines of carbon, silicon, and aluminium and discuss the origin of these lines and the nature of their variations. The maximum abundances of helium and silicon are not high enough to cause significant light variations. The detected variability of many ultraviolet lines is most likely of atmospheric origin and reflects the inhomogeneous elemental surface distribution. The variability of the Civ resonance lines is stronger and it probably results from the dependence of the wind mass-loss rate on the chemical composition and magnetic field orientation. We have not been able to detect a clear signature of the matter trapped in the circumstellar clouds.

Even though it is broadly accepted that single Be stars are rapidly rotating stars surrounded by a flat rotating circumstellar disk, there is still a debate about how fast these stars rotate and also about the mechanisms involved in the angular-momentum and mass input in the disk. We study the properties of stars that rotate near their critical-rotation rate and investigate the properties of the disks formed by equatorial mass ejections. We used the most recent Geneva stellar evolutionary tracks for rapidly rotating stars that reach the critical limit and used a simple model for the disk structure. We obtain that for a 9 Msun star at solar metallicity, the minimum average velocity during the Main Sequence phase to reach the critical velocity is around 330 km/s, whereas it would be 390 km/s at the metallicity of the Small Magellanic Cloud (SMC). Red giants or supergiants originating from very rapid rotators rotate six times faster and show N/C ratios three times higher than those originating from slowly rotating stars. This difference becomes stronger at lower metallicity. It might therefore be very interesting to study the red giants in clusters that show a large number of Be stars on the MS band. On the basis of our single-star models, we show that the observed Be-star fraction with cluster age is compatible with the existence of a temperature-dependent lower limit in the velocity rate required for a star to become a Be star. The mass, extension, and diffusion time of the disks produced when the star is losing mass at the critical velocity, obtained from simple parametrized expressions, are not too far from those estimated for disks around Be-type stars. At a given metallicity, the mass and the extension of the disk increase with the initial mass and with age on the MS phase. Denser disks are expected in low-metallicity regions.

Different diagnostics of hot star wind mass-loss rates provide results that are difficult to reconcile with each other. The widely accepted presence of clumping in hot star winds implies a significant reduction of observational mass-loss rate estimates from diagnostics that depend on the square of the density. Moreover, the ultraviolet Pv resonance lines indicate a possible need for even stronger reduction of hot star mass-loss rates, provided that Pv is a dominant ionization stage of phosphorus at least in some hot stars. The latter assumption is challenged by a possible presence of the XUV radiation. Here we study the influence of the XUV radiation on the Pv ionization fraction in the hot star winds. By a detailed solution of the hydrodynamical, radiative transfer, and statistical equilibrium equations we confirm that sufficiently strong XUV radiation source may decrease the Pv ionization fraction, possibly depreciating the Pv lines as a reliable mass-loss rate indicator. On the other hand, the XUV radiation influences also the ionization fraction of heavier ions that drive the wind, leading to a decrease of the wind terminal velocity. Consequently, we conclude that the XUV radiation alone can not bring theory and observations in accord. We fit our predicted wind mass-loss rates by a suitable formula and compare the results with the observational mass-loss rate diagnostics. We show that for supergiants and giants the theoretical predictions do not contradict the mass-loss rate estimates based on X-ray line profiles or density squared diagnostics. On the other hand, for main-sequence stars the predicted mass-loss rates are still significantly higher than that inferred from Pv or X-ray lines. This indicates that the "weak wind problem" recently detected in low-luminosity main-sequence stars may occur to some extent also for the stars with higher luminosity.

Vela X-1 is the archetype of high-mass X-ray binaries, composed of a neutron star and a massive B supergiant. The supergiant is a source of a strong radiatively-driven stellar wind. The neutron star sweeps up this wind, and creates a huge amount of X-rays as a result of energy release during the process of wind accretion. Here we provide detailed NLTE models of the Vela X-1 envelope. We study how the X-rays photoionize the wind and destroy the ions responsible for the wind acceleration. The resulting decrease of the radiative force explains the observed reduction of the wind terminal velocity in a direction to the neutron star. The X-rays create a distinct photoionized region around the neutron star filled with a stagnating flow. The existence of such photoionized bubbles is a general property of high-mass X-ray binaries. We unveiled a new principle governing these complex objects, according to which there is an upper limit to the X-ray luminosity the compact star can have without suspending the wind due to inefficient line driving

The spectral energy distribution (SED) in chemically peculiar stars may be significantly affected by their abundance anomalies. The observed SED variations are usually assumed to be a result of inhomogeneous surface distribution of chemical elements, flux redistribution and stellar rotation. However, the direct evidence for this is still only scarce. We aim to identify the processes that determine the SED and its variability in the UV and visual spectral domains of the helium-weak star CU Vir. We used the model atmospheres to obtain the emergent flux and predict the rotationally modulated flux variability of the star. We show that most of the light variations in the vby filters of the Stromgren photometric system are a result of the uneven surface distribution of silicon, chromium, and iron. Our models are only able to explain a part of the variability in the u filter, however. The observed UV flux distribution is very well reproduced, and the models are able to explain most of the observed features in the UV light curve. The variability observed in the visible is merely a faint gleam of that in the UV. While the amplitude of the light curves reaches only several hundredths of magnitude in the visual domain, it reaches about 1 mag in the UV. The visual and UV light variability of CU Vir is caused by the flux redistribution from the far UV to near UV and visible regions, inhomogeneous distribution of the elements and stellar rotation. Bound-free transitions of silicon and bound-bound transitions of iron and chromium contribute the most to the flux redistribution. This mechanism can explain most of the rotationally modulated light variations in the filters centred on the Paschen continuum and on the UV continuum of the star CU Vir. However, another mechanism(s) has to be invoked to fully explain the observed light variations in the u filter and in the region 2000-2500 A.

CU Vir and V901 Ori belong among these few magnetic chemically peculiar stars whose rotation periods vary on timescales of decades. We aim to study the stability of the periods in CU Vir and V901 Ori using all accessible observational data containing phase information. We collected all available relevant archived observations supplemented with our new measurements of these stars and analysed the period variations of the stars using a novel method that allows for the combination of data of diverse sorts. We found that the shapes of their phase curves were constant, while the periods were changing. Both stars exhibit alternating intervals of rotational braking and acceleration. The rotation period of CU Vir was gradually shortening until the year 1968, when it reached its local minimum of 0.52067198 d. The period then started increasing, reaching its local maximum of 0.5207163 d in the year 2005. Since that time the rotation has begun to accelerate again. We also found much smaller period changes in CU Vir on a timescale of several years. The rotation period of V901 Ori was increasing for the past quarter-century, reaching a maximum of 1.538771 d in the year 2003, when the rotation period began to decrease. A theoretically unexpected alternating variability of rotation periods in these stars would remove the spin-down time paradox and brings a new insight into structure and evolution of magnetic upper-main-sequence stars.

The velocity law of hot star winds is usually parameterized via the so-called beta velocity law. Although this parameterization stems from theoretical considerations, it is not the most accurate description of the wind velocity law that follows from hydrodynamical calculations. We show that the velocity profile of our hydrodynamical wind models is described much better by polynomial approximation. This approximation provides a better fit than the beta velocity law already for the same number of free parameters.

Photometrically and spectroscopically variable chemically peculiar (CP) stars are the optimum laboratories for testing the rotational evolution of main sequence stars. A vast majority of well-studied CP stars have quite steady rotational periods (e.g. SrCrEu star CQ UMa). However, there are several CP stars that exhibit apparent period variations. The origin of the period variations is unclear in many cases. We describe the observed period variations of several individual CP stars, especially V901 Ori, sigma Ori E, HR 7355, CU Vir, SX Ari, and EE Dra. The CP stars with unsteady periods now represent a very diverse group with dissimilar O-C diagrams and time scales. We also discuss the causes of the period changes found and a possible cyclicity or chaoticity of them.

Context. TW Draconis is one of the best studied Algol-type eclipsing binaries. There is significant evidence for miscellaneous physical processes between interacting binary components manifesting themselves by period and light curve changes. Aims. Obtaining new set of photometric and spectroscopic observations, we analysed them together with the older spectroscopic and photometric data to build model of this eclipsing system with respect to observed changes of O-C diagram and light curve. Methods. Reduction of new spectra was carried out in the IRAF and SPEFO programs. Radial velocities were determined manually using SPEFO, by CCF using the Zverko's code and from the program KOREL. Orbital elements were derived with the FOTEL program and via disentangling with KOREL. The final combined solution was obtained with the programs PHOEBE and FOTEL. Results. Photometry shows small irregularities in light curves as a results of pulsating of one component and spot activity. Using net of KOREL outputs we found the mass ratio q=0.405(3).We confirm the presence of stellar matter around the primary. Even after subtraction of ADS 9706B influence, light curve solutions show third light in the system. Models in FOTEL and Phoebe are presented. Conclusions. TW Dra is an astrophysically very interesting eclipsing binary. Future combination of interferometry, spectroscopy, and photometry is promising. It could definitely confirm the hypothesis of quadruple system for TW Dra and explain behaviour of this system in complex.

We examine the nature and role of mass loss via an equatorial decretion disk in massive stars with near-critical rotation induced by evolution of the stellar interior. In contrast to the usual stellar wind mass loss set by exterior driving from the stellar luminosity, such decretion-disk mass loss stems from the angular momentum loss needed to keep the star near and below critical rotation, given the interior evolution and decline in the star's moment of inertia. Because the specific angular momentum in a Keplerian disk increases with the square root of the radius, the decretion mass loss associated with a required level of angular momentum loss depends crucially on the outer radius for viscous coupling of the disk, and can be significantly less than the spherical mass loss the spherical, wind-like mass loss commonly assumed in evolutionary calculations. We discuss the physical processes that affect the outer disk radius, including thermal disk outflow, and ablation of the disk material via a line-driven wind induced by the star's radiation. We present parameterized scaling laws for taking account of decretion-disk mass loss in stellar evolution codes, including how these are affected by metallicity, or by presence within a close binary and/or a dense cluster. Effects similar to those discussed here should also be present in accretion disks during star formation, and may play an important role in shaping the distribution of rotation speeds on the ZAMS.

We simulate the light variability of the Ap star epsUMa using the observed surface distributions of Fe, Cr, Ca, Mn, Mg, Sr and Ti obtained with the help of Doppler Imaging technique. Using all photometric data available we specified light variations of epsUMa modulated by its rotation from far UV to IR. We employed the LLmodels stellar model atmosphere code to predict the light variability in different photometric systems. The rotational period of epsUMa is refined to 5d088631(18). It is shown that the observed light variability can be explained as a result of the redistribution of radiative flux from the UV spectral region to the visual caused by the inhomogeneous surface distribution of chemical elements. Among seven mapped elements, only Fe and Cr significantly contribute to the amplitude of the observed light variability. In general, we find a very good agreement between theory and observations. We confirm the important role of Fe and Cr to the magnitude of the well-known depression around 5200 \AA\ through the analysis of the peculiar $a$-parameter. Finally, we show that the abundance spots of considered elements cannot explain the observed variability in near UV and $\beta$ index which are likely due to some other causes. The inhomogeneous surface distribution of chemical elements can explain most of the observed light variability of the A-type CP star epsUMa.

We provide hot star wind models with radiative force calculated using the solution of comoving frame (CMF) radiative transfer equation. The wind models are calculated for the first stars, O stars, and the central stars of planetary nebulae. We show that without line overlaps and with solely thermal line broadening the pure Sobolev approximation provides a reliable estimate of the radiative force even close to the wind sonic point. Consequently, models with the Sobolev line force provide good approximations to solutions obtained with non-Sobolev transfer. Taking line overlaps into account, the radiative force becomes slightly lower, leading to a decrease in the wind mass-loss rate by roughly 40%. Below the sonic point, the CMF line force is significantly lower than the Sobolev one. In the case of pure thermal broadening, this does not influence the mass-loss rate, as the wind mass-loss rate is set in the supersonic part of the wind. However, when additional line broadening is present (e.g., the turbulent one) the region of low CMF line force may extend outwards to the regions where the mass-loss rate is set. This results in a decrease in the wind mass-loss rate. This effect can at least partly explain the low wind mass-loss rates derived from some observational analyses of luminous O stars.

During the evolution of rotating first stars, which initially consisted of only hydrogen and helium, CNO elements may emerge to their surface. These stars may therefore have winds that are driven only by CNO elements. We study weak wind effects (Gayley-Owocki heating and multicomponent effects) in stellar winds of first generation stars driven purely by CNO elements. We apply our NLTE multicomponent models and hydrodynamical simulations. The multicomponent effects (frictional heating and decoupling) are important particularly for low metallicity winds, but they influence mass loss rate only if they cause decoupling for velocities lower than the escape velocity. The multicomponent effects also modify the feedback from first stars. As a result of the decoupling of radiatively accelerated metals from hydrogen and helium, the first low-energy cosmic ray particles are generated. We study the interaction of these particles with the interstellar medium concluding that these particles easily penetrate the interstellar medium of a given minihalo. We discuss the charging of the first stars by means of their winds. Gayley-Owocki heating, frictional heating, and the decoupling of wind components occur in the winds of evolved low-metallicity stars and the solar metallicity main-sequence stars.

We calculate NLTE models of stellar winds of hot compact stars (central stars of planetary nebulae and subdwarf stars). The studied range of subdwarf parameters is selected to cover a large part of these stars. The models predict the wind hydrodynamical structure and provide mass-loss rates for different abundances. Our models show that CNO elements are important drivers of subdwarf winds, especially for low-luminosity stars. We study the effect of X-rays and instabilities on these winds. Due to the line-driven wind instability, a significant part of the wind could be very hot.

Strong meridional mixing induced by rapid rotation is one reason why all hot main-sequence stars are not chemically peculiar. However, the finding that the He-strong CP star HR 7355 is a rapid rotator complicates this concept. Our goal is to explain the observed behaviour of HR 7355 based on period analysis of all available photometry. Over two years, we acquired 114 new BV observations of HR 7355 at observatories in Arizona, U.S.A and Cape Town, South Africa. We performed period analyses of the new observations along with new analyses of 732 archival measurements from the Hipparcos and ASAS projects. We find that the light curves of HR 7355 in various filters are quite similar, with amplitudes 0.035(4), 0.036(4), and 0.038(3) mag in B, Hp and V, respectively. The light curves are double-peaked, with unevenly deep minima. We substantially refine the rotational period to be P=0.5214410(4) d, indicating that HR 7355 is the most rapidly rotating CP star known. Our period analyses reveal a possible lengthening of the rotational period with dP/dt/P = 2.4(8)x 10^-6/yr. We conclude that the shape and amplitude of HR 7355 light curves are typical of magnetic He-strong CP stars, for which light variations are the result of photospheric spots on the surface of a rotating star. We hypothesise that the light variations are caused mainly by an uneven distribution of overabundant helium on the star's surface. We briefly describe and discuss the cause of the rapid rotational braking of the star.

Hot stars are sources of X-ray emission originating in their winds. Although hydrodynamical simulations that are able to predict this X-ray emission are available, the inclusion of X-rays in stationary wind models is usually based on simplifying approximations. To improve this, we use results from time-dependent hydrodynamical simulations of the line-driven wind instability (seeded by the base perturbation) to derive the analytical approximation of X-ray emission in the stellar wind. We use this approximation in our non-LTE wind models and find that an improved inclusion of X-rays leads to a better agreement between model ionization fractions and those derived from servations. Furthermore, the slope of the L_x-L relation is in better agreement with observations, however the X-ray luminosity is underestimated by a factor of three. We propose a possible solution for this discrepancy.

Magnetic chemically peculiar (mCP) stars of the upper part of the main sequence usually exhibit strictly periodic variations of their light, spectrum and magnetic field. These changes can be well explained using the oblique rotator model of the rigidly rotating star with persistent photometric spots, uneven horizontal abundance of chemical elements and dipole-like magnetic field. Long-term observations of mCP stars variability enable to determine their rotational periods with an extraordinary accuracy and reliability. We compare rotational periods of mCP stars with those of normal main sequence stars of the same spectral type and discuss the found discrepancies. We also discuss the cases of mCP stars V901 Ori, SX Ari, CU Vir, (and possibly also EE Dra), which have been reported to display an increase of their rotational periods.

Although photometric variations of chemically peculiar (CP) stars are frequently used to determine their rotational periods, the detailed mechanism of their light variability remains poorly understood. We simulate the light variability of the star HR 7224 using the observed surface distribution of silicon and iron. We used the TLUSTY model atmospheres calculated for the appropriate silicon and iron abundances to obtain the emergent flux and to predict the rotationally modulated light curve of the star. We also obtained additional photometric measurements and employed our own regression procedure to derive a more precise estimate of the light elements. We show that the light variation of the star can be explained as a result of i) the uneven surface distribution of the elements, ii) the flux redistribution from the ultraviolet to the visible part of the spectrum, and iii) rotation of the star. We show that the silicon bound-free transitions and iron bound-bound transitions provide the main contribution to the flux redistribution, although an additional source of opacity is needed. We confirm that numerous iron lines significantly contribute to the well-known depression at 5200 A and discuss the connection between iron abundance and the value of peculiarity index a. The uneven surface distribution of silicon and iron is able to explain most of the rotationally modulated light variation in the star HR 7224.

We study the influence of X-rays on the wind structure of selected O stars. For this purpose we use our non-local thermodynamic equilibrium (NLTE) wind code with inclusion of additional artificial source of X-rays, assumed to originate in the wind shocks. We show that the influence of shock X-ray emission on wind mass-loss rate is relatively small. Wind terminal velocity may be slightly influenced by the presence of strong X-ray sources, especially for stars cooler than Teff < 35 000 K. We discuss the origin of the Lx/L \sim 10^-7 relation. For stars with thick wind this relation can be explained assuming that the cooling time depends on wind density. Stars with optically thin winds exhibiting the "weak wind problem" display enhanced X-ray emission which may be connected with large shock cooling length. We propose that this effect can explain the "weak wind problem". Inclusion of X-rays leads to a better agreement of the model ionization structure with observations. However, we do not found any significant influence of X-rays on Pv ionization fraction implying that the presence of X-rays cannot explain the Pv problem. We study the implications of modified ionization equilibrium due to shock emission on the line transfer in the X-ray region. We conclude that the X-ray line profiles of helium-like ions may be affected by the line absorption within the cool wind.

During the evolution of first stars, the CNO elements may emerge on their surfaces due to the mixing processes. Consequently, these stars may have winds driven purely by CNO elements. We study the properties of such stellar winds and discuss their influence on the surrounding environment. For this purpose, we used our own NLTE models and tested which stellar parameters of the first stars at different evolutionary stages result in CNO winds. If such winds are possible, we calculate their hydrodynamic structure and predict their parameters. We show that, while the studied stars do not have any wind driven purely by hydrogen and helium, CNO driven winds exist in more luminous stars. On the other hand, for very hot stars, CNO elements are too ionized to drive a wind. In most cases the derived mass-loss rate is much less than calculated with solar mixture of elements. This is because wind mass-loss rate in present hot stars is dominated by elements heavier than CNO. We conclude that, until a sufficient amount of these elements is created, the influence of line-driven winds is relatively small on the evolution of hot stars (which are not close to the Eddington limit).

Thermal wind emission in the form of free-free and free-bound emission is known to show up in the infrared and radio continuum of hot and massive stars. For OB supergiants with moderate mass loss rates and a wind velocity distribution with \beta = 0.8...1.0, no influence of the wind to the optical continuum, i.e. for \lambda < 1 micron, is expected. Investigations of stellar and wind parameters of OB supergiants over the last few years suggest, however, that for many objects \beta is much higher than 1.0, reaching values up to 3.5. We investigate the influence of the free-free and free-bound emission on the emerging radiation, especially at optical wavelengths, from OB supergiants having wind velocity distributions with \beta > 1. For the case of a spherically symmetric, isothermal wind in local thermodynamical equilibrium (LTE) we calculate the free-free and free-bound processes and the emerging wind and total continuum spectra. We localize the generation region of the optical wind continuum and especially focus on the influence of a \beta-type wind velocity distribution with \beta > 1 on the formation of the wind continuum at optical wavelengths. The optical wind continuum is found to be generated within about 2 R_* which is exactly the wind region where \beta strongly influences the density distribution. We find that for \beta > 1, the continuum of a typical OB supergiant can indeed be contaminated with thermal wind emission, even at optical wavelengths. The strong increase in the optical wind emission is dominantly produced by free-bound processes.

We analyzed uvby and Hp light curves of 19 well observed magnetic CP stars selected from the "On-line database of photometric observations of mCP stars" of which light curves in all the five colours were similar. We assumed that among these photometrically simply behaving (PSB) stars could be found such ones which have a single photometric spot. The insight into such simple situations would help us to comprehend more complicated cases. Light curves of the 19 PSB mCP stars proved to be generally nearly symmetric but surprisingly diverse. The analysis shows that only in the case of HD 110956B, HD 188041, and perhaps HD 193722 we are able to explain their photometric behaviour by a simple one-spot model. Consequently, occurrence of more than one photometric spot on an mCP star is typical.

We study the long-term light and spectral variations in the He-strong magnetic chemically peculiar star HD 37776 (V901 Ori) to search for changes of its 1.5387 d period in 1976-2007. We analyze all published photometric observations and spectrophotometry in the HeI 4026 A line. The data were supplmented with 506 new (U)VB observations obtained during the last 2 observing seasons, 66 estimates of HeI equivalent widths on 23 CFHT spectrograms and 35 of the 6-m Zeeman spectrograms. All the 1895 particular observations heve been processed simultaneously. We confirm the previously suspected increase of the period in HD 37776 which is a record-breaking among CP stars. The mean rate of the period increase during the last 31 years is 0.541+-0.020 s per year. We interpret this ongoing period increase as the slowing down of the star's surface rotation due to momentum loss through events and processes in its magnetosphere.

We study the influence of clumping on the predicted wind structure of O-type stars. For this purpose we artificially include clumping into our stationary wind models. When the clumps are assumed to be optically thin, the radiative line force increases compared to corresponding unclumped models, with a similar effect on either the mass-loss rate or the terminal velocity (depending on the onset of clumping). Optically thick clumps, alternatively, might be able to decrease the radiative force.

The influence of the wind to the total continuum of OB supergiants is discussed. For wind velocity distributions with \beta > 1.0, the wind can have strong influence to the total continuum emission, even at optical wavelengths. Comparing the continuum emission of clumped and unclumped winds, especially for stars with high \beta values, delivers flux differences of up to 30% with maximum in the near-IR. Continuum observations at these wavelengths are therefore an ideal tool to discriminate between clumped and unclumped winds of OB supergiants.

We simulate light curves of the helium strong chemically peculiar star HD 37776 assuming that the observed periodic light variations originate as a result of inhomogeneous horizontal distribution of chemical elements on the surface of a rotating star. We show that chemical peculiarity influences the monochromatic radiative flux, mainly due to bound-free processes. Using the model of the distribution of silicon and helium on HD 37776 surface, derived from spectroscopy, we calculate a photometric map of the surface and consequently the uvby light curves of this star. Basically, the predicted light curves agree in shape and amplitude with the observed ones. We conclude that the basic properties of variability of this helium strong chemically peculiar star can be understood in terms of the model of spots with peculiar chemical composition.

We present preliminary results of the generalized Principal Component Analysis (PCA) of light curves of 82 magnetic chemically peculiar (further mCP) stars applied to 54 thousand individual photometric observations in the uvby and Hp colours taken from the "On-line database of photometric observations of mCP stars". We show that every of the observed light curves (LCs) can be, with a sufficient accuracy, represented using five parameters of a harmonic polynomial of the second order, and that the third order and higher harmonics reflect only a noise. We found that a prevailing majority of the uvbyHp LCs can be satisfactorily well represented by a linear combination of a constant term and one or two (extraordinarily three) basic, mutually orthonormal functions, which reduces the number of parameters necessarily needed to describe the set of the LCs of a particular star. While the shape of an individual LC depends on the individual distribution of photometric spots on the surface of a rotating star, its amplitude depends on the mechanism of variability. To describe the amplitude of LCs we introduce a robust quantity effective amplitude. Applying the PCA to the set of the effective amplitudes of the LCs in all five colours we revealed at least three different sources of the light variability showing different wavelength dependence. The amplitudes of all the LCs harbor a component monotonously decreasing with increasing wavelengths. The second component reaches a remarkable extremum in the $v$ colour, and is in antiphase relative to the remaining four colours. It occurs particularly at cooler mCP-s having a high degree of the chemical peculiarity. The third component reflects a remarkable diversity of the light curves in the $u$ colour found in a few mCP stars.

We calculate the uvby light curves of the helium-strong chemically peculiar (CP) star HD 37776. Assuming that the chemical peculiarity influences the monochromatic radiative flux due to mainly the bound-free processes and using the model of the surface distribution of helium and silicon on HD 37776, as derived from spectroscopy, we calculate the emergent fluxes for a series of its rotational phases. We show that it is possible to consistently simulate light curves of a chemically peculiar star and that the basic properties of variability of helium CP stars can be understood in terms of the model of spots with a peculiar chemical composition.

We compare the hot star wind models calculated assuming older solar abundance determination with models calculated using the recently published values derived from 3D hydrodynamical model atmospheres. We show that the use of new abundances with lower metallicity improves the agreement between wind observation and theory in several aspects: (1) The predicted wind mass-loss rates are lower by a factor of 0.76. This leads to a better agreement with mass-loss rate determinations derived from observations with account of clumping. (2) As a result of the lowering of mass-loss rates, there is a better agreement between predicted modified wind momentum-luminosity relationship and that derived from observations with account of clumping. (3) Both the lower mass fraction of heavier elements and lower mass-loss rates lead to a decrease of the opacity in the X-ray region. This has influence on the prediction of the X-ray line profile shapes. (4) There is a better agreement between predicted PV ionization fractions and those derived from observations.

We study the influence of the frictional heating on the wind line-profiles of SMC stars. For this purpose we use our NLTE wind code to obtain consistent occupation numbers of studied levels and our radiative transfer code to solve the radiative transfer equation in moving media. We compare predicted wind line profiles calculated with and without frictional heating for a low-luminosity SMC star with a weak wind and discuss the relevance of frictional heating as a solution of a "weak-wind problem" for this star.

We study the influence of X-ray radiation on the wind parameters of O stars. For this purpose we use our own NLTE wind code. The X-ray emission (assumed to be generated in wind shocks) is treated as an input quantity. We study its influence on the mass-loss rate, terminal velocity and ionization state of the stellar wind of Galactic O stars.

We study the possibility of the helium decoupling in the stellar wind of sigma Ori E. To obtain reliable wind parameters for this star we first calculate an NLTE wind model and derive wind mass-loss rate and terminal velocity. Using corresponding force multipliers we study the possibility of helium decoupling. We find that helium decoupling is not possible for realistic values of helium charge (calculated from NLTE wind models). Helium decoupling seems only possible for a very low helium charge. The reason for this behaviour is the strong coupling between helium and hydrogen. We also find that frictional heating becomes important in the outer parts of the wind of sigma Ori E due to the collisions between some heavier elements and the passive components -- hydrogen and helium. For a metallicity ten times lower than the solar one both hydrogen and helium decouple from the metals and may fall back onto the stellar surface. However, this does not explain the observed chemical peculiarity since both these components decouple together from the absorbing ions. Although we do not include the effects of the magnetic field into our models, we argue that the presence of a magnetic field will likely not significantly modify the derived results because in such case model equations describe the motion parallel to the magnetic field.

We calculate NLTE line-driven wind models of selected O stars in the spectral range of O4 to O9 in the Small Magellanic Cloud (SMC). We compare predicted basic wind properties, i.e. the terminal velocity and the mass-loss rate with values derived from observation. We found relatively good agreement between theoretical and observed terminal velocities. On the other hand, predicted mass-loss rates and mass-loss rates derived from observation are in a good agreement only for higher mass-loss rates. Theoretical mass-loss rates lower than approximately 10^{-7} M_o/year are significantly higher than those derived from observation. These results confirm the previously reported problem of weak winds, since our calculated mass-loss rates are in a fair agreement with predictions of Vink et al. (2001). We study multicomponent models for these winds. For this purpose we develop a more detailed description of wind decoupling. We show that the instability connected with the decoupling of individual wind elements may occur for low-density winds. In the case of winds with very low observed mass-loss rate the multicomponent effects are important for the wind structure, however this is not able to consistently explain the difference between predicted mass-loss rate and mass-loss rate derived from observation for these stars. Similar to previous studies, we found the dependence of wind parameters on the metallicity. We conclude that the wind mass-loss rate significantly increases with metallicity as dM/dt ~ Z^0.67, whereas the wind terminal velocity on average depends on metallicity only slightly, namely v_infty ~ Z^0.06 (for studied stars).

We review basic physics of line-driven stellar winds of OB stars. We discuss elementary processes due to which stellar winds are accelerated on a microscopic level. We show how these microscopic processes may enable the outflow and how they determine wind properties on a macroscopic level. We discuss shortcomings of present wind theories and future wind model improvements.

We present a principle component analysis of multicolour light variations of the He-strong magnetic chemically peculiar star HD 37776 and its (O-C) diagram. The period increase rate found, (dP/dt) / P = 5.5 x 10^{-6}/year, is consistent with the conception of angular momentum loss through magnetically confined stellar wind.

We calculated spectral energy distribution of the helium strong chemically peculiar star HD 37776 using adequate model atmospheres. We show that the chemical peculiarity influences the stellar energy distribution. Consequently, spots of peculiar chemical composition on the surface of a rotating star may cause detectable light variations of the star. However, the observed light variations are not likely caused by an uneven surface distribution of helium, but may be due to spots of metals (mainly carbon).

We study stellar wind properties of selected late O stars in the Small Magellanic Cloud (SMC). We calculate NLTE line-driven wind models for these stars and compare predicted wind parameters with observed values. We found satisfactory agreement between theoretical and observed terminal velocities. On the other hand, predicted and observed mass-loss rates are in a good agreement only for higher mass-loss rates. For mass-loss rates lower than approximately 10^{-7} M_sun / year we found large discrepancy between theoretical and observed values. We propose a new explanation of this effect based on dynamical decoupling of some atoms. Finally, we study the dependence of wind terminal velocities and mass-loss rates on metallicity.

We study dynamical aspects of circumstellar environment around massive zero-metallicity first stars. For this purpose we apply our NLTE wind models. We show that the hydrogen-helium stellar wind from stationary massive first generation (Population III) stars (driven either by the line (bound-bound) or continuum (bound-free and free-free) transitions) is unlikely. The possibility of expulsion of chemically homogeneous wind and the role of minor isotopes are also discussed. Finally, we estimate the importance of hydrogen and helium lines for shutting off the initial accretion onto first stars and its influence on initial mass function of first stars.

We calculate multicomponent radiatively driven stellar wind models suitable for central stars of planetary nebulae. Some of these stellar winds may be adequately modelled using one-component models, however for some of them multicomponent models are necessary. We obtain a range of stellar parameters for different types of mass-loss.

We present an analysis of new spectroscopic observations of the bright Be star kappa Dra obtained at the Ondrejov observatory during 1992 - 2003 and UBV photometric observations secured at several observatories. The fundamental stellar parameters have been obtained from a comparison with a grid of NLTE model atmospheres. The best fit was found for Teff=14000K, log g = 3.5, and v sin i = 170km/s. These values together with a Hipparcos parallax lead to a stellar mass M=4.8+-0.8 M_solar and radius R=6.4+-0.5 R_solar. Long-term variations of kappa Dra were analysed using measurements of equivalent widths, central intensities, peak intensities of emission lines and emission peak velocity differences for H-alpha, H-beta, H-gamma, H-delta, and some helium, silicon, and iron lines. It turned out that the previously reported period of 23 years in the variation of the emission strength is probably a cyclic, not a strictly periodic phenomenon. An attempt to find out a period from all available records of the H-beta emission strength led to a value of 8044+-167 days (22.0 years) but the phase plots show that each cycle has a different shape and length. For the moment, the nature and origin of the disk around kappa Dra still remains unknown. From the comparison of the electronic spectra obtained at different phases of the long-term cycle and synthetic spectra it appears that there are no detectable changes in the photospheric part of the Balmer lines related to variations in the Balmer emission strength which could be attributed to an extended photosphere corresponding to inner parts of the disk, optically thick in continuum.

The Discrete Absorption Components (DACs) commonly observed in the ultraviolet lines of hot stars have previously been modelled by dynamical simulations of Corotating Interaction Regions (CIRs) in their line-driven stellar winds. Here we apply the kinematic DAC inversion method of Brown et al. to the hydrodynamical CIR models and test the reliability of the results obtained. We conclude that the inversion method is able to recover valuable information on the velocity structure of the mean wind and to trace movement of velocity plateaux in the hydrodynamical data, though the recovered density profile of the stream is correct only very near to the stellar surface.

New numerical models of line-driven stellar winds of late O stars are presented. Statistical equilibrium (NLTE) equations of the most abundant elements are solved. Properly obtained occupation numbers are used to calculate consistent radiative force and radiative heating terms. Wind density, velocity and temperature are calculated as a solution of model hydrodynamical equations. Contrary to other published models we account for a multicomponent wind nature and do not simplify the calculation of the radiative force (e.g. using force multipliers). We discuss the convergence behaviour of our models. The ability of our models to predict correct values of mass-loss rates and terminal velocities of selected late O stars (mainly giants and supergiants) is demonstrated. The systematic difference between predicted and observed terminal velocities reported in the literature has been removed. Moreover, we found good agreement between the theoretical wind momentum-luminosity relationship and the observed one for Cyg OB2 supergiants.

We calculate multicomponent stellar wind models with inclusion of a helium component applicable to He-rich and He-poor stars. We show that helium does not decouple from the stellar wind of He-rich stars due to its coupling to hydrogen and its ionized state. For He-poor stars helium may decouple from the stellar wind, however this effect is not able to explain the chemical peculiarity of these stars. We conclude that the explanation of chemical peculiarity of these stars based purely on helium or hydrogen decoupling from the stellar wind is unlikely.

We present a simple method for the solution of one-component and multicomponent hydrodynamic equations based on the Newton-Raphson method. We show that this method can be used for the solution of stationary hydrodynamic equations. This method has been used for the calculation of the low density stellar wind models for which the multicomponent nature of the wind influences the overall wind structure.

We calculate multicomponent line-driven wind models of stars at extremely low metallicity suitable for massive first generation stars. For most of the models we find that the multicomponent wind nature is not important for either wind dynamics or for wind temperature stratification. However, for stars with the lowest metallicities we find that multicomponent effects influence the wind structure. These effects range from pure heating to possible fallback of the nonabsorbing wind component. We present a simple formula for the calculation of metallicity for which the multicomponent effects become important. We show that the importance of the multicomponent nature of winds of low metallicity stars is characterised not only by the low density of driving ions, but also by lower mass-loss rate.

We study stability of isothermal two-component radiatively driven stellar winds against one-dimensional perturbations larger than the Sobolev length, and radiative-acoustic waves in such stellar winds. We perform linear perturbation analysis in comoving fluid-frames of individual components and obtain dispersion relation in the common fluid-frame. For high density winds the velocity difference between velocities of both components is relatively small and the wind is stable for radiative-acoustic waves discovered originally by Abbott, in accordance with the previous studies of the one-component wind. However, for such high density winds we found new types of waves including a special case of "frozen-in" wavy patterns. On the other hand, if the velocity difference between wind components is sufficiently large (for low density winds) then the multicomponent stellar wind is unstable even for large scale perturbations and ion runaway occurs. Thus, isothermal two-component stationary solutions of the radiatively line driven stellar wind with an abrupt lowering of the velocity gradient are unstable.

We show that the so-called Gayley-Owocki (Doppler) heating is important for the temperature structure of the wind of main sequence stars cooler than the spectral type O6. The formula for Gayley-Owocki heating is derived directly from the Boltzmann equation as a direct consequence of the dependence of the driving force on the velocity gradient. Since Gayley-Owocki heating deposits heat directly to the absorbing ions, we also investigated the possibility that individual components of the radiatively driven stellar wind have different temperatures. This effect is negligible in the wind of O stars, whereas a significant temperature difference takes place in the winds of main sequence B stars for stars cooler than B2. Typical temperature difference between absorbing ions and other flow components for such stars is of the order 10^3 K. However, in the case when passive component falls back onto the star the absorbing component reaches temperatures of order 10^6 K, which allows for emission of X-rays. Moreover, we compare our computed terminal velocities with the observed ones. We found quite good agreement between predicted and observed terminal velocities. The systematic difference coming from the using of the so called "cooking formula" has been removed.

An analysis of accurate radial velocities (RVs) of the Be star gamma Cas from 295 Reticon spectrograms secured between October 1993 and May 2000 allowed us to prewhiten the RVs for the long-term changes and to obtain the first orbital RV curve of this star. The orbital period is 203d59 and the orbit has an eccentricity of 0.26. The orbital motion is detectable even in the published velocities, based on photographic spectra. This implies that gamma Cas is a primary component of a spectroscopic binary. The secondary has a mass of about 1 solar mass, appropriate for a white dwarf or a neutron star, but it could also be a normal late-type dwarf. The ultimate solution of the dispute whether the observed X-ray emission is associated with the secondary or with the primary will need further dedicated studies.