• ### Gaia Data Release 2: Summary of the variability processing & analysis results(1804.09373)

May 8, 2018 astro-ph.SR, astro-ph.IM
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.
• ### Gaia Data Release 2: The first Gaia catalog of Long Period Variable candidates(1805.02035)

May 5, 2018 astro-ph.SR, astro-ph.IM
Gaia Data Release 2 (DR2) provides a unique all-sky catalog of 550'737 variable stars, of which 151'761 are Long Period Variable (LPV) candidates with G variability amplitudes larger than 0.2 mag. About one fifth of the LPV candidates are Mira candidates, the majority of the rest being semi-regular variable candidates. For each source, G, BP and RP photometric time series are published, together with some LPV-specific attributes for the subset of 89'617 candidates with periods in G larger than 60 days. We describe this first Gaia catalog of LPV candidates, give an overview of its content, and present various validation checks. Various samples of LPVs are used to validate the catalog: a sample of well-studied very bright LPVs having AAVSO light curves concomitant with Gaia light curves, a sample of Gaia LPV candidates with good parallaxes, the ASAS_SN all-sky catalog of LPVs, and the OGLE-III catalogs of LPVs towards the Magellanic Clouds and the Galactic Bulge. The analyses of these samples show a good agreement between Gaia DR2 and literature periods. The same is globally true for bolometric corrections of M-type stars. The main contaminant to our DR2 catalog comes from Young Stellar Objects (YSOs) in the Solar vicinity (within ~1 kpc), their number in the whole catalog being only at the percent level. A cautious note is provided about parallax-dependent LPV attributes published in the catalog. This first Gaia catalog of LPVs about doubles the number of known LPVs with amplitudes larger than 0.2 mag, despite the conservative candidate selection criteria in order to prioritize low contamination over high completeness, and despite the limited DR2 time coverage compared to the long periods characteristic of LPVs. It also contains a small set of YSO candidates, which offers the serendipitous opportunity to study these objects at an early stage of the Gaia data releases.
• ### Gaia Data Release 2: Specific characterisation and validation of all-sky Cepheids and RR Lyrae stars(1805.02079)

May 5, 2018 astro-ph.SR
Gaia second Data Release (DR2) presents a first mapping of full-sky RR Lyrae stars and Cepheids observed by the spacecraft during the initial 22 months of observations and publishes characteristic parameters derived for these sources by the Specific Objects Study (SOS) pipeline developed to validate and fully characterise Cepheids and RR Lyrae stars (SOS Cep&RRL) observed by Gaia. The SOS Cep&RRL processing uses tools such as the period-amplitude (PA) and the period-luminosity (PL) relations in the G-band. For the Gaia DR2 data processing we also used tools based on the G_BP and G_RP photometry, such as the period-Wesenheit (PW) relation in G,G_RP. Furthermore, we implemented the use of parallaxes working directly in parallax space and applied different PL, PW relations depending on the source position on sky, whether in the LMC, in the SMC or outside them. G, G_BP and G_RP time series photometry and characterisation by the SOS Cep&RRL pipeline (mean magnitudes and pulsation characteristics) are published in Gaia DR2 for a total of 150,359 sources distributed all over the sky: 9,575 are classified as Cepheids and 140,784 as RR Lyrae stars. These samples include also variables in 87 globular clusters and 12 dwarf galaxies. To the best of our knowledge, as of 25 April 2018, 50,570 of these variables (about 350 Cepheids and 50,220 RR Lyrae stars) do not have a known counterpart in the literature, hence they are likely new discoveries by Gaia. Furthermore, an estimate of the interstellar absorption is published for 54,272 fundamental mode RR Lyrae stars from a relation based on the amplitude of the light variation in the G-band and the star period. Photometric metal abundances ([Fe/H]) derived from the Fourier parameters of the light curves are also released for 64,932 RR Lyrae stars and for 3,738 fundamental-mode classical Cepheids with period shorter than 6.3 days.
• ### Gaia Data Release 2: The Short Timescale Variability Processing and Analysis(1805.00747)

May 3, 2018 astro-ph.IM
The Gaia DR2 short timescale variable candidates sample results from the investigation of the first 22 months of Gaia $G$ per-CCD, $G_{BP}$ and $G_{RP}$ photometry, for a subsample of sources at the Gaia faint end ($G \sim 16.5 - 20\,$mag). For this first Gaia short timescale variability search, we limit ourselves to the case of rapid, suspected periodic variability. Our study combines fast variability detection through variogram analysis, Least-Square high frequency search, and empirical selection criterion based on various statistics and built from the investigation of specific sources seen through Gaia eyes (e.g. known variables or visually identified objects with peculiar features in their light-curves). The progressive selection criterion definition, improvement and validation also make use of supplementary ground-based photometric monitoring, performed at the Flemish Mercator telescope in La Palma (Canary Islands, Spain) between August and November 2017. We publish a list of 3018 bona fide, suspected periodic, short timescale variable candidates, spread all over the sky, with a contamination level from false positives and non-periodic variables up to 10-20\% in the Magellanic Clouds. Though its completeness is around 0.05\%, the Gaia DR2 short timescale variables sample recover very interesting known short period variables, such as Post Common Envelope Binaries or Cataclysmic Variables, and points fascinating newly discovered variables sources. Several improvements in the short timescale variability processing are considered for the future Gaia Data Releases, by enhancing the existing variogram and period search algorithms or going one step beyond with the classification of the identified candidates. The encouraging outcome of our analysis demonstrates the power of the Gaia mission for such fast variability studies and opens great perspectives for this domain of astrophysics.
• ### Gaia Data Release 2: Rotational modulation in late-type dwarfs(1805.00421)

