Several dozen optical echelle spectra demonstrate that HR 6819 is a hierarchical triple. A classical Be star is in a wide orbit with an unconstrained period around an inner 40 d binary consisting of a B3 III star and an unseen companion in a circular orbit. The radial-velocity semi-amplitude of 61.3 km/s of the inner star and its minimum (probable) mass of 5.0 Msun (6.3 +- 0.7 Msun) imply a mass of the unseen object of >= 4.2 Msun (>= 5.0 +- 0.4 Msun), that is, a black hole (BH). The spectroscopic time series is stunningly similar to observations of LB-1. A similar triple-star architecture of LB-1 would reduce the mass of the BH in LB-1 from ~70 Msun to a level more typical of Galactic stellar remnant BHs. The BH in HR 6819 probably is the closest known BH to the Sun, and together with LB-1, suggests a population of quiet BHs. Its embedment in a hierarchical triple structure may be of interest for models of merging double BHs or BH + neutron star binaries. Other triple stars with an outer Be star but without BH are identified; through stripping, such systems may become a source of single Be stars.

OGLE-LMC-DPV-065 is an interacting binary whose double-hump long photometric cycle remains hitherto unexplained. We analyze photometric time series available in archive datasets spanning 124 years and present the analysis of new high-resolution spectra. A refined orbital period is found of 10\fd0316267 $\pm$ 0\fd0000056 without any evidence of variability. In spite of this constancy, small but significant changes in timings of the secondary eclipse are detected. We show that the long period continuously decreases from 350 to 218 days during 13 years, then remains almost constant for about 10 years. Our study of radial velocities indicates a circular orbit for the binary and yields a mass ratio of 0.203 $\pm$ 0.001. From the analysis of the orbital light curve we find that the system contains 13.8 and 2.81 \msun\ stars of radii 8.8 and 12.6 \rsun\ and absolute bolometric magnitudes -6.4 and -3.0, respectively. The orbit semi-major axis is 49.9 \rsun\ and the stellar temperatures are 25460 K and 9825 K. We find evidence for an optically and geometrically thick disk around the hotter star. According to our model, the disk has a radius of 25 \rsun, central and outer vertical thickness of 1.6 \rsun\ and 3.5 \rsun, and temperature of 9380 K at its outer edge. Two shock regions located at roughly opposite parts of the outer disk rim can explain the light curves asymmetries. The system is a member of the double periodic variables and its relatively high-mass and long photometric cycle make it similar in some aspects to $\beta$ Lyrae.

We present a novel method for interpreting observations of high-mass X-ray binaries (HMXBs) based on a combination of spectroscopic data and numerical results from a radiation hydrodynamic model of stellar winds. We calculate synthetic Doppler tomograms of predicted emission in low/hard and high/soft X-ray states and compare them with Doppler tomograms produced using spectra of Cygnus X-1, a prototype of HMXBs. Emission from HMXBs is determined by local conditions within the circumstellar medium, namely density, temperature, and ionization state. These quantities depend strongly on the X-ray state of the systems. By increasing intensity of an X-ray emission produced by the compact companion in the HMXB model, we achieved a complete redistribution of the circumstellar medium in the vicinity of the modelled system. These changes (which simulate the transitions between two major spectral states) are also apparent in the synthetic Doppler tomograms which are in good agreement with the observations.

Aims: We have developed a new code for the three-dimensional time-dependent raditation hydrodynamic simulation of the stellar wind in interacting binaries to improve models of accretion in high-mass X-ray binaries and to quantitatively clarify the observed variability of these objects. We used the code to test the influence of various parameters on the structure and properties of circumstellar matter. Methods: Our code takes into account acceleration of the wind due to the Roche effective potential, Coriolis force, gas pressure, and (CAK-) radiative pressure in the lines and continuum of the supergiant radiation field that is modulated by its gravity darkening and by the photo-ionization caused by X-ray radiation from the compact companion. The parameters of Cygnus X-1 were used to test the properties of our model. Results: Both two- and three-dimensional numerical simulations show that the Coriolis force substantially influences the mass loss and consequently the accretion rate onto the compact companion. The gravitational field of the compact companion focuses the stellar wind, which leads to the formation of a curved cone-like gaseous tail behind the companion. The changes of X-ray photo-ionization of the wind material during X-ray spectral-state transitions significantly influence the wind structure and offer an explanation of the variability of Cygnus X-1 in optical observations (the H$\alpha$ emission).

