We aim to investigate variations in the arrival time of coherent stellar pulsations due to the light-travel time effect to test for the presence of sub-stellar companions. Those companions are the key to one possible formation scenario of apparently single sub-dwarf B stars. We made use of an extensive set of ground-based observations of the four large amplitude p-mode pulsators DW Lyn, V1636 Ori, QQ Vir, and V541 Hya. Observations of the TESS space telescope are available on two of the targets. The timing method compares the phase of sinusoidal fits to the full multi-epoch light curves with phases from the fit of a number of subsets of the original time series. Observations of the TESS mission do not sample the pulsations well enough to be useful due to the (currently) fixed two-minute cadence. From the ground-based observations, we infer evolutionary parameters from the arrival times. The residual signals show many statistically significant periodic signals, but no clear evidence for changes in arrival time induced by sub-stellar companions. The signals can be explained partly by mode beating effects. We derive upper limits on companion masses set by the observational campaign.

We present rotation period measurements for 107 M dwarfs in the mass range $0.15-0.70 M_\odot$ observed within the context of the APACHE photometric survey. We measure rotation periods in the range 0.5-190 days, with the distribution peaking at $\sim$ 30 days. We revise the stellar masses and radii for our sample of rotators by exploiting the Gaia DR2 data. For $\sim 20\%$ of the sample, we compare the photometric rotation periods with those derived from different spectroscopic indicators, finding good correspondence in most cases. We compare our rotation periods distribution to the one obtained by the Kepler survey in the same mass range, and to that derived by the MEarth survey for stars in the mass range $0.07-0.25 M_\odot$. The APACHE and Kepler periods distributions are in good agreement, confirming the reliability of our results, while the APACHE distribution is consistent with the MEarth result only for the older/slow rotators, and in the overlapping mass range of the two surveys. Combining the APACHE/Kepler distribution with the MEarth distribution, we highlight that the rotation period increases with decreasing stellar mass, in agreement with previous work. Our findings also suggest that the spin-down time scale, from fast to slow rotators, changes crossing the fully convective limit at $\approx0.3 M_\odot$ for M dwarfs. The catalogue of 107 rotating M dwarfs presented here is particularly timely, as the stars are prime targets for the potential identification of transiting small planets with TESS and amenable to high-precision mass determination and further atmospheric characterization measurements.

V391 Peg (alias HS2201+2610) is a subdwarf B (sdB) pulsating star that shows both p- and g-modes. By studying the arrival times of the p-mode maxima and minima through the O-C method, in a previous article the presence of a planet was inferred with an orbital period of 3.2 yr and a minimum mass of 3.2 M_Jup. Here we present an updated O-C analysis using a larger data set of 1066 hours of photometric time series (~2.5x larger in terms of the number of data points), which covers the period between 1999 and 2012 (compared with 1999-2006 of the previous analysis). Up to the end of 2008, the new O-C diagram of the main pulsation frequency (f1) is compatible with (and improves) the previous two-component solution representing the long-term variation of the pulsation period (parabolic component) and the giant planet (sine wave component). Since 2009, the O-C trend of f1 changes, and the time derivative of the pulsation period (p_dot) passes from positive to negative; the reason of this change of regime is not clear and could be related to nonlinear interactions between different pulsation modes. With the new data, the O-C diagram of the secondary pulsation frequency (f2) continues to show two components (parabola and sine wave), like in the previous analysis. Various solutions are proposed to fit the O-C diagrams of f1 and f2, but in all of them, the sinusoidal components of f1 and f2 differ or at least agree less well than before. The nice agreement found previously was a coincidence due to various small effects that are carefully analysed. Now, with a larger dataset, the presence of a planet is more uncertain and would require confirmation with an independent method. The new data allow us to improve the measurement of p_dot for f1 and f2: using only the data up to the end of 2008, we obtain p_dot_1=(1.34+-0.04)x10**-12 and p_dot_2=(1.62+-0.22)x10**-12 ...

Routinely operating since July 2012, the APACHE survey has celebrated its second birthday. While the main goal of the Project is the detection of transiting planets around a large sample of bright, nearby M dwarfs in the northern hemisphere, the APACHE large photometric database for hundreds of different fields represents a relevant resource to search for and provide a first characterization of new variable stars. We celebrate here the conclusion of the second year of observations by reporting the discovery of 14 new variables.

We present more than 80 new variable stars discovered during the first observing season of the APACHE survey. APACHE is a project aimed at detecting extrasolar planets transiting nearby, bright M dwarfs by using an array of small-aperture telescopes. Despite the fact that the survey is targeted to a well-defined sample of cool stars, we also reduce and analyze data for all the detected field stars. Since July 2012 dozens of different stellar fields have been monitored, leading to the detection of several variables for which we propose a classification and estimate a period, when a periodicity is evident in the data. Thanks to the SuperWASP public archive, we have also retrieved and analyzed photometric data collected by the SWASP survey, which helped us to refine the classification and the period estimation of many variables found in the APACHE database. Some of the variables present peculiarities and thus are discussed separately.

First, we summarize the four-year long efforts undertaken to build the final setup of the APACHE Project, a photometric transit search for small-size planets orbiting bright, low-mass M dwarfs. Next, we describe the present status of the APACHE survey, officially started in July 2012 at the site of the Astronomical Observatory of the Autonomous Region of the Aosta Valley, in the Western Italian Alps. Finally, we briefly discuss the potentially far-reaching consequences of a multi-technique characterization program of the (potentially planet-bearing) APACHE targets.

We have used 605 days of photometric data from the Kepler spacecraft to study KIC 6614501, a close binary system with an orbital period of 0.15749747(25) days (3.779939 hours), that consists of a low-mass subdwarf B (sdB) star and a white dwarf. As seen in many other similar systems, the gravitational field of the white dwarf produces an ellipsoidal deformation of the sdB which appears in the light curve as a modulation at two times the orbital frequency. The ellipsoidal deformation of the sdB implies that the system has a maximum inclination of ~40 degrees, with i \approx 20 degrees being the most likely. The orbital radial velocity of the sdB star is high enough to produce a Doppler beaming effect with an amplitude of 432 \pm 5 ppm, clearly visible in the folded light curve. The photometric amplitude that we obtain, K1 = 85.8 km/s, is ~12 per cent less than the spectroscopic RV amplitude of 97.2 \pm 2.0 km/s. The discrepancy is due to the photometric contamination from a close object at about 5 arcsec North West of KIC 6614501, which is difficult to remove. The atmospheric parameters of the sdB star, Teff = 23 700 \pm 500 K and log g = 5.70 \pm 0.10, imply that it is a rare object below the Extreme Horizontal Branch (EHB), similar to HD 188112 (Heber et al. 2003). The comparison with different evolutionary tracks suggests a mass between ~0.18 and ~0.25 Msun, too low to sustain core helium burning. If the mass was close to 0.18-0.19 Msun, the star could be already on the final He-core WD cooling track. A higher mass, up to ~0.25 Msun, would be compatible with a He-core WD progenitor undergoing a cooling phase in a H-shell flash loop. A third possibility, with a mass between ~0.32 and ~0.40 Msun, can not be excluded and would imply that the sdB is a "normal" (but with an unusually low mass) EHB star burning He...

New data and results on the optical behaviour of the blazar PKS 0735+178 are presented. In addition the whole historical light curve, and a new photometric calibration of comparison stars are reported. Optical spectral indexes are calculated and studied on years scales, while several methods for time-series analysis are applied to the whole historical series and to each observing season of our data set. This allowed to search and identify optical variability modes, characteristic timescales and the signal power spectrum over 3 decades in time. In the last 10 years the optical flux of PKS 0735+178 exhibited a rather achromatic long-term behaviour and a variability mode resembling the shot-noise. The brightness level was in an intermediate/low state with a mild flaring activity and a superimposition/succession of rapid and slower flares, with no extraordinary/isolated outbursts but, at any rate, characterized by 1 major active phase in 2001. Several mid-term scales (days, weeks) were found, the more common falling into values of about 27-28 days, 50-56 days and 76-79 days. The rapid variability in the historical curve appear to be modulated by a general, slower and rather oscillating trend, where typical timescales of about 4.5, 8.5 and 11-13 years can be identified. This spectral and temporal analysis, accompanying our data publication, suggests the occurrence of distinctive variability signatures at days/weeks scales, that can likely be of transitory nature. On the other hand the possible pseudo-cyclical or multi-component modulations at long times could be more stable, recurrent and correlated to the bimodal radio flux behaviour and the twisted radio structure observed by many years in this blazar.

Context: Being dominated by non-thermal emission from aligned relativistic jets, blazars allow us to elucidate the physics of extragalactic jets, and, ltimately, how the energy is extracted from the central black hole in radio-loud active galactic nuclei. Aims: Crucial information is provided by broad-band spectral energy distributions (SEDs), their trends with luminosity and correlated multi-frequency variability. With this study we plan to obtain a database of contemporaneous radio-to-optical spectra of a sample of blazars, which are and will be observed by current and future high-energy satellites. Methods: Since December 2004 we are performing a monthly multi-frequency radio monitoring of a sample of 35 blazars at the antennas in Medicina and Noto. Contemporaneous near-IR and optical observations for all our observing epochs are organised. Results: Until June 2006 about 4000 radio measurements and 5500 near-IR and optical measurements were obtained. Most of the sources show significant variability in all observing bands. Here we present the multi-frequency data acquired during the first eighteen months of the project, and construct the SEDs for the best-sampled sources.

We present parallax determinations for six white dwarf candidates in the Torino Observatory Parallax Program. The absolute parallaxes are found with precisions at the 2-3 milliarcsecond level. For WD 1126+185 we find a distance incompatible with being a white dwarf, implying an incorrect classification. For WD 2216+484 we find our distance is consistent with a simple DA white dwarf rather than a composite system as previously proposed in the literature. In general it is found that the published photometric distance is an overestimate of the distance found here.

Eight optical and four radio observatories have been intensively monitoring the BL Lac object 0716+714 in the last years: 4854 data points have been collected in the UBVRI bands since 1994, while radio light curves extend back to 1978. Many of these data are presented here for the first time. The long-term trend shown by the optical light curves seems to vary with a characteristic time scale of about 3.3 years, while a longer period of 5.5-6 years seems to characterize the radio long-term variations. In general, optical colour indices are only weakly correlated with brightness. The radio flux behaviour at different frequencies is similar, but the flux variation amplitude decreases with increasing wavelength. The radio spectral index varies with brightness (harder when brighter), but the radio fluxes seem to be the sum of two different-spectrum contributions: a steady base level and a harder-spectrum variable component. Once the base level is removed, the radio variations appear as essentially achromatic, similarly to the optical behaviour. Flux variations at the higher radio frequencies lead the lower-frequency ones with week-month time scales. The behaviour of the optical and radio light curves is quite different, the broad radio outbursts not corresponding in time to the faster optical ones and the cross-correlation analysis indicating only weak correlation with long time lags. However, minor radio flux enhancements simultaneous with the major optical flares can be recognized, which may imply that the mechanism producing the strong flux increases in the optical band also marginally affects the radio one.

We present the most continuous data base of optical $BVR_{c}I_{c}$ observations ever published on the BL Lacertae object GC 0109+224, collected mainly by the robotic telescope of the Perugia University Observatory in the period November 1994-February 2002. These observations have been complemented by data from the Torino Observatory, collected in the period July 1995-January 1999, and Mt. Maidanak Observatory (December 2000). GC 0109+224 showed rapid optical variations and six major outbursts were observed at the beginning and end of 1996, in fall 1998, at the beginning and at the end of 2000, and at the beginning of 2002. Fast and large-amplitude drops characterized its flux behaviour. The $R_c$ magnitude ranged from 13.3 (16.16 mJy) to 16.46 (0.8 mJy), with a mean value of 14.9 (3.38 mJy). In the periods where we collected multi-filter observations, we analyzed colour and spectral indexes, and the variability patterns during some flares. The long-term behaviour seems approximatively achromatic, but during some isolated outbursts we found evidence of the typical loop-like hysteresis behaviour, suggesting that rapid optical variability is dominated by non-thermal cooling of a single emitting particle population. We performed also a statistical analysis of the data, through the discrete correlation function (DCF), the structure function (SF), and the Lomb-Scargle periodogram, to identify characteristic times scales, from days to months, in the light curves, and to quantify the mode of variability. We also include the reconstruction of the historical light curve and a photometric calibration of comparison stars, to favour further extensive optical monitoring of this interesting blazar.

New optical and radio data on the BL Lacertae object AO 0235+16 have been collected in the last four years by a wide international collaboration, which confirm the intense activity of this source. The optical data also include the results of the Whole Earth Blazar Telescope (WEBT) first-light campaign organized in November 1997. The optical spectrum is observed to basically steepen when the source gets fainter. We have investigated the existence of typical variability time scales and of possible correlations between the optical and radio emissions by means of visual inspection, Discrete Correlation Function analysis, and Discrete Fourier Transform technique. The major radio outbursts are found to repeat quasi-regularly with a periodicity of about 5.7 years; this period is also in agreement with the occurrence of some of the major optical outbursts, but not all of them.

Light curves of 3C 279 are presented in optical (R-band), X-rays (RXTE/PCA), and gamma rays (CGRO/EGRET) for 1999 Jan-Feb and 2000 Jan-Mar. During both of those epochs the gamma-ray levels were high, and all three observed bands demonstrated substantial variation, on time scales as short as one day. Correlation analyses provided no consistent pattern, although a rather significant optical/gamma-ray correlation was seen in 1999, with a gamma-ray lag of ~2.5 days, and there are other suggestions of correlations in the light curves. For comparison, correlation analysis is also presented for the gamma-ray and X-ray light curves during the large gamma ray flare in 1996 Feb and the two gamma-bright weeks leading up to it; the correlation at that time was strong, with a gamma-ray/X-ray offset of no more than 1 day.

BeppoSAX observed 8 times the gamma-ray bright blazar PKS 0528+134 in Feb. and Mar. 1997, during a multiwavelength campaign involving EGRET, onboard CGRO, and ground based telescopes. The source was in its faintest X-ray state observed so far, with an unabsorbed [2-10] keV flux of 2.5E-12 erg/cm2/s. The spectrum can be fitted with a power law, with energy index alpha_X=0.49+-0.07 between 0.1 and 10 keV. The high energy instrument PDS detected a flux in the 15-100 keV band which appears disconnected from the lower X-ray energies, most likely due to a contaminating source. We consider our findings in the context of the overall spectral energy distribution and discuss their implications for synchrotron and inverse Compton models.

We report on the BeppoSAX observations of the gamma-ray blazar PKS 0528+134 performed in Feb and Mar 1997, during a multiwavelength campaign involving EGRET and ground based telescopes. The source was in a faint and hard state, with energy spectral index alpha=0.48+-0.04 between 0.1 and 10 keV, and [2-10] keV flux of 2.7E-12 erg/cm2/s. No significant variability was observed. The source was detected in the 20-120 keV band by the PDS, with a flux lying slightly above the extrapolation from lower X-ray energies. Comparing this low state with previous higher states of the source, there is an indication that the X-ray spectrum hardens and the gamma-ray spectrum steepens when the source is fainter.