May 1, 2018 astro-ph.SR
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.
• 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.
• ### Gaia Data Release 2: Variable stars in the colour-absolute magnitude diagram(1804.09382)

April 25, 2018 astro-ph.SR
The ESA Gaia mission provides a unique time-domain survey for more than 1.6 billion sources with G ~ 21 mag. We showcase stellar variability across the Galactic colour-absolute magnitude diagram (CaMD), focusing on pulsating, eruptive, and cataclysmic variables, as well as on stars exhibiting variability due to rotation and eclipses. We illustrate the locations of variable star classes, variable object fractions, and typical variability amplitudes throughout the CaMD and illustrate how variability-related changes in colour and brightness induce `motions' using 22 months worth of calibrated photometric, spectro-photometric, and astrometric Gaia data of stars with significant parallax. To ensure a large variety of variable star classes to populate the CaMD, we crossmatch Gaia sources with known variable stars. We also used the statistics and variability detection modules of the Gaia variability pipeline. Corrections for interstellar extinction are not implemented in this article. Gaia enables the first investigation of Galactic variable star populations across the CaMD on a similar, if not larger, scale than previously done in the Magellanic Clouds. Despite observed colours not being reddening corrected, we clearly see distinct regions where variable stars occur and determine variable star fractions to within Gaia's current detection thresholds. Finally, we show the most complete description of variability-induced motion within the CaMD to date. Gaia enables novel insights into variability phenomena for an unprecedented number of stars, which will benefit the understanding of stellar astrophysics. The CaMD of Galactic variable stars provides crucial information on physical origins of variability in a way previously accessible only for Galactic star clusters or external galaxies.
• 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.
• ### Gaia eclipsing binary and multiple systems. Two-Gaussian models applied to OGLE-III eclipsing binary light curves in the Large Magellanic Cloud(1703.10597)

March 30, 2017 astro-ph.SR, astro-ph.IM
The advent of large scale multi-epoch surveys raises the need for automated light curve (LC) processing. This is particularly true for eclipsing binaries (EBs), which form one of the most populated types of variable objects. The Gaia mission, launched at the end of 2013, is expected to detect of the order of few million EBs over a 5-year mission. We present an automated procedure to characterize EBs based on the geometric morphology of their LCs with two aims: first to study an ensemble of EBs on a statistical ground without the need to model the binary system, and second to enable the automated identification of EBs that display atypical LCs. We model the folded LC geometry of EBs using up to two Gaussian functions for the eclipses and a cosine function for any ellipsoidal-like variability that may be present between the eclipses. The procedure is applied to the OGLE-III data set of EBs in the Large Magellanic Cloud (LMC) as a proof of concept. The bayesian information criterion is used to select the best model among models containing various combinations of those components, as well as to estimate the significance of the components. Based on the two-Gaussian models, EBs with atypical LC geometries are successfully identified in two diagrams, using the Abbe values of the original and residual folded LCs, and the reduced $\chi^2$. Cleaning the data set from the atypical cases and further filtering out LCs that contain non-significant eclipse candidates, the ensemble of EBs can be studied on a statistical ground using the two-Gaussian model parameters. For illustration purposes, we present the distribution of projected eccentricities as a function of orbital period for the OGLE-III set of EBs in the LMC, as well as the distribution of their primary versus secondary eclipse widths.
• 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.
• ### Gaia Data Release 1: The variability processing & analysis and its application to the south ecliptic pole region(1702.03295)

Feb. 10, 2017 astro-ph.SR, astro-ph.IM
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]
• ### Gaia Data Release 1 - The Cepheid & RR Lyrae star pipeline and its application to the south ecliptic pole region(1609.04269)

Oct. 14, 2016 astro-ph.GA, astro-ph.SR
We present an overview of the Specific Objects Study (SOS) pipeline developed within the Coordination Unit 7 (CU7) of the Gaia Data Processing and Analysis Consortium (DPAC), the coordination unit charged with the processing and analysis of variable sources observed by Gaia, to validate and fully characterise Cepheids and RR Lyrae stars observed by the spacecraft. We describe how the SOS for Cepheids and RR Lyrae stars (SOS Cep&RRL) was specifically tailored to analyse Gaia's G-band photometric time-series with a South Ecliptic Pole (SEP) footprint, which covers an external region of the Large Magellanic Cloud (LMC). G-band time-series photometry and characterization by the SOS Cep&RRL pipeline (mean magnitude and pulsation characteristics) are published in Gaia Data Release 1 (Gaia DR1) for a total sample of 3,194 variable stars, 599 Cepheids and 2,595 RR Lyrae stars, of which 386 (43 Cepheids and 343 RR Lyrae stars) are new discoveries by Gaia. All 3,194 stars are distributed over an area extending 38 degrees on either side from a point offset from the centre of the LMC by about 3 degrees to the north and 4 degrees to the east. The vast majority, but not all, are located within the LMC. The published sample also includes a few bright RR Lyrae stars that trace the outer halo of the Milky Way in front of the LMC.
• ### Automated classification of Hipparcos unsolved variables(1301.1545)

Jan. 8, 2013 astro-ph.SR, astro-ph.IM
We present an automated classification of stars exhibiting periodic, non-periodic and irregular light variations. The Hipparcos catalogue of unsolved variables is employed to complement the training set of periodic variables of Dubath et al. with irregular and non-periodic representatives, leading to 3881 sources in total which describe 24 variability types. The attributes employed to characterize light-curve features are selected according to their relevance for classification. Classifier models are produced with random forests and a multistage methodology based on Bayesian networks, achieving overall misclassification rates under 12 per cent. Both classifiers are applied to predict variability types for 6051 Hipparcos variables associated with uncertain or missing types in the literature.
• ### Search for high-amplitude Delta Scuti and RR Lyrae stars in Sloan Digital Sky Survey Stripe 82 using principal component analysis(1203.6196)

March 28, 2012 astro-ph.IM
We propose a robust principal component analysis (PCA) framework for the exploitation of multi-band photometric measurements in large surveys. Period search results are improved using the time series of the first principal component due to its optimized signal-to-noise ratio.The presence of correlated excess variations in the multivariate time series enables the detection of weaker variability. Furthermore, the direction of the largest variance differs for certain types of variable stars. This can be used as an efficient attribute for classification. The application of the method to a subsample of Sloan Digital Sky Survey Stripe 82 data yielded 132 high-amplitude Delta Scuti variables. We found also 129 new RR Lyrae variables, complementary to the catalogue of Sesar et al., 2010, extending the halo area mapped by Stripe 82 RR Lyrae stars towards the Galactic bulge. The sample comprises also 25 multiperiodic or Blazhko RR Lyrae stars.
• ### Hipparcos Variable Star Detection and Classification Efficiency(1107.3638)

July 21, 2011 astro-ph.SR
A complete periodic star extraction and classification scheme is set up and tested with the Hipparcos catalogue. The efficiency of each step is derived by comparing the results with prior knowledge coming from the catalogue or from the literature. A combination of two variability criteria is applied in the first step to select 17 006 variability candidates from a complete sample of 115 152 stars. Our candidate sample turns out to include 10 406 known variables (i.e., 90% of the total of 11 597) and 6600 contaminating constant stars. A random forest classification is used in the second step to extract 1881 (82%) of the known periodic objects while removing entirely constant stars from the sample and limiting the contamination of non-periodic variables to 152 stars (7.5%). The confusion introduced by these 152 non-periodic variables is evaluated in the third step using the results of the Hipparcos periodic star classification presented in a previous study (Dubath et al. [1]).
• ### Random forest automated supervised classification of Hipparcos periodic variable stars(1101.2406)

July 19, 2011 astro-ph.SR
We present an evaluation of the performance of an automated classification of the Hipparcos periodic variable stars into 26 types. The sub-sample with the most reliable variability types available in the literature is used to train supervised algorithms to characterize the type dependencies on a number of attributes. The most useful attributes evaluated with the random forest methodology include, in decreasing order of importance, the period, the amplitude, the V-I colour index, the absolute magnitude, the residual around the folded light-curve model, the magnitude distribution skewness and the amplitude of the second harmonic of the Fourier series model relative to that of the fundamental frequency. Random forests and a multi-stage scheme involving Bayesian network and Gaussian mixture methods lead to statistically equivalent results. In standard 10-fold cross-validation experiments, the rate of correct classification is between 90 and 100%, depending on the variability type. The main mis-classification cases, up to a rate of about 10%, arise due to confusion between SPB and ACV blue variables and between eclipsing binaries, ellipsoidal variables and other variability types. Our training set and the predicted types for the other Hipparcos periodic stars are available online.
• ### Variability type classification of multi-epoch surveys(0901.2835)

Jan. 19, 2009 astro-ph.IM
The classification of time series from photometric large scale surveys into variability types and the description of their properties is difficult for various reasons including but not limited to the irregular sampling, the usually few available photometric bands, and the diversity of variable objects. Furthermore, it can be seen that different physical processes may sometimes produce similar behavior which may end up to be represented as same models. In this article we will also be presenting our approach for processing the data resulting from the Gaia space mission. The approach may be classified into following three broader categories: supervised classification, unsupervised classifications, and "so-called" extractor methods i.e. algorithms that are specialized for particular type of sources. The whole process of classification- from classification attribute extraction to actual classification- is done in an automated manner.