The upper main sequence stars CU Virginis is the most enigmatic object among magnetic chemically peculiar (mCP) stars. It is an unusually fast rotator showing strictly periodic light variations in all regions of the electromagnetic spectrum, as well as spectroscopic and spectropolarimetric changes. At same time, it is also the first radio main-sequence pulsar. Exploiting information hidden in phase variations, we monitored the secular oscillation of the rotational period during the last 53 years. Applying own phenomenological approach, we analyzed 37\,975 individual photometric and spectroscopic measurements from 72 data sources and improved the O-C model. All the relevant observations indicate that the secular period variations can be well approximated by the fifth degree polynomial.

The inhomogeneous surface distribution of heavy elements is known to cause periodic light variability of magnetic chemically peculiar stars. It is unclear to what extent the same paradigm applies to mercury-manganese (HgMn) stars. We aim to model the photometric variability of the HgMn star {\phi}Phe using abundance maps obtained from high-resolution spectroscopy and to study how this variability evolves with time. We compute a grid of ATLAS12 model atmospheres and the corresponding SYNSPEC synthetic spectra. Interpolating within this grid and integrating the specific intensity over the visible stellar surface at different rotational phases, we obtain theoretical light curves of the star. We predict the variability of {\phi} Phe in the ultraviolet and in the visible spectral regions with amplitude of the order of millimagnitudes, mainly caused by absorption in lines of yttrium, chromium, and titanium. We also show how this variability is affected by changes of the distribution of the heavy elements over time. The main characteristics of the predicted light variability of {\phi} Phe correspond roughly to the variability of the star observed with the TESS satellite.

HD\,37776 (V901\,Ori, B2\,Vp), also known as Landstreet's Star, is possibly the most remarkable magnetic chemically peculiar (mCP) star known. Zeeman Doppler Imaging revealed this young, rapidly rotating star's surface magnetic field to be not only the strongest ($\sim 30$ kG) of the He-strong class of hot mCP stars but also by far the most topologically complex. In contrast to the overwhelming majority of mCP stars, which are well described by tilted dipoles, Landstreet's Star's non-axisymmetric surface magnetic field is entirely dominated by high-order spherical harmonics. It is one of the handful of stars for which rotational period change has been measured, and over the past two decades of monitoring, the object has demonstrated an unexpected acceleration in its rotation that so far defies explanation. Recently acquired TESS data have provided a photometric data set of unprecedented precision. These data have revealed a highly stable yet multi-featured light curve, making Landstreet's Star the prototype of hot mCP stars whose light curves are difficult to reproduce using the standard model of chemical/photometric spots modulated by solid-body rotation.

We present the discovery of a unique object, a chemically peculiar Ap-type star showing $\delta$ Scuti pulsations which is bound in an eclipsing binary system with an orbital period shorter than 3 days. HD 99458 is, therefore, a complex astrophysical laboratory opening doors for studying various, often contradictory, physical phenomena at the same time. It is the first Ap star ever discovered in an eclipsing binary. The orbital period of 2.722 days is the second shortest among all known chemically peculiar (CP2) binary stars. Pulsations of $\delta$ Scuti type are also extremely rare among CP2 stars and no unambiguously proven candidate has been reported. HD 99458 was formerly thought to be a star hosting an exoplanet, but we definitely reject this hypothesis by using photometric observations from the K2 mission and new radial velocity measurements. The companion is a low-mass red dwarf star ($M_{2}=0.45(2)$ M$_{\odot}$) on an inclined orbit ($i=73.2(6)$ degrees) that shows only grazing eclipses. The rotation and orbital periods are synchronized, while the rotation and orbital axes are misaligned. HD 99458 is an interesting system deserving of more intense investigations.

