The VRI passbands light curves of V344 Lac were presented and analyzed by using the latest version of the W-D code. The observed spectrum reveals that V344 Lac is not an A3 type but would be a later F type star according to the yielded temperature. The results of solution show that V344 Lac is an A-subtype contact binary, with a mediate photometric mass ratio of 0.387(0.003) and a mediate contact factor of 44.6(3.0)%. Based on the parallax given by Gaia, the parameters of the components are estimated as: M1 = 1.16 Ms, M2 = 0.45 Ms, R1 = 1.31 Rs, R2 = 0.88 Rs, L1 = 2.512 Ls, L2 = 1.057 Ls. The period investigation indicates that V344 Lac may have an eccentric orbital oscillation, with P3 = 12.4(0.5) yr, A3 = 0.0020(0.0002) d, and e = 0.38(0.16). Analysis shows such oscillation would be caused by a magnetic activity which can be explained by the Applegate mechanism. Meanwhile, according to the value of l3 and the estimated physical parameters of V344 Lac, the mass of the third companion may be 0.79 Ms. This third body could be a wide company.

As a characteristic feature of generic blazars the polarization behavior of the quasi-periodic optical outbursts observed in OJ287 is investigated. The optical light-curves of the December/2015 outburst are also simulated in terms of the precessing jet nozzle model previously proposed. The polarization behavior of three primary quasi-periodic optical outbursts peaking in 1983.0, 2007.8 and 2015.8 are analyzed in order to understand the nature of their optical radiation. A two-component model has been applied,showing that the variations in flux density, polarization degree and polarization position angle can be consistently interpreted with two polarized components: one steady-component with constant polarization and one burst-component with varying polarization (e.g., relativistic shock propagating along the jet-beam axis). The flux light curves of the December/2015 outburst (including its first flare and second flare) are well model-simulated in terms of 14 elementary synchrotron sub-flares, each having a symmetric profile. The model simulations of polarization behavior for the three major outbursts (in 1983.0, 2007.8 and 2015.8) demonstrate that they all exhibit rapid and large rotations in polarization position angle, implying that they are synchrotron flares produced in the jet. Combining with the results previously obtained for interpreting the optical light curves in terms of lighthouse effect for both quasi-periodic and non-periodic outbursts, we suggest that relativistic jet models may be the most appropriate models for understanding the nature of the optical flaring radiation in blazar OJ287: its optical outbursts may comprise a number of blended "elementary synchrotron flares" , each produced by the helical motion of individual superluminal optical knots via lighhouse effect.

ASAS J174406+2446.8 was originally found as a $\delta$ Scuti-type pulsating star with the period P=0.189068 $days$ by ASAS survey. However, the LAMOST stellar parameters reveal that it is far beyond the red edge of pulsational instability strip on the $\log g-T$ diagram of $\delta$ Scuti pulsating stars. To understand the physical properties of the variable star, we observed it by the 1.0-m Cassegrain reflecting telescope at Yunnan Observatories. Multi-color light curves in B, V, R$_{c}$ and I$_{c}$ bands were obtained and are analyzed by using the W-D program. It is found that this variable star is a shallow-contact binary with an EB-type light curve and an orbital period of 0.3781\,days rather than a $\delta$ Scuti star. It is a W-subtype contact binary with a mass ratio of $1.135(\pm0.019)$ and a fill-out factor of $10.4(\pm5.6)\,\%$. The situation of ASAS J174406+2446.8 resembles those of other EB-type marginal-contact binaries such as UU Lyn, II Per and GW Tau. All of them are at a key evolutionary phase from a semi-detached configuration to a contact system predicted by the thermal relaxation oscillation theory. The linear ephemeris was corrected by using 303 new determined times of light minimum. It is detected that the O - C curve shows a sinusoidal variationthat could be explained by the light-travel-time effect via the presence of a cool red dwarf. The present investigation reveals that some of the $\delta$ Scuti-type stars beyond the red edge of pulsating instability strip on the $\log g-T$ diagram are misclassified eclipsing binaries. To understand their structures and evolutionary states, more studies are required in the future.

V752 Cen is a triple-lined spectroscopic contact binary. Its multi-color light curves were obtained in the years 1971 and 2018, independently. Photometric analyses reveal that the two sets of light curves produce almost consistent results. It contains a W-subtype totally eclipsing binary, and its mass ratio and fill-out factor are $q = 3.35(1)$ and $f = 29(2)\,\%$. The absolute elements of its two component stars were determined to be $M_{1} = 0.39(2)M_\odot$, $M_{2} = 1.31(7)M_\odot$, $R_{1} = 0.77(1)R_\odot$, $R_{2} = 1.30(2)R_\odot$, $L_{1} = 0.75(3)L_\odot$ and $L_{2} = 2.00(7)L_\odot$. The period of V752 Cen is 0.37023198 day. The 0.37-d period remained constant from its first measurement in 1971 until the year 2000. However, it changed suddenly around the year 2000 and has been increasing continuously at a rate of $dP/dt=+5.05\times{10^{-7}}day\cdot year^{-1}$ since then, which can be explained by mass transfer from the less massive component star to the more massive one with a rate of $\frac{dM_{2}}{dt}=2.52\times{10^{-7}}M_\odot/year$. The period variation of V752 Cen over the 48 years in which the period has been monitored is really unusual, and is potentially related to effects from the possible presence of a nearby third star or of a pair of stars in a second binary.

Blazar OJ287 is a unique source in which optical outbursts with double-peak structure have been observed quasi-periodically with a cycle of 12yr. It may be one of the best candidates for searching supermassive black hole binaries. We investigate the connection between its optical and radio/millimeter variations and interpret the emissions in terms of relativistic jet models. Specifically, we make a detailed analysis and model simulation of the optical and radio/mm light curves for the outburst during the period of 1995.8--1996.1. It is shown that the multi-wavelength light curves at optical V-band and radio/mm wavelengths (37, 22, 14.5 and 8 GHz) can be decomposed into 36 individual elementary flares, each of which has a symmetric profile. The elementary flares can be understood to be produced through lighthouse effect due to the helical motion of corresponding superluminal optical/radio knots. Helical motion of superluminal knots should be prevailing in the inner regions of its relativistic jet formed in the magnetosphere of the putative supermassive black hole/accretion disk system. A comprehensive and compatible framework for understanding the entire phenomena in OJ287 is described.

YZ Phe is a very short-period contact binary (Sp.= $K2\,V$) with an orbital period of 0.2347 days near the short period limit (0.22 d). We present the complete light curves in $VRI$ bands, which photometric data were obtained with the 0.61-m telescope of PROMPT-8 at CTIO in Chile during June to October 2016 and August 2017. The photometric solutions were determined by using the W-D method and the results reveal that YZ Phe is a W-subtype shallow contact binary ($f\sim$ 10%, $q$ = 2.635 or $1/q$ = 0.379 for W subtype) with rotational motion of a large hot spot on the more massive component, showing a strong O'Connell effect with variation of maxima in photometric time series at period of 4.20 yr and stellar cycle at period of 1.28 yr. By compiling all available eclipse times, the result shows a long-term period decrease at a rate of $\mathrm{d}P/\mathrm{d}t = -2.64(\pm 0.02)\times 10^{-8}$ d $yr^{-1}$, superimposed on a cyclic variation ($A_3$ = 0.0081 days and $P_3$ = 40.76 years). This variation cannot be explained by Applegate mechanism. Thus, the cyclic change may be interpreted as light-travel time effect via the presence of a cool third body. Based on photometric solutions, the third light was detected with 2% contribution of total light in $V$ and $I$ bands. Those support the existence of a third body. For the long-term period decrease, it can be explained by mass transfer from the more massive component ($M_2 \sim 0.74 M_{\odot}$) to the less massive one ($M_1 \sim 0.28 M_{\odot}$) or plus AML via magnetic braking. With $1/q$ $<$ 0.4 and long-term period decrease, all suggest that YZ Phe is on the AML-controlled state and its fill-out factor will increase, as well as the system will evolve into a deeper normal contact binary.

Six mid-eclipse times of the eclipsing binary XZ And are obtained, which are analysed together with others collected from the literature. Two sets of cyclic variations with periods of 33.43 and 100.4 yr are found if a double-Keplerian model is used to fit the data. The 1:3 ratio of the periods suggests that both cyclic variations arise from dynamic motions of two companions rather than magnetic activity of the eclipsing pair. According to the double-Keplerian model, the companions have the masses of about 1.32 Msun and 1.33 Msun, respectively. Comparing the total masses of the eclipsing pair of 3.12 Msun, it is obvious that XZ And is a general N-body system. The strong gravitational perturbation between two companions invalidates the double-Keplerian model. It is strange that two Keplerian periods with a 1:3 ratio are derived from the best fits with the inappropriate model. The illogical, but interesting phenomena also appear in other two Algol systems, suggesting that our discoveries deserve attention from astronomers.

