We present analysis of medium resolution optical spectra and long term V band photometry of four cool stars, BD+13 5000, BD+11 3024, TYC 3557-919-1 and TYC 5163-1764-1. Our spectroscopic analysis reveals that the stars are giant or sub-giant from K0 or K1 spectral type, and all of them exhibit emission features in their Ca ii H& K lines. These features appear to be modulated with the rotation of the stars. Except BD+11 3024, we observe that the radial velocities of the target stars are not stable, which suggests that each of them might be a member of a binary system. Global analysis of photometric data indicates clear cyclic variation for BD+13 5000 and TYC 5163-1764-1 with a period of 8.0$\pm$0.3 and 5.04$\pm$0.04 year, respectively. Besides that, we observe a dramatic increase ($\sim$0\fm7) in the mean brightness of BD+11 3024, accompanied with a 2.87$\pm$0.12 cyclic variation, embedded into the global brightening trend, which indicates possible multiple cycles on this star.

We present spectroscopic and photometric analysis of KIC\,9451096, where the latter is based on very high precision long cadence photometry obtained by $Kepler$ space craft. Combined spectroscopic and photometric modeling show that the system is a detached eclipsing binary in a circular orbit and composed of F5V + K2V components. Subtracting the best--fit light curve model from whole long cadence data reveals additional low (mmag) amplitude light variation in time and occasional flares, suggesting low, but still remarkable level of magnetic spot activity on the K2V component. Analyzing rotational modulation of light curve residuals enables us to estimate differential rotation coefficient of the K2V component as $k = 0.069 \pm 0.008$, which is 3 times weaker compared with the solar value of $k = 0.19$, assuming a solar type differential rotation. We find stellar flare activity frequency for K2V component as 0.000368411 $h^{-1}$ indicating low magnetic activity level.

We present spectroscopic and photometric analysis of the eclipsing binary KIC\,7385478. We find that the system is formed by F1V + K4III-IV components. Combining results from analysis of spectroscopic data and $Kepler$ photometry, we calculate masses and radii of the primary and the secondary components as M$_{1}$ = 1.71 $\pm$ 0.08 \Msun, M$_{2}$ = 0.37 $\pm$ 0.04 \Msun~ and R$_{1}$ = 1.59 $\pm$ 0.03\Rsun, R$_{2}$ = 1.90 $\pm$ 0.03\Rsun, respectively. Position of the primary component in HR diagram is in the region of $\gamma$ Doradus type pulsators and residuals from light curve modeling exhibit additional light variation with a dominant period of $\sim$0.5 day. These are clear evidences of the $\gamma$ Doradus type pulsations on the primary component. We also observe occasional increase in amplitude of the residuals, where the orbital period becomes the most dominant period. These may be attributed to the cool star activity originating from the secondary component.

HD 208472 is among the most active RS~CVn binaries with cool starspots. Decade-long photometry has shown that the spots seem to change their longitudinal appearance with a period of about six years, coherent with brightness variations. Our aim is to spatially resolve the stellar surface of HD 208472 and relate the photometric results to the true longitudinal and latitudinal spot appearance. Furthermore, we investigate the surface differential rotation pattern of the star. We employed three years of high-resolution spectroscopic data with a high signal-to-noise ratio (S/N) from the STELLA robotic observatory and determined new and more precise stellar physical parameters. Precalculated synthetic spectra were fit to each of these spectra, and we provide new spot-corrected orbital elements. A sample of 34 absorption lines per spectrum was used to calculate mean line profiles with a S/N of several hundred. A total of 13 temperature Doppler images were reconstructed from these line profiles with the inversion code iMap. Differential rotation was investigated by cross-correlating successive Doppler images in each observing season. Spots on HD 208472 are distributed preferably at high latitudes and less frequently around mid-to-low latitudes. No polar-cap like structure is seen at any epoch. We observed a flip-flop event between 2009 and 2010, manifested as a flip of the spot activity from phase 0.0 to phase 0.5, while the overall brightness of the star continued to increase and reached an all-time maximum in 2014. Cross-correlation of successive Doppler images suggests a solar-like differential rotation that is ~15 times weaker than that of the Sun.

