The second Gaia data release contains the identification of 147 535 low-mass ($\le 1.4 M_{\odot}$) rotational modulation variable candidates on (or close to) the main sequence, together with their rotation period and modulation amplitude. The richness, the period and amplitude range, and the photometric precision of this sample make it possible to unveil, for the first time, signatures of different surface inhomogeneity regimes in the amplitude-period density diagram. The modulation amplitude distribution shows a clear bimodality, with an evident gap at periods $P \le 2$ d. The low amplitude branch, in turn, shows a period bimodality with a main clustering at periods $P \approx$ 5 - 10 d and a secondary clustering of ultra-fast rotators at $P \le 0.5$ d. The amplitude-period multimodality is correlated with the position in the period-absolute magnitude (or period-color) diagram, with the low- and high-amplitude stars occupying different preferential locations. Here we argue that such a multimodality represents a further evidence of the existence of different regimes of surface inhomogeneities in young and middle-age low-mass stars and we lay out possible scenarios for their evolution, which manifestly include rapid transitions from one regime to another. In particular, the data indicate that stars spinning up close to break-up velocity undergo a very rapid change in their surface inhomogeneities configuration, which is revealed here for the first time. The multimodality can be exploited to identify field stars of age $\sim$ 100 -- 600 Myr belonging to the slow-rotator low-amplitude sequence, for which age can be estimated from the rotation period via gyrochronology relationships.

We highlight the power of the Gaia DR2 in studying many fine structures of the Hertzsprung-Russell diagram (HRD). Gaia allows us to present many different HRDs, depending in particular on stellar population selections. We do not aim here for completeness in terms of types of stars or stellar evolutionary aspects. Instead, we have chosen several illustrative examples. We describe some of the selections that can be made in Gaia DR2 to highlight the main structures of the Gaia HRDs. We select both field and cluster (open and globular) stars, compare the observations with previous classifications and with stellar evolutionary tracks, and we present variations of the Gaia HRD with age, metallicity, and kinematics. Late stages of stellar evolution such as hot subdwarfs, post-AGB stars, planetary nebulae, and white dwarfs are also analysed, as well as low-mass brown dwarf objects. The Gaia HRDs are unprecedented in both precision and coverage of the various Milky Way stellar populations and stellar evolutionary phases. Many fine structures of the HRDs are presented. The clear split of the white dwarf sequence into hydrogen and helium white dwarfs is presented for the first time in an HRD. The relation between kinematics and the HRD is nicely illustrated. Two different populations in a classical kinematic selection of the halo are unambiguously identified in the HRD. Membership and mean parameters for a selected list of open clusters are provided. They allow drawing very detailed cluster sequences, highlighting fine structures, and providing extremely precise empirical isochrones that will lead to more insight in stellar physics. Gaia DR2 demonstrates the potential of combining precise astrometry and photometry for large samples for studies in stellar evolution and stellar population and opens an entire new area for HRD-based studies.

The Gaia Data Release 2 (DR2): we summarise the processing and results of the identification of variable source candidates of RR Lyrae stars, Cepheids, long period variables (LPVs), rotation modulation (BY Dra-type) stars, delta Scuti & SX Phoenicis stars, and short-timescale variables. In this release we aim to provide useful but not necessarily complete samples of candidates. The processed Gaia data consist of the G, BP, and RP photometry during the first 22 months of operations as well as positions and parallaxes. Various methods from classical statistics, data mining and time series analysis were applied and tailored to the specific properties of Gaia data, as well as various visualisation tools. The DR2 variability release contains: 228'904 RR Lyrae stars, 11'438 Cepheids, 151'761 LPVs, 147'535 stars with rotation modulation, 8'882 delta Scuti & SX Phoenicis stars, and 3'018 short-timescale variables. These results are distributed over a classification and various Specific Object Studies (SOS) tables in the Gaia archive, along with the three-band time series and associated statistics for the underlying 550'737 unique sources. We estimate that about half of them are newly identified variables. The variability type completeness varies strongly as function of sky position due to the non-uniform sky coverage and intermediate calibration level of this data. The probabilistic and automated nature of this work implies certain completeness and contamination rates which are quantified so that users can anticipate their effects. This means that even well-known variable sources can be missed or misidentified in the published data. The DR2 variability release only represents a small subset of the processed data. Future releases will include more variable sources and data products; however, DR2 shows the (already) very high quality of the data and great promise for variability studies.