The light curves of optical outbursts observed in blazar OJ287 during 1983-2015 are analyzed and model-simulated to investigate the nature of its optical radiation. It is shown that the December/2015 outburst has its multi-wavelength variability behavior very similar to that of the synchrotron outburst in March/2016, indicating that the 2015 outburst may originate from synchrotron process. In combination with helical motion of superluminal components, the precessing jet nozzle scenario previously proposed is used to model-simulate the lightcurves of all the optical outbursts discussed. The optical light curves for both periodic and non-periodic outbursts observed in blazar OJ287 can be well interpreted in terms of lighthouse effect due to the helical motion of superluminal optical knots, showing their common origin in synchrotron process. A coherent and compatible framework is tentatively suggested to understand the entire phenomena in OJ287. The double-peak structure of the periodic outbursts might be explained by invoking the cavity-accretion flare models for comparable-mass binary systems in eccentric motion.

TY Pup is a well-known bright eclipsing binary in southern hemisphere with an orbital period of 0.8192 days. New light curves in $B, V, (RI)_C$ bands were obtained with the 0.61-m reflector robotic telescope (PROMPT-8) at CTIO in Chile from January to February 2015 and from March to April 2017. By analyzing those photometric data with the W-D method, it is found that TY Pup is a low-mass-ratio ($q \sim$ 0.184) and deep contact binary with a high fill-out factor ($84.3\,\%$). An investigation of all available times of minimum light including three new ones obtained with the 60-cm and the 1.0-m telescopes at Yunnan Observatories in China indicates that the period change of TY Pup is complex. An upward parabolic variation in the $O-C$ diagram is detected to be superimposed on a cyclic oscillation. The upward parabolic change reveals a long-term continuous increase in the orbital period at a rate of $\mathrm{d}P/\mathrm{d}t = 5.57(\pm 0.08)\times10^{-8}$ d $\textrm{yr}^{-1}$. The period increase can be explained by mass transfer from the less massive component ($M_2 \sim 0.3 M_{\odot}$) to the more massive one ($M_1 \sim 1.65 M_{\odot}$) and the mass ratio of the binary system will become more smaller. In this way, the binary will be merging when it meets the criterion that the orbital angular momentum is less than 3 times the total spin angular momentum, i.e., $J_{orb} < 3J_{rot}$. This suggests that the system will finally merge into a rapid-rotating single star and may produce a luminous red nova. The cyclic oscillation in the $O-C$ diagram can be interpreted by the light-travel time effect (LITE) via the presence of an additional companion.

New CCD photometric light curves of short period (P=0.285d) eclipsing binary RW Dor are presented. The observations performed with the PROMPT-8 robotic telescope at CTIO in Chile from March 2015 to March 2017. The other eclipse timings were obtained from the 2.15-m JS telescope at CASLEO, San Juan, Argentina in December 2011. By light-curve analysis, it is found that RW Dor is a W-type shallow contact binary with a fill-out factor $f \sim 11\%$ and high mass ratio $q \sim 1.587$ (1/q = 0.63), where the hotter component is the less massive one ($M_1 \sim 0.52M_{\odot}$ and $M_2 \sim 0.82M_{\odot}$). For orbital period investigation, the new fifteen eclipse times and those in previous published were compiled. Based on $O-C$ analysis with very weak evidence suggests that a long-term period decrease with a rate of $\mathrm{d}P/\mathrm{d}t = -9.61\times10^{-9}$ d $\textrm{yr}^{-1}$ is superimposed on a cyclic variation ($A_3$ = 0.0054 days and $P_3$ = 49.9 yrs). The long-term period decrease can be interpreted as mass transfer from the more massive component to the less massive one or combine with the angular momentum loss (AML) via magnetic braking. In addition, with the marginal contact phase, high mass ratio (1/q $>$ 0.4) and the long-term period decrease, all suggest that RW Dor is a newly formed contact binary via a Case A mass transfer and it will evolve into a deeper normal contact binary. If the cyclic change is correct, the light-travel time effect via the presence of a cool third body will be more plausible to explain for this.

GU Ori was observed with the 1m telescope at Yunnan Observatories in 2005. To determine its physical properties, the Wilson-Devinney program is used. The results reveal that GU Ori is a W-subtype shallow contact binary with a more massive but cooler star 2. The mass of its two component stars are estimated to be $M_1 = 0.45M_\odot$, $M_2 = 1.05M_\odot$. The O'Connell effect was reported to be negative on the light curves observed in 2005. However, it changed to a positive one on the light curves observed from 2011 to 2012. The mean surface temperatures of star 2 ($T_2$) determined by the two sets of light curves were different, which may result from stellar activities. The O - C diagram shows that the period of GU Ori is decreasing at a rate of $dP/dt=-6.24\times{10^{-8}}day\cdot year^{-1}$, which may be caused by mass transfer from star 2 to star 1 with a rate of $\frac{dM_{2}}{dt}= - 2.98\times{10^{-8}}M_\odot/year$. GU Ori is a contact binary with quite high metallicity.

We try to propose a tentative relativistic jet model for explaining the entire radio-optical phenomena observed in blazar OJ287, which has been observed quasi-periodically with a cycle of about 12yr in its optical light curve. We investigate the currently available theoretical and observational studies on the phenomena observed in OJ287 and try to find clues to its double-jet structure. It is found that the kinematic features of its superluminal components observed at 43GHz and 15GHz could be well interpreted in terms of a precessing double-jet nozzle model. And the light curves of a few optical double-peaked outbursts could be interpreted in terms of relativistic jet models. Both jets precess with the same period of 12yr, equal to the optical period and the precession of the jets could be originated from the orbital motion of the binary. Thus the masses of the binary black holes could be estimated. We also tentatively suggest a comprehensive framework for understanding the entire phenomena in OJ287 in terms of relativistic jet models.

We present a timing analysis of the eclipsing post-common envelope binary (PCEB) DE CVn. Based on new CCD photometric observations and the published data, we found that the orbital period in DE CVn has a cyclic period oscillation with an amplitude of $28.08$ s and a period of $11.22$ years plus a rapid period decrease at a rate of $\dot{P}=-3.35\times10^{-11}ss^{-1}$. According to the evolutionary theory, secular period decreases in PCEBs arise from angular momentum losses (AMLs) driven by gravitational radiation (GR) and magnetic braking (MB). However, the observed orbital decay is too fast to be produced by AMLs via GR and MB, indicating that there could be other AML mechanism. We suggest that a circumbinary disk around DE CVn may be responsible for the additional AML. The disk mass was derived as a few$\times$$10^{-4}$-$10^{-3}$$M_{\odot}$ , which is in agreement with that inferred from previous studies in the order of magnitude. The cyclic change is most likely result of the gravitational perturbation by a circumbinary object due to the Applegate's mechanism fails to explain such a large period oscillation. The mass of the potential third body is calculated as $M_{3}\sin{i'}=0.011(\pm0.003)M_{\odot}$. Supposing the circumbinary companion and the eclipsing binary is coplanar, its mass would correspond to a giant planet. This hypothetical giant planet is moving in a circular orbit of radius $\sim5.75(\pm2.02)$ AU around its host star.

Four-color charge-coupled device (CCD) light curves in $B$, $V$, $Rc$ and $Ic$ bands of the total-eclipsing binary system, V1853 Ori, are presented. By comparing our light curves with those published by previous investigators, it is detected that the O'Connell effect on the light curves is disappeared. By analyzing those multi-color light curves with the Wilson-Devinney code (W-D code), it is discovered that V1853 Ori is an A-type intermediate-contact binary with a degree of contact factor of $f=33.3\%(3.7\%)$ and a mass ratio of $q=0.1896(0.0013)$. Combining our 10 new determined times of light minima together with the others published in the literature, the period changes of the system is investigated. We found that the general trend of the observed-calculated $(O-C)$ curve shows a downward parabolic variation that corresponds to a long-term decrease in the orbital period with a rate of $dP/dt=-1.96(0.46)\times{10^{-7}}$ d yr$^{-1}$. The long-term period decrease could be explained by mass transfer from the more-massive component to the less-massive one. By combining our photometric solutions with the Gaia DR 2 data, absolute parameters were derived as $M_{1}$ = 1.20 M$_{\odot}$, $M_{2}$ = 0.23 M$_{\odot}$, $R_{1}$ = 1.36 R$_{\odot}$, and $R_{2}$ = 0.66 R$_{\odot}$. The long-term period decrease and the intermediate-contact configuration suggest that V1853 Ori will evolve into a high fill-out overcontact binary.

New photometric observations of NSVS 01286630 were performed and two sets of four-color (B, V , Rc, Ic) light curves (LCs) were obtained. Using the 2013 version of the Wilson-Devinney (W-D) code, we analyzed these data. The photometric solutions reveal that NSVS 01286630 is an active detached eclipsing binary (EB) with a high orbit inclination. Remarkably, the temperature of the primary component (the hotter star) is higher than the secondary one, but, the value of mass ratio q(M2/M1) for NSVS 0128663 is more than 1, which can be explained that the surface of the secondary component of NSVS 01286630 is covered with big cool-star-spots. Based on our new CCD mid-eclipse times and the data published until now, the variations of the mid-eclipse times were reanalyzed in detail using a weighted least-squares method. It is discovered that the (O-C) diagram of the system shows a cyclic oscillation with a period of 3.61 yr and an amplitude of 0.001 days. The cyclic variation may be caused by the light travel time effect (LTTE) due to the presence of a third companion, whose mass we calculated is about 0.11 solar mass. The third body may affect the orbital evolution of the central binary system by transferring the angular momentum.

