Since the end of 2018, the Transiting Exoplanet Survey Satellite (TESS) provides high-quality space data on stellar photometry to the astronomical community. We present the results of an analysis of TESS photometric data for known slowly rotating magnetic chemically peculiar (mCP) stars. In general, mCP stars show an inhomogeneous distribution of elements in their stellar atmospheres that leads to spectroscopic (line profile) and photometric (light curve) variations over the rotation period. In the frame of the oblique magnetic rotator (OMR) model, patches of enhanced chemical abundance on the stellar surface reveal the frequency of stellar rotation. Using this approach, we have compiled a list of slowly rotating mCP stars with rotation periods longer than two days from the analysis of the photometric data provided by TESS for the first eight sectors of observations. Slowly rotating mCP stars usually possess a hydrodynamically stable stellar atmosphere where a magnetic field can amplify the process of atomic diffusion and leads to the horizontal and vertical stratification of chemical abundances.

HD~63401 is known magnetic chemically peculiar (mCP) star that shows a slow rotation and probably possesses a hydrodynamically stable stellar atmosphere. In the latter case the atomic diffusion mechanism enforced by the magnetic field can lead to stratification of elemental abundance with optical depth. HD~63401 was recently observed with the space telescope \textit{TESS}, and its light curve shows distinct variability, which is usually detected in the $\alpha^2$ CVn type stars. Based on the analysis of the light curve we derived its rotational period $P$ = 2.414 $\pm$ 0.018~d and studied variability of the effective temperature with this period. The estimates of the effective temperature and surface gravity were obtained from the best fit of Balmer line profiles observed in seven high-resolution spectra acquired with the spectropolarimeter ESPaDOnS for HD~63401. The same spectra were used to perform an abundance analysis employing the modified ZEEMAN2 code. We have found that He, C, P, V, Y, and Dy are in deficit in stellar atmosphere of HD~63401, while Na, Al, Si, Fe, Zn, and Sr are significantly overabundant.

This research is our forth article related to the topic of cyclic O-C changes, third components as the physical process that cause these changes and elliptical orbits of the third components. Here five more eclipsing binary stars were investigated: EQ Tau, IR Cas, IV Cas, RY Aqr and RZ Com. All of them have cyclic O-C curve with superposition of parabolic trend. We computed the mass transfer rate, minima possible mass of the third component and their errors for each of the researched stars.

In our research, we studied 9 eclipsing binary stars: AR Lac, U CrB, S Equ, SU Boo, VV UMa, WW Gem, YY Eri, V0404 Lyr, HP Aur. We collected large sets of moments of minima from BRNO and observational data from AAVSO databases. Then, we obtained moments of minima for AAVSO observations (totally - 397 minima) using method of approximation with symmetric polynomial, realized in software MAVKA. This software was provided by Kateryna D. Andrych and Ivan L. Andronov. After we combined obtained moments of minima with data taken from BRNO and plotted O-C diagrams. For all stars these diagrams represented sinusoidal-like oscillations with superposition of linear (for SU Boo, VV UMa, WW Gem, V0404 Lyr and HP Aur) or parabolic trend (for AR Lac, U CrB, S Equ and YY Eri). The oscillations could be described as presence of the third component, which does not take part in eclipses, but causes the well-known light-time effect (LTE). That effect could be easily detected using long data series of observations. However, the oscillations have clear asymmetry, which we interpreted as elliptical shape of the third components orbit. Parabolic trend we explained as mass transfer between components of binary system. For all stars we computed minimal possible mass of the third component. In addition, we developed our own code in computed language Python and using it we computed orbital elements of the third component. Moreover, for stars with parabolic trend we calculated rate of the mass transfer. Finally, for all computed values we estimated errors.

In this research, five eclipsing binary stars were studied: AB Cas, AF Gem, AR Boo, BF Vir and CL Aur. The large sets of moments of minima were used: from the international BRNO database and amateur observations from the database AAVSO. Firstly, moments of minima for AAVSO observations were obtained (totally - 222 minima). The software MAVKA was used. It was kindly provided by Kateryna D. Andrych and Ivan L. Andronov (2019OEJV..197...65A) and approximation with various methods in order to find the best fit. Then all obtained moments of minima were combined and O-C diagrams were plotted. For all stars these diagrams represented sinusoidal-like oscillations with superposition of parabolic trend. One of the possible reasons for such oscillations could be presence of well-known light-time effect (LTE) caused by third component with elliptical orbit. Parabolic trend was explained as mass transfer between components of binary system. For all these stars we computed possible mass of the third component, orbital elements, mass transfer rate and errors for all computed values.

We investigated two variable stars AH Tau and ZZ Cas. They are eclipsing binaries of W UMa and $\beta$ Lyr types. The period changes with time. The reason for steady changes of the period could be the mass transfer (the flow of matter) between components of these stellar systems. For ZZ Cas the changes of the period are cyclic caused by the presence of the hypothetical third component (either a small star or a large planet). The cyclic changes of the period for AH Tau superimpose on steady ones (the period decrease) indicating for a third component and mass transfer. The third components do not take part in the eclipses. However, due to their gravity, they make the visible close binary systems rotate and become closer or further to an observer, causing the Light-Time Effect (LTE). Generally, an orbit of a third component is not a circle, but an ellipse and it is inclined relatively to the observer's line of sight. Using the O-C curve, we estimated the parameters of a third component's orbit, such as a semi-major axis, an eccentricity, angles of orientation and a period of a third component's rotation. We elaborated a modeling program using the computer language Python. We calculated errors of third component's orbit parameters and even estimated its mass. The values of masses of the third components within errors of calculations show that the third components are probably stars. All these calculations were made using all available data from international database BRNO (Brno Regional Network for Observers). Moreover, we used moments of minima which we calculated as the result of observation processing from the AAVSO database.

In our research, we have studied 6 close binary stellar systems which are eclipsing variables of $\beta$ Lyrae and W Ursae Majoris types of variability. We have studied their O-C curves. To build them, we used moments of minima, listed in the BRNO database. Also, we used ones, which we obtained as the result of processing of observational data taken from AAVSO database. As the result of the O-C analysis, we detected that all of these stars have parabolic O-C curves, which might be caused by mass transfer from one component to another. In an addition, 3 from researched stars (U Peg, V0523 Cas and WZ Cyg) have superimposed parabolic and cyclic O-C changes that could be caused by presence of the third components in the systems. Also, we calculated minimal possible masses of third components, rates of the mass transfer for these systems and corrected the ephemerides.