We present the methods devised to identify the BY Dra variables candidates in Gaia DR2 and infer their variability parameters. BY Dra candidates are pre-selected from their position in the HR diagram, built from Gaia parallaxes, $G$ magnitudes, and $(G_{BP} - G_{RP})$ colours. Since the time evolution of the stellar active region can disrupt the coherence of the signal, segments not much longer than their expected evolution timescale are extracted from the entire photometric time-series and period search algorithms are applied to each segment. For the Gaia DR2, we select sources having similar period in at least two segments as candidates BY Dra. Results are further filtered considering the time series phase coverage and the expected approximate light curve shape. Gaia DR2 includes rotational periods and modulation amplitudes of 147 535 BY Dra candidates. The data unveil the existence of two populations with distinctive period and amplitude distributions. The sample covers 38% of the whole sky when divided in bins (HEALPix) of $\approx$0.84 square degrees and we estimate that represents 0.7 -- 5 % of all BY Dra stars potentially detectable by Gaia. The preliminary data contained in Gaia DR2 illustrate the vast and unique information that the mission is going to provide on stellar rotation and magnetic activity. This information, complemented by Gaia exquisite parallaxes, proper motions, and astrophysical parameter, is opening new and unique perspectives for our understanding of the evolution of stellar angular momentum and dynamo action.

Magnetic cycles have been detected in tens of solar-like stars. The relationship between the cycle properties and global stellar parameters is not fully understood yet. We searched for activity cycles in 90 solar-like stars with ages between 4 and 95 Myr aiming to investigate the properties of activity cycles in this age range. We measured the length $P_{ cyc}$ of a given cycle by analyzing the long-term time-series of three activity indexes. For each star, we computed also the global magnetic activity index <IQR> that is proportional to the amplitude of the rotational modulation and is a proxy of the mean level of the surface magnetic activity. We detected activity cycles in 67 stars. Secondary cycles were also detected in 32 stars. The lack of correlation between $P_{ cyc}$ and $P_{ rot}$ suggest that these stars belong to the Transitional Branch and that the dynamo acting in these stars is different from the solar one. This statement is also supported by the analysis of the butterfly diagrams. We computed the Spearman correlation coefficient $r_{ S}$ between $P_{ cyc}$, <IQR> and different stellar parameters. We found that $P_{ cyc}$ is uncorrelated with all the investigated parameters. The <IQR> index is positively correlated with the convective turn-over time-scale, the magnetic diffusivity time-scale $\tau_{ diff}$, and the dynamo number $D_{ N}$, whereas it is anti-correlated with the effective temperature $T_{ eff}$, the photometric shear $\Delta\Omega_{\rm phot}$ and the radius $R_{ C}$ at which the convective zone is located. We found that $P_{ cyc}$ is about constant and that <IQR> decreases with the stellare age in the range 4-95 Myr. We investigated the magnetic activity of AB Dor A by merging ASAS time-series with previous long-term photometric data. We estimated the length of the AB Dor A primary cycle as $P_{ cyc} = 16.78 \pm 2 \rm yr$.

Parallaxes for 331 classical Cepheids, 31 Type II Cepheids and 364 RR Lyrae stars in common between Gaia and the Hipparcos and Tycho-2 catalogues are published in Gaia Data Release 1 (DR1) as part of the Tycho-Gaia Astrometric Solution (TGAS). In order to test these first parallax measurements of the primary standard candles of the cosmological distance ladder, that involve astrometry collected by Gaia during the initial 14 months of science operation, we compared them with literature estimates and derived new period-luminosity ($PL$), period-Wesenheit ($PW$) relations for classical and Type II Cepheids and infrared $PL$, $PL$-metallicity ($PLZ$) and optical luminosity-metallicity ($M_V$-[Fe/H]) relations for the RR Lyrae stars, with zero points based on TGAS. The new relations were computed using multi-band ($V,I,J,K_{\mathrm{s}},W_{1}$) photometry and spectroscopic metal abundances available in the literature, and applying three alternative approaches: (i) by linear least squares fitting the absolute magnitudes inferred from direct transformation of the TGAS parallaxes, (ii) by adopting astrometric-based luminosities, and (iii) using a Bayesian fitting approach. TGAS parallaxes bring a significant added value to the previous Hipparcos estimates. The relations presented in this paper represent first Gaia-calibrated relations and form a "work-in-progress" milestone report in the wait for Gaia-only parallaxes of which a first solution will become available with Gaia's Data Release 2 (DR2) in 2018.