We investigate the anisotropy of stellar winds in binaries to improve the models of accretion in high-mass X-ray binaries. We model numerically the stellar wind from a supergiant component of a binary in radial and three-dimensional radiation hydrodynamic approximation taking into account the Roche potential, Coriolis force, and radiative pressure in the continuum and spectral lines. The Coriolis force influences substantially the mass loss and thus also the accretion rate. The focusing of the stellar wind by the gravitational field of the compact companion leads to the formation of a gaseous tail behind the companion.

VO-KOREL is a web service exploiting the technology of Virtual Observatory for providing the astronomers with the intuitive graphical front-end and distributed computing back-end running the most recent version of Fourier disentangling code KOREL. The system integrates the ideas of the e-shop basket, conserving the privacy of every user by transfer encryption and access authentication, with features of laboratory notebook, allowing the easy housekeeping of both input parameters and final results, as well as it explores a newly emerging technology of cloud computing. While the web-based front-end allows the user to submit data and parameter files, edit parameters, manage a job list, resubmit or cancel running jobs and mainly watching the text and graphical results of a disentangling process, the main part of the back-end is a simple job queue submission system executing in parallel multiple instances of FORTRAN code KOREL. This may be easily extended for GRID-based deployment on massively parallel computing clusters. The short introduction into underlying technologies is given, briefly mentioning advantages as well as bottlenecks of the design used.

The first stars are assumed to be predominantly massive. Although, due to the low initial abundances of heavy elements the line-driven stellar winds are supposed to be inefficient in the first stars, these stars may loose a significant amount of their initial mass by other mechanisms. In this work, we study the evolution with a prescribed mass loss rate of very massive, galactic and pregalactic, Population III stars, with initial metallicities $Z=10^{-6}$ and $Z=10^{-9}$, respectively, and initial masses 100, 120, 150, 200, and 250$\,M_{\odot}$ during the hydrogen and helium burning phases. The evolution of these stars depends on their initial mass, metallicity and the mass loss rate. Low metallicity stars are hotter, compact and luminous, and they are shifted to the blue upper part in the Hertzprung-Russell diagram. With mass loss these stars provide an efficient mixing of nucleosynthetic products, and depending on the He-core mass their final fate could be either pair-instability supernovae or energetic hypernovae. These stars contributed to the reionization of the universe and its enrichment with heavy elements, which influences the subsequent star formation properties.

We studied both components of a slightly overlooked visual binary HR 1847 spectroscopically to determine its basic physical and orbital parameters. Basic stellar parameters were determined by comparing synthetic spectra to the observed echelle spectra, which cover both the optical and near-IR regions. New observations of this system used the Ond\v{r}ejov and Rozhen 2-m telescopes and their coud\'e spectrographs. Radial velocities from individual spectra were measured and then analysed with the code {\FOTEL} to determine orbital parameters. The spectroscopic orbit of HR 1847A is presented for the first time. It is a single-lined spectroscopic binary with a B-type primary, a period of 719.79 days, and a highly eccentric orbit with e=0.7. We confirmed that HR 1847B is a Be star. Its H\alpha emission significantly decreased from 2003 to 2008. Both components have a spectral type B7-8 and luminosity class IV-V.

