We present spectroscopic observations of the Be/X-ray binary X Per obtained during the period December 2017 - March 2019 (MJD58095 - MJD58568). In December 2017 the H-alpha, H-beta, and HeI6678 emission lines were symmetric with violet-to-red peak ratio V/R=1. During the first half of the period (MJD58095 - MJD58360) the V/R-ratio decreased to 0.5 and the asymmetry developed simultaneously in all three lines. During the second half of the period (MJD58360 - MJD58570) the V/R-ratio increased, however with the HeI6678 line showing a markedly different behaviour than the Balmer lines. A third component with velocity >250 km/s appeared in the red side of HeI line profile. Possibly, during the first half of the period the entire disc become asymmetric and/or elliptical. The third component in HeI can be an indication for inner elliptical disc with eccentricity 0.4.

High-resolution observations in the region, centered at 4400 \AA\ and those of the lines HeII 4686, Hbeta and HeI 6678 of the spectrum of the symbiotic binary Z And were performed during its outburst in 2006. The line Hbeta had additional satellite high-velocity emission components situated on either side of its central peak. The lines of neutral helium presented two components, consisting of a nebular emission situated close to the reference wavelength and a highly variable P Cyg absorption. Close to the optical maximum the line HeII 4686 was weak emission feature, but with the fading of the light it changed into an intensive emission consisting of a central narrow component and a broad component with a low intensity. The lines of NIII and CIII were very broadened. We demonstrate that all of these groups of lines with very different profiles can be interpreted in the light of the same model, where a disc-shaped material surrounding the compact object collimates its stellar wind and gives rise to bipolar outflow.

The optical light of the symbiotic binary BF Cygni during its last eruption after 2006 shows orbital variations because of an eclipse of the outbursting compact object. The first orbital minimum is deeper than the following ones. Moreover, the Balmer profiles of this system acquired additional satellite components indicating bipolar collimated outflow at one time between the first and second orbital minima. This behaviour is interpreted in the framework of the model of collimated stellar wind from the outbursting object. It is supposed that one extended disc-like envelope covering the accretion disc of the compact object and collimating its stellar wind forms in the period between the first and second minima. The uneclipsed part of this envelope is responsible for the decrease of the depth of the orbital minimum. The calculated $UBVR_{C}I_{C}$ fluxes of this uneclipsed part are in agreement with the observed residual of the depths of the first and second orbital minima. The parameters of the envelope require that it is the main emitting region of the line H$\alpha$ but the H$\alpha$ profile is less determined from its rotation and mostly from other mechanisms. It is concluded that the envelope is a transient nebular region and its destruction determines the increase of the depth of the orbital minimum with fading of the optical light.

We present new radial velocities of the high-mass X-ray binary star 4U 2206+54 based on optical spectra obtained with the Coude spectrograph at the 2m RCC telescope at the Rozhen National Astronomical Observatory, Bulgaria in the period November 2011 -- July 2013. The radial velocity curve of the HeI $\lambda$6678 \AA\ line is modeled with an orbital period P$_{orb}$ = 9.568~d and an eccentricity of $e$ = 0.3. These new measurements of the radial velocity resolve the disagreements of the orbital period discussions.

We report high resolution spectral observations of the yellow symbiotic star V417 Cen obtained in 2004. We find that the equivalent widths of the emission lines decreased, while the brightness increased. The FWHMs and wavelengths of the emission lines do not change. We estimated the interstellar extinction towards V417 Cen as $E_{B-V}$ = 0.95 $\pm$0.10, using the KI interstellar line. Using the [O III] lines, we obtain a rough estimation of the density and the temperature in the forbidden lines region $\log N_e \approx 4.5 \pm$ 0.5 and $T_e=100000 \pm 25000$~K. Tidal interaction in this binary is also discussed.

We studied the average orbital modulation of various parameters (gamma-ray flux, H-alpha emission line, optical V band brightness) of the radio- and gamma-ray emitting Be/X-ray binary LSI+61303. Using the Spearman rank correlation test, we found highly significant correlations between the orbital variability of the equivalent width of the blue hump of the H-alpha and Fermi-LAT flux with a Spearman p-value 2e-5, and the equivalent widths ratio EW_B/EW_R and Fermi-LAT flux with p-value 9e-5. We also found a significant anti-correlation between Fermi-LAT flux and V band magnitude with p-value 7.10^{-4}. All these correlations refer to the average orbital variability, and we conclude that the H-alpha and gamma-ray emission processes in LSI+61303 are connected. The possible physical scenario is briefly discussed.