Context. The first Gaia Data Release contains the Tycho-Gaia Astrometric Solution (TGAS). This is a subset of about 2 million stars for which, besides the position and photometry, the proper motion and parallax are calculated using Hipparcos and Tycho-2 positions in 1991.25 as prior information. Aims. We investigate the scientific potential and limitations of the TGAS component by means of the astrometric data for open clusters. Methods. Mean cluster parallax and proper motion values are derived taking into account the error correlations within the astrometric solutions for individual stars, an estimate of the internal velocity dispersion in the cluster, and, where relevant, the effects of the depth of the cluster along the line of sight. Internal consistency of the TGAS data is assessed. Results. Values given for standard uncertainties are still inaccurate and may lead to unrealistic unit-weight standard deviations of least squares solutions for cluster parameters. Reconstructed mean cluster parallax and proper motion values are generally in very good agreement with earlier Hipparcos-based determination, although the Gaia mean parallax for the Pleiades is a significant exception. We have no current explanation for that discrepancy. Most clusters are observed to extend to nearly 15 pc from the cluster centre, and it will be up to future Gaia releases to establish whether those potential cluster-member stars are still dynamically bound to the clusters. Conclusions. The Gaia DR1 provides the means to examine open clusters far beyond their more easily visible cores, and can provide membership assessments based on proper motions and parallaxes. A combined HR diagram shows the same features as observed before using the Hipparcos data, with clearly increased luminosities for older A and F dwarfs.

The ESA Gaia mission provides a unique time-domain survey for more than one billion sources brighter than G=20.7 mag. Gaia offers the unprecedented opportunity to study variability phenomena in the Universe thanks to multi-epoch G-magnitude photometry in addition to astrometry, blue and red spectro-photometry, and spectroscopy. Within the Gaia Consortium, Coordination Unit 7 has the responsibility to detect variable objects, classify them, derive characteristic parameters for specific variability classes, and provide global descriptions of variable phenomena. We describe the variability processing and analysis that we plan to apply to the successive data releases, and we present its application to the G-band photometry results of the first 14 months of Gaia operations that comprises 28 days of Ecliptic Pole Scanning Law and 13 months of Nominal Scanning Law. Out of the 694 million, all-sky, sources that have calibrated G-band photometry in this first stage of the mission, about 2.3 million sources that have at least 20 observations are located within 38 degrees from the South Ecliptic Pole. We detect about 14% of them as variable candidates, among which the automated classification identified 9347 Cepheid and RR Lyrae candidates. Additional visual inspections and selection criteria led to the publication of 3194 Cepheid and RR Lyrae stars, described in Clementini et al. (2016). Under the restrictive conditions for DR1, the completenesses of Cepheids and RR Lyrae stars are estimated at 67% and 58%, respectively, numbers that will significantly increase with subsequent Gaia data releases. Data processing within the Gaia Consortium is iterative, the quality of the data and the results being improved at each iteration. The results presented in this article show a glimpse of the exceptional harvest that is to be expected from the Gaia mission for variability phenomena. [abridged]

In case of pre-main sequence objects, the only way to determine age and mass is by fitting theoretical isochrones on color-magnitude (alternatively luminosity-temperature) diagrams. Since young stellar objects exhibit photometric variability over wide range in magnitude and colors, the age and mass determined by fitting isochrones is expected to be inaccurate, if not erroneous. These in turn will badly affect any study carried out on age spread and process of star formation. Since we have carried out very extensive photometric observations of the Orion Nebula Cluster (ONC), we decided to use our multi-band data to explore the influence of variability in determining mass and age of cluster members. In this study, we get the amplitudes of the photometric variability in V, R, and I optical bands of a sample of 346 ONC members and use it to investigate how the variability affects the inferred masses and ages and if it alone can take account for the age spread among the ONC members reported by earlier studies. We find that members that show periodic and smooth photometric rotational modulation have their masses and ages unaffected by variability. On other hand, we found that members with periodic but very scattered photometric rotational modulation and members with irregular variability have their masses and ages significantly affected. Moreover, using Hertzsprung-Russell (HR) diagrams we find that the observed I band photometric variability can take account of only a fraction (about 50%) of the inferred age spread, whereas the V band photometric variability is large enough to mask any age spread.

Average stellar radii in open clusters can be estimated from rotation periods and projected rotational velocities under the assumption of random orientation of the spin axis. Such estimates are independent of distance, interstellar absorption, and models, but their validity can be limited by missing data (truncation) or data that only represent upper/lower limits (censoring). We present a new statistical analysis method to estimate average stellar radii in the presence of censoring and truncation. We use theoretical distribution functions of the projected stellar radius $R \sin i$ to define a likelihood function in the presence of censoring and truncation. Average stellar radii in magnitude bins are then obtained by a maximum likelihood parametric estimation procedure. This method is capable of recovering the average stellar radius within a few percent with as few as $\approx$ 10 measurements. Here it is applied for the first time to the dataset available for the Pleiades. We find an agreement better than $\approx$ 10 percent between the observed $R$ vs $M_K$ relationship and current standard stellar models for 1.2 $\ge M/M_{\odot} \ge$ 0.85 with no evident bias. Evidence of a systematic deviation at $2\sigma$ level are found for stars with 0.8 $\ge M/M_{\odot} \ge$ 0.6 approaching the slow-rotator sequence. Fast-rotators ($P$ < 2 d) agree with standard models within 15 percent with no systematic deviations in the whole 1.2 $\ge M/M_{\odot} \ge$ 0.5 range. The evidence found of a possible radius inflation just below the lower mass limit of the slow-rotator sequence indicates a possible connection with the transition from the fast to the slow-rotator sequence.

Surface differential rotation (SDR) plays a key role in dynamo models. SDR estimates are therefore essential for constraining theoretical models. We measure a lower limit to SDR in a sample of solar-like stars belonging to young associations with the aim of investigating how SDR depends on global stellar parameters in the age range (4-95 Myr). The rotation period of a solar-like star can be recovered by analyzing the flux modulation caused by dark spots and stellar rotation. The SDR and the latitude migration of dark-spots induce a modulation of the detected rotation period. We employ long-term photometry to measure the amplitude of such a modulation and to compute the quantity DeltaOmega_phot =2p/P_min -2pi/P_max that is a lower limit to SDR. We find that DeltaOmega_phot increases with the stellar effective temperature and with the global convective turn-over time-scale tau_c. We find that DeltaOmega_phot is proportional to Teff^2.18pm 0.65 in stars recently settled on the ZAMS. This power law is less steep than those found by previous authors, but closest to recent theoretical models. We find that DeltaOmega_phot steeply increases between 4 and 30 Myr and that itis almost constant between 30 and 95 Myr in a 1 M_sun star. We find also that the relative shear increases with the Rossby number Ro. Although our results are qualitatively in agreement with hydrodynamical mean-field models, our measurements are systematically higher than the values predicted by these models. The discrepancy between DeltaOmega_phot measurements and theoretical models is particularly large in stars with periods between 0.7 and 2 d. Such a discrepancy, together with the anomalous SDR measured by other authors for HD 171488 (rotating in 1.31 d), suggests that the rotation period could influence SDR more than predicted by the models.

We present the results of our analysis on V1481 Ori (JW 239), a young SB2 in the Orion Nebula Cluster with a circumbinary disc accreting on the lower-mass component. The analysis is based on high-resolution spectroscopic data and high-quality photometric time series about 20-yr long. Thanks to the spectroscopy, we confirm the binary nature of this system consisting of M3 + M4 components and derive the mass ratio M_B/M_A = 0.54, a variable luminosity ratio L_B/L_A = 0.68--0.94, and an orbital period P_orb = 4.433d. The photometric data allowed us to measure the rotation periods of the two components P_phot = 4.4351d and they are found to be synchronized with the orbital period. The simultaneous modeling of V-, I-band, and radial velocity curves in the 2005 season suggests that the variability is dominated by one hot spot on the secondary component covering at least about 3.5% of the stellar surface and about 420K hotter than the unperturbed photosphere. Such a spot may originate from the material of the circumbinary disc accreting onto the secondary component. We also detect an apparent 6-yr periodic variation in the position of this hot spot, which is inferred from the phase migration of the light curve maximum, which we interpret as due to either the presence of surface differential rotation as large as 0.065%, a value compatible with the fully convective components, or to a periodic exchange of angular momentum between the disc and the star, which implies a minimum magnetic field strength of 650G at the stellar surface.

We present a study on the determination of rotation periods (P) of solar-like stars from the photometric irregular time-sampling of the ESA Gaia mission, currently scheduled for launch in 2013, taking into account its dependence on ecliptic coordinates. We examine the case of solar-twins as well as thousands of synthetic time-series of solar-like stars rotating faster than the Sun. In the case of solar twins we assume that the Gaia unfiltered photometric passband G will mimic the variability of the total solar irradiance (TSI) as measured by the VIRGO experiment. For stars rotating faster than the Sun, light-curves are simulated using synthetic spectra for the quiet atmosphere, the spots, and the faculae combined by applying semi-empirical relationships relating the level of photospheric magnetic activity to the stellar rotation and the Gaia instrumental response. The capabilities of the Deeming, Lomb-Scargle, and Phase Dispersion Minimisation methods in recovering the correct rotation periods are tested and compared. The false alarm probability (FAP) is computed using Monte Carlo simulations and compared with analytical formulae. The Gaia scanning law makes the rate of correct detection of rotation periods strongly dependent on the ecliptic latitude (beta). We find that for P ~ 1 d, the rate of correct detection increases with ecliptic latitude from 20-30 per cent at beta ~ 0{\deg} to a peak of 70 per cent at beta=45{\deg}, then it abruptly falls below 10 per cent at beta > 45{\deg}. For P > 5 d, the rate of correct detection is quite low and for solar twins is only 5 per cent on average.

We present some preliminary results from our program of intensive near-infrared photometric monitoring of a sample of confirmed and candidate Luminous Blue Variables (LBVs) conducted from 2008 to 2010. Clear long-term variability has been observed for Wray 17-96 and V481 Sct, with overall brightness variation greater than 1 mag in the J band. Other sources, such as LBV 1806-20 showed detectable variability with amplitudes of few tenths of a magnitude with time-scale of about 60 days.

Since 2004, we have been engaged in a long-term observing program to monitor young stellar objects in the Orion Nebula Cluster. We have collected about two thousands frames in V, R, and I broad-band filters on more than two hundred nights distributed over five consecutive observing seasons. The high-quality and time-extended photometric data give us an opportunity to address various phenomena associated with young stars. The prime motivations of this project are i) to explore various manifestations of stellar magnetic activity in very young low-mass stars; ii) to search for new pre-main sequence eclipsing binaries; and iii) to look for any EXor and FUor like transient activities associated with YSOs. Since this is the first paper on this program, we give a detailed description of the science drivers, the observation and the data reduction strategies as well. In addition to these, we also present a large number of new periodic variables detected from our first five years of time-series photometric data. Our study reveals that about 72% of CTTS in our FoV are periodic, whereas, the percentage of periodic WTTS is just 32%. This indicates that inhomogeneities patterns on the surface of CTTS of the ONC stars are much more stable than on WTTS. From our multi-year monitoring campaign we found that the photometric surveys based on single-season are incapable of identifying all periodic variables. And any study on evolution of angular momentum based on single-season surveys must be carried out with caution.

We report the outcome of a deep multi-wavelength study of the IC2391 young open cluster. We aim at uncovering new low-mass and sub-stellar members of the cluster and identifying new debris disk objects. A 30*30 square arcmin area in IC 2391 was observed using the wide-field imager at the ESO 2.2m telescope. The completeness limits of the photometry at 3 sigma level are V=24.7, Rc=23.7 and Ic=23.0, faint enough to reveal sub-stellar members down to about 0.03 solar masses. Our membership criteria are based on the use of our optical data, in combination with JHKs magnitudes from the 2MASS catalog. We also estimate the physical parameters of the selected candidates. Debris disk candidates are identified on the basis of their infrared excess emission using near- and mid-infrared photometry from the Spitzer Space Telescope. Our optical survey, which has a limiting magnitude at 3 sigma level 1-2 mag fainter than previous optical surveys conducted in IC2391, revealed 29 new low-mass member candidates of the cluster. We estimate the contamination to be at least 50%. We constrain the fraction of sub-stellar objects in the range 8-15% and discuss possible explanations for the deficit of brown dwarfs in this cluster. We also identified 10 candidates in the cluster showing IR excess emission consistent with the presence of debris disks.

Beating the Earth's day-night cycle is mandatory for long and continuous time-series photometry and had been achieved with either large ground-based networks of observatories at different geographic longitudes or when conducted from space. A third possibility is offered by a polar location with astronomically-qualified site characteristics. Aims. In this paper, we present the first scientific stellar time-series optical photometry from Dome C in Antarctica and analyze approximately 13,000 CCD frames taken in July 2007. We conclude that high-precision CCD photometry with exceptional time coverage and cadence can be obtained at Dome C in Antarctica and be successfully used for time-series astrophysics.

Rotation and solar-type magnetic activity are closely related to each other in stars of G or later spectral types. Presence and level of magnetic activity depend on star's rotation and rotation itself is strongly influenced by strength and topology of the magnetic fields. Open clusters represent especially useful targets to investigate the connection between rotation and activity. The open cluster NGC2099 has been studied as a part of the RACE-OC project, which is aimed at exploring the evolution of rotation and magnetic activity in the late-type members of open clusters of different ages. Time series CCD photometric observations of this cluster were collected during January 2004. The relations between activity manifestations, such as the light curve amplitude, and global stellar parameters are investigated. We have discovered 135 periodic variables, 122 of which are candidate cluster members. Determination of rotation periods of G- and K-type stars has allowed us to better explore evolution of angular momentum at an age of about 500 Myr. A comparison with the older Hyades cluster (~625 Myr) shows that the newly determined distribution of rotation periods is consistent with the scenario of rotational braking of main-sequence spotted stars as they age. However, a comparison with the younger M34 cluster (~200 Myr) shows that the G8-K5 members of these clusters have the same rotation period distribution, that is G8-K5 members in NGC2099 seem to have experienced no significant braking in the age range from ~200 to ~500 Myr. Finally, NGC2099 members have a level of photospheric magnetic activity, as measured by light curve amplitude, smaller than in younger stars of same mass and rotation, suggesting that the activity level also depends on some other age-dependent parameters.

Spectral variability is the main tool for constraining emission models of BL Lac objects. By means of systematic observations of the BL Lac prototype PKS 2155-304 in the infrared-optical band, we explore variability on the scales of months, days and hours. We made our observations with the robotic 60 cm telescope REM located at La Silla, Chile. VRIJHK filters were used. PKS 2155-304 was observed from May to December 2005. The wavelength interval explored, the total number of photometric points and the short integration time render our photometry substantially superior to previous ones for this source. On the basis of the intensity and colour we distinguish three different states of the source, each of duration of months, which include all those described in the literature. In particular, we report the highest state ever detected in the H band. The source varied by a factor of 4 in this band, much more than in the V band (a factor ~2). The source softened with increasing intensity, contrary to the general pattern observed in the UV-X-ray bands. On five nights of November we had nearly continuous monitoring for 2-3 hours. A variability episode with a time scale of ~24 h is well documented, a much more rapid flare with t=1-2 h, is also apparent, but is supported by relatively few points.

The successful launch of the CoRoT (COnvection, ROtation and planetary Transits) satellite opens a new era in asteroseismology. The space photometry is complemented by high-resolution spectroscopy and multicolour photometry from ground, to disclose the pulsational content of the asteroseismic targets in the most complete way. Some preliminary results obtained with both types of data are presented. The paper is based on observations collected at S. Pedro Martir, Serra La Nave, La Silla, and Telescopio Nazionale Galileo Observatories.