Preparing for the expected wealth of Gaia detections, we consider here a simple algorithm for classifying unresolved astrometric binaries with main-sequence (MS) primary into three classes: binaries with a probable MS secondary, with two possible values for the mass ratio; probable hierarchical triple MS systems with an astrometric secondary as a close binary, with a limited range of mass-ratio values; and binaries with a compact-object secondary, with a minimal value of the mass ratio. This is done by defining a unit-less observational parameter 'Astrometric Mass-Ratio Function' (AMRF), $\mathcal{A}$, of a binary, based on primary-mass estimation, in addition to the astrometric parameters - the angular semi-major axis, the period and the parallax. We derive the $\mathcal{A}$ value that differentiates the three classes by forward modeling representative binaries of each class, assuming some mass-luminosity relation.To demonstrate the potential of the algorithm, we consider the orbits of 98 Hipparcos astrometric binaries with main-sequence primaries, using the Hipparcos parallaxes and the primary-mass estimates. For systems with known spectroscopic orbital solution, our results are consistent with the spectroscopic elements, validating the suggested approach. The algorithm will be able to identify hierarchical triple systems and dormant neutron-star and black-hole companions in the Gaia astrometric binaries.

We report here on a detailed study of a c-type RR Lyrae variable (RRc variable), SB 796, serendipitously discovered in a search of the WASP public data for stars that display large photometric periodic modulation. SB 796 displays a period of $P = 0.26585$ d and semi-amplitude of $\sim$ 0.1 mag. Comparison of the modulation shape and period with the detailed analysis of LMC variables indicates that SB 796 is an RRc variable. Gaia DR2 classification corroborated our result. Radial-velocity (RV) follow-up observations revealed a periodic variation consistent with a sine modulation, with a semi-amplitude of $5.6\pm0.2$ km/s, and a minimum at phase of maximum brightness. Similar amplitude and phase were previously seen in other RRc variables. The stellar averaged RV is $\sim 250$ km/s, turning SB 796 to be a high-velocity star, while its present position, as derived from the Gaia astrometry, is at $\sim 3.5$ kpc below the Galactic plane. Integration of the stellar Galactic motion shows that SB 796 oscillates at a range of 0.5--20 kpc Galacto-centric distance, passing near the Galactic center about three times in 1 Gyr. The Galactic radial motion takes SB 796 up and down the plane to a scale height of $\sim 10$ kpc. During its $\sim10$ Gyrs estimated life time, SB 796 therefore passed $\sim 30$ times near the Galactic center.

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.