As an eclipsing polar with 3.39 hrs orbital period, MN Hya was going through state change when we observed it during 2009-2016. 10 new mid-eclipse times, along with others obtained from literature, allow us to give a new ephemeris. The residuals of linear fit show that period decreased during the phase of state change. It means angular momentum was lost during this phase. The X-ray observation indicates the mass accretion rate as about $3.6\times10^{-9}M_{\odot}yr^{-1}$. The period decrease gives that at least 60 percent of mass being transfered from secondary was lost, maybe in form of the spherically symmetric stellar wind. In high state, the data shows the intensity of the flickering reduced when system had higher accretion rate, and that flickering sticks out with primary timescale about 2 minutes, which implies the position of the threading point as about 30 radius of the white dwarf above the surface of it. The trend of light curves of the system in high state follows that of low state for a large fraction of phase interval from phase 0 to phase 0.4 since which the cyclotron feature is visible, and the primary intensity hump of light curves near phase 0.7 when the system is in high state did not appear on the curve when it is in low state. Those facts contradict the predictions of the two-pole model.

About 766 $\delta$ Scuti stars were observed by LAMOST by June 16, 2017. Stellar atmospheric parameters of 525 variables were determined, while spectral types were obtained for all samples. In the paper, those spectroscopic data are catalogued. We detect a group of 131 unusual and cool variable stars (UCVs) that are distinguished from the normal $\delta$ Scuti stars (NDSTs). On the H-R diagram and the $\log g-T$ diagram, the UCVs are far beyond the red edge of pulsational instability trip. Their effective temperatures are lower than 6700\,K with periods in the range from 0.09 to 0.22\,days. NDSTs have metallicity close to that of the Sun as expected, while UCVs are slightly metal poor than NDSTs. The two peaks on the distributions of the period and stellar atmospheric parameters are all caused by the existence of UCVs. When those UCVs are excluded, it is discovered that the effective temperature, the gravitational acceleration and the metallicity all are correlated with the pulsating period for NDSTs and their physical properties and evolutionary states are discussed. The detecting of those UCVs indicates that about 25\% of the known $\delta$ Scuti stars may be misclassified. However, if some of them are confirmed to be pulsating stars, they will be a new-type pulsator and their investigations will shed light on theoretical instability domains and on the theories of interacting between the pulsation and the convection of solar-type stars. Meanwhile, 88 $\delta$ Scuti stars are detected to be the candidates of binary or multiple systems.

About 3196 EA-type binaries (EAs) were observed by LAMOST by June 16, 2017 and their spectral types were derived. Meanwhile stellar atmospheric parameters of 2020 EAs were determined. In the paper, those EAs are catalogued and their physical properties and evolutionary states are investigated. The period distribution of EAs suggests that the period limit of tidal locking for the close binaries is about 6 days. It is found that the metallicity of EAs is higher than that of EWs indicating that EAs are generally younger than EWs and they are the progenitors of EWs. The metallicities of long-period EWs (0.4 < P < 1 days) are the same as those of EAs with the same periods, while their values of Log (g) are usually smaller than those of EAs. These support the evolutionary process that EAs evolve into long-period EWs through the combination of angular momentum loss (AML) via magnetic braking and case A mass transfer. For short-period EWs, their metallicities are lower than those of EAs, while their gravitational accelerations are higher. These reveal that they may be formed from cool short-period EAs through AML via magnetic braking with little mass transfer. For some EWs with high metallicities, they may be contaminated by material from the evolution of unseen neutron stars and black holes or they have third bodies that may help them to form rapidly through a short timescale of pre-contact evolution. The present investigation suggests that the modern EW populations may be formed through the combination of aforementioned mechanisms.

In this paper, CCD photometric light curves for the short-period eclipsing binary 1SWASP J140533.33+114639.1 (hereafter J1405) in the $BVR$ bands are presented and analyzed using the 2013 version of the Wilson-Devinney (W-D) code. It is discovered that the J1405 is a W-subtype shallow contact binary with a contact degree of $f$ = 7.9$\pm$ 0.5\% and a mass ratio of $q$ = 1.55 $\pm$ 0.02. In order to explain the asymmetric light curves of the system, a cool star-spot on the more massive component was employed. This shallow contact eclipsing binary may be formed from a short-period detached system through the orbital shrinkage due to angular momentum loss. Based on $(O-C)$ method, the variation of the orbital period was studied using all the available times of the minimum light. The $(O-C)$ diagram reveals that the period is increasing continuously at a rate of $dP/dt=+2.09\times{10^{-7}}$days yr$^{-1}$, which can be explained by mass transfer from the less massive component to the more massive one.

W UMa-type contact binaries belong to close binary systems whose components exactly overflow their Roche lobes and share a common convective envelope (CCE). In the last twenty years, the long-term variations of their orbital periods have been thought to depend on several mechanisms. Now, we suggest a new mechanism: CCE-dominated mechanism. The CCE-dominated mechanism is found based on our numerical result, especially at high mass ratios, that the orbital periods (P) of contact binaries change very much with their fill-out factors (f). Because f is taken as a measurement of the thickness of CCE, the physical cause for the variation of P is a mass transfer between CCE and components. Further, an f-dominated simplification model for this mechanism is introduced. According to it, P may change in a long-term oscillation way with a similar time scale of the thermal modulation, meanwhile q is decreasing slowly till the two components merge. It could be also applied to explain the presence of extremely short period, high mass ratio and deep contact binaries. Moreover, the CCE-dominated mechanism should always work due to mass transfer and mass loss both occurring via CCE. Therefor, the effect of CCE on the variations of orbital periods may have been underestimated before.