Massive stars are characterized by a significant loss of mass either via spherically symmetric stellar winds or pre-explosion pulses, or by aspherical forms of circumstellar matter (CSM) such as bipolar lobes or outflowing circumstellar equatorial disks. A significant fraction of most massive stars end their lives by a core collapse; supernovae (SNe) are always located inside circumstellar envelopes created by their progenitors. We study the dynamics and thermal effects of collision between expanding ejecta of SNe and CSM that may be formed during, for example, a sgB[e] star phase, a luminous blue variable phase, around PopIII stars, or by various forms of accretion. For time-dependent hydrodynamic modeling we used our own code built with a finite volumes method. The code is highly efficient for calculations of shocks and physical flows with large discontinuities. The initial geometry of the disks corresponds to a density structure of a material that orbits in Keplerian trajectories. We examine the behavior of the density, pressure, velocity of expansion, and temperature structure in the interaction zone under various geometrical configurations of dense equatorial disks. Our `low density' model shows significant asphericity in the case of the disk mass-loss rate $\dot{M}_\text{csd}=10^{-6}\,M_\odot\,\text{yr}^{-1}$. The models also show the zones of overdensity in the SN - disk contact region and indicate the development of Kelvin-Helmholtz instabilities within the zones of shear between the disk and the more freely expanding material outside the disk.

We present a comprehensive overview of the properties of a sample of 41 magnetic chemically peculiar stars that have been recently identified in the Kepler field by our team (H\"{u}mmerich et al. 2018). The stars populate the whole age range from zero-age to terminal-age main sequence in the mass interval from 1.5 to 4 Ms. Several of the studied objects exhibit a hitherto unobserved wealth of detail in their light curves indicative of persisting complex surface structures. Monoperiodic variability and light curve stability were identified as cardinal criteria for selecting mCP star candidates among early-type objects in photometric surveys. Subsequent studies will be concerned with an exhaustive follow-up analysis of the new mCP stars, which we expect to lead to new insights on the physics of the CP star phenomenon.

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.

The presence of ultra-high excitation (UHE) absorption lines (e.g., O VIII) in the optical spectra of several of the hottest white dwarfs poses a decades-long mystery and is something that has never been observed in any other astrophysical object. The occurrence of such features requires a dense environment with temperatures near $10^6$K, by far exceeding the stellar effective temperature. Here we report the discovery of a new hot wind white dwarf, GALEXJ014636.8+323615. Astonishingly, we found for the first time rapid changes of the equivalent widths of the UHE features, which are correlated to the rotational period of the star ($P=0.242035$d). We explain this with the presence of a wind-fed circumstellar magnetosphere in which magnetically confined wind shocks heat up the material to the high temperatures required for the creation of the UHE lines. The photometric and spectroscopic variability of GALEXJ014636.8+323615 can then be understood as consequence of the obliquity of the magnetic axis with respect to the rotation axis of the white dwarf. This is the first time a wind-fed circumstellar magnetosphere around an apparently isolated white dwarf has been discovered and finally offers a plausible explanation of the ultra hot wind phenomenon.

Context. Evolution of massive stars is affected by a significant loss of mass either via (nearly) spherically symmetric stellar winds or by aspherical mass-loss mechanisms, namely the outflowing equatorial disks. However, the scenario that leads to the formation of a disk or rings of gas and dust around massive stars is still under debate. Aims. We study the hydrodynamic and thermal structure of optically thick, dense parts of outflowing circumstellar disks that may be formed around various types of critically rotating massive stars, for example, Be stars, B[e] supergiant (sgB[e]) stars or Pop III stars. Methods. We specify the optical depth of the disk along the line-of-sight from stellar poles. Within the optically thick dense region we calculate the vertical disk thermal structure using the diffusion approximation while for the optically thin outer layers we assume a local thermodynamic equilibrium with the impinging stellar irradiation. We use two of our own types of hydrodynamic codes: two-dimensional operator-split numerical code and unsplit code based on the Roe's method. Results. Our models show the geometric distribution and contribution of viscous heating that begins to dominate in the central part of the disk. In the models of dense viscous disks the viscosity increases the central temperature up to several tens of thousands of Kelvins. The high mass-loss rates and high viscosity lead to instabilities with significant waves or bumps in density and temperature in the very inner disk region. Conclusions. The two-dimensional radial-vertical models of dense outflowing disks including the full Navier-Stokes viscosity terms show very high temperatures that are however limited to only the central disk cores inside the optically thick area, while near the edge of the optically thick region the temperature may be low enough for the existence of neutral hydrogen.

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.

Context: The presence of heavier chemical elements in stellar atmospheres influences the spectral energy distribution (SED) of stars. An uneven surface distribution of these elements, together with flux redistribution and stellar rotation, are commonly believed to be the primary causes of the variability of chemically peculiar (CP) stars. Aims: We aim to model the photometric variability of the CP star PHI Dra based on the assumption of inhomogeneous surface distribution of heavier elements and compare it to the observed variability of the star. We also intend to identify the processes that contribute most significantly to its photometric variability. Methods: We use a grid of TLUSTY model atmospheres and the SYNSPEC code to model the radiative flux emerging from the individual surface elements of PHI Dra with different chemical compositions. We integrate the emerging flux over the visible surface of the star at different phases throughout the entire rotational period to synthesise theoretical light curves of the star in several spectral bands. Results: The synthetic light curves in the visible and in the near-UV regions are in very good agreement with the observed variability of the star. The lack of usable far-UV measurements of the star precludes making any conclusions about the correctness of our model in this spectral region. We also obtained 194 new BVRI observations of PHI Dra and improved its rotational period to P=1.716500(2). Conclusions: We show that the inhomogeneous distribution of elements, flux redistribution, and rotation of the star are fully capable of explaining the stellar variability in the visible and the near-UV regions. The flux redistribution is mainly caused by bound--free transitions of silicon and bound--bound transitions of iron.

Sigma Ori E (HD 37479) is the prototypical helium-strong star shown to harbor a strong magnetic field, as well as a magnetosphere consisting of two clouds of plasma. The observed optical (ubvy) light curve of sigma Ori E is dominated by eclipse features due to circumstellar material, however, there remain additional features unexplained by the Rigidly Rotating Magnetosphere (RRM) model of Townsend & Owocki. Using the technique of magnetic Doppler imaging (MDI), spectropolarimetric observations of sigma Ori E are used to produce maps of both the magnetic field topology and various elemental abundance distributions. We also present an analysis utilizing these computed MDI maps in conjunction with NLTE TLUSTY models to study the optical brightness variability of this star arising from surface inhomogeneities. It has been suggested that this physical phenomena may be responsible for the light curve inconsistencies between the model and observations.

During their evolution massive stars can reach the phase of critical rotation when a further increase in rotational speed is no longer possible. Direct centrifugal ejection from a critically or near-critically rotating surface forms a gaseous equatorial decretion disk. Anomalous viscosity provides the efficient mechanism for transporting the angular momentum outwards. The outer part of the disk can extend up to a very large distance from the parent star. We study the evolution of density, radial and azimuthal velocity, and angular momentum loss rate of equatorial decretion disks out to very distant regions. We investigate how the physical characteristics of the disk depend on the distribution of temperature and viscosity. We calculated stationary models using the Newton-Raphson method. For time-dependent hydrodynamic modeling we developed the numerical code based on an explicit finite difference scheme on an Eulerian grid including full Navier-Stokes shear viscosity. The sonic point distance and the maximum angular momentum loss rate strongly depend on the temperature profile and are almost independent of viscosity. The rotational velocity at large radii rapidly drops accordingly to temperature and viscosity distribution. The total amount of disk mass and the disk angular momentum increase with decreasing temperature and viscosity. The time-dependent one-dimensional models basically confirm the results obtained in the stationary models as well as the assumptions of the analytical approximations. Including full Navier-Stokes viscosity we systematically avoid the rotational velocity sign change at large radii. The unphysical drop of the rotational velocity and angular momentum loss at large radii (present in some models) can be avoided in the models with decreasing temperature and viscosity.

The majority of magnetic chemically peculiar (mCP) stars exhibit periodic light, magnetic, radio, and spectroscopic variations that can be adequately modelled as a rigidly-rotating main-sequence star with persistent surface structures. Nevertheless, there is a small sample of diverse mCP stars whose rotation periods vary on timescales of decades while the shapes of their phase curves remain unchanged. Alternating period increases and decreases have been suspected in the hot CP stars CU Vir and V901 Ori, while rotation in the moderately cool star BS Cir has been decelerating. These examples bring new insight into this theoretically unpredicted phenomenon. We discuss possible causes of such behaviour and propose that dynamic interactions between a thin, outer, magnetically-confined envelope braked by the stellar wind, and an inner faster-rotating stellar body are able to explain the observed rotational variability

Massive stars may during their evolution reach the phase of critical rotation when the further increase in rotational speed is no longer possible. The ejection of matter in the equatorial region forms the gaseous outflowing disk, which allows the star to remove the excess of its angular momentum. The outer part of the disk can extend up to large distance from the parent star. We study the evolution of density, radial and azimuthal velocity and angular momentum loss rate of equatorial decretion disk up to quite distant outer regions. We calculate the evolution of density, radial and azimuthal velocity from the initial Keplerian state until the disk approaches the final stationary state.