In this study, we present BVR photometry of V369 Gem in 2006 - 2008. We combined our V data with the data taken from the literature. We analysed all the data together and calculated the new light elements. General features of the light curves were discussed. The light curves were modelled with the Fourier method and the Fourier coefficients calculated. The cos(2theta) term was found to be dominant. The pre-whitened V light curves demonstrated that stellar spot activity also exist. There are generally two minima in the pre-whitened light curves. Thus, there must be two spotted areas, which are separated about 180 deg from each other. As conclusions, we have remarked for the first time that V369 Gem is an ellipsoidal variable with the high level chromospheric activity.

We present multi-colour CCD photometry of GSC4589-2999 obtained in 2008 and 2009. The observations indicate that the system is an active Algol binary. Based on the new data, the mean brightness of the system is decreasing through the years 2007-2009. The light curves obtained in 2008-2009 are modelled using the Wilson-Devinney code. We also discussed the light and colour variations of the system at different orbital phases. Evidence suggests that these brightness and colour variations are due to the rotation of unevenly distributed starspots on two components of the system.

In this paper, we present standard Johnson UBV photometry of the eclipsing binary BD+36 3317 which is known as a member of Delta Lyrae (Stephenson 1) cluster. We determined colors and brightness of the system, calculated E(B - V) color excess. We discovered that the system shows total eclipse in secondary minimum. Using this advantage, we found that the primary component of the system has B9 - A0 spectral type. Although there is no published orbital solution, we tried to estimate the physical properties of the system from simultaneous analysis of UBV light curves with 2003 version of Wilson-Devinney code. Then we considered photometric solution results together with evolutionary models and estimated the masses of the components as M1 = 2.5 Msun and M2 = 1.6 Msun. Those estimations gave the distance of the system as 353 pc. Considering the uncertainties in distance estimation, resulting distance is in agreement with the distance of Delta Lyrae cluster.

We present BVR observations of DK CVn from 2007 and 2008. We analysed the BVR light curves of the system and obtained the system's parameters. Using the `q-search' method, we measured the mass ratio of the system (q) as 0.55. Taking the temperature of the primary component as 4040K, the temperature of the secondary was found to be 3123K. Several flares were detected, and the distributions of flare equivalent duration versus flare total duration were modelled using the One-Phase Exponential Association Function for these flares. The parameters of the model demonstrated that the flares are the same as those detected from UV Ceti stars. We also demonstrate that the variation at out-of-eclipse must be caused by some cool spot(s) on one of the components. The star is found to show two active longitudes in which the spots are mainly formed. Consequently, this study reveals that DK CVn should be a chromospherically active binary star.

Taking into account results obtained from light-curve analysis and out-of-eclipse analyses, we discuss the nature of GSC 02038-00293 and also its magnetic activity behaviour. We obtained light curves of the system during observing seasons 2007, 2008 and 2011. We obtained its secondary minimum clearly in I-band observations in 2008 for the first time. Analysing this light curve, we found the physical parameters of the components. The light-curve analysis indicates that the possible mass ratio of the system is 0.35. We obtained the remaining V-band light curves, extracting the eclipses. We modelled these remaining curves using the spotmodel program and found possible spot configurations of the magnetically active component for each observing season. The models demonstrated that there are two active longitudes for the active component. The models reveal that both active longitudes migrate in the direction of decreasing longitude. We also examined the light curves in out-of-eclipse phases with respect to minimum and maximum brightness, amplitude, etc. The amplitude of the curves during out-of-eclipse phases varies in a sinusoidal way with a period of ~8.9yr the mean brightness of the system is dramatically decreasing. The phases of the deeper minimum during out-of-eclipse periods exhibit a migration toward decreasing phase.