New Multi-color $B$ $V$ $R_c$ $I_c$ photometric observation are presented for W UMa type eclipsing binary V1073 Cyg. The multi-color light curves analysis with the Wilson-Devinney(W-D) procedure acquired the absolute parameters of this system, showing that V1073Cyg is a shallow contact binary system with fill-out factor $f=0.124(\pm0.011)$. We collected all available times of light minima spanning 119 years including CCD data to construct the O-C curve and made detailed O-C analysis. The O-C diagram shows that the period change is complex. There exist a long-term continuous decrease and a cyclic variation. The period is decreasing at a rate of $\dot P=-1.04(\pm0.18)\times 10^{-10}$ $days\cdot{cycle}^{-1}$, and with the period decrease, V1073 Cyg will evolve to deep contact stage. The cyclic variation with a period of $P_3=82.7(\pm3.6) years$ and an amplitude of $A=0.028(\pm0.002) day$ may be explained by the magnetic activity of one or both components or the light travel time effect (LTTE) caused by a distant third companion with $M_3({i'}=90^{\circ})=0.511M_\odot$.

Aims. We aim to find more eclipsing multiple systems and obtain their parameters, thus increasing our understanding of multiple systems. Methods. The extraneous eclipses on the \textit{kepler} binary light curves indicating extraneous bodies were searched. The binary light curves were analyzed using the binary model, and the extraneous eclipses were studied on their periodicity and shape changes. Results. Three binaries with extraneous eclipses on the binary light curves were found and studied based on the \textit{Kepler} observations. The object KIC 5255552 is an eclipsing triple system with a fast changing inner binary and an outer companion uncovered by three groups of extraneous eclipses of $862.1(\pm0.1)$ d period. The KIC 10091110 is suggested to be a double eclipsing binary system with several possible extraordinary coincidences: the two binaries share similar extremely small mass ratios ($0.060(13)$ and $0.0564(18)$), similar mean primary densities ($0.3264(42)\;\rho_\odot$ and $0.3019(28)\;\rho_\odot$), and, most notably, the ratio of the two binaries' periods is very close to integer 2 (8.5303353/4.2185174 = 2.022). The KIC 11495766 is a probable triple system with a $\sim120.73$ d period binary and (at least) one non-eclipse companion. Furthermore, very close to it in the celestial sphere, there is a blended background stellar binary of 8.3404432 d period. A first list of 25 eclipsing multiple candidates is presented, with the hope that it will be beneficial for study of eclipsing multiples.

Two twin binaries, KIC 4826439 and KIC 6045264, with very similar component stars were found photometrically based on $\textit{Kepler}$ eclipsing binary light curves. The absolute parameters of the massive components are 1.156(0.03)$M_\odot$, 1.881(0.02)$R_\odot$, 6065K for KIC 4826439, and 0.874(0.3)$M_\odot$, 1.206(0.02)$R_\odot$, 6169(30)K for KIC 6045264. The differences between the components are less than two percents for all the parameters. A very low proportion of the twin binaries ($2/1592\approx0.13\%$) was found, which does not support the previous findings of the excesses of twins on binary mass ratio distribution, but support a deficiently low proportion of twins. A new method is practiced to work out the absolute parameters of the two twins without the radial velocities. This method requires the solution of the light curves, the spectra and the evolutionary isochrones of covering the complete stellar parameter space, simultaneously. We also studied their evolution tracks that: KIC 4826439 will experience an unstable mass transfer stage followed by an unclear ending, and KIC 6045264 will become a single star via an over-contact phase. It seems highly unlikely that the two twin binaries will produce twin degenerate binaries, although they have quite similar components.

1SWASP\,J162117.36$+$441254.2 was originally classified as an EW-type binary with a period of 0.20785\,days. However, it was detected to have undergone a stellar outburst on June 3, 2016. Although the system was latter classified as a cataclysmic variable (CV) and the event was attributed as a dwarf-nova outburst, the physical reason is still unknown. This binary has been monitored photometrically since April 19, 2016 and many light curves were obtained before, during and after the outburst. Those light and color curves observed before the outburst indicate that the system is a special CV. The white dwarf is not accreting material from the secondary and there are no accretion disks surrounding the white dwarf. By comparing the light curves obtained from April 19 to September 14, 2016, it was found that magnetic activity of the secondary is associated with the outburst. We show strong evidence that the $L_1$ region on the secondary was heavily spotted before and after the outburst and thus quench the mass transfer, while the outburst is produced by a sudden mass accretion of the white dwarf. These results suggest that J162117 is a good astrophysical laboratory to study stellar magnetic activity and its influences on CV mass transfer and mass accretion.

Numerous EWs were discovered by several deep photometric survey and there are about 40785 EW-type binary systems listed in the international variable star index (VSX) by March 13, 2017. 7938 of them were observed by LAMOST by November 30, 2016 and their spectral types were given. Stellar atmospheric parameters of 5363 EW-type binary stars were determined based on good spectroscopic observations. In the paper, those EWs were catalogued and their properties are analyzed. The distributions of the orbital period (P), the effect temperature (T), the gravitational acceleration (Log(g)), the metallicity ([Fe/H]) and the radial velocity (RV) are presented for those observed EW-type systems. It is shown that about 80.6\% sample stars have metallicity below zero indicating that EW-type systems are old stellar population. This is in agreement with the the conclusion that the EW binaries are formed from moderately close binaries through angular momentum loss via magnetic braking that takes a few hundred million to a few billion years. The unusual high metallicities of a few percent of EWs may be caused by contaminating of material from the evolution of unseen neutron stars and black holes in the systems. The correlations between the orbital period and the effect temperature, the gravitational acceleration and the metallicity are presented and their scatters are mainly caused by (i) the presence of the third bodies and (ii) the wrong determined periods sometimes. It is shown that some EW contain evolved component stars and the physical properties of EWs are mainly depending on their orbital periods. It is found that the extremely short-period EWs may be older than their long-period cousins because they have lower metallicities. This reveals that they have a longer timescale of pre-contact evolution and their formation and evolution are mainly driven by angular momentum loss via magnetic braking.

Multi-color photometric light curves of NO Cam in $V$ $R_C$ and $I_C$ bands are obtained and analyzed simultaneously using the Wilson-Devinney (W-D) program. The solutions suggest that NO Cam is an A-subtype overcontact binary with its mass ratio to be $q = 0.439$ and a contact degree of $f = 55.5\,\%$. The small temperature difference ($\Delta T = 44K$) between its two components indicates that the system is under thermal contact. The high orbital inclination ($i = 84.5^{\circ}$) strengthens our confidence in the parameters determined from the light curves. All available times of minimum light are collected and period variations are analyzed for the first time. The $O - C$ curve reveals that its period is increasing continuously at a rate of $dP/dt=+1.46\times{10^{-9}}$, which can be explained by mass transfer from the less massive component to the more massive one. After the upward parabolic variation is subtracted, the residuals suggest that there may be a cyclic variation with a period of 2.23 years and an amplitude of $A_3 = 0.00153$ days, which may due to the light-travel time effect (LTTE) arising from the gravitational influence of a close-in tertiary component. The close-in companion reveals that early dynamic interaction among triple system may have played very important role in the formation of the W UMa type binaries.

Multi-color ($B$ $V$ $R_c$ $I_c$) CCD photometric light curves of the contact binary V502 Oph are analyzed using the Wilson-Devinney (W-D) program. The solutions reveal that V502 Oph is a W-subtype contact ($f = 35.3\,\%$) binary system. The temperature difference between its two components is $240K$ and the more massive star has a lower surface temperature. A cool spot is added in our model to account for the light curves' asymmetry (O'Connell effect) and third light is detected for the first time in the light curves' modeling. Combining the orbital inclination ($i = 76.4^{\circ}$) with the published mass function of V502 Oph, the absolute physical parameters of the two components are determined, which are $M_1= 0.46(\pm0.02)M_\odot, M_2=1.37(\pm0.02)M_\odot, R_1=0.94(\pm0.01)R_\odot, R_2=1.51(\pm0.01)R_\odot, L_1=1.13(\pm0.02)L_\odot$ and $L_2=2.49(\pm0.03)L_\odot$. The formation and evolutionary status of V502 Oph are discussed. All photoelectric and CCD times of light minimum about V502 Oph are gathered and its orbital period variations are analyzed. The results show that the orbital period of V502 Oph is decreasing continuously at a rate of $dP/dt=-1.69\times{10^{-7}}day\cdot year^{-1}$, which corresponds to a conservative mass transfer rate of $\frac{dM_{2}}{dt}=-3.01\times{10^{-8}}M_\odot/year$. The light-travel time effect (LTTE) is due to the presence of a close-in tertiary component with a period of $P_3=18.7$ years and an amplitute of 0.00402days. V502 Oph is an ideal target to test the formation and evolution theories of binary and multiple systems in which the light curves, $O - C$ curve and spectroscopic observations are comprehensively researched.

New eclipse timings of the Z Cam-type dwarf nova AY Psc were measured and the orbital ephemeris was revised. Based on the long-term AAVSO data, moreover, the outburst behaviors were also explored. Our analysis suggests that the normal outbursts are quasi-periodic, with an amplitude of $\sim2.5(\pm0.1)$ mag and a period of $\sim18.3(\pm0.7)$ days. The amplitude vs. recurrence-time relation of AY Psc is discussed, and we concluded that this relation may represents general properties of dwarf nova (DN) outbursts. The observed standstill ends with an outburst, which is inconsistent with the general picture of Z Cam-type stars. This unusual behavior was considered to be related to the mass-transfer outbursts. Moreover, the average luminosity is brighter during standstills than during outburst cycles. The changes in brightness marks the variations in $\dot{M}_{2}$ due to the disc of AY Psc is nearly steady state. $\dot{M}_{2}$ value was limited to the range from $6.35\times10^{-9}$ to $1.18\times10^{-8}M_{\odot}yr^{-1}$. More detailed examination shows that there are a few small outbursts presence during standstills. These events with amplitudes of $\sim0.5-0.9$ mag are very similar to the stunted outbursts reported in some NLs. We discussed several possible mechanisms and suggested that the most reasonable mechanism for these stunted outbursts is a changing mass-transfer rate.

GSC 4560-02157 is a new eclipsing cataclysmic variable with an orbital period of $0.265359$ days. By using the published $V-$ and $R-$band data together with our observations, we discovered that the $O-C$ curve of GSC 4560-02157 may shows a cyclic variation with the period of $3.51$ years and an amplitude of $1.40$ min. If this variation is caused by a light travel-time effect via the existence of a third body, its mass can be derived as $M_{3}sini'\approx91.08M_{Jup}$, it should be a low-mass star. In addition, several physical parameters were measured. The colour of the secondary star was determined as $V-R=0.77(\pm0.03)$ which corresponds to a spectral type of K2-3. The secondary star's mass was estimated as $M_{2}=0.73(\pm0.02)M_{\odot}$ by combing the derived $V-R$ value around phase 0 with the assumption that it obeys the mass-luminosity relation of the main sequence stars. This mass is consistent with the mass$-$period relation of CV donor stars. For the white dwarf, the eclipse durations and contacts of the white dwarf yield an upper limit of the white dwarf's radius corresponding to a lower limit mass of $M_{1}\approx0.501M_{\odot}$. The overestimated radius and previously published spectral data indicate that the boundary layer may has a very high temperature.

EX Dra is a long-period eclipsing dwarf nova with $\sim2-3$ mag amplitude outbursts. This star has been monitored photometrically from November, 2009 to March, 2016 and 29 new mid-eclipse times were obtained. By using new data together with the published data, the best fit to the $O-C$ curve indicate that the orbital period of EX Dra have an upward parabolic change while undergoing double-cyclic variations with the periods of 21.4 and 3.99 years, respectively. The upward parabolic change reveals a long-term increase at a rate of $\dot{P}={+7.46}\times10^{-11}{s} {s^{-1}}$. The evolutionary theory of cataclysmic variables (CVs) predicts that, as a CV evolves, the orbital period should be decreasing rather than increasing. Secular increase can be explained as the mass transfer between the secondary and primary or may be just an observed part of a longer cyclic change. Most plausible explanation for the double-cyclic variations is a pair of light travel-time effect via the presence of two companions. Their masses are determined to be $M_{A}sini'_{A}=29.3(\pm0.6) M_{Jup}$ and $M_{B}sini'_{B}=50.8(\pm0.2) M_{Jup}$. When the two companions are coplanar to the orbital plane of the central eclipsing pair, their masses would match to brown dwarfs.

DV UMa is an eclipsing dwarf nova with an orbital period of $\sim2.06$ h, which lies just at the bottom edge of the period gap. To detect its orbital period changes we present 12 new mid-eclipse times by using our CCD photometric data and archival data. Combining with the published mid-eclipse times in quiescence, spanning $\sim30$ yr, the latest version of the $O-C$ diagram was obtained and analyzed. The best fit to those available eclipse timings shows that the orbital period of DV UMa is undergoing a cyclic oscillation with a period of $17.58(\pm0.52)$ yr and an amplitude of $71.1(\pm6.7)$ s. The periodic variation most likely arises from the light-travel-time effect via the presence of a circumbinary object because the required energy to drive the Applegate mechanism is too high in this system. The mass of the unseen companion was derived as $M_{3}\sin{i'}=0.025(\pm0.004)M_{\odot}$. If the third body is in the orbital plane (i.e. $i'=i=82.9^{\circ}$) of the eclipsing pair, it would match to a brown dwarf. This hypothetical brown dwarf is orbiting its host star at a separation of $\sim8.6$ AU in an eccentric orbit ($e=0.44$).

We present the photometric results of the eclipsing cataclysmic variable (CV) WZ Sge near the period minimum ($P_{min}$). Eight new mid-eclipse times were determined and the orbital ephemeris was updated. Our result shows that the orbital period of WZ Sge is decreasing at a rate of $\dot{P}=-2.72(\pm0.23)\times{10^{-13}}\,s s^{-1}$. This secular decrease, coupled with previous detection of its donor, suggest that WZ Sge is a pre-bounce system. Further analysis indicates that the observed period decrease rate is about $1.53$ times higher than pure gravitational radiation (GR) driving. We constructed the evolutionary track of WZ Sge, which predicts that $P_{min}$ of WZ Sge is $\sim77.98 (\pm0.90)$ min. If the orbital period decreases at the current rate, WZ Sge will evolve past its $P_{min}$ after $\sim25.3$ Myr. Based on the period evolution equation we find $\dot{M}_{2}\simeq4.04(\pm0.10)\times10^{-11}M_{\odot}yr^{-1}$, which is compatible with the current concept of CV evolution at ultrashort orbital periods.

Four candidates of eclipsing multiples, based on new extraneous eclipses found on Kepler binary light curves, are presented and studied. KIC 7622486 is a double eclipsing binary candidate with orbital period of 2.2799960 days and 40.246503 days. The two binary systems do not eclipse each other in the line of sight, but there is mutual gravitational influence between them which leads to the small but definite eccentricity 0.0035(0.0022) on the short 2.2799960 days period orbit. KIC 7668648 is hierarchical quadruple system candidate, with two sets of solid 203(+-5) days period extraneous eclipses and another independent set of extraneous eclipses. A clear and credible extraneous eclipse is found on the binary light curve of KIC 7670485 which made it a triple system candidates. Two sets of extraneous eclipse of about 390 days and 220 days period are found on KIC 8938628 binary curves, which not only confirms the previous conclusion of $388.5(+-0.3) triple system, but also proposed a new additional objects that make KIC 8938628 a hierarchical quadruple system candidates. It is wished that the four candidates here will make some contributions to the field of eclipsing multiples.

The $B$ $V$ $R_c$ $I_c$ bands light curves of the newly discovered binary system \astrobj{GSC 03122-02426} are obtained and analyzed using the Wilson-Devinney (W-D) code. The solutions suggest that the mass ratio of the binary system is $q = 2.70$ and the less massive component is $422K$ hotter than the more massive one. We conclude that \astrobj{GSC 03122-02426} is a W-subtype shallow contact (with a contact degree of $f = 15.3\,\%$) binary system. It may be a newly formed contact binary system which is just under geometrical contact and will evolve to be a thermal contact binary system. The high orbital inclination ($i = 81.6^{\circ}$) implies that \astrobj{GSC 03122-02426} is a total eclipsing binary system and the photometric parameters obtained by us are quite reliable. We also estimate the absolute physical parameters of the two components in \astrobj{GSC 03122-02426}, which will provide fundamental information for the research of contact binary systems. The formation and evolutionary scenario of \astrobj{GSC 03122-02426} is discussed.

The four-color ($B$ $V$ $R_c$ $I_c$) light curves of V776 Cas are presented and analyzed using the Wilson-Devinney (W-D) method. It is discovered that V776 Cas is an early F-type (F2V) overcontact binary with a very high contact degree ($ f=64.6\,\%$) and an extremely low mass ratio ($q=0.130$), which indicate that it is at the final evolutionary stage of cool short-period binaries. The mass of the primary and secondary stars are calculated to be $M_1 = 1.55(\pm0.04)M_\odot$, $M_2 = 0.20(\pm0.01)M_\odot$. V776 Cas is supposed to be formed from an initially detached binary system via the loss of angular momentum due to the magnetic wind. The initial mass of the present primary and secondary components are calculated to be $M_{1i} = 0.86(\pm0.10)M_\odot$ and $M_{2i} = 2.13(\pm0.04)M_\odot$. The observed-calculated ($O$-$C$) curve exhibits a cyclic period variation, which is due to the light-travel time effect (LTTE) caused by the presence of a third component with a period of 23.7 years. The mass of the third component is estimated to be $M_3 = 1.04(\pm0.03)M_\odot$ and the orbital inclination of the third component is calculated to be $i' = 33.1^{\circ}$. The distance of the binary system to the mass center of the triple system is calculated to be $a'_{12} = 3.45AU$. The presence of the close-in tertiary component may play an important role in the formation and evolution of this binary system by drawing angular momentum from the central system.

II UMa is a late F-type (F5) contact binary with a close-in tertiary and a distant visual companion. According to the four-color ($B$ $V$ $R_c$ $I_c$) light curves' solutions of II UMa, it is a high fill-out (f=$86.6\,\%$) and low mass ratio ($q = 0.172$) contact binary system, which indicate that it is at the late evolutionary stage of late-type tidal-locked binary stars. The mass of the primary star and secondary one are calculated to be $M_1 = 1.99M_\odot$, $M_2 = 0.34M_\odot$. The primary star has evolved from ZAMS, but it is still before TAMS, and the secondary star is even more evolved. Considering the mass ratio ($M_3/M_1 = 0.67$) obtained by spectroscopic observations, the mass of the close-in tertiary is estimated to be $M_3 = 1.34M_\odot$. The period variations of the binary system is investigated for the first time. According to the observed-calculated ($O$-$C$) curve analysis, a continuous period increase at a rate of $dP/dt=4.88\times{10^{-7}}day\cdot year^{-1}$ is determined. It may be just a part of a cyclic period change, or the combinational period change of a parabolic variation and a cyclic one. More times of minimum light are needed to confirm this. The presence of the tertiary component may play an important role in the formation and evolution of this binary system by drawing angular momentum from the central system during the pre-contact stage.

Two sets of light curves in $V$ $R_c$ $I_c$ bands for a newly discovered binary system UCAC4 436-062932 are obtained and analyzed using the Wilson-Devinney (W-D) code. The two sets of light curves get almost consistent results. The determined mass ratio is about $q = 2.7$ and the less massive component is nearly $250K$ hotter than the more massive one. The solutions conclude that UCAC4 436-062932 is a W-subtype shallow contact (with a contact degree of $f = 20\,\%$) binary system. Since the O'Connell effect appears on one set of the light curves, theories proposed to explain the effect are discussed. We assume that spot model may be the more plausible one to the O'Connell effect appeared on the asymmetric light curves of the binary system UCAC4 436-062932. Therefore, we add a cool spot on the surface of the more massive star (component with lower effective temperature) and get a quite approving results for the light curve fitting. It will provide evidence to support the spot model in the explanatory mechanism of O'Connell effect.

The light curve of V921 Her in ultraviolet band observed by the Lunar-based Ultraviolet Telescope (LUT) is analyzed by the Wilson-Devinney code. Our solutions conclude that V921 Her is an early type marginal contact binary system with an additional close-in component. The binary system is under poor thermal contact with a temperature difference of nearly $700K$ between the two components. The close-in component contributes about $19\,\%$ of the total luminosity in the triple system. Combining the radial velocity study together with our photometric solutions, the mass of the primary star and secondary one are calculated to be $M_1 = 1.784(\pm0.055)M_\odot$, $M_2 = 0.403(\pm0.012)M_\odot$. The evolutionary scenario of V921 Her is discussed. All times of light minimum of V921 Her available in the bibliography are taken into account and the $O - C$ curve is analyzed for the first time. The most probable fitting results are discussed in the paper, which also confirm the existence of a third component ($P_3=10.2$ year) around the binary system. The period of V921 Her is also undergoing a continuously rapid increase at a rate of $dP/dt=+2.79\times{10^{-7}}day\cdot year^{-1}$, which may due to mass transfer from the less massive component to the more massive one.

The first V Rc Ic bands light curves of MQ UMa are presented and analyzed using the Wilson-Devinney (W-D) program. It is discovered that MQ UMa is an A-subtype contact binary with a high fill-out (f=82 %) and a low mass ratio (q = 0:195), which indicates that it is in the late evolutionary stage of latetype tidal-locked binary stars. The mass of the primary and secondary stars are estimated and the evolutionary status of the two components are placed in the H-R diagram. The W-D solutions also indicates that MQ UMa may be a triple system with an F5V type additional component. A sample of sixteen high fillout, low mass ratio overcontact binaries are collected and their possible evolution scenarios are discussed. Using the five times of light minimum recently observed together with those collected from literatures, the authors find that the observedcalculate (O-C) curve exhibits a cyclic period variation. The cyclic period change also reveals the presence of a tertiary component, which may play an important role in the formation and evolution of this binary system by drawing angular momentum from the central system.

We report the detection of orbital modulation, a model solution, and X-ray properties of a newly discovered contact binary, 2MASS J11201034$-$2201340. We serendipitously found this X-ray point source outside the error ellipse when searching for possible X-ray counterparts of $\gamma$-ray millisecond pulsars among the unidentified objects detected by the {\it Fermi Gamma-ray Space Telescope}. The optical counterpart of the X-ray source (unrelated to the $\gamma$-ray source) was then identified using archival databases. The long-term CRTS survey detected a precise signal with a period of $P=0.28876208(56)$ days. A follow-up observation made by the SLT telescope of Lulin Observatory revealed the binary nature of the object. Utilizing archived photometric data of multi-band surveys, we construct the spectral energy distribution, which is well fitted by a K2V spectral template. The fitting result of the orbital profile using the Wilson--Devinney code suggests that 2MASS J11201034-2201340 is a short-period A-type contact binary and the more massive component has a cool spot. The X-ray emission was first noted in observations made by \swift\, then further confirmed and characterized by an \xmm\ observation. The X-ray spectrum can be described by a power law or thermal Bremsstrahlung. Unfortunately, we could not observe significant X-ray orbital modulation. Finally, according to the spectral energy distribution, this system is estimated to be 690 pc from Earth with a calculated X-ray intensity of $(0.7-1.5)\times 10^{30}$ erg s$^{-1}$, which is in the expected range of an X-ray emitting contact binary.

All available mid-eclipse times of the eclipsing binary Z Draconis are analysed, and three sets of cyclic variations with periods of 20.1, 29.96 and 59.88 yr are found. The low-amplitude variations with the period of 20.1 yr may be attributed to the unavoidable and slight imperfection in the double-Keplerian model, which gives the periods of 29.96 and 59.88 yr. Interestingly, the Z Draconis system is close to a 2:1 mean-motion resonance, or 6:3:2 mean-motion resonance if the period of 20.1 yr is true. We also find that the best solutions tend to give the minimum eccentricities. Based on Kepler's third law, the outermost companion has the minimum mass of $\sim0.77M_{\bigodot}$, whereas the middle companion is an M dwarf star with a mass of $\sim0.40M_{\bigodot}$, suggesting that Z Draconis is a general N-body system.

We present a new few-parameter phenomenological model of light curves of eclipsing binaries and stars with transiting planets that is able to fit the observed light curves with the accuracy better than 1\% of their amplitudes. The model can be used namely for appropriate descriptions of light curve shapes, classification, mid-eclipse time determination, and fine period analyses.

New multi-color light curves of the very short period K-type eclipsing binary V1799 Ori were obtained and analyzed with the W-D code. The photometric solutions reveal that the system is a W-type shallow-contact binary with a mass ratio of $q=1.335(\pm0.005)$ and a degree of contact about $f = 3.5(\pm1.1)\%$. In general, the results are in good agreement with which is reported by Samec. The remarkable O'Connell effects in the light curves are well explained by employing star spots on the binary surface, which confirms that the system is active at present. Several new times of light minimum were obtained. All the available times of light minimum were collected, along with the recalculated and new obtained. Applying a least-squares method to the constructed O-C diagram, a new ephemeris was derived for V1799 Ori. The orbital period is found to show a continuous weak increase at a rate of $1.8(\pm0.6)\times10^{-8}$ days$\cdot$yr$^{-1}$. The extreme-shallow contact, together with the period increase, suggests that the binary may be at a critical stage predicted by the TRO theory. \keywords{binaries : close -- binaries : eclipsing -- stars: individual (V1799 Ori)

NGC 7789 is an intermediate-age open cluster with an age that is similar to the mean age of contact binary stars. V12 is a bright UMa-type binary star with an orbital period of 0.397 day. The first complete light curves of V12 in V, R, and I bands are presented and analyzed with the W-D (Wilson and Devinney) method. The results shows that V12 is an intermediate-contact binary (f = 43%) with a mass ratio of 3.848, and it is a W-type contact binary where the less massive component is slightly hotter than the more massive one. The asymmetries of the light curves are explained by the presence of a dark spot on the more massive component. The derived orbital inclination (i = 83.6o) indicates that it is a total eclipsing binary, which suggests that the determined parameters are reliable. The orbital period may show a long-term increase at a rate of P = +2.48 x 10-6 days/year that reveals a rapid mass transfer from the less massive component to the more massive one. However, more observations are needed to confirm this conclusion. The presence of an intermediate-contact binary in an intermediate-age open cluster may suggest that some contact binaries have a very short pre-contact timescale. The presence of a third body or/and stellar collision may help to shorten the pre-contact evolution.

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.

This work reports new photometric results of eclipsing cataclysmic variable V617 Sagittarii (V617 Sgr). We analyzed the orbital period change of V617 Sgr, by employing three new CCD eclipse timings since 2010 along with all the available data from the literature. It was found that the orbital period of V617 Sgr undergoes an obvious long-term increase, which confirms the result revealed by Steiner et al. (2006). The rate of orbital period increase was calculated to be ${\dot{P}}$ = +2.14(0.05) $\times$ 10$^{-7}$ day/year. This suggests the lifetime of the secondary star will attain to the end in a timescale of 0.97 $\times$ 10$^6$ years faster than that predicted previously. In particular, a cyclic variation with a period of 4.5 year and an amplitude of 2.3 minutes may present in the O-C diagram. Dominated by the wind-accretion mechanism, high mass transfer from the low mass secondary to the white dwarf is expected to sustain in the V Sge-type star V617 Sgr during its long-term evolution. The mass transfer rate $|\dot{M}_{tr}|$ was estimated to be in the range of about 2.2 $\times$ 10$^{-7}$ to 5.2 $\times$ 10$^{-7}$ M$_{\odot}$ yr$^{-1}$. Accordingly, the already massive ($\geq$ 1.2 M$_{\odot}$) white dwarf primary will process stable nuclear burning, accrete a fraction of mass from its companion to reach the standard Chandrasekhar mass limit ($\simeq$ 1.38 M$_{\odot}$), and ultimately produce a type Ia supernova (SN Ia) within about 4 $\sim$ 8 $\times$ 10$^{5}$ years or earlier.

Twenty-one new optical light curves, including five curves obtained in 2009 and sixteen curves detected from the AAVSO International Database spanning from 1977 to 2011, demonstrate 16 new primary minimum light times in the high state. Furthermore, seven newly found low-state transient events from 2006 to 2009 were discovered, consisting of five Gaussian-shaped events and two events with an exponential form with decay timescales of $<$0.005 days; these timescales are one order of magnitude shorter than those of previous X-ray flare events. In the state transition, two special events were detected: a "disrupted event" with an amplitude of $\sim$ 2 mag and a duration of $\sim$ 72 minutes and continuing R-band twin events larger than all known R-band flares detected in M-type red dwarfs. All 45 available high-state data points spanning over 35 yr were used to construct an updated O-C diagram of \objectname{AM Herculis}, which clearly shows a significant sine-like variation with a period of 12-15 yr and an amplitude of 6-9 minutes. Using the inspected physical parameters of the donor star, the secular variation in the O-C diagram cannot be interpreted by any decided angular momentum loss mechanism, but can satisfy the condition $\tau_{\dot{\rm M}_{2}}\simeq\tau_{\rm KH}>>\tau_{\dot{\rm R}_{\rm 2}}$, which is required by numerical calculations of the secular evolution of cataclysmic variables. In order to explain the prominent periodic modulation, three plausible mechanisms - spot motion, the light travel-time effect, and magnetic active cycles - are discussed in detail.

We analyzed the B and V light curves of two EA-type binaries V211 and NV358 using the WD code for the first time. Our analysis shows that V211 is a typical Algol-type binary and NV358 is a well detached binary system. As the two binaries are definite proper motion members of omega Centauri, we estimated their physical parameters, obtaining M_1=1.13\pm0.03M_\odot, R_1=0.98\pm0.01R_\odot and M_2=0.33\pm0.01M_\odot, R_2=0.92\pm0.01R_\odot for V211; M_1=1.30\pm0.05M_\odot, R_1=1.03\pm0.01R_\odot and M_2=0.58\pm0.02M_\odot, R_2=0.78\pm0.01R_\odot for NV358. On the color-magnitude diagram of omega Centauri, V211 is located in the faint blue straggler region and its primary component is more massive than a star at the main-sequence turnoff. Therefore, V211 is a blue straggler and should be formed by mass transfer from the secondary component to the primary. The age of NV358 is less than 1.93 Gyr, indicating that it is much younger than first-generation stars in omega Centauri. As NV364 in omega Centauri, NV358 should be a second-generation binary.

Models of period variations are basic tools for period analyzes of variable stars. We introduce phase function and instant period and formulate basic relations and equations among them. Some simple period models are also presented.

Plenty of mid-eclipse timings of short-periodic eclipsing binaries derived from series of visual observations appear to be an acceptable source of information for orbital period analyses, namely if they were done in time-intervals not covered by other types of observations. However, our thorough period analysis of the nearly contact eclipsing binary BS Vulpeculae proves that visually determined timings done in 1979--2003 were undoubtedly biased to accommodate the existing linear ephemeris. The heavily subjective character of visual observations disqualifies them as a source of true phase information apt for fine eclipsing binary period analyses. Consequently we warn against the use of visual timings without a preceding careful verification.

By using six new determined mid-eclipse times together with those collected from the literature, we found that the Observed-Calculated (O-C) curve of RR Cae shows a cyclic change with a period of 11.9 years and an amplitude of 14.3s, while it undergoes an upward parabolic variation (revealing a long-term period increase at a rate of dP/dt =+4.18(+-0.20)x10^(-12). The cyclic change was analyzed for the light-travel time effect that arises from the gravitational influence of a third companion. The mass of the third body was determined to be M_3*sin i' = 4.2(+-0.4) M_{Jup} suggesting that it is a circumbinary giant planet when its orbital inclination is larger than 17.6 degree. The orbital separation of the circumbinary planet from the central eclipsing binary is about 5.3(+-0.6)AU. The period increase is opposite to the changes caused by angular momentum loss via magnetic braking or/and gravitational radiation, nor can it be explained by the mass transfer between both components because of its detached configuration. These indicate that the observed upward parabolic change is only a part of a long-period (longer than 26.3 years) cyclic variation, which may reveal the presence of another giant circumbinary planet in a wide orbit.

We report here the tentative discovery of a Jovian planet in orbit around the rapidly pulsating subdwarf B-type (sdB-type) eclipsing binary NY Vir. By using new determined eclipse times together with those collected from the literature, we detect that the observed-calculated (O-C) curve of NY Vir shows a small-amplitude cyclic variation with a period of 7.9\,years and a semiamplitude of 6.1\,s, while it undergoes a downward parabolic change (revealing a period decrease at a rate of $\dot{P}=-9.2\times{10^{-12}}$). The periodic variation was analyzed for the light-travel time effect via the presence of a third body. The mass of the tertiary companion was determined to be $M_3\sin{i^{\prime}}=2.3(\pm0.3)$\,$M_{Jupiter}$ when a total mass of 0.60\,$M_{\odot}$ for NY Vir is adopted. This suggests that it is most probably a giant circumbinary planet orbiting NY Vir at a distance of about 3.3 astronomical units (AU). Since the rate of period decrease can not be explained by true angular momentum loss caused by gravitational radiation or/and magnetic braking, the observed downward parabolic change in the O-C diagram may be only a part of a long-period (longer than 15 years) cyclic variation, which may reveal the presence of another Jovian planet ($\sim2.5$$M_{Jupiter}$) in the system.

The new discovered short-period close binary star, XY LMi, was monitored photometrically since 2006. It is shown that the light curves are typical EW-type and show complete eclipses with an eclipse duration of about 80 minutes. By analyzing the complete B, V, R, and I light curves with the 2003 version of the W-D code, photometric solutions were determined. It is discovered that XY LMi is a high fill-out, extreme mass ratio overcontact binary system with a mass ratio of q=0.148 and a fill-out factor of f=74.1%, suggesting that it is on the late evolutionary stage of late-type tidal-locked binary stars. As observed in other overcontact binary stars, evidence for the presence of two dark spots on both components are given. Based on our 19 epoches of eclipse times, it is found that the orbital period of the overcontact binary is decreasing continuously at a rate of dP/dt=-1.67\times10^{-7}\,days/year, which may be caused by the mass transfer from the primary to the secondary or/and angular momentum loss via magnetic stellar wind. The decrease of the orbital period may result in the increase of the fill-out, and finally, it will evolve into a single rapid-rotation star when the fluid surface reaching the outer critical Roche Lobe.

Using the precise times of mid-egress of the eclipsing polar HU Aqr, we discovered that this polar is orbited by two or more giant planets. The two planets detected so far have masses of at least 5.9 and 4.5\,M_{Jup}. Their respective distances from the polar are 3.6 AU and 5.4 AU with periods of 6.54 and 11.96 years, respectively. The observed rate of period decrease derived from the downward parabolic change in O-C curve is a factor 15 larger than the value expected for gravitational radiation. This indicates that it may be only a part of a long-period cyclic variation, revealing the presence of one more planet. It is interesting to note that the two detected circumbinary planets follow the Titus-Bode law of solar planets with n=5 and 6. We estimate that another 10 years of observations will reveal the presence of the predicted third planet.

Four newest CCD eclipse timings of the white dwarf for polar UZ Fornacis and Six updated CCD mid-eclipse times for SW Sex type nova-like V348 Puppis are obtained. The detailed O-C analyses for both CVs inside period gap are made. Orbital period increases at a rate of $2.63(\pm0.58)\times10^{-11} s\;s^{-1}$ for UZ Fornacis and of $5.8(\pm1.9)\times10^{-12} s\;s^{-1}$ for V348 Puppis, respectively, are discovered in their new O-C diagrams. However, the conservative mass transfer from the secondary to the massive white dwarf cannot explain the observed orbital period increases for both CVs, which are regarded as part of modulations at longer periods. Moreover, the O-C diagram of UZ Fornacis shows a possible cyclical change with a period of $\sim23.4(\pm5.1)yr$. For explaining the observed cyclical period changes in UZ Fornacis, both mechanisms of magnetic activity cycles in the late-type secondary and the light travel-time effect are regarded as two probable causes. Not only does the modulation period 23.4yr obey the empirical correlation derived by \cite{lan99}, but also the estimated fractional period change $\Delta P/P\sim7.3\times10^{-7}$ displays a behavior similar to that of the CVs below the period gap. On the other hand, a calculation for the light travel-time effect implies that the tertiary component in UZ Fornacis may be a brown dwarf with a high confidence level, when the orbital inclination of the third body is larger than $16^{\circ}$.

Combining with our newest CCD times of light minimum of EM Cygni, all 45 available times of light minimum including 7 data with large scatters are compiled and the updated O-C analysis is made. The bestfit for the O-C diagram of EM Cygni is a quadratic-plus-sinusoidal fit. The secular orbital period decrease rate -2.5(\pm 0.3)x10^{-11} s s^{-1} means that magnetic braking effect with a rate of mass loss via stellar wind, 2.3x10^{-10}Msunyr^{-1}, is needed for explaining the observed orbital period decrease. Moreover, for explaining the significant cyclical period change with a period of \sim 17.74(\pm 0.01)yr shown in the O-C diagram, magnetic activity cycles and light travel-time effect are discussed in detail. The O-C diagram of EM Cygni cannot totally rule the possibility of multi-periodic modulation out due to the gaps presented after 25000 cycles. Based on the hypothesis of a K-type third star in literature, light trave-time effect may be a more plausible explanation. However, the low orbital inclination of the third body (\sim 7.4 degree) suggests that the hypothetic K-type third star may be captured by EM Cygni. But assuming the spectral contamination from a block of circumbinary material instead of a K-type third star, the third star may be a brown dwarf in case of the coplanar orbit with parent binary.

We searched the orbital period changes in 182 EA-type (including the 101 Algol systems used by \cite{hal89}), 43 EB-type and 53 EW-type binaries with known both the mass ratio and the spectral type of their secondary components. We reproduced and improved the same diagram as Hall's (1989) according to the new collected data. Our plots do not support the conclusion derived by \cite{hal89} that all cases of cyclical period changes are restricted to binaries having the secondary component with spectral types later than F5. The presence of period changes also among stars with secondary component of early type indicates that the magnetic activity is one cause, but not the only one, for the period variation. It is discovered that cyclic period changes, likely due to the presence of a third body are more frequent in EW-type binaries among close binaries. Therefore, the most plausible explanation of the cyclical period changes is the LTTE via the presence of a third body. By using the century-long historical record of the times of light minimum, we analyzed the cyclical period change in the Algol binary WW Dra. It is found that the orbital period of the binary shows a $\sim112.2 \textbf{\textrm{yr}}$ cyclic variation with an amplitude of $\sim0.1977\textbf{\textrm{days}}$. The cyclic oscillation can be attributed to the LTTE via a third body with a mass no less than $6.43 M_{\odot}$. However, no spectral lines of the third body were discovered indicating that it may be a candidate black hole. The third body is orbiting the binary at a distance shorter than 14.4 AU and it may play an important role in the evolution of this system.

The two CCD photometries of the intermediate polar TV Columbae are made for obtaining the two updated eclipse timings with high precision. There is an interval time \sim 17yr since the last mid-eclipse time observed in 1991. Thus, the new mid-eclipse times can offer an opportunity to check the previous orbital ephemerides. A calculation indicates that the orbital ephemeris derived by Augusteijn et al. (1994) should be corrected. Based on the proper linear ephemeris (Hellier, 1993), the new orbital period analysis suggests a cyclical period variation in the O-C diagram of TV Columbae. Using Applegate's mechanism to explain the periodic oscillation in O-C diagram, the required energy is larger than that a M0-type star can afford over a complete variation period \sim 31.0(\pm 3.0)yr. Thus, the light travel-time effect indicates that the tertiary component in TV Columbae may be a dwarf with a low mass, which is near the mass lower limit \sim 0.08Msun as long as the inclination of the third body high enough.

HS 0705+6700 is a short-period (P=2.3 hours), close binary containing a hot sdB-type primary and a fully convective secondary. We have monitored this eclipsing binary for more than 2 years and as a result, 32 times of light minimum were obtained. Based on our new eclipse times together with these compiled from the literature, it is discovered that the O-C curve of HS 0705+6700 shows a cyclic variation with a period of 7.15 years and a semiamplitude of 92.4 s. The periodic change was analyzed for the light-travel time effect that may be due to the presence of a tertiary companion. The mass of the third body is determined to be M3 sin i = 0.0377 (+/-0.0043) Msun when a total mass of 0.617 Msun for HS 0705+6700 is adopted. For orbital inclinations i >= 32.8, the mass of the tertiary component would be below the stable hydrogen-burning limit of M3~0.072 Msun, and thus it would be a brown dwarf. The third body is orbiting the sdB-type binary at a distance shorter than 3.6 astronomical units (AU). HS 0705+6700 was formed through the evolution of a common envelope after the primary becomes a red giant. The detection of a sub-stellar companion in HS 0705+6700 system at this distance from the binary could give some constraints on stellar evolution in such systems and the interactions between red giants and their companions.

Up to now, most stellar-mass black holes were discovered in X-ray emitting binaries, in which the black holes are formed through a common-envelope evolu tion. Here we give evidence for the presence of a massive black hole candidate as a tertiary companion in the massive eclipsing binary V Puppis. We found that the orbital period of this short-period binary (P=1.45 days) shows a periodic variation while it undergoes a long-term increase. The cyclic period oscillation can be interpreted by the light-travel time effect via the presence of a third body with a mass no less than 10.4 solar mass. However, no spectral lines of the third body were discovered indicating that it is a massive black hole candidate. The black hole candidate may correspond to the weak X-ray source close to V Puppis discovered by Uhuru, Copernicus, and ROSAT satellites produced by accreting materials from the massive binary via stellar wind. The circumstellar matter with many heavy elements around this binary may be formed by the supernova explosion of the progenitor of the massive black hole. All of the observations suggest that a massive black hole is orbiting the massive close binary V Puppis with a period of 5.47 years. Meanwhile, we found the central close binary is undergoing slow mass transfer from the secondary to the primary star on a nuclear time scale of the secondary component, revealing that the system has passed through a rapid mass-transfer stage.

During the period 1966.5 - 2006.2 the 15GHz and 8GHz lightcurves of 3C454.3 (z=0.859) show a qsasi-periodicity of ~12.8 yr (~6.9 yr in the rest frame of the source) with a double-bump structure. This periodic behaviour is interpreted in terms of a rotating double-jet model in which the two jets are created from the black holes in a binary system and rotate with the period of the orbital motion. The periodic variations in the radio fluxes of 3C454.3 are suggested to be mainly due to the lighthouse effects (or the variation in Doppler boosting) of the precessing jets which are caused by the orbital motion. In addition, variations in the mass-flow rates accreting onto the black holes may be also involved.

It is generally thought that conservative mass transfer in Algol binaries causes their orbits to be wider, in which the less massive star overflows its Roche-lobe. The observed decrease in the orbital periods of some Algol binaries suggests orbital angular momentum loss during the binary evolution, and the magnetic braking mechanism is often invoked to explain the observed orbital shrinkage. Here we suggest an alternative explanation, assuming that a small fraction of the transferred mass forms a circumbinary disk, which extracts orbital angular momentum from the binary through tidal torques. We also perform numerical calculations of the evolution of Algol binaries with typical initial masses and orbital periods. The results indicate that, for reasonable input parameters, the circumbinary disk can significantly influence the orbital evolution, and cause the orbit to shrink on a sufficiently long timescale. Rapid mass transfer in Algol binaries with low mass ratios can also be accounted for in this scenario.

There is new evidence which identifies seasonal changes of the variability time scale in intraday variable sources with refractive interstellar scintillation effects. Such a RISS-model takes the annual change of the Earth velocity relative to the scattering medium into account. In September 1998 we found a remarkable prolongation of the variability time scale in the IDV-source 0917+624 with only small variations in flux density during a period of 5 days. This was explained as a seasonal effect, in which the velocity vector of the Earth and the interstellar medium nearly canceled. In order to further investigate the applicability of the model for 0917+624, we performed an Effelsberg 6cm-flux monitoring program over the course of one year. Since September 2000, the source appears to be remarkably inactive and yet (May 2001), no return to its normal, faster and stronger variability pattern is observed. Here, our observational results and a possible explanation for the current quiescence are presented.

The quasar 0917+624 has been one of the best studied intraday variable (IDV) radio sources. However, debate continues as to whether the underlying cause is intrinsic or extrinsic. Much previous work has assumed the IDV to be intrinsic and which implies an extraordinarily compact source for the radio emission; in contrast, an extrinsic variation due to interstellar scintillation (ISS) implies a relatively larger source diameter, though at the smaller end of the range expected for relativistic jet models. Kraus et al. (1999) reported a marked slowing of the IDV at 6cm wavelength in September 1998, and suggested a change in the source was responsible. However, here we show that the slowing is consistent with the annual modulation in scintillation time-scale expected for ISS, under the assumption that the scattering medium moves with the local standard of rest (LSR). The ISS time scale is governed by the ISS spatial scale divided by the Earth's velocity relative to the scattering plasma. It happens that in the direction of 0917+624 the transverse velocity of the Earth with respect to the LSR varies widely with a deep minimum in the months of September to November. Hence the slowing of the IDV in September 1998 strongly suggests that ISS, rather than intrinsic variation of the source is the dominant cause of the IDV.

We present radio observations at three frequencies and contemporaneous optical monitoring of the peculiar BL Lac object AO 0235+164. During a three-week campaign with the VLA we observed intraday variability in this source and found a distinct peak which can be identified throughout the radio frequencies and tentatively connected to the R-band variations. This event is characterized by unusual properties: its strength increases, and its duration decreases with wavelength, and it peaks earlier at 20cm than at 3.6 and 6cm. We discuss several generic models (a `standard' shock-in-jet model, a precessing beam, free-free-absorption in a foreground screen, interstellar scattering, and gravitational microlensing), and explore whether they can account for our observations. Most attempts at explaining the data on 0235+164 require an extremely small source size, which can be reconciled with the $10^{12}$ K inverse Compton limit only when the Doppler factor of the bulk flow is of order 100. However, none of the models is completely satisfactory, and we suggest that the observed variability is due to a superposition of intrinsic and propagation effects.

We present a simple model for the apparent superluminal motion along curved trajectories of features observed in the compact radio jet of the quasar 3C345. The model is inspired by the magneto-hydrodynamic approach of Camenzind (1986), and assumes a conical geometry of the jet within about three milliarcseconds of the radio core. Simultaneous conservation of three out of four physical parameters is used to derive analytic solutions for the equation of motion of the idealized jet features. These are the jet's kinetic energy, angular momentum, and momentum along the jet axis. Conservation of angular momentum is required to fit the observed trajectories at distances larger than 3 milliarcsec, in agreement with the kinematic predictions of Camenzind's model. The best-fit model variant preserves angular momentum, kinetic energy, and opening angle, with an inclination of the jet axis to the line of sight 6.8 degrees, an opening angle 0.95 degrees, and Lorentz factors 5.8 and 4.6 for components C4 and C5, respectively (H_o=100 km/(s Mpc)). The model ignores turbulent small-scale motion of the jet components, and it does not take into account the overall bending of the jet beyond about 4 milliarcseconds, restricting it's applicability to the inner few milliarcseconds from the core.