We developed a novel direct algorithm to derive the mass-ratio distribution (MRD) of short-period binaries from an observed sample of single-lined spectroscopic binaries (SB1). The algorithm considers a class of parameterized MRDs and finds the set of parameters that best fits the observed sample. The algorithm consists of four parts. First, we define a new observable, the `modified mass function', that can be calculated for each binary in the sample. We show that the distribution of the modified mass function follows the shape of the underlying MRD, turning it more advantageous than the previously used mass function, reduced mass function or reduced mass function logarithm. Second, we derive the likelihood of the sample of modified mass functions to be observed given an assumed MRD. An MCMC search enables the algorithm to find the parameters that best fit the observations. Third, we suggest to express the unknown MRD by a linear combination of a basis of functions that spans the possible MRDs. We suggest two such bases. Fourth, we show how to account for the undetected systems that have an RV amplitude below a certain threshold. Without the correction, this observational bias suppresses the derived MRD for low mass ratios. Numerous simulations show that the algorithm works well with either of the two suggested bases. The four parts of the algorithm are independent, but the combination of the four turn the algorithm to be highly effective in deriving the MRD of the binary population.

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.

This thesis consists of seven scientific papers that cover the proof-of-concept, the development, and discoveries made through the use of the BEER (BEaming, Ellipsoidal, and Reflection) algorithm for searching for companions in the light curves from the Kepler and CoRoT space telescopes. Paper I presents the detection of the ellipsoidal and the beaming effects in the CoRoT light curve of CoRoT-3, a system of a $22$$M_{\rm Jup}$ brown dwarf orbiting an F star with an orbital period of $4.3$ days. This work served as a proof-of-concept that these effects are detectable in the space light curves of systems with brown-dwarf or planetary secondaries, thus indicating that similar modulations may be detected in the light curves of non-transiting systems. Abridged... The last study, Paper VII, demonstrates the different strengths and utility of the BEER search algorithm. It presents the discovery of four short-period eclipsing binaries in the Kepler light curves, consisting of an A-star primary and a low-mass WD secondary (dA+WD). The systems show BEER phase modulations together with primary and secondary eclipses. These add to the 6 Kepler, and 18 WASP, previously known short-period eclipsing dA+WD binaries. The paper shows that three of the new systems harbor the smallest WDs detected so far in such binaries. These three binaries extend the previously known population to older systems with cooler and smaller WD secondaries, allowing to test binary evolution theories in a parameter region not observed before. The seven papers illustrate the effectiveness of the BEER algorithm in finding both common stellar binaries and rare astrophysical objects. As such, the BEER tool can serve as an important component in the virtual astronomy toolbox for mining the vast astronomical data produced by current and future photometric surveys. Abridged.

A few studies have reported a significant dearth of exoplanets with Neptune mass and radius with orbital periods below $2$--$4$ d. This cannot be explained by observational biases because many Neptunian planets with longer orbital periods have been detected. The existence of this desert is similar to the appearance of the so-called brown-dwarf desert that suggests different formation mechanisms of planets and stellar companions with short orbital periods. Similarly, the Neptunian desert might indicate different mechanisms of formation and evolution for hot Jupiters and short-period super-Earths. We here follow a previous study and examine the location and shape of the desert in both the period-mass and period-radius planes, using the currently available large samples of planets. The desert in the period-mass plane has a relatively sharp upper edge, with a planetary mass that is inversely proportional to the planetary orbital period, while the lower, somewhat blurred, boundary is located along masses that are apparently linearly proportional to the period. The desert in the period-radius plane of the transiting planets is less clear. It seems as if the radius along the upper boundary is inversely proportional to the period to the power of one-third, while the lower boundary shows a radius that is proportional to the period to the power of two-thirds. The combination of the two upper bounds of the desert, in the period-mass and period-radius planes, yields a planetary mass-radius relation of $R_{\rm p}/R_{\rm Jup}\simeq (1.2\pm0.3)(M_{\rm p}/M_{\rm Jup})^{0.27\pm0.11}\,$ for $ 0.1\lesssim M_{\rm p}/M_{\rm Jup}\lesssim 1$. The derived shape of the desert, which might extend up to periods of $5$--$10$ d, could shed some light on the formation and evolution of close-in planets.

We report the discovery of the sixth known eclipsing double white dwarf (WD) system, SDSS J1152+0248, with a 2.3968 +/- 0.0003 h orbital period, in data from the Kepler Mission's K2 continuation. Analysing and modelling the K2 data together with ground-based fast photometry, spectroscopy, and radial-velocity measurements, we determine that the primary is a DA-type WD with mass M1 = 0.47 +/- 0.11 Msun, radius R1 = 0.0197 +/- 0.0035 Rsun, and cooling age t1 = 52 +/- 36 Myr. No lines are detected, to within our sensitivity, from the secondary WD, but it is likely also of type DA. Its central surface brightness, as measured from the secondary eclipse, is 0.31 of the primary's surface brightness. Its mass, radius, and cooling age, respectively, are M2 = 0.44 +/- 0.09 Msun, R2 = 0.0223 +0.0064 -0.0050 Rsun, and t2 = 230 +/- 100 Myr. SDSS J1152+0248 is a near twin of the double-lined eclipsing WD system CSS 41177.

We report the discovery of four short-period eclipsing systems in the Kepler light curves, consisting of an A-star primary and a low-mass white dwarf (WD) secondary (dA+WD) - KIC 4169521, KOI-3818, KIC 2851474, and KIC 9285587. The systems show BEaming, Ellipsoidal and Reflection (BEER) phase modulations together with primary and secondary eclipses. These add to the 6 Kepler and 18 WASP short-period eclipsing dA+WD binaries that were previously known. The light curves, together with follow-up spectroscopic observations, allow us to derive the masses, radii, and effective temperatures of the two components of the four systems. The orbital periods, of 1.17-3.82 days, and WD masses, of 0.19-0.22 Msun, are similar to those of the previously known systems. The WD radii of KOI-3818, KIC 2851474, and KIC 9285587 are 0.026, 0.035, and 0.026 Rsun, respectively, the smallest WD radii derived so far for short-period eclipsing dA+WD binaries. These three binaries extend the previously known population to older systems with cooler and smaller WD secondaries. KOI-3818 displays evidence for a fast-rotating primary and a minute but significant eccentricity of ~0.0015. These features are probably the outcome of the mass-transfer process.

The CoRoT space mission was operating for almost 6 years, producing thousands of continuous photometric light curves. The temporal series of exposures are processed by the production pipeline, correcting the data for known instrumental effects. But even after these model-based corrections, some collective trends are still visible in the light curves. We propose here a simple exposure-based algorithm to remove instrumental effects. The effect of each exposure is a function of only two instrumental stellar parameters, position on the CCD and photometric aperture. The effect is not a function of the stellar flux, and therefore much more robust. As an example, we show that the $\sim2\%$ long-term variation of the early run LRc01 is nicely detrended on average. This systematics removal process is part of the CoRoT legacy data pipeline.

(abridged for arXiv) The BEER algorithm searches stellar light curves for the BEaming, Ellipsoidal, and Reflection photometric modulations that are caused by a short-period companion. Applying the search to the first five long-run center CoRoT fields, we identified $481$ non-eclipsing candidates with periodic flux amplitudes of $0.5-87$ mmag. Optimizing the Anglo-Australian-Telescope pointing coordinates and the AAOmega fiber-allocations with dedicated softwares, we acquired six spectra for $231$ candidates and seven spectra for another $50$ candidates in a seven-night campaign. Analysis of the red-arm AAOmega spectra, which covered the range of $8342-8842\AA{}$, yielded a radial-velocity precision of $\sim1$ km/s. Spectra containing lines of more than one star were analyzed with the two-dimensional correlation algorithm TODCOR. The measured radial velocities confirmed the binarity of seventy of the BEER candidates$-45$ single-line binaries, $18$ double-line binaries, and $7$ diluted binaries. We show that red giants introduce a major source of false candidates and demonstrate a way to improve BEER's performance in extracting higher fidelity samples from future searches of CoRoT light curves. The periods of the confirmed binaries span a range of $0.3-10$ days and show a rise in the number of binaries per $\Delta$log$P$ toward longer periods. The estimated mass ratios of the double-line binaries and the mass ratios assigned to the single-line binaries, assuming an isotropic inclination distribution, span a range of $0.03-1$. On the low-mass end, we have detected two brown-dwarf candidates on a $\sim1$ day period orbit. This is the first time non-eclipsing beaming binaries are detected in CoRoT data, and we estimate that $\sim300$ such binaries can be detected in the CoRoT long-run light curves.

We applied the BEER algorithm to the CoRoT lightcurves from the first five LRc fields and identified $481$ non-eclipsing BEER candidates with periodic lightcurve modulations and amplitudes of $0.5-87$ mmag. Medium-resolution spectra of $281$ candidates were obtained in a seven-night AAOmega radial-velocity (RV) campaign, with a precision of $\sim1$ km/s. The RVs confirmed the binarity of $70$ of the BEER candidates, with periods of $0.3-10$ days.

We analyzed the Kepler light curves of four transiting hot Jupiter systems --- KOI-13, HAT-P-7, TrES-2, and Kepler-76, which show BEaming, Ellipsoidal and Reflection (BEER) phase modulations. The mass of the four planets can be estimated from either the beaming or the ellipsoidal amplitude, given the mass and radius of their parent stars. For KOI-13, HAT-P-7, and Kepler-76 we find that the beaming-based planetary mass estimate is larger than the mass estimated from the ellipsoidal amplitude, consistent with previous studies. This apparent discrepancy may be explained by equatorial superrotation of the planet atmosphere, which induces an angle shift of the planet reflection/emission phase modulation, as was suggested for Kepler-76 in the first paper of this series. We propose a modified BEER model that supports superrotation, assuming either a Lambertian or geometric reflection/emission phase function, and provides a photometry-consistent estimate of the planetary mass. Our analysis shows that for Kepler-76 and HAT-P-7, the Lambertian superrotation BEER model is highly preferable over an unshifted null model, while for KOI-13 it is preferable only at a 1.4 sigma level. For TrES-2 we do not find such preference. For all four systems the Lambertian superrotation model mass estimates are in excellent agreement with the planetary masses derived from, or constrained by, radial velocity measurements. This makes the Lambertian superrotation BEER model a viable tool for estimating the masses of hot Jupiters from photometry alone. We conclude that hot Jupiter superrotation may be a common phenomenon that can be detected in the visual light curves of Kepler.

We present the first case in which the BEER algorithm identified a hot Jupiter in the Kepler light curve, and its reality was confirmed by orbital solutions based on follow-up spectroscopy. The companion Kepler-76b was identified by the BEER algorithm, which detected the BEaming (sometimes called Doppler boosting) effect together with the Ellipsoidal and Reflection/emission modulations (BEER), at an orbital period of 1.54 days, suggesting a planetary companion orbiting the 13.3 mag F star. Further investigation revealed that this star appeared in the Kepler eclipsing binary catalog with estimated primary and secondary eclipse depths of 5e-3 and 1e-4 respectively. Spectroscopic radial-velocity follow-up observations with TRES and SOPHIE confirmed Kepler-76b as a transiting 2.0+/-0.26 Mjup hot Jupiter. The mass of a transiting planet can be estimated from either the beaming or the ellipsoidal amplitude. The ellipsoidal-based mass estimate of Kepler-76b is consistent with the spectroscopically measured mass while the beaming-based estimate is significantly inflated. We explain this apparent discrepancy as evidence for the superrotation phenomenon, which involves eastward displacement of the hottest atmospheric spot of a tidally-locked planet by an equatorial super-rotating jet stream. This phenomenon was previously observed only for HD 189733b in the infrared. We show that a phase shift of 10.3+/-2.0 degrees of the planet reflection/emission modulation, due to superrotation, explains the apparently inflated beaming modulation, resolving the ellipsoidal/beaming amplitude discrepancy. Kepler-76b is one of very few confirmed planets in the Kepler light curves that show BEER modulations and the first to show superrotation evidence in the Kepler band. Its discovery illustrates for the first time the ability of the BEER algorithm to detect short-period planets and brown dwarfs.

We present the study of the CoRoT transiting planet candidate 101186644, also named LRc01_E1_4780. Analysis of the CoRoT lightcurve and the HARPS spectroscopic follow-up observations of this faint (m_V = 16) candidate revealed an eclipsing binary composed of a late F-type primary (T_eff = 6090 +/- 200 K) and a low-mass, dense late M-dwarf secondary on an eccentric (e = 0.4) orbit with a period of ~20.7 days. The M-dwarf has a mass of 0.096 +/- 0.011 M_Sun, and a radius of 0.104 +0.026/-0.006 R_Sun, which possibly makes it the smallest and densest late M-dwarf reported so far. Unlike the claim that theoretical models predict radii that are 5%-15% smaller than measured for low-mass stars, this one seems to have a radius that is consistent and might even be below the radius predicted by theoretical models.

KOI-13 was presented by the Kepler team as a candidate for having a giant planet - KOI-13.01, with orbital period of 1.7 d and transit depth of ~0.8%. We have analyzed the Kepler Q2 data of KOI-13, which was publicly available at the time of the submission of this paper, and derived the amplitudes of the beaming, ellipsoidal and reflection modulations: 8.6 +/- 1.1, 66.8 +/- 1.6 and 72.0 +/- 1.5 ppm (parts per million), respectively. After the paper was submitted, Q3 data were released, so we repeated the analysis with the newly available light curve. The results of the two quarters were quite similar. From the amplitude of the beaming modulation we derived a mass of 10 +/- 2 M_Jup for the secondary, suggesting that KOI-13.01 was a massive planet, with one of the largest known radii. We also found in the data a periodicity of unknown origin with a period of 1.0595 d and a peak-to-peak modulation of ~60 ppm. The light curve of Q3 revealed a few more small-amplitude periodicities with similar frequencies. It seemed as if the secondary occultation of KOI-13 was slightly deeper than the reflection peak-to-peak modulation by 16.8 +/- 4.5 ppm. If real, this small difference was a measure of the thermal emission from the night side of KOI-13.01. We estimated the effective temperature to be 2600 +/- 150 K, using a simplistic black-body emissivity approximation. We then derived the planetary geometrical and Bond albedos as a function of the day-side temperature. Our analysis suggested that the Bond albedo of KOI-13.01 might be substantially larger than the geometrical albedo.

We present seven newly discovered non-eclipsing short-period binary systems with low-mass companions, identified by the recently introduced BEER algorithm, applied to the publicly available 138-day photometric light curves obtained by the Kepler mission. The detection is based on the beaming effect (sometimes called Doppler boosting), which increases (decreases) the brightness of any light source approaching (receding from) the observer, enabling a prediction of the stellar Doppler radial-velocity modulation from its precise photometry. The BEER algorithm identifies the BEaming periodic modulation, with a combination of the well known Ellipsoidal and Reflection/heating periodic effects, induced by short-period companions. The seven detections were confirmed by spectroscopic radial-velocity follow-up observations, indicating minimum secondary masses in the range of 0.07-0.4 Msun. The discovered binaries establish for the first time the feasibility of the BEER algorithm as a new detection method for short-period non-eclipsing binaries, with the potential to detect in the near future non-transiting brown dwarfs secondaries, or even massive planets.

We present a simple algorithm, BEER, to search for a combination of the BEaming, Ellipsoidal and the Reflection/heating periodic modulations, induced by short-period non-transiting low-mass companions. The beaming effect is due to the increase (decrease) of the brightness of any light source approaching (receding from) the observer. To first order, the beaming and the reflection/heating effects modulate the stellar brightness at the orbital period, with phases separated by a quarter of a period, whereas the ellipsoidal effect is modulated with the orbital first harmonic. The phase and harmonic differences between the three modulations allow the algorithm to search for a combination of the three effects and identify stellar candidates for low-mass companions. The paper presents the algorithm, including an assignment of a likelihood factor to any possible detection, based on the expected ratio of the beaming and ellipsoidal effects, given an order-of-magnitude estimate of the three effects. As predicted by Loeb & Gaudi (2003) and Zucker, Mazeh & Alexander (2007), the Kepler and the CoRoT lightcurves are precise enough to allow detection of massive planets and brown-dwarf/low-mass-stellar companions with orbital period up to 10-30 days. To demonstrate the feasibility of the algorithm, we bring two examples of candidates found in the first 33 days of the Q1 Kepler lightcurves. Although we used relatively short timespan, the lightcurves were precise enough to enable the detection of periodic effects with amplitudes as small as one part in 10,000 of the stellar flux.

CoRoT-3b is a 22 Jupiter-mass massive-planet/brown-dwarf object, orbiting an F3-star with a period of 4.3 days. We analyzed the out-of-transit CoRoT-3 red-channel lightcurve obtained by the CoRoT mission and detected the ellipsoidal modulation, with half the orbital period and amplitude of 59+/-9 ppm (parts per million) and the relativistic beaming effect, with the orbital period and an amplitude of 27+/-9 ppm. Phases and amplitudes of both modulations were consistent with our theoretical approximation.