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.

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.

Context: Multiple stellar systems are unique laboratories for astrophysics. Analysis of their orbital dynamics may reveal invaluable information about the physical properties of the participating stars. Unfortunately, there are only a few known and well described multiple systems, this is even more so for systems located outside the Milky Way galaxy. A particularly interesting situation occurs when the inner binary in a compact triple system is eclipsing. This is because the stellar interaction, typically resulting in precession of orbital planes, may be observable as a variation of depth of the eclipses on a long timescale. Aims: We aim to present a novel method to determine compact triples using publicly available photometric data from large surveys. Here we apply it to eclipsing binaries (EBs) in Magellanic Clouds from OGLE III database. Methods: We analyzed light curves (LCs) of 26121 LMC and 6138 SMC EBs with the goal to identify those for which the orbital inclination varies in time. Archival LCs of the selected systems, when complemented by our own observations with Danish 1.54m telescope, were thoroughly analyzed using the PHOEBE program. Time dependence of the EB's inclination was described using the theory of orbital-plane precession. By observing the parameter-dependence of the precession rate, we were able to constrain the third companion mass and its orbital period around EB. Results: We identified 58 candidates of new compact triples in Magellanic Clouds. This is the largest published sample of such systems so far. Eight of them were analyzed thoroughly and physical parameters of inner binary were determined together with an estimation of basic characteristics of the third star. These data may provide important clues about stellar formation mechanisms for objects with different metalicity than found in our galactic neighborhood.

The study of chemically peculiar (CP) stars in open clusters provides valuable information about their evolutionary status. Their detection can be performed using the $\Delta a$ photometric system, which maps a characteristic flux depression at $\lambda \sim 5200$ \AA. This paper aims at studying the occurrence of CP stars in the earliest stages of evolution of a stellar population by applying this technique to Hogg 16, a very young Galactic open cluster (about 25 Myr). We identified several peculiar candidates: two B-type stars with a negative $\Delta a$ index (CD-60 4701, CPD-60 4706) are likely emission-line (Be) stars, even though spectral measurements are necessary for a proper classification of the second one; a third object (CD-60 4703), identified as a Be candidate in literature, appears to be a background B-type supergiant with no significant $\Delta a$ index, which does not rule out the possibility that it is indeed peculiar as the normality line of $\Delta a$ for supergiants has not been studied in detail yet. A fourth object (CD-60 4699) appears to be a magnetic CP star of 8 M$_\odot$, but obtained spectral data seem to rule out this hypothesis. Three more magnetic CP star candidates are found in the domain of early F-type stars. One is a probable nonmember and close to the border of significance, but the other two are probably pre-main sequence cluster objects. This is very promising, as it can lead to very strong constraints to the diffusion theory. Finally, we derived the fundamental parameters of Hogg 16 and provide for the first time an estimate of its metal content.

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.

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.