The Virtual Observatory is a new technology of the astronomical research allowing the seamless processing and analysis of a heterogeneous data obtained from a number of distributed data archives. It may also provide astronomical community with powerful computational and data processing on-line services replacing the custom scientific code run on user's computers. Despite its benefits the VO technology has been still little exploited in stellar spectroscopy. As an example of possible evolution in this field we present an experimental web-based service for disentangling of spectra based on code KOREL. This code developed by P. Hadrava enables Fourier disentangling and line-strength photometry, i.e. simultaneous decomposition of spectra of multiple stars and solving for orbital parameters, line-profile variability or other physical parameters of observed objects. We discuss the benefits of the service-oriented approach from the point of view of both developers and users and give examples of possible user-friendly implementation of spectra disentangling methods as a standard tools of Virtual Observatory.

HD172189 is a spectroscopic eclipsing binary system with a rapidly-rotating pulsating delta Scuti component. It is also a member of the open cluster IC4756. These combined characteristics make it an excellent laboratory for asteroseismic studies. To date, HD172189 has been analysed in detail photometrically but not spectroscopically. For this reason we have compiled a set of spectroscopic data to determine the absolute and atmospheric parameters of the components. We determined the radial velocities (RV) of both components using four different techniques. We disentangled the binary spectra using KOREL, and performed the first abundance analysis on both disentangled spectra. By combining the spectroscopic results and the photometric data, we obtained the component masses, 1.8 and 1.7 solar masses, and radii, 4.0 and 2.4 solar radii, for inclination i = 73.2 degrees, eccentricity e = 0.28, and orbital period 5.70198 days. Effective temperatures of 7600 K and 8100 K were also determined. The measured vsini are 78 and 74 km/s, respectively, giving rotational periods of 2.50 and 1.55 days for the components. The abundance analysis shows [Fe/H] = -0.28 for the primary (pulsating) star, consistent with observations of IC4756. We also present an assessment of the different analysis techniques used to obtain the RVs and the global parameters.

The determination of pulsation velocities from observed spectra of Cepheids is needed for the Baade-Wesselink calibration of these primary distance markers. The applicability of the Fourier-disentangling technique for the determination of pulsation velocities of Cepheids and other pulsating stars is studied. The KOREL-code was modified to enable fitting of free parameters of a prescribed line-profile broadening function corresponding to the radial pulsations of the stellar atmosphere. It was applied to spectra of delta Cep in the H-alpha region observed with the Ondrejov 2-m telescope. The telluric lines were removed using template-constrained disentangling, phase-locked variations of line-strengths were measured and the curves of pulsational velocities obtained for several spectral lines. It is shown that the amplitude and phase of the velocities and line-strength variations depend on the depth of line formation and the excitation potential. The disentangling of pulsations in the Cepheid spectra may be used for distance determination.

In this document a review of the author's method of Fourier disentangling of spectra of binary and multiple stars is presented for the purpose of the summer school organized at Ondrejov observatory in September 2008. Related methods are also discussed and some practical hints for the use of the author's code KOREL and related auxiliary codes with examples are given.

We present an examination of high resolution, ultraviolet spectroscopy from Hubble Space Telescope of the photospheric spectrum of the O-supergiant in the massive X-ray binary HD 226868 = Cyg X-1. We analyzed this and ground-based optical spectra to determine the effective temperature and gravity of the O9.7 Iab supergiant. Using non-local thermodynamic equilibrium (non-LTE), line blanketed, plane parallel models from the TLUSTY grid, we obtain T_eff = 28.0 +/- 2.5kK and log g > 3.00 +/- 0.25, both lower than in previous studies. The optical spectrum is best fit with models that have enriched He and N abundances. We fit the model spectral energy distribution for this temperature and gravity to the UV, optical, and IR fluxes to determine the angular size of and extinction towards the binary. The angular size then yields relations for the stellar radius and luminosity as a function of distance. By assuming that the supergiant rotates synchronously with the orbit, we can use the radius - distance relation to find mass estimates for both the supergiant and black hole as a function of the distance and the ratio of stellar to Roche radius. Fits of the orbital light curve yield an additional constraint that limits the solutions in the mass plane. Our results indicate masses of 23^{+8}_{-6} M_sun for the supergiant and 11^{+5}_{-3} M_sun for the black hole.

Context: The technique of disentangling has been applied to numerous high-precision studies of spectroscopic binaries and multiple stars. Although, its possibilities have not yet been fully understood and exploited. Aims: Theoretical background aspects of the method, its latest improvements and hints for its use in practice are explained in this series of papers. Methods: In this first paper of the series, we discuss spectral-resolution limitations due to a discrete representation of the observed spectra and introduce a new method how to achieve a precision higher than the step of input-data binning. Results: Based on this principle, the latest version of the KOREL code for Fourier disentangling achieves an increase in precision for an order of magnitude.

We present the results of the spectroscopic observations of HDE 226868, the optical counterpart to the black hole X-ray binary Cyg X-1, from 2001 to 2006. We analyze the variabilities of the two components in the complex H$\alpha$ line: one P-Cygni shaped component which follows the motion of the supergiant and another emission component moving with an antiphase orbital motion relative to the supergiant, which is attributed to a focused-stellar wind. The results of KOREL disentangling of our spectra indicate that the focused stellar wind is responsible for the major part of the variability of the H$\alpha$ emission line. The emission of the supergiant component had a small difference between the low/hard and high/soft states, while the focused wind component became strong in the low/hard state and weak in the high/soft state. The wind is nearly undisturbed by the X-ray photoionization during the low/hard state. However, during the high/soft state, the X-rays from the compact object could decelerate the line-driven wind and result in a high mass accretion rate, due to the effect of the X-ray photoionization. The X-ray illuminating could also change the temperature profile of the stellar wind and increase its temperature, and thus decrease the H$\alpha$ emissivity of the wind, which could explain the H$\alpha$ variabilities of Cyg X-1 during different X-ray states.

We present results from Hubble Space Telescope UV spectroscopy of the massive X-ray binary system, HD226868 = Cyg X-1. The spectra were obtained at both orbital conjunction phases in two separate runs in 2002 and 2003 when the system was in the X-ray high/soft state. The stellar wind lines suffer large reductions in strength when the black hole is in the foreground due to the X-ray ionization of the wind ions. A comparison of HST and archival IUE spectra shows that similar photoionization effects occur in both the X-ray states. We constructed model UV wind line profiles assuming that X-ray ionization occurs everywhere in the wind except the zone where the supergiant blocks the X-ray flux. The good match between the observed and model profiles indicates that the wind ionization extends to near to the hemisphere of the supergiant facing the X-ray source. The H-alpha emission strength is generally lower in the high/soft state compared to the low/hard state, but the He II 4686 emission is relatively constant between states. The results suggest that mass transfer in Cyg X-1 is dominated by a focused wind flow that peaks along the axis joining the stars and that the stellar wind contribution is shut down by X-ray photoionization effects. The strong stellar wind from the shadowed side of the supergiant will stall when Coriolis deflection brings the gas into the region of X-ray illumination. This stalled gas component may be overtaken by the orbital motion of the black hole and act to inhibit accretion from the focused wind. The variations in the strength of the shadow wind component may then lead to accretion rate changes that ultimately determine the X-ray state.

The star HDE 226868 known as an optical counterpart of the black hole candidate Cyg X-1 has been observed in H_alpha region using spectrograph at Ondrejov 2-m telescope. The orbital parameters are determined from HeI-line by means of the author's method of Fourier disentangling. Preliminary results are also presented of disentangling the H_alpha-line into a P-Cyg profile of the (optical) primary and an emission profile of the circumstellar matter (and a telluric component).

The use of a method of spectra disentangling for telluric lines is explained in detail, with a particular emphasis on high-precision radial-velocity measurements for the search for extrasolar planets. New improvements to the method are introduced.

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.

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.