We report 138 spectral observations of the H-alpha emission line of the radio- and gamma-ray-emitting Be/X-ray binary LSI+61303 obtained during the period of September 1998 -- January 2013. From measuring various H-alpha parameters, we found that the orbital modulation of the H-alpha is best visible in the equivalent width ratio EW(B)/EW(R), the equivalent width of the blue hump, and in the radial velocity of the central dip. The periodogram analysis confirmed that the H-alpha emission is modulated with the orbital and superorbital periods. For the past 20 years the radius of the circumstellar disk is similar to the Roche lobe size at the periastron. It is probably truncated by a 6:1 resonance. The orbital maximum of the equivalent width of H-alpha emission peaks after the periastron and coincides on average with the X-ray and gamma-ray maxima. All the spectra are available upon request from the authors and through the CDS.

Line profiles containing indication of bipolar collimated outflow along with P Cyg absorption during phase of activity of several symbiotic systems are considered. The Hgamma profile of Z And during its 2006 outburst consisted of four groups of components. The profile of the HeI {\lambda} 5876 line of the Hen 3-1341 system during its 1998-2004 outburst had high-velocity satellite components on one hand and a broad P Cyg absorption on other hand. The Halpha and Hbeta profiles of the BF Cyg system during its 2006-2012 outburst had satellite components, observed for the first time, along with P Cyg absorptions. These profiles are interpreted in the light of a model related to the strong recurrent outbursts of Z And. The AG Dra system during its 2006 outburst did not contain indication of collimated outflow. Its profiles are interpreted in the light of the model related to the first outburst of 2000-2012 active phase of Z And.

Two observational campaigns were carried out during the eclipses of EE Cep in 2003 and 2008/9 to verify whether the eclipsing body in the system is indeed a dark disk and to understand the observed changes in the depth and durations of the eclipses. Multicolour photometric data and spectroscopic observations at both low and high resolution were collected. We numerically modelled the variations in brightness and colour during the eclipses. We tested models with different disk structure. We considered the possibility of disk precession. The complete set of observational data collected during the last three eclipses are made available to the astronomical community. Two blue maxima in the colour indices were detected during these two eclipses, one before and one after the photometric minimum. The first (stronger) blue maximum is simultaneous with a "bump" that is very clear in all the UBVRI light curves. Variations in the spectral line profiles seem to be recurrent during each cycle. NaI lines always show at least three absorption components during the eclipse minimum and strong absorption is superimposed on the H_alpha emission. These observations confirm that the eclipsing object in EE Cep system is indeed a dark, dusty disk around a low luminosity object. The primary appears to be a rapidly rotating Be star that is strongly darkened at the equator and brightened at the poles. Some of the conclusions of this work require verification in future studies: (i) a complex, possibly multi-ring structure of the disk in EE Cep; (ii) our explanation of the "bump" observed during the last two eclipses in terms of the different times of obscuration of the hot polar regions of the Be star by the disk; and (iii) our suggested period of the disk precession (~11-12 P_orb) and predicted depth of about 2 mag the forthcoming eclipse in 2014.

Results of the study of the symbiotic binary Z And during its recent active phase 2000 -- 2010 when it experienced a series of six optical outbursts are presented. High-resolution spectra obtained during the first and fourth outburst, which was the strongest one, have been analyzed. These data are compared with results of theoretical computations. The comparison provides information about the behaviour of the system during the entire active phase rather than during an individual outburst. In particular it was found fundamental difference between the first outburst, which opened the active phase, and the recurrent outbursts - namely, the presence of bipolar collimated optical outflow during some of the recurrent outbursts. A scenario that can explain all the spectroscopic phenomena observed during this active phase as well as previous active phases of Z And is proposed. The possibility to use this scenario for explanation of the line spectrum of other classical symbiotic stars during their active phases is motivated.

High resolution observations in the region of the lines Halpha, He II 4686 and Hgamma of the spectrum of the symbiotic binary Z And were performed during its small-amplitude brightening at the end of 2002. The profiles of the hydrogen lines were double-peaked. These profiles give a reason to suppose that the lines can be emitted mainly by an optically thin accretion disc. The Halpha line is strongly contaminated by the emission of the envelope, therefore for consideration of accretion disc properties we use the Hgamma line. The Halpha line had broad wings which are supposed to be determined mostly from radiation damping but high velocity stellar wind from the compact object in the system can also contribute to their appearance. The Hgamma line had a broad emission component which is assumed to be emitted mainly from the inner part of the accretion disc. The line He II 4686 had a broad emission component too, but it is supposed to appear in a region of a high velocity stellar wind. The outer radius of the accretion disc can be calculated from the shift between the peaks. Assuming, that the orbit inclination can ranges from 47$^\circ$ to 76$^\circ$, we estimate the outer radius as 20 - 50 R_sun. The behaviour of the observed lines can be considered in the framework of the model proposed for interpretation of the line spectrum during the major 2000 - 2002 brightening of this binary.

The photometric and spectroscopic observational campaign organized for the 2008/9 eclipse of EE Cep revealed features, which indicate that the eclipsing disk in the EE Cep system has a multi-ring structure. We suggest that the gaps in the disk can be related to the possible planet formation.

High resolution data in the region of the line Halpha have been obtained at the time of the light maximum and after it of the 2006 optical outburst of the symbiotic binary Z And. A blue-shifted absorption component indicating outflow velocity of about 1400 km/s as well as additional emission components with similar velocities, situated on the two sides of the main peak of the line were observed during that time. It is suggested that all of them are spectral signature of bipolar outflow, observed for the first time in the optical spectrum of this binary. The emission measure and the mass of the nebular part of the streams were estimated approximately.

Two-dimensional gas-dynamical modeling of the mass-flow structure is used to study the outburst development in the classical symbiotic star Z And. The stage-by-stage rise of the light during the outburst can be explained in the framework of the colliding winds model. We suggest a scenario for the development of the outburst and study the possible influence of the changes of the flow structure on the light of the system. The model variations of the luminosity due to the formation of a system of shocks are in good agreement with the observed light variations.

We report that the H_alpha emission line of RS Oph was strongly variable during our 2004 observations on a time scale of 1 month. The line consisted of both a double peaked central narrow component (FWHM ~ 220 km/s) and a strongly variable broad one (FWHM ~ 2000 km/s). The base of the H_alpha line was very broad with FWZI 4600 km/s on all spectra from 1986 to 2004. The variability of the broad component extends from -2000 to +2000 km/s. Most probably this is due to either blobs ejected from the white dwarf (with a typical blob mass estimated to be ~10^{-10} M_solar or a variable accretion disk wind. We also detected variability of the HeII4686 line on a time scale shorter than 1 day. The possible origin is discussed.

We report our search for variability in the H-alpha emission line of the recurrent novae T CrB with time resolution 10-15 minutes. This is comparable with the time scale of the photometric flickering observed in this object. This is the first time that observations of the short time scale variation in emission lines have been made in this object. On two nights (990106 and 990107) we detected statistically significant variability (at 0.99 confidence level) in the line profile of H-alpha. This variability is confined in the central part of the emission line (\pm 100 km/s) although FWZI(H-alpha) is ~800 km/s. The variability in the line profile is accompanied with variability of the total equivalent width, EW(H-alpha), \pm 8% for 990106, and \pm 6% for 990107 (calculated from the mean value of EW). Assuming Keplerian motion, the variability is generated at a distance ~20-30 R_Solar from the white dwarf, which is approximately the radius of the ring the stream of gas forms as it flows away from L_1. On three other nights we are only able to put upper limits on the variability, \Delta EW(H-alpha) \pm 2% for 980415, \pm 4% for 980802, and \pm 3% for 980803.

We analyze H-alpha observations of the recurrent nova T CrB obtained during the last decade. For the first time the H-alpha emission profile is analyzed after subtraction of the red giant contribution. Based on our new radial velocity measurements of the H-alpha emission line we estimate the component masses of T CrB. It is found that the hot component is most likely a massive white dwarf. We estimate the inclination and the component masses to be i~67 deg, Mwd = 1.37 +/-0.13 Msun and Msec=1.12 +/-0.23 Msun, respectively. The radial velocity of the central dip in the H-alpha profile changes nearly in phase with that of the red giant's absorption lines. This suggests that the dip is most likely produced by absorption in the giant's wind. Our observations cover an interval when the H-alpha and the U-band flux vary by a factor of ~6, while the variability in B and V is much smaller. Based on our observations, and archival ultraviolet and optical data we show that the optical, ultraviolet and H-alpha fluxes strongly correlate. We argue that the presence of an accretion disc can account for most of the observed properties of T CrB.

We present a comparative study of the circumstellar disks in Be/X-ray binaries and isolated Be stars based upon the H-alpha emission line. From this comparison it follows that the overall structure of the disks in the Be/X-ray binaries is similar to the disks of other Be stars, i.e. they are axisymmetric and rotationally supported. The factors for the line broadening (rotation and temperature) in the disks of the Be stars and the Be/X-ray binaries seem to be identical. However, we do detect some intriguing differences between the envelopes. On average, the circumstellar disks of the Be/X-ray binaries are twice as dense as the disks of the isolated Be stars. The different distribution of the Be/X-ray binaries and the Be stars seen in the full width at half maximum (FWHM) versus peak separation diagram indicates that the disks in Be/X-ray binaries have on average a smaller size, probably truncated by the compact object.