We investigate in detail the LiTE of the most probable binary candidate TU UMa. We used a new code applied on previously available and newly determined maxima timings to confirm binarity and refine parameters of the orbit of the RRab component in the binary system. The binary hypothesis is also tested using RV measurements. We used new approach to determine brightness maxima timings based on template fitting. This can also be used on sparse or scattered data. This approach was successfully applied on measurements from different sources. To determine the orbital parameters of the double star TU UMa, we developed a new code to analyse LiTE that also includes secular variation in the pulsation period. Its usability was successfully tested on CL Aur, an eclipsing binary with mass-transfer in a triple system that shows similar changes in the O-C diagram. Since orbital motion would cause systematic shifts in mean RVs (dominated by pulsations), we computed and compared our model with centre-of-mass velocities. They were determined using high-quality templates of RV curves of RRab stars. Maxima timings adopted from the GEOS database (168) together with those newly determined from sky surveys and new measurements (85) were used to construct an O-C diagram spanning almost five proposed orbital cycles. This data set is three times larger than data sets used by previous authors. Modelling of the O-C dependence resulted in 23.3-year orbital period, which translates into a minimum mass of the second component of about 0.33 M_s. Secular changes in the pulsation period of TU UMa over the whole O-C diagram were satisfactorily approximated by a parabolic trend with a rate of -2.2 ms/yr. To confirm binarity, we used RV measurements from nine independent sources. Although our results are convincing, additional long-term monitoring is necessary to unambiguously confirm the binarity of TU UMa.

The trustworthiness of orbital period analyses of eclipsing binaries strictly depends on the correctness of the observed mid-eclipse time determination, as well the reliability of its uncertainty estimation. The majority of them has been determined by means of the Kwee-van Woerden method (KWM). There are also other possibilities: to use physical models of eclipsing binaries or light curve templates and to determine mid-eclipse times using the least square method (LSM). We compared results yielded by both methods by means of a computer simulation on the synthetic model of AR\,Aur primary minimum. Minima times determined by the KWM and the exact LSM approach are nearly the same, while the scatter of LSM times is always smaller than the scatter of KWM times. KWM uncertainties are systematically underestimated. We think that the time is ripe for the Kwee-van Woerden method to retire.

In-depth analysis of eclipsing binary (EB) observational data collected for several decades can inform us about a lot of astrophysically interesting processes taking place in the systems. We have developed a wide-ranging method for the phenomenological modelling of eclipsing binary phase curves that enables us to combine even very disparate sources of phase information. This approach is appropriate for the processing of both standard photometric series of eclipses and data from photometric surveys of all kind. We conclude that mid-eclipse times, determined using the latest version of our 'hi-fi' phenomenological light curve models, as well as their accuracy, are nearly the same as the values obtained using much more complex standard physical EB models.

We present development of the collimated bipolar jets from the symbiotic prototype Z And that appeared and disappeared during its 2006 outburst. In 2006 July Z And reached its historical maximum at U ~ 8.0. During this period, rapid photometric variations with Dm ~ 0.06 mag on the timescale of hours developed. Simultaneously, high-velocity satellite components appeared on both sides of the H-alpha and H-beta emission line profiles. They were launched asymmetrically with the red/blue velocity ratio of 1.2 - 1.3. From about mid-August they became symmetric. Their spectral properties indicated ejection of bipolar jets collimated within an average opening angle of 6.1 degrees. We estimated average outflow rate via jets to dM(jet)/dt ~ 2xE10-6(R(jet)/1AU)**(1/2) M(Sun)/year, during their August - September maximum, which corresponds to the emitting mass in jets, M(jet, emitting) ~ 6xE-10(Rjet)/1AU)^{3/2} M(Sun). During their lifetime, the jets released the total mass of M(jet, total) approx 7.4x1E-7 M(Sun). Evolution in the rapid photometric variability and asymmetric ejection of jets around the optical maximum can be explained by a disruption of the inner parts of the disk caused by radiation-induced warping of the disk.