• ### 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.
• 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.
• ### KIC 8164262: a Heartbeat Star Showing Tidally Induced Pulsations with Resonant Locking(1706.05051)

Oct. 12, 2017 astro-ph.SR
We present the analysis of KIC 8164262, a heartbeat star with a high-amplitude (~1 mmag), tidally resonant pulsation (a mode in resonance with the orbit) at 229 times the orbital frequency and a plethora of tidally induced g-mode pulsations (modes excited by the orbit). The analysis combines Kepler light curves with follow-up spectroscopic data from the Keck telescope, KPNO (Kitt Peak National Observatory) 4-m Mayal telescope and the 2.7-m telescope at the McDonald observatory. We apply the binary modelling software, PHOEBE, to the Kepler light curve and radial velocity data to determine a detailed binary star model that includes the prominent pulsation and Doppler boosting, alongside the usual attributes of a binary star model (including tidal distortion and reflection). The results show that the system contains a slightly evolved F star with an M secondary companion in a highly eccentric orbit (e = 0.886). We use the results of the binary star model in a companion paper (Fuller et al., 2017) where we show that the prominent pulsation can be explained by a tidally excited oscillation mode held near resonance by a resonance locking mechanism.
• 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 eclipsing binary and multiple systems. Supervised classification and self-organizing maps(1702.06296)

Feb. 21, 2017 astro-ph.SR, astro-ph.IM
Large surveys producing tera- and petabyte-scale databases require machine-learning and knowledge discovery methods to deal with the overwhelming quantity of data and the difficulties of extracting concise, meaningful information with reliable assessment of its uncertainty. This study investigates the potential of a few machine-learning methods for the automated analysis of eclipsing binaries in the data of such surveys. We aim to aid the extraction of samples of eclipsing binaries from such databases and to provide basic information about the objects. We estimate class labels according to two classification systems, one based on the light curve morphology (EA/EB/EW classes) and the other based on the physical characteristics of the binary system (system morphology classes; detached through overcontact systems). Furthermore, we explore low-dimensional surfaces along which the light curves of eclipsing binaries are concentrated, to use in the characterization of the binary systems and in the exploration of biases of the full unknown Gaia data with respect to the training sets. We explore the performance of principal component analysis (PCA), linear discriminant analysis (LDA), random forest classification and self-organizing maps (SOM). We pre-process the photometric time series by combining a double Gaussian profile fit and a smoothing spline, in order to de-noise and interpolate the observed light curves. We achieve further denoising, and selected the most important variability elements from the light curves using PCA. We perform supervised classification using random forest and LDA based on the PC decomposition, while SOM gives a continuous 2-dimensional manifold of the light curves arranged by a few important features. We estimate the uncertainty of the supervised methods due to the specific finite training set using ensembles of models constructed on randomized training sets.
• ### 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]
• ### IAU 2015 Resolution B3 on Recommended Nominal Conversion Constants for Selected Solar and Planetary Properties(1510.07674)

Astronomers commonly quote the properties of celestial objects in units of parameters for the Sun, Jupiter, or the Earth. The resolution presented here was proposed by the IAU Inter-Division Working Group on Nominal Units for Stellar and Planetary Astronomy and passed by the XXIXth IAU General Assembly in Honolulu. IAU 2015 Resolution B3 adopts a set of nominal solar, terrestrial, and jovian conversion constants for stellar and (exo)planetary astronomy which are defined to be exact SI values. While the nominal constants are based on current best estimates (CBEs; which have uncertainties, are not secularly constant, and are updated regularly using new observations), they should be interpreted as standard values and not as CBEs. IAU 2015 Resolution B3 adopts five solar conversion constants (nominal solar radius, nominal total solar irradiance, nominal solar luminosity, nominal solar effective temperature, and nominal solar mass parameter) and six planetary conversion constants (nominal terrestrial equatorial radius, nominal terrestrial polar radius, nominal jovian equatorial radius, nominal jovian polar radius, nominal terrestrial mass parameter, and nominal jovian mass parameter).
• ### IAU 2015 Resolution B2 on Recommended Zero Points for the Absolute and Apparent Bolometric Magnitude Scales(1510.06262)

The XXIXth IAU General Assembly in Honolulu adopted IAU 2015 Resolution B2 on recommended zero points for the absolute and apparent bolometric magnitude scales. The resolution was proposed by the IAU Inter-Division A-G Working Group on Nominal Units for Stellar and Planetary Astronomy after consulting with a broad spectrum of researchers from the astronomical community. Resolution B2 resolves the long-standing absence of an internationally-adopted zero point for the absolute and apparent bolometric magnitude scales. Resolution B2 defines the zero point of the absolute bolometric magnitude scale such that a radiation source with $M_{\rm Bol}$ = 0 has luminosity L$_{\circ}$ = 3.0128e28 W. The zero point of the apparent bolometric magnitude scale ($m_{\rm Bol}$ = 0) corresponds to irradiance $f_{\circ}$ = 2.518021002e-8 W/m$^2$. The zero points were chosen so that the nominal solar luminosity (3.828e26 W) adopted by IAU 2015 Resolution B3 corresponds approximately to $M_{\rm Bol}$(Sun) = 4.74, the value most commonly adopted in recent literature. The nominal total solar irradiance (1361 W/m$^2$) adopted in IAU 2015 Resolution B3 corresponds approximately to apparent bolometric magnitude $m_{\rm bol}$(Sun) = -26.832. Implicit in the IAU 2015 Resolution B2 definition of the apparent bolometric magnitude scale is an exact definition for the parsec (648000/$\pi$ au) based on the IAU 2012 Resolution B2 definition of the astronomical unit.
• ### KIC 10080943: An eccentric binary system containing two pressure- and gravity-mode hybrid pulsators(1509.00781)

Sept. 24, 2015 astro-ph.SR
Gamma Doradus and delta Scuti pulsators cover the transition region between low mass and massive main-sequence stars, and as such, are critical for testing stellar models. When they reside in binary systems, we can combine two independent methods to derive critical information, such as precise fundamental parameters to aid asteroseismic modelling. In the Kepler light curve of KIC10080943, clear signatures of gravity- and pressure-mode pulsations have been found. Ground-based spectroscopy revealed this target to be a double-lined binary system. We present the analysis of four years of Kepler photometry and high-resolution spectroscopy to derive observational constraints with which to evaluate theoretical predictions of the stellar structure and evolution for intermediate-mass stars. We used the method of spectral disentangling to determine atmospheric parameters for both components and derive the orbital elements. With PHOEBE, we modelled the ellipsoidal variation and reflection signal of the binary in the light curve and used classical Fourier techniques to analyse the pulsation modes. We show that the eccentric binary system KIC10080943 contains two hybrid pulsators with masses $M_1=2.0\pm0.1~M_\odot$ and $M_2=1.9\pm0.1~M_\odot$, with radii $R_1=2.9\pm0.1~R_\odot$ and $R_2=2.1\pm0.2~R_\odot$. We detect rotational splitting in the g modes and p modes for both stars and use them to determine a first rough estimate of the core-to-surface rotation rates for the two components, which will be improved by future detailed seismic modelling.
• ### Discovery of Two New Thermally Bloated Low-Mass White Dwarfs Among the Kepler Binaries(1502.02303)

Feb. 12, 2015 astro-ph.SR
We report the discovery of two new low-mass, thermally bloated, hot white dwarfs among the Kepler sample of eclipsing binaries. These are KIC 9164561 and KIC 10727668 with orbital periods of 1.2670 and 2.3058 days, respectively. The current primary in both systems is an A star of about 2 Msun. This brings the number of similar binaries among the Kepler sample to six, and the two new systems have the shortest orbital periods among them. The white dwarf in KIC 9164561 has the largest thermal bloating, compared to its cold degenerate radius, of about a factor of 14. We utilize RV measurements of the A star in KIC 9164561 to determine the white dwarf mass rather accurately: 0.197 +/- 0.005 Msun. The mass of the white dwarf in KIC 10727668 is based on the Doppler boosting signal in the Kepler photometry, and is less accurately determined to be 0.266 +/- 0.035 Msun. Based on the inferred radii and effective temperatures of these two white dwarfs we are able to make an independent theoretical estimate of their masses to within ~0.01 Msun based on evolutionary models of their cooling history after they lose their hydrogen-rich envelopes. We also present evidence that there is a third body in the KIC 9164561 system with an orbital period of 8-14 years.
• ### HD 183648: a Kepler eclipsing binary with anomalous ellipsoidal variations and a pulsating component(1407.1788)

July 7, 2014 astro-ph.SR
KIC 8560861 (HD 183648) is a marginally eccentric (e=0.05) eclipsing binary with an orbital period of P_orb=31.973d, exhibiting mmag amplitude pulsations on time scales of a few days. We present the results of the complex analysis of high and medium-resolution spectroscopic data and Kepler Q0 -- Q16 long cadence photometry. The iterative combination of spectral disentangling, atmospheric analysis, radial velocity and eclipse timing variation studies, separation of pulsational features of the light curve, and binary light curve analysis led to the accurate determination of the fundamental stellar parameters. We found that the binary is composed of two main sequence stars with an age of 0.9\+-0.2 Gyr, having masses, radii and temperatures of M_1=1.93+-0.12 M_sun, R_1=3.30+-0.07 R_sun, T_eff1=7650+-100 K for the primary, and M_2=1.06+-0.08 M_sun, R_2=1.11+-0.03 R_sun, T_eff2=6450+-100 K for the secondary. After subtracting the binary model, we found three independent frequencies, two of which are separated by twice the orbital frequency. We also found an enigmatic half orbital period sinusoidal variation that we attribute to an anomalous ellipsoidal effect. Both of these observations indicate that tidal effects are strongly influencing the luminosity variations of HD 183648. The analysis of the eclipse timing variations revealed both a parabolic trend, and apsidal motion with a period of (P_apse)_obs=10,400+-3,000 y, which is three times faster than what is theoretically expected. These findings might indicate the presence of a distant, unseen companion.
• ### Stellar statistics along the ecliptic and the impact on the K2 mission concept(1403.5888)

March 24, 2014 astro-ph.SR
K2 is the mission concept for a repurposed Kepler mission that uses two reaction wheels to maintain the satellite attitude and provide ~85 days of coverage for ten 105-deg^2 fields along the ecliptic in the first 2.5 years of operation. We examine stellar populations based on the improved Besancon model of the Galaxy, comment on the general properties for the entire ecliptic plane, and provide stellar occurrence rates in the first 6 tentative K2 campaigns grouped by spectral type and luminosity class. For each campaign we distinguish between main sequence stars and giants and provide their density profile as a function of galactic latitude. We introduce the crowding metric that serves for optimized target selection across the campaigns. For all main sequence stars we compute the expected planetary occurrence rates for three planet sizes: 2-4 REarth, 4-8 REarth, and 8-32 REarth with orbital periods up to 50 days. In conjunction with Gaia and the upcoming TESS and Plato missions, K2 will become a gold mine for stellar and planetary astrophysics.
• ### Pulsating red giant stars in eccentric binary systems discovered from Kepler space-based photometry(1312.4500)

Dec. 16, 2013 astro-ph.SR
The unparalleled photometric data obtained by NASA's Kepler space telescope led to an improved understanding of red giant stars and binary stars. Seismology allows us to constrain the properties of red giants. In addition to eclipsing binaries, eccentric non-eclipsing binaries, exhibiting ellipsoidal modulations, have been detected with Kepler. We aim to study the properties of eccentric binary systems containing a red giant star and derive the parameters of the primary giant component. We apply asteroseismic techniques to determine masses and radii of the primary component of each system. For a selected target, light and radial velocity curve modelling techniques are applied to extract the parameters of the system. The effects of stellar on the binary system are studied. The paper presents the asteroseismic analysis of 18 pulsating red giants in eccentric binary systems, for which masses and radii were constrained. The orbital periods of these systems range from 20 to 440days. From radial velocity measurements we find eccentricities between e=0.2 to 0.76. As a case study we present a detailed analysis of KIC5006817. From seismology we constrain the rotational period of the envelope to be at least 165 d, roughly twice the orbital period. The stellar core rotates 13 times faster than the surface. From the spectrum and radial velocities we expect that the Doppler beaming signal should have a maximum amplitude of 300ppm in the light curve. Through binary modelling, we determine the mass of the secondary component to be 0.29$\pm$0.03\,$M_\odot$. For KIC5006817 we exclude pseudo-synchronous rotation of the red giant with the orbit. The comparison of the results from seismology and modelling of the light curve shows a possible alignment of the rotational and orbital axis at the 2$\sigma$ level. Red giant eccentric systems could be progenitors of cataclysmic variables and hot subdwarf B stars.
• ### KIC 4544587: an Eccentric, Short Period Binary System with delta Sct Pulsations and Tidally Excited Modes(1306.1819)

June 7, 2013 astro-ph.SR
We present Kepler photometry and ground based spectroscopy of KIC 4544587, a short-period eccentric eclipsing binary system with self-excited pressure and gravity modes, tidally excited modes, tidally influenced p modes, and rapid apsidal motion of 182 y per cycle. The primary and secondary components of KIC 4544587 reside within the delta Scuti and gamma Dor instability region of the Hurtzsprung-Russell diagram, respectively. By applying the binary modelling software PHOEBE to prewhitened Kepler photometric data and radial velocity data obtained using the William Herschel Telescope and 4-m Mayall telescope at KPNO, the fundamental parameters of this important system have been determined, including the stellar masses, 1.98+/-0.07 MSun and 1.60+/-0.06 MSun, and radii, 1.76+/-0.03 RSun and 1.42+/-0.02 RSun, for the primary and secondary components, respectively. Frequency analysis of the residual data revealed 31 modes, 14 in the gravity mode region and 17 in the pressure mode region. Of the 14 gravity modes 8 are orbital harmonics: a signature of tidal resonance. While the measured amplitude of these modes may be partially attributed to residual signal from binary model subtraction, we demonstrate through consideration of the folded light curve that these frequencies do in fact correspond to tidally excited pulsations. Furthermore, we present an echelle diagram of the pressure mode frequency region (modulo the orbital frequency) and demonstrate that the tides are also influencing the p modes. A first look at asteroseismology hints that the secondary component is responsible for the p modes, which is contrary to our expectation that the hotter star should pulsate in higher radial overtone, higher frequency p modes.
• ### Luminosity Discrepancy in the Equal-Mass, Pre--Main Sequence Eclipsing Binary Par 1802: Non-Coevality or Tidal Heating?(1111.2322)

Nov. 9, 2011 astro-ph.SR
Parenago 1802, a member of the ~1 Myr Orion Nebula Cluster, is a double-lined, detached eclipsing binary in a 4.674 d orbit, with equal-mass components (M_2/M_1 = 0.985 \pm 0.029). Here we present extensive VIcJHKs light curves spanning ~15 yr, as well as a Keck/HIRES optical spectrum. The light curves evince a third light source that is variable with a period of 0.73 d, and is also manifested in the high-resolution spectrum, strongly indicating the presence of a third star in the system, probably a rapidly rotating classical T Tauri star. We incorporate this third light into our radial velocity and light curve modeling of the eclipsing pair, measuring accurate masses (M_1 = 0.391 \pm 0.032, M_2 = 0.385 \pm 0.032 M\odot), radii (R_1 = 1.73 \pm 0.02, R_2 = 1.62 \pm 0.02 R\odot), and temperature ratio (T_1/T_2 = 1.0924 \pm 0.0017). Thus the radii of the eclipsing stars differ by 6.9 \pm 0.8%, the temperatures differ by 9.2 \pm 0.2%, and consequently the luminosities differ by 62 \pm 3%, despite having masses equal to within 3%. This could be indicative of an age difference of ~3x10^5 yr between the two eclipsing stars, perhaps a vestige of the binary formation history. We find that the eclipsing pair is in an orbit that has not yet fully circularized, e = 0.0166 \pm 0.003. In addition, we measure the rotation rate of the eclipsing stars to be 4.629 \pm 0.006 d; they rotate slightly faster than their 4.674 d orbit. The non-zero eccentricity and super-synchronous rotation suggest that the eclipsing pair should be tidally interacting, so we calculate the tidal history of the system according to different tidal evolution theories. We find that tidal heating effects can explain the observed luminosity difference of the eclipsing pair, providing an alternative to the previously suggested age difference.
• ### Kepler photometry of KIC 10661783: a binary star with total eclipses and delta Scuti pulsations(1102.3599)

Feb. 17, 2011 astro-ph.SR
We present Kepler satellite photometry of KIC 10661783, a short-period binary star system which shows total eclipses and multi-periodic delta Scuti pulsations. A frequency analysis of the eclipse-subtracted light curve reveals at least 68 frequencies of which 55 or more can be attributed to pulsation modes. The main limitation on this analysis is the frequency resolution within the 27-day short-cadence light curve. Most of the variability signal lies in the frequency range 18 to 31 c/d, with amplitudes between 0.1 and 4 mmag. One harmonic term (2.f) and a few combination frequencies (f_i+f_j) have been detected. From a plot of the residuals versus orbital phase we assign the pulsations to the primary star in the system. The pulsations were removed from the short-cadence data and the light curve was modelled using the Wilson-Devinney code. We are unable to get a perfect fit due to the residual effects of pulsations and also to the treatment of reflection and reprocessing in the light curve model. A model where the secondary star fills its Roche lobe is favoured, which means that KIC 10661783 can be classified as an oEA system. Further photometric and spectroscopic observations will allow the masses and radii of the two stars to be measured to high precision and hundreds of delta Scuti pulsation frequencies to be resolved. This could lead to unique constraints on theoretical models of delta Scuti stars, if the evolutionary history of KIC 10661783 can be accounted for.
• ### Pre-Spectroscopic False Positive Elimination of Kepler Planet Candidates(1001.0392)

Jan. 3, 2010 astro-ph.EP
Ten days of commissioning data (Quarter 0) and thirty-three days of science data (Quarter 1) yield instrumental flux timeseries of ~150,000 stars that were combed for transit events, termed Threshold Crossing Events (TCE), each having a total detection statistic above 7.1-sigma. TCE light curves are modeled as star+planet systems. Those returning a companion radius smaller than 2R_J are assigned a KOI (Kepler Object of Interest) number. The raw flux, pixel flux, and flux-weighted centroids of every KOI are scrutinized to assess the likelihood of being an astrophysical false-positive versus the likelihood of a being a planetary companion. This vetting using Kepler data is referred to as data validation. Herein, we describe the data validation metrics and graphics used to identify viable planet candidates amongst the KOIs. Light curve modeling tests for a) the difference in depth of the odd- versus even-numbered transits, b) evidence of ellipsoidal variations, and c) evidence of a secondary eclipse event at phase=0.5. Flux-weighted centroids are used to test for signals correlated with transit events with a magnitude and direction indicative of a background eclipsing binary. Centroid timeseries are complimented by analysis of images taken in-transit versus out-of-transit, the difference often revealing the pixel contributing the most to the flux change during transit. Examples are shown to illustrate each test. Candidates passing data validation are submitted to ground-based observers for further false-positive elimination or confirmation/characterization.
• ### HD 174884: a strongly eccentric, short-period early-type binary system discovered by CoRoT(0910.3513)

Oct. 19, 2009 astro-ph.SR
Accurate photometric CoRoT space observations of a secondary seismological target, HD 174884, led to the discovery that this star is an astrophysically important double-lined eclipsing spectroscopic binary in an eccentric orbit (e of about 0.3), unusual for its short (3.65705d) orbital period. The high eccentricity, coupled with the orientation of the binary orbit in space, explains the very unusual observed light curve with strongly unequal primary and secondary eclipses having the depth ratio of 1-to-100 in the CoRoT 'seismo' passband. Without the high accuracy of the CoRoT photometry, the secondary eclipse, 1.5 mmag deep, would have gone unnoticed. A spectroscopic follow-up program provided 45 high dispersion spectra. The analysis of the CoRoT light curve was performed with an adapted version of PHOEBE that supports CoRoT passbands. The final solution was obtained by simultaneous fitting of the light and the radial velocity curves. Individual star spectra were derived by spectrum disentangling. The uncertainties of the fit were derived by bootstrap resampling and the solution uniqueness was tested by heuristic scanning. The results provide a consistent picture of the system composed of two late B stars. The Fourier analysis of the light curve fit residuals yields two components, with orbital frequency multiples and an amplitude of about 0.1 mmag, which are tentatively interpreted as tidally induced pulsations. An extensive comparison with theoretical models is carried out by means of the Levenberg-Marquardt minimization technique and the discrepancy between models and the derived parameters is discussed. The best fitting models yield a young system age of 125 million years which is consistent with the eccentric orbit and synchronous component rotation at periastron.
• ### CoRoT photometry and high-resolution spectroscopy of the interacting eclipsing binary AU Mon(0909.1546)

Sept. 8, 2009 astro-ph.SR
Analyses of very accurate CoRoT space photometry, past Johnson V photoelectric photometry and high-resolution \'echelle spectra led to the determination of improved and consistent fundamental stellar properties of both components of AU Mon. We derived new, accurate ephemerides for both the orbital motion (with a period of 11.113d) and the long-term, overall brightness variation (with a period of 416.9d) of this strongly interacting Be + G semi-detached binary. It is shown that this long-term variation must be due to attenuation of the total light by some variable circumbinary material. We derived the binary mass ratio $M_{\rm G}/M_{\rm B}$ = 0.17\p0.03 based on the assumption that the G-type secondary fills its Roche lobe and rotates synchronously. Using this value of the mass ratio as well as the radial velocities of the G-star, we obtained a consistent light curve model and improved estimates of the stellar masses, radii, luminosities and effective temperatures. We demonstrate that the observed lines of the B-type primary may not be of photospheric origin. We also discover rapid and periodic light changes visible in the high-quality residual CoRoT light curves. AU Mon is put into perspective by a comparison with known binaries exhibiting long-term cyclic light changes.
• ### Fully Automated Approaches to Analyze Large-Scale Astronomy Survey Data(0904.0739)

Observational astronomy has changed drastically in the last decade: manually driven target-by-target instruments have been replaced by fully automated robotic telescopes. Data acquisition methods have advanced to the point that terabytes of data are flowing in and being stored on a daily basis. At the same time, the vast majority of analysis tools in stellar astrophysics still rely on manual expert interaction. To bridge this gap, we foresee that the next decade will witness a fundamental shift in the approaches to data analysis: case-by-case methods will be replaced by fully automated pipelines that will process the data from their reduction stage, through analysis, to storage. While major effort has been invested in data reduction automation, automated data analysis has mostly been neglected despite the urgent need. Scientific data mining will face serious challenges to identify, understand and eliminate the sources of systematic errors that will arise from this automation. As a special case, we present an artificial intelligence (AI) driven pipeline that is prototyped in the domain of stellar astrophysics (eclipsing binaries in particular), current results and the challenges still ahead.
• ### In search of RR Lyrae type stars in eclipsing binary systems. OGLE052218.07-692827.4: an optical blend(0808.3560)

Aug. 26, 2008 astro-ph
During the OGLE-2 operation, Soszynski et al. (2003) found 3 LMC candidates for an RR Lyr-type component in an eclipsing binary system. Two of those have orbital periods that are too short to be physically plausible and hence have to be optical blends. For the third, OGLE052218.07-692827.4, we developed a model of the binary that could host the observed RR Lyr star. After being granted HST/WFPC2 time, however, we were able to resolve 5 distinct sources within a 1.3" region that is typical of OGLE resolution, proving that OGLE052218.07-692827.4 is also an optical blend. Moreover, the putative eclipsing binary signature found in the OGLE data does not seem to correspond to a physically plausible system; the source is likely another background RR Lyr star. There are still no RR Lyr stars discovered so far in an eclipsing binary system.
• ### Artificial Intelligence Approach to the Determination of Physical Properties of Eclipsing Binaries. I. The EBAI Project(0807.1724)

July 10, 2008 astro-ph
Achieving maximum scientific results from the overwhelming volume of astronomical data to be acquired over the next few decades will demand novel, fully automatic methods of data analysis. Artificial intelligence approaches hold great promise in contributing to this goal. Here we apply neural network learning technology to the specific domain of eclipsing binary (EB) stars, of which only some hundreds have been rigorously analyzed, but whose numbers will reach millions in a decade. Well-analyzed EBs are a prime source of astrophysical information whose growth rate is at present limited by the need for human interaction with each EB data-set, principally in determining a starting solution for subsequent rigorous analysis. We describe the artificial neural network (ANN) approach which is able to surmount this human bottleneck and permit EB-based astrophysical information to keep pace with future data rates. The ANN, following training on a sample of 33,235 model light curves, outputs a set of approximate model parameters (T2/T1, (R1+R2)/a, e sin(omega), e cos(omega), and sin i) for each input light curve data-set. The whole sample is processed in just a few seconds on a single 2GHz CPU. The obtained parameters can then be readily passed to sophisticated modeling engines. We also describe a novel method polyfit for pre-processing observational light curves before inputting their data to the ANN and present the results and analysis of testing the approach on synthetic data and on real data including fifty binaries from the Catalog and Atlas of Eclipsing Binaries (CALEB) database and 2580 light curves from OGLE survey data. [abridged]
• ### Introducing Powell's Direction Set Method to a Fully Automated Analysis of Eclipsing Binary Stars(astro-ph/0607340)

July 14, 2006 astro-ph
With recent observational advancements, substantial amounts of photometric and spectroscopic eclipsing binary data have been acquired. As part of an ongoing effort to assemble a reliable pipeline for fully automatic data analysis, we put Powell's direction set method to the test. The method does not depend on numerical derivatives, only on function evaluations, and as such it cannot diverge. Compared to differential corrections (DC) and Nelder & Mead's downhill simplex (NMS) method, Powell's method proves to be more efficient in terms of solution determination and the required number of iterations. However, its application is still not optimal in terms of time cost. Causes for this deficiency are identified and two steps toward the solution are proposed: non-ortogonality of the parameter set should be removed and better initial directions should be determined before the minimization is initiated. Once these setbacks are worked out, Powell's method will probably replace DC and NMS as the default minimizing algorithm in PHOEBE modeling package.