EB Science with Extrasolar Planet Surveys

Eclipsing Binary Science Via the Merging of Transit and Doppler Exoplanet Survey Data - A Case Study With the MARVELS Pilot Project and SuperWASP

Scott W. Fleming11affiliation: Department of Astronomy, University of Florida, 211 Bryant Space Science Center, Gainesville, FL 32611, USA , Pierre F. L. Maxted22affiliation: Astrophysics Group, Keele University, Staffordshire, ST5 5BG, UK , Leslie Hebb33affiliation: Department of Physics and Astronomy, Vanderbilt University, Nashville, TN 37235, USA , Keivan G. Stassun33affiliation: Department of Physics and Astronomy, Vanderbilt University, Nashville, TN 37235, USA 44affiliation: Department of Physics, Fisk University, 1000 17th Ave. N., Nashville, TN 37208, USA , Jian Ge11affiliation: Department of Astronomy, University of Florida, 211 Bryant Space Science Center, Gainesville, FL 32611, USA , Phillip A. Cargile33affiliation: Department of Physics and Astronomy, Vanderbilt University, Nashville, TN 37235, USA , Luan Ghezzi55affiliation: Observatório Nacional, Rua General José Cristino, 77, 20921-400 São Cristóvão, Rio de Janeiro, RJ, Brazil 1717affiliation: Laboratório Interinstitucional de e-Astronomia, - LIneA, Rua Gal. José Cristino 77, Rio de Janeiro, RJ - 20921-400, Brazil , Nathan M. De Lee11affiliation: Department of Astronomy, University of Florida, 211 Bryant Space Science Center, Gainesville, FL 32611, USA , John Wisniewski66affiliation: Department of Astronomy, University of Washington, P.O. Box 351580, Seattle, WA 98195, USA , Bruce Gary33affiliation: Department of Physics and Astronomy, Vanderbilt University, Nashville, TN 37235, USA , G. F. Porto de Mello77affiliation: Universidade Federal do Rio de Janeiro, Observatório do Valongo, Ladeira do Pedro Antônio, 43, CEP: 20080-090, Rio de Janeiro, RJ, Brazil 1717affiliation: Laboratório Interinstitucional de e-Astronomia, - LIneA, Rua Gal. José Cristino 77, Rio de Janeiro, RJ - 20921-400, Brazil , Leticia Ferreira77affiliation: Universidade Federal do Rio de Janeiro, Observatório do Valongo, Ladeira do Pedro Antônio, 43, CEP: 20080-090, Rio de Janeiro, RJ, Brazil 1717affiliation: Laboratório Interinstitucional de e-Astronomia, - LIneA, Rua Gal. José Cristino 77, Rio de Janeiro, RJ - 20921-400, Brazil , Bo Zhao11affiliation: Department of Astronomy, University of Florida, 211 Bryant Space Science Center, Gainesville, FL 32611, USA , David R. Anderson22affiliation: Astrophysics Group, Keele University, Staffordshire, ST5 5BG, UK , Xiaoke Wan11affiliation: Department of Astronomy, University of Florida, 211 Bryant Space Science Center, Gainesville, FL 32611, USA , Coel Hellier22affiliation: Astrophysics Group, Keele University, Staffordshire, ST5 5BG, UK , Pengcheng Guo11affiliation: Department of Astronomy, University of Florida, 211 Bryant Space Science Center, Gainesville, FL 32611, USA , Richard G. West88affiliation: Department of Physics and Astronomy, University of Leicester, Leicester, LE1 7RH, UK , Suvrath Mahadevan99affiliation: Department of Astronomy and Astrophysics, The Pennsylvania State University, 525 Davey Laboratory, University Park, PA 16802, USA. 1010affiliation: Center for Exoplanets and Habitable Worlds, The Pennsylvania State University, University Park, PA 16802, USA. , Don Pollacco1111affiliation: Astrophysics Research Centre, School of Mathematics & Physics, Queen’s University, University Road, Belfast, BT7 1NN, UK , Brian Lee11affiliation: Department of Astronomy, University of Florida, 211 Bryant Space Science Center, Gainesville, FL 32611, USA , Andrew Collier Cameron1212affiliation: SUPA, School of Physics & Astronomy, University of St Andrews, North Haugh, KY16 9SS, St Andrews, Fife, Scotland, UK , Julian C. van Eyken11affiliation: Department of Astronomy, University of Florida, 211 Bryant Space Science Center, Gainesville, FL 32611, USA 1313affiliation: NASA Exoplanet Science Institute, Caltech, MS 100-22, 770 South Wilson Avenue, Pasadena, CA 91125, USA , Ian Skillen1414affiliation: Isaac Newton Group of Telescopes, Apartado de correos 321, E-38700 Santa Cruz de la Palma, Canary Islands, Spain , Justin R. Crepp11affiliation: Department of Astronomy, University of Florida, 211 Bryant Space Science Center, Gainesville, FL 32611, USA 1515affiliation: Department of Astronomy, California Institute of Technology, 1200 E. California Blvd., Pasadena, CA 91125, USA , Duy Cuong Nguyen11affiliation: Department of Astronomy, University of Florida, 211 Bryant Space Science Center, Gainesville, FL 32611, USA , Stephen R. Kane11affiliation: Department of Astronomy, University of Florida, 211 Bryant Space Science Center, Gainesville, FL 32611, USA 1313affiliation: NASA Exoplanet Science Institute, Caltech, MS 100-22, 770 South Wilson Avenue, Pasadena, CA 91125, USA , Martin Paegert33affiliation: Department of Physics and Astronomy, Vanderbilt University, Nashville, TN 37235, USA , Luiz Nicolaci da Costa55affiliation: Observatório Nacional, Rua General José Cristino, 77, 20921-400 São Cristóvão, Rio de Janeiro, RJ, Brazil 1717affiliation: Laboratório Interinstitucional de e-Astronomia, - LIneA, Rua Gal. José Cristino 77, Rio de Janeiro, RJ - 20921-400, Brazil , Marcio A. G. Maia55affiliation: Observatório Nacional, Rua General José Cristino, 77, 20921-400 São Cristóvão, Rio de Janeiro, RJ, Brazil 1717affiliation: Laboratório Interinstitucional de e-Astronomia, - LIneA, Rua Gal. José Cristino 77, Rio de Janeiro, RJ - 20921-400, Brazil , Basilio X. Santiago1616affiliation: Instituto de Física, UFRGS, Caixa Postal 15051, Porto Alegre, RS - 91501-970, Brazil 1717affiliation: Laboratório Interinstitucional de e-Astronomia, - LIneA, Rua Gal. José Cristino 77, Rio de Janeiro, RJ - 20921-400, Brazil scfleming@astro.ufl.edu

Exoplanet transit and Doppler surveys discover many binary stars during their operation that can be used to conduct a variety of ancillary science. Specifically, eclipsing binary stars can be used to study the stellar mass-radius relationship and to test predictions of theoretical stellar evolution models. By cross-referencing 24 binary stars found in the MARVELS Pilot Project with SuperWASP photometry, we find two new eclipsing binaries, TYC 0272-00458-1 and TYC 1422-01328-1, which we use as case studies to develop a general approach to eclipsing binaries in survey data. TYC 0272-00458-1 is a single-lined spectroscopic binary for which we calculate a mass of the secondary and radii for both components using reasonable constraints on the primary mass through several different techniques. For a primary mass of , we find , and , and find that both stars have masses and radii consistent with model predictions. TYC 1422-01328-1 is a triple-component system for which we can directly measure the masses and radii of the eclipsing pair. We find that the eclipsing pair consists of an evolved primary star (, ) and a G-type dwarf secondary (, ). We provide the framework necessary to apply this analysis to much larger datasets.

1 Introduction

Photometric and Doppler radial velocity (RV) surveys have generated comprehensive data sets in their quest for extrasolar planets. These surveys have been successful in discovering several hundred planets111http://exoplanet.eu/ and furthering our understanding of how planets form, their dynamically evolution, and their distribution of physical properties. These data sets also provide a wealth of ancillary science studies involving variable stars, binary stellar systems and brown dwarf companions. The amplitudes of the photometric or Doppler variability from such objects are often several orders of magnitude larger than extrasolar planet signals, making them comparatively simpler to identify and characterize.

Cross-referencing stars in both photometric and Doppler surveys can minimize the amount of follow-up data required to conduct these ancillary projects (e.g., Kane et al., 2009). For binary systems, there are several advantages in combining data from both types of surveys. One disadvantage of using spectroscopic binaries (SBs) as a source of detached eclipsing binary (EB) candidates is that the line-of-sight inclination of the orbital plane is unconstrained. The a priori probability that a given SB eclipses can be estimated as , where and are the radii of the primary and secondary, respectively, and is the orbital semimajor axis. Typical a priori values range from a few percent up to %, therefore conducting photometric follow-up results in null detections a majority of the time.

In addition, the predicted ephemerides for eclipses based on modeling the Doppler RV data depend sensitively on the orbital phase coverage, the precise values of the eccentricity and argument of periastron, and the amount of time elapsed since the last measurement. Conversely, archival data from photometric surveys consist of tens of thousands of epochs over several years, resulting in nearly continuous phase coverage for short orbital periods. Given the orbital period from the spectroscopic orbit, it is more efficient to check archival photometry from transit surveys than to attempt follow-up photometric observations based on the transit ephemeris from the radial velocity data.

Attacking the problem from the other direction results in a different set of drawbacks. The advantage of starting with an EB detected from a photometric survey is that the favorable geometry for an eclipse is already known to exist. However, obtaining high-precision, time-series RV follow-up is a resource-intensive enterprise. Such observations are conducted using heavily-subscribed instruments (generally Echelle spectrographs) on moderately large telescopes (usually larger than two meters in diameter) for which access to observing time is quite competitive. Obtaining the best orbital solution nominally requires a minimum of seven epochs that sample the orbital phase well, requiring many nights of follow-up observations.

Doppler RV surveys, on the other hand, already have the necessary number of epochs and sampling rate because an identical observing strategy is required to detect planets. The observing cadences are selected with great care in order to provide synoptic coverage of all possible exoplanet orbital periods and phases. Some surveys cease observations if a star exhibits Doppler variation of many over short timescales, indicative of an SB, to maximize planet yield. Other surveys, such as the multi-object MARVELS survey (Ge et al., 2009), continue to observe the SB to fully characterize the orbit. For the specific case of MARVELS, the multiplicity of a given target is usually not known prior to the start of observations. Continuing to observe SBs to fully characterize their orbit works particularly well in multi-object surveys, where keeping a few SBs per field maximizes science results while not significantly decreasing the capability of detecting extrasolar planets.

Another disadvantage of EBs that only have photometric survey data is that the masses of the two components can only be estimated. Although the mass of the primary can be estimated using colors (e.g., the infrared flux method, Casagrande et al., 2010), deriving the mass of the secondary usually requires an observation of the reflex motion of the primary. Given the large number of EBs present in transit surveys, the ability to select a subset that have the precise masses of interest for the project is more efficient.

2 The Mass-Radius Relationship of K and M Dwarfs

One example of binary science that benefits greatly from the merging of Doppler and photometric survey data is the study of the stellar mass-radius relationship. Observations of EBs have been in excellent agreement with theoretical models with the exception of K and M dwarfs (, López-Morales, 2007; Torres et al., 2010). Stars in this mass range have partially convective outer atmospheres, and can have radii that are 10-20% larger than models predict. Stars above this mass range have primarily radiative atmospheres and agree well with current models (e.g., Popper, 1997; Imbert, 2002; Torres et al., 2008; Meibom et al., 2009, and additional references in Torres et al. (2010)).

One possible cause of this inflation is the stronger magnetic fields for K and M dwarfs in short-period binaries as compared to isolated stars of the same mass (Chabrier et al., 2007). Tidal synchronization increases the rotational period of the stars to match the orbital period, and the stronger magnetic dynamo acts to suppress convection in the outer layers of the star. This results in an increase in cool starspots, reducing the effective temperature, and therefore the radius of the star expands to maintain the luminosity of the star. Recent theoretical work has shown that starspots predominantly distributed at polar latitudes, and/or suppression of convective transport in the outer layers of the atmosphere, are both effects caused by magnetic fields that can reproduce the observations (Morales et al., 2010).

In order to place constraints on the interaction of magnetic fields and the stellar radius, a large sample of EBs with masses and radii measured to precisions of a few percent is required. Specifically, the sample should have at least one component in each EB that is and should have a wide range of orbital periods in order to sample a range of magnetic field strengths. However, despite several decades of surveys for photometric variability, there remains a paucity of such EB systems with precisely determined masses and radii. For example, there were only 13 EB systems in the compilation of López-Morales (2007) that fit the mass criterion above. Combining photometric and Doppler survey data can help increase the sample size by factors of several over the next few years, particularly at longer orbital periods (days). Ground-based surveys for exoplanets identify many SBs and EBs that are, by design, bright stars (V 12). This translates to higher signal-to-noise, requiring less telescope time for follow-up, and therefore a larger number of these binaries can be measured with high precision. The Kepler mission is uncovering thousands of EBs with exquisite photometric precision (Coughlin et al., 2010; Prsǎ et al., 2010), but a majority of those binaries are several magnitudes fainter than EBs found from ground-based exoplanet surveys, resulting in resource-intensive spectroscopic follow-up. In addition, any such ground-based observations are in direct competition with programs dedicated to exoplanet follow-up.

Although several groups are working on studying the mass-radius relationship of double-lined spectroscopic binaries (SB2s) from ground-based surveys (e.g., Devor et al., 2008), less attention is placed on binaries consisting of F/G primaries and K/M secondaries (but see Fernandez et al., 2009). These binaries are generally brighter than K/M pairs, and avoid any biases introduced by studying binaries with mass ratios close to unity. The difficulty with studying binaries that have small flux ratios () is that the mass of the fainter component usually cannot be measured using the survey data. In order to measure the Doppler motion of the fainter companion, large-aperture telescopes equipped with high resolution spectrographs must be used in conjunction with two-dimensional correlation techniques (TODCOR, Zucker & Mazeh, 1994; Zucker et al., 2003). Only a minimum number of epochs per target using these large telescopes are required, because other parameters such as the orbital period, semiamplitude of the primary, and epoch of transit are known from the survey data. Observations can be taken near the RV extrema to measure the semiamplitude of the fainter component and thereby obtain both dynamical masses.

As a case study, we present initial results from a cross-referencing of 24 SBs with days from the MARVELS Pilot Project (hereafter MPP) with photometry from the SuperWASP transit survey (Pollacco et al., 2006). This cross-referencing was performed after the RV observations from the MPP were completed, and the SuperWASP lightcurves were visually examined to identify which of the SBs were eclipsing. Out of those 24 SBs, five did not have a sufficient number of photometric observations and two are blended with other stars. Out of the remaining 17 SBs we have found two EBs: TYC 0272-00458-1 and TYC 1422-01328-1, which we refer to as TYC 272 and TYC 1422 throughout the rest of the paper. TYC 272 is a single-lined spectroscopic binary (SB1), but we constrain the mass of the primary using several different techniques to estimate the properties of the host star. We then calculate the mass of the secondary, as well as the radii of the primary and secondary, from the survey data. TYC 1422 is found to be a triple system with all three components visible in high-resolution spectra. We note that although this work analyzed the photometry and spectroscopic data separately, many software packages exist to simultaneously model all observational data to derive the best-fit binary solutions, (e.g., Prsǎ & Zwitter, 2005). We demonstrate the success of cross-referencing Doppler and transit survey data and lay the foundations for more comprehensive studies using significantly larger data sets.

3 Survey Data

3.1 MPP Doppler Observations

The MPP was a trial survey conducted in 2007 that used the Keck ET instrument (Ge et al., 2006a) on the Sloan Digital Sky Survey (SDSS) 2.5m telescope (Gunn et al., 2006) at the Apache Point Observatory. Although the primary objective was to serve as a test study for the MARVELS extrasolar planet survey, the MPP obtained between 5 and 38 RV measurements of 708 stars over 1-5 month baselines. The 59-object, fiber-fed instrument observed a total of twelve fields, and successfully demonstrated the ability to conduct a multi-object, DFDI (dispersed fixed-delay interferometry, Ge et al., 2002; Ge, 2002; Erskine, 2002; Erskine et al., 2003) survey for extrasolar planets. The spectral resolution of the Keck ET was and had a wavelength range of nm. Despite its relatively small scale, the MPP was also able to demonstrate the ancillary science that can be done with such a survey, including the discovery of a short-period, brown dwarf candidate with minimum mass (Fleming et al., 2010).

3.2 SuperWASP Photometry

The WASP instruments provide flux measurements for millions of stars using wide-angle images of the night sky over a band pass of 400-700 nm defined by a broad-band filter. The survey has produced more than 35 transiting exoplanets to-date, orbiting stars of A through K spectral types and with orbital periods days (e.g., Collier Cameron et al., 2007a, 2010; Hebb et al., 2010; Queloz et al., 2010), in addition to studies of variable stars (Norton et al., 2007; Antipin et al., 2010; Wils, 2010). Eight cameras on each instrument provide images covering approximately 7.8x7.8 degrees using Canon 200mm f/1.8 camera lenses and e2v 2048x2048 CCDs. Synthetic aperture photometry using an aperture radius of 49 arcseconds at the position of cataloged stars is performed on the images (Pollacco et al., 2006). We extracted 1003 observations from the WASP archive for TYC 272 from 10 nights with good quality data covering either the primary or secondary eclipse. For TYC 1422, 649 data points from 12 nights of good quality data were used.

4 Follow-up Observations

4.1 Hereford Arizona Observatory Absolute Photometry

We used the Hereford Arizona Observatory (HAO), a private facility in southern Arizona (observatory code G95 in the IAU Minor Planet Center) to measure multi-band, absolute photometry of TYC 272 and TYC 1422. The facility has been used to measure transiting exoplanets with high photometric precision222http://brucegary.net/book_EOA/x.htm, and has recently been used to observe EBs and variable star lightcurves with similar levels of precision. HAO employs an 11-inch Celestron Schmidt-Cassegrain (model CPC 1100) telescope fork-mounted on an equatorial wedge and located in a dome. The 11-inch telescope’s CCD is an SBIG ST-8XE CCD with a KAF 1602E detector. HAO also includes a 14-inch Meade Schmidt-Cassegrain (model LX200GPS) telescope, fork-mounted on an equatorial wedge located in a dome. That telescope’s CCD is an SBIG ST-10XME with a KAF-3200ME detector. A 10-position filter wheel accommodates SDSS and Johnson/Cousins filter sets.

The stars were observed with Sloan filters with the 11-inch telescope. To derive the target star’s magnitude in band we use a generic photometry equation:


where is the desired magnitude for observed band , is an instrumental zero-shift for band , is the measured intensity in the observed band, is the integration time, is the zenith extinction for the observed band, is the airmass, is the star color sensitivity for the observed band and is a star color defined using two bands (e.g., or ). For the Johnson-Kron-Cousins bands, standard stars are taken from the list published by Landolt & Uomoto (2007) and (Landolt, 2009). For the SDSS bands, standard stars are taken from the list published by Smith et al. (2002). and magnitudes are determined using the conversion equations given in Smith et al. (2002) that convert , and magnitudes to , and , while is determined using the fluxes.

4.2 SMARTS Echelle Spectroscopy

We obtained spectra from the SMARTS 1.5m telescope at CTIO using the R 42000 Echelle spectrograph that covers nm. These observations were used to resolve the components of TYC 1422 and get individual RV measurements for each star in the system. Seven epochs were observed using ten minute integrations between February and April 2010 with typical signal-to-noise ratios of 20 per resolution element at nm.

The spectra were processed in the standard way for cross-dispersed Echelle spectra, using a pipeline written specifically for spectra taken with the Echelle spectrograph on the CTIO-1.5m. The routine processes the data using biases and quartz lamp observations taken at the beginning of the night, median combines three individual images while performing cosmic ray rejection, extracts the individual orders from the combined image and performs the wavelength solution on each order using a ThAr arc lamp taken either before or after each set of science exposures.

A total of 22 orders are used each night to derive radial velocities via cross-correlation with a standard template. We use the bright G0 star HD 108510 as the heliocentric RV standard star. Each spectral order is cross-correlated separately, then an iterative 3- clipping is performed prior to performing a weighted average to obtain a final RV measurement for each night. The components are identified each night via the peak and width of each feature in the CCF. The typical RV precisions ranged from 0.3 to 1.7 for the various components.

4.3 APO Echelle Spectroscopy

We obtained ARCES (Wang et al., 2003) spectroscopy from the APO 3.5m telescope to check for additional spectra from unresolved companions. ARCES is an Echelle spectrograph with spectral resolution 31500 that covers the entire wavelength range from nm. on a single 2048x2048 SITe CCD. Three exposures for a total integration time of 18 minutes were obtained for TYC 1422 on 2009 March 23, while two 20-minute exposures were taken of TYC 272 on 2010 June 19. The spectra were reduced using an IRAF script. Spectra are corrected for bias and dark subtraction, cosmic rays and bad pixels. Flatfielding is done using a combination of two different sets of quartz lamp exposures; one with a blue filter in place and another set without the filter. These two different sets of flats allow for optimal extraction of both the blue and red orders. A ThAr lamp is used for wavelength calibration. The co-added spectrum for TYC 1422 had a signal-to-noise of per resolution element at 600 nm, while TYC 272 had a signal-to-noise of per resolution element.

5 Case Studies

5.1 Tyc 272

5.1.1 Orbital Parameters

We first identified the binaries spectroscopically using the MPP RV data. Following Fleming et al. (2010), the RV uncertainties for each star are scaled by a “quality factor” that is a multiplicative scaling that modifies the formal RV uncertainties to better account for systematics defined as:


where represents the RV measurements and RMS is the root-mean-square residual. The quality factor is derived from the RMS of the RVs for the other 58 objects in the same field-of-view as the target, under the assumption that a majority of the objects per field should not have significant RV variability. The field containing TYC 272 has a = 3.44, suggesting the formal uncertainties are underestimated. Table 1 contains the times of observations, RVs, and both the formal and scaled RV uncertainties for TYC 272. A total of 14 epochs were obtained with an average, scaled RV uncertainty of 0.345 for TYC 272.

TYC 272 was fit using the RVSIM software program (Kane et al., 2007). Fig. 1 shows the best-fit orbital solution to the MPP RV data. The , after scaling by , is 2.62 with eight degrees of freedom. Because it is a short-period binary, there is additional RV jitter expected due to chromospheric activity, however, chromospheric effects are typically on the order of tens of (Santos et al., 2000). The eccentricity is consistent with a circularized orbit, as expected for the best-fit orbital period days (Mazeh, 2008). The semi-amplitude corresponds to a minimum mass of the secondary assuming the mass of the primary is , making this EB an excellent system to further study in the context of the mass-radius relationship. The best-fit orbital properties are listed in Table 2.

5.1.2 Stellar Properties

We utilize a spectral synthesis technique as part of the process to estimate basic stellar parameters such as , and for the primary. In Valenti & Fischer (2005), they outline a general procedure for deriving stellar properties using Spectroscopy Made Easy (SME; Valenti & Piskunov, 1996). We follow the exact methodology described in Valenti & Fischer (2005), but include some changes that allow for a more in-depth search of space. Specifically, we use the same spectral regions, line lists and continuum regions as in Valenti & Fischer (2005), but do not allow any of the individual abundances to be free parameters.

Our SME spectral analysis pipeline is divided into two steps, because the spectral synthesis parameters are highly correlated, and therefore a range of initial guesses is required to adequately explore the space. First, a forward modeling procedure is used to calculate values for model spectra over a large grid of parameter space. This allows us to focus our analysis on the region of parameter space with the most likely, global, best-fit stellar parameter solution. Second, we use the SME Levenberg-Marquardt minimization algorithm as was used in Valenti & Fischer (2005) to fit observed spectra with synthetic spectra for 100 different initial input parameters. We select these different combinations of input values from the 100 grid points with the lowest chi-square values determined from the first step of our analysis pipeline. Accessing Vanderbilt University’s ACCRE parallel computer facility allows us to efficiently explore the dependence of these initial parameters on our output parameters, as well as providing a good measure of the internal precision of SME’s minimization procedure.

Using our implementation of SME on the ARCES spectrum of TYC 272, we derive best-fit values of K, , metallicity . The value for is larger than expected for stars of this temperature, but we note that is a very challenging parameter to measure spectroscopically, even at high spectral resolution, and is particularly challenging at the ARCES spectral resolution of 31500.

In addition to SME spectral synthesis, we analyzed the spectrum following a standard spectroscopic method based on the requirements of excitation and ionization equilibria. Our automated analysis merges the approaches described in Porto de Mello et al. (2008) and Ghezzi et al. (2010), and will be described in a forthcoming paper of the MARVELS series Fleming et al. (2011). A total of 57 Fe I and 3 Fe II lines were used, after a 2- clipping was applied to remove those lines with too high or too low iron abundances. The equivalent widths (EWs) were automatically measured using the task ”bplot” in IRAF. We have derived a K, and metallicity . We investigated the effects of a contaminating spectrum by modifying the measured EWs by a correction factor equal to the inverse of the expected flux contribution from the primary. In the case of a solar-type dwarf with a mid-K companion, the expected flux ratio is 10% in V band, and the EW correction factor is 1.1. The general effect is to increase the derived effective temperature, surface gravity and metallicity. However, the increase is small and well within the internal errors of our technique. We therefore conclude that a flux ratio of 10% or lower would not significantly affect our measured values based on the EW method.

We also examine the H line profile, which serves as a and chromospheric activity indicator, following Lyra & Porto de Mello (2005) and derive a best-fit K. This temperature is lower than the value determined via the Fe line EW method, but is affected by contamination from the secondary that fills in the H wings, resulting in an effective temperature determination that is systematically cooler. An analysis of the H core filling suggests that TYC 272 is much more chromospherically active than the Sun, however, there is a degeneracy between flux contamination from a companion and core emission due to activity, both of which can fill in the core of the H line. We further note that it is often the case that H-based and color-based are systematically lower than temperatures determined spectroscopically, one possible cause being chromospheric activity (Porto de Mello et al., 2008). Given the relatively short orbital period of 5.728 days, it is likely that some amount of tidal spin-up has occurred, resulting in stronger magnetic fields and increased starspot activity.

We complement the spectroscopic analysis with a variety of additional techniques to estimate stellar parameters based on other observables. To further estimate the luminosity class, we use a reduced proper motion (RPM) diagram (Collier Cameron et al., 2007b) to determine whether the host star is a likely giant or a dwarf/subgiant. The RPM diagram is excellent at distinguishing giants from dwarfs/subgiants. We tested the RPM technique on the Kepler Input Catalog (KIC)333http://archive.stsci.edu/kepler/kic10/search.php, which has had its determinations of luminosity class verified using Kepler’s high-precision photometry (Koch et al., 2010). We find that only 2% of KIC giants were classified as dwarfs/subgiants using our RPM analysis. We use the proper motion values from the Tycho-2 catalog and NIR photometry from 2MASS (Table 3) to place TYC 272 on the RPM-J diagram (Fig. 2, bottom). The solid, curved line in the RPM diagram is the dividing line between giants and dwarfs. The labels show the approximate regions where dwarfs, subgiants and giants are located. Given how precise the RPM diagram is at identifying giant stars, the fact that TYC 272 lies in the dwarf/subgiant region of the RPM diagram means it is highly unlikely that TYC 272 is a giant star.

We use the infrared flux method (IRFM) (Casagrande et al., 2010) to estimate the effective temperature of the primary star via a Monte Carlo approach. We sample the observed , and magnitudes assuming Gaussian errors, sample uniformly in [Fe/H] from and uniformly in extinction from , where the upper bound is determined from the Schlegel et al. dust maps (Schlegel et al., 1998). Fig. 2 shows the results of the Monte Carlo IR flux method (top). We estimate the effective temperature based on the (-) color to be K, where the uncertainties are 1- equivalents centered on the median of the asymmetrical distribution. The difference in temperature between the IRFM and spectroscopic results can be caused by flux contribution from the secondary. Indeed, the peak of the (-)-based temperature distribution is slightly cooler, which can occur because the flux ratio in the NIR bands is larger than in the optical bands. This IRFM is therefore a lower bound on the true of the primary due to reddening by the companion.

To further study the effective temperature of the primary, we make use of the spectral energy distribution (SED). We use the optical band measurements from HAO, near-IR measurements via the 2MASS (Skrutskie et al., 2006) Point Source Catalog and UV measurements from GALEX (Martin et al., 2005), which are compiled in Table 3. NextGen models from Hauschildt et al. (1999) are used to construct the theoretical SED. Based on the results from the spectroscopic and IRFM analysis, we construct a model with , and and compare with the observed fluxes. These stellar parameters were selected as representative values between the spectral synthesis, EW measurements, H profile analysis and IRFM temperature calculations. Distance and extinction are free parameters

Fig. 3 (top) shows the result of the SED fit and the model spectrum. The blue points are the passband-integrated model fluxes for each filter, while the red crosses are the bandpasses of each filter (horizontal bars) and the measured error in the fluxes (vertical bars). A clear excess in the near-IR is seen when we do not consider flux from the secondary. Adopting a primary mass of and a radius of , we calculate a secondary mass of (K5 spectral type dwarf). We then add flux contribution from a second star with K and a radius of , and find that the inclusion of flux from a K-type companion can explain the near-IR excess seen in the SED Fig. 3 (bottom). This SED analysis does not provide an independent measurement of the stellar parameters, but rather serves as a plausibility experiment that demonstrates flux contribution from a K-type companion can explain the apparent NIR excess if the primary has a K.

Because the secondary is too faint to have its spectral features detected in the MPP spectra, the mass of the primary cannot be directly measured via RV. Instead, we make use of Padova isochrones (Marigo et al., 2008) to estimate the mass of the TYC 272 primary. Fig. 4 plots radii and masses as a function of effective temperature for two assumed ages (100 Myr and 9 Gyr) and a metalicity 0.01 . This metallicity is chosen based on the SME and EW determinations, which are in good agreement with each other. The uncertainty in does not significantly affect the derived range of probable masses for TYC 272, therefore we elect to use a single value for the tracks’ metallicities, for clarity. The derived values from the SME, EW, H and IRFM are plotted as horizontal bars representing the 1- confidence intervals. We elect to use a range equal to the EW 3- confidence interval as our best estimate for the primary because it is expected to be minimally affected by flux contribution from the secondary. We note however that all four are consistent to within 2-, and in fact the SME, EW and IRFM temperatures all agree to within 1-. Based on this range, we find the Padova models constrain for and .

5.1.3 Light Curve Parameters

The lightcurve for TYC 272 is shown in Fig. 5, with the best-fit model overplotted as the solid line. The top panel contains the primary eclipse, while the shallowness of the secondary eclipse can be seen in the bottom panel. The secondary eclipse of TYC 272 appears to be grazing, leading to model degeneracies between the ratio of the radii and the surface brightness ratio. This is a well-known problem for grazing EBs, and the solution is to constrain the flux ratio when fitting the lightcurve (Andersen et al., 1983). The flux ratio is usually constrained spectroscopically (Clausen et al., 2003; Southworth et al., 2004; Blake et al., 2008), and we utilized the ARCES spectrum to search for any evidence of a secondary spectrum. The ARCES spectrum was obtained at an orbital phase corresponding to an expected velocity shift of -27.25 , however, the signal-to-noise was not sufficient to detect a secondary spectrum for all flux ratios. Fig. 6 compares the ARCES spectrum with the (lower-resolution, R=12,000) FOE (Montes et al., 1997; Montes & Martin, 1998) spectrum of 9 Cet, a G2 V standard. Two different wavelength regions are examined that contain absorption features present in stars of mid-F through mid-M spectral types. No set of secondary spectra are seen, and we estimate an upper limit of % for the flux ratio at these wavelengths. In the absence of sufficiently high resolution spectra, we estimate the flux ratio using the mass range expected for in §5.1.2 and the measured RV semiamplitude of the primary . The flux ratio is constrained to be %, the estimated flux ratio in of a G-type primary and K-type companion. As shown in the SED analysis, such a flux ratio reproduces the observed SED well, particularly in the 2MASS NIR bands.

To model the lightcurves, we used the jktebop444http://www.astro.keele.ac.uk/~jkt/codes/jktebop.html implementation of the EBOP lightcurve model (Nelson & Davis, 1972; Popper & Etzel, 1981) to produce a least-squares fit to the data. We modified jktebop to include the normalization for each night of data as additional free parameters. We consulted various tabulations of linear limb darkening coefficients (van Hamme, 1993; Claret, 2000, 2004) and decided to adopt a limb darkening coefficient of 0.6 in the WASP band for the G-type primary of TYC 272. For the K-type secondary star we used a linear limb darkening coefficient of 0.7. An uncertainty of 0.1 is assumed for these parameters, and the small effect of this uncertainty is accounted for in the standard errors quoted for other parameters of the lightcurve solution. The values and uncertainties of and from the RV orbital solution are included as constraints to the least-squares fit to the lightcurve. Without including such a constraint, the best-fit value of , however, because this is determined via the widths of the eclipses, we do not get a robust estimate from the lightcurve data alone.

Table 4 summarizes the geometric parameters derived from fitting the photometry. To account for correlated errors (“red noise”) in the photometry, uncertainties in the parameters are derived using a “prayer-bead” bootstrap Monte Carlo method. This method uses the residuals from the best fit to create synthetic lightcurves with noise characteristics that are similar to the real data via circular permutation. Comparison of the parameters that are common to both the RV analysis and the photometric analysis show excellent agreement. The orbital period agrees to within 1-, while the and values are consistent with a circular orbit.

5.1.4 Determination of Mass and Radius

The mass of the secondary, as well as the radii of the two binary components, can be determined in an SB1 if a primary and secondary eclipse are both observed and the mass of the primary star is known (or assumed). This is accomplished via the mass function:


Where is the gravitational constant, is the orbital period, and is determined observationally using the definition of the RV semiamplitude from the radial velocity equation:


where is the orbital eccentricity, is the line-of-sight orbital inclination and the subscript represents the primary and secondary, respectively. The radii of the two stars can be calculated directly from the geometric lightcurve parameters and found in Table 4 and the semimajor axis .

The uncertainties for , , and are determined as follows. The uncertainty of is taken from the independent measurement of the primary’s mass, if such a determination exists, or from the range of estimates for based on the star’s color and/or spectra. Using implicit partial differentiation, the propagation of error equation, and removing the cross terms since and are independent, the uncertainty in from Eqn. 3 is given by:


where is the right-hand-side of Eqn. 3:


and is found by using the propagation of error equation and assuming , , and are all independent. Likewise, the uncertainties in the stellar radii are calculated via the propagation of error equation using the measured uncertainties in , (Table 4) and the calculated found via propagating the uncertainties in and .

We adopt a primary mass of based on the stellar characterization presented in §5.1.2 and calculate , and (Fig. 7). Both the primary and secondary are consistent (1-) with the model predictions for mass and radii. Metalicity is a minor effect on the mass and radius of the isochrones, so only the tracks corresponding to are shown.

Although the relative errors for the masses and radii in the specific case of TYC 272 are larger than typically used to study the mass-radius relationship (% for , 8% and 18% for and , respectively), relative errors of a few percent are possible if better constraints on the primary mass could be determined and more precise photometry of the grazing secondary eclipse could be obtained. Specifically, large-aperture telescopes and multi-dimensional correlation software (Zucker & Mazeh, 1994; Zucker et al., 2003) can be used to extract RV solutions for faint companions, such that the masses of both components can be measured without relying on models or empirical relations to constrain . In addition, high precision photometric follow-up of the shallow secondary eclipse will improve the precision in the radius. TYC 272 is not an ideal benchmark system because of the shallow, grazing secondary eclipse, however, the analysis presented here demonstrates the required steps to derive precise masses and radii for EBs using a combination of transit and RV exoplanet survey data. A summary of these analysis steps for EBs with low flux ratios, such as TYC 272, are described in §6.

5.2 Tyc 1422

5.2.1 Orbital Parameters

TYC 1422 had a total of 16 epochs observed during the MPP, with a = 3.07 for its field and an average, scaled RV uncertainty of 0.304 . When fitting the MPP RV measurements for TYC 1422, the residuals were found to be exceedingly large (RMS ). Such a large scatter is often caused by contamination of an additional spectrum from an unresolved source, a likely scenario for short-period SBs with similar fluxes.

The MPP DFDI pipeline was not designed to treat multiple sources within a single spectrum, and can produce large systematic errors in those cases. It is important to note that the inability to obtain precise RV measurements for SB2s applies to the MPP instrument only. Other survey instruments, including the MARVELS survey instrument itself, will be able to treat the case of double-lined binaries. The MARVELS survey pipeline stores information about the fringe fitting that can be used to remove the fringes from the stellar DFDI spectrum and re-create a ”traditional”, de-fringed stellar spectrum. This spectrum can then be used to perform cross-correlation analysis and identify and measure multiple CCF components, thereby measuring RVs of double-lined binary stars.

Because the MPP instrument passed starlight through an Iodine cell, a similar reconstruction of the stellar spectrum is not possible due to the Iodine lines that are present. Attempting to model the Iodine lines and remove them from the spectrum prior to de-fringing is also difficult, because the Iodine lines are numerous and very narrow. In fact, most of the Iodine lines are unresolved at the resolution of the MPP instrument. For this reason, independent RV follow-up observations (SMARTS data) were obtained to measure the RVs of the multiple components. In general, such additional follow-up is not required and RV measurements for both components of double-lined binaries can be obtained directly from the survey data itself.

Despite the large scatter in the residuals, the derived orbital period from the MPP data is in good agreement with the period derived from the SuperWASP photometry ( 6.2 days). To derive an improved orbital solution, we use the higher resolution SMARTS spectra to resolve the components and measure the Doppler reflex velocities of the individual components. We do not use the MPP data in our final analysis because it is so heavily affected by systematic errors. However, for completeness, we present the MPP RV data in Table 5.

The TYC 1422 system consists of three sources, as demonstrated in the triple set of Na D doublets from the ARCES spectrum in Fig. 8. We do not attempt an SED analysis due to the complexity of the system, however, the basic stellar properties, including the HAO absolute photometry, are compiled in Table 6. The RVs for each component are phase-folded to the orbital period and ephemeris as determined by the SuperWASP photometry in Fig. 9. The RV measurements for each component are connected by different line types for visualization. Components 2 and 3 are identified as the eclipsing binary pair. Their RVs are in antiphase, as expected for an orbiting pair, and have relative velocities that cross the zero-line near the phases of primary and secondary eclipse. Component 1 is found to have a gradual linear trend of . Fig. 10 shows the RVs from Component 1 and the best-fit linear relation. The large uncertainty for the third epoch is caused by blending. Table 7 contains the times and measured RVs for each of the three components.

The SMARTS RVs were fit using MPRVFIT2, an IDL procedure written by N. De Lee which combines a periodogram searching algorithm and a non-linear least squares fitter to find the optimal Keplerian orbit for a given set of radial velocity points. The period search algorithm consists of a modified Lomb-Scargle periodogram based on Cumming (2004) with an analytical False Alarm Probability (FAP) derived from Baluev (2008). The software allows for simultaneous fitting of RVs from both components of a binary system. The SMARTS spectra consist of seven epochs, resulting in only one degree of freedom for a full Keplerian fit. However, because Components 2 and 3 are a physically bound pair, we can exploit the fact that several parameters are identical for the two components. Specifically, the period (), eccentricity (), epoch of periastron () and systemic RV () are identical, while the argument of periastron () differs by radians. This results in seven degrees of freedom between the 14 combined data points, the seven parameters being , where and are the RV semi-amplitudes of the two components. Due to the linear trend observed in Component 1, we also include a linear term that reduces the degrees of freedom to six. Fig. 11 shows the best-fit orbital solution for Components 2 and 3, and the orbital parameters are presented in Table 8. We find a best-fit period of 6.2005 days and a slightly eccentric orbit with and a mass ratio of 0.78. The orbital period is in agreement with the photometric period to within . The slope of the linear term is , which is of comparable magnitude and direction as the slope measured on Component 1, and therefore may be residual instrument drift.

5.2.2 Light Curve Parameters

We performed the same analysis using the jktebop software as done for TYC 272, however, third light (flux contributed from Component 1, in Table 9) is included as a free parameter in the least-squares fit to the lightcurve of TYC 1422. In principle one can compare a spectroscopic flux ratio based on the EWs of spectral lines with the best-fit flux ratio derived from the lightcurve solution. However, in the case of TYC 1422 and the ARCES spectrum, most of the lines are blended at the resolution of the spectrograph. The Na D lines shown in Fig. 8 are quite visible, however, the wings of the features are still blended, and it is challenging to find an appropriate pseudo-continuum from which to measure the EWs. Furthermore, the EWs of the Na D doublets are known to depend on (Tripicchio et al., 1997; Díaz et al., 2007), and the of the third component is not constrained in our analysis.

We use limb darkening coefficients of 0.55 for the primary and secondary. The third light is accounted for within the WASP band, which ranges from 400-700 nm. and is similar to . The lightcurve with the best-fit model overplotted is shown in Fig. 12. Table 9 summarizes the geometric parameters derived from photometry fitting. Correlated errors are accounted for using the “prayer-bead” bootstrap Monte Carlo approach, as described in §5.1.3. Timing of the primary and secondary eclipse events confirm the system is in an eccentric orbit.

5.2.3 Determination of Mass and Radius

Since both components of the EB pair have Doppler measurements, we can measure the masses and radii of Components 2 and 3 directly from the RV semi-amplitudes and the lightcurve. We find the semimajor axes and using the RV semiamplitudes from Eqn. 4. The masses can then be determined using the center of mass equation and the complete version of Kepler’s Third Law:


where and are the masses of the primary and secondary, respectively.

The radii of the two stars can be calculated directly from the geometric parameters measured from the lightcurves in Table 9 and the value of the semimajor axis . The primary star is evolved off the main sequence and has a mass of and a radius of . The secondary has a mass of and a radius of and is consistent with a main-sequence G-type dwarf. The relative uncertainties for and are 2.8% and 4.2%, very close to the desired precision goal of %. Additional, high precision photometry of the eclipses will be helpful to further improve the precisions. The relative uncertainties for and are 2.9% and 7.4%. The larger uncertainty for is in part due to the larger uncertainties in the RV measurements for that component. The SMARTS spectra had relatively low signal-to-noise ( per resolution element), and obtaining higher signal-to-noise spectra will reduce the uncertainty in to better than 3%. Fig. 13 compares the masses and radii with Padova isochrones at an age of and metalicities of . TYC 1422 is consistent with solar metalicity (solid line), whereas the dashed lines are additional metalicity tracks for comparison purposes.

6 Framework For EBs From Surveys

When the RV semiamplitudes for both components can be measured spectroscopically (SB2s), the masses and radii for both stars can be determined directly from observable quantities in the survey data. TYC 1422 is such an example, although it is not a good benchmark system due to the presence of a third spectrum. Additional spectra can blend with the other components and affect their derived RVs, and the additional flux contribution complicates the photometry modeling. Despite these challenges, we were able to obtain 2.9% relative precision on , 2.8% relative precision on and 4.2% relative precision on . The MPP data were unable to be de-fringed due to the presence of the Iodine lines from the gas cell, and therefore follow-up observations with SMARTS data were required. Obtaining better signal-to-noise spectra would reduce the uncertainties in the masses further. We note that the MARVELS survey instrument does not pass starlight through an Iodine cell, and spectra of SBs with similar flux ratios can be de-fringed and analyzed using cross-correlation techniques like the ones used on the SMARTS spectra of TYC 1422 without additional follow-up.

EBs with smaller flux ratios require model-dependent constraints on the primary mass to derive the mass of the secondary and both radii, until observations with larger telescopes can be obtained. In this paper, we provide a framework for identifying EBs of interest from SB1s using survey data and catalog information. We summarize the steps taken in a suggested “analysis flow diagram” presented in Fig. 14. This framework applies only to those EBs that have both a measurable primary and secondary eclipse. Once the mass of the primary and an uncertainty are adopted via comparison with isochrones and/or using empirical relations, the mass of the secondary, both stellar radii, and their uncertainties can be calculated using the equations presented in §5.1.4.

We apply this technique to TYC 272, however, the grazing secondary eclipse makes it a poor benchmark EB as well. The 5.9% relative error on is dominated by the large uncertainty for the estimate of . Observations with high resolution spectrographs at high signal-to-noise can be used to derive and using TODCOR (see Bender & Simon, 2008, for an example in the near-IR). Such a technique can exploit the fact that the other orbital parameters are well-determined from the survey data, and only a measurement of is required to directly measure and thereby complete the analysis. In addition, obtaining a high signal-to-noise, high resolution spectrum at a phase where the spectral lines are well-separated can yield spectroscopically-determined flux ratios, which are then used as additional constraints in the lightcurve modeling of eclipses that are grazing.

In favorable cases where there are deep, total eclipses and no third light contamination, transit survey photometry itself can yield radii to better than 2% relative error (Southworth et al., 2011). However, in both SB2 and SB1 situations, the survey photometry may sometimes be insufficient to obtain such high precisions. Because the epoch of central transit is well-determined from the survey data, follow-up observations can be planned when a primary or secondary eclipse is known to occur. In addition, obtaining transits in multiple filters can aid in the EB analysis. We note that this can be a problem for any study of EBs using transiting exoplanet survey data, and is not unique to the methodology presented here.

7 Conclusion

We have presented a framework for eclipsing binary science via cross-referencing of photometric and Doppler exoplanet survey data. We have found two new EBs around F/G primaries; TYC 272, an SB1 consisting of a K dwarf secondary, and TYC 1422, a triple-lined SB that includes an evolved primary star and a G dwarf secondary. Although neither of these two specific binaries are ideal benchmark systems for studying the mass-radius relationship, we have shown that the combination of survey data can be a rich source of EBs, for which many physical parameters can be directly measured or constrained from the survey data itself. In particular, the analysis of TYC 272 provides a framework for estimating the masses and radii of SB1s from survey data. Obtaining additional, high resolution, high signal-to-noise spectra can then directly yield the semiamplitudes, and thus masses, of both components. Such an approach can be used in particular to study the mass-radius relationship of F/G + K/M EBs.

We thank Roger Cohen for useful discussions on theoretical isochrones. L.G. thanks Dr. Simone Daflon and Dr. Herman Hensberge for the helpful discussions. L.G. acknowledges financial support provided by the PAPDRJ - CAPES/FAPERJ Fellowship. L. Dutra-Ferreira acknowledges financial support provided by CAPES fellowship. G.F.P.M. acknowledges financial support from CNPq grants no. 476909/2006-6 and 474972/2009-7, plus a FAPERJ grant no. APQ1/26/170.687/2004. Funding for the multi-object Doppler instrument was provided by the W.M. Keck Foundation. The pilot survey was funded by NSF with grant AST-0705139, NASA with grant NNX07AP14G and the University of Florida. The SuperWASP Consortium consists of astronomers primarily from the Queen’s University Belfast, St. Andrews, Keele, Leicester, The Open University, Isaac Newton Group La Palma and Instituto de Astrofísica de Canarias. The SuperWASP Cameras were constructed and operated with funds made available from Consortium Universities and the UK’s Science and Technology Facilities Council. This research has made use of the SIMBAD database, operated at CDS, Strasbourg, France. Based on observations with the SDSS 2.5-meter telescope. Funding for the SDSS and SDSS-II has been provided by the Alfred P. Sloan Foundation, the Participating Institutions, the National Science Foundation, the U.S. Department of Energy, the National Aeronautics and Space Administration, the Japanese Monbukagakusho, the Max Planck Society, and the Higher Education Funding Council for England. The SDSS Web Site is http://www.sdss.org/. The SDSS is managed by the Astrophysical Research Consortium for the Participating Institutions. The Participating Institutions are the American Museum of Natural History, Astrophysical Institute Potsdam, University of Basel, University of Cambridge, Case Western Reserve University, University of Chicago, Drexel University, Fermilab, the Institute for Advanced Study, the Japan Participation Group, Johns Hopkins University, the Joint Institute for Nuclear Astrophysics, the Kavli Institute for Particle Astrophysics and Cosmology, the Korean Scientist Group, the Chinese Academy of Sciences (LAMOST), Los Alamos National Laboratory, the Max-Planck-Institute for Astronomy (MPIA), the Max-Planck-Institute for Astrophysics (MPA), New Mexico State University, Ohio State University, University of Pittsburgh, University of Portsmouth, Princeton University, the United States Naval Observatory, and the University of Washington. Based on observations obtained with the Apache Point Observatory 3.5-meter telescope, which is owned and operated by the Astrophysical Research Consortium. This publication makes use of data products from the Two Micron All Sky Survey, which is a joint project of the University of Massachusetts and the Infrared Processing and Analysis Center/California Institute of Technology, funded by the National Aeronautics and Space Administration and the National Science Foundation. The Center for Exoplanets and Habitable Worlds is supported by the Pennsylvania State University, the Eberly College of Science, and the Pennsylvania Space Grant Consortium. This work was conducted in part using the resources of the Advanced Computing Center for Research and Education at Vanderbilt University, Nashville, TN.


  • Andersen et al. (1983) Andersen, J., Clausen, J. V., Nordstroem, B., Reipurth, B. 1983, A&A, 121, 271
  • Antipin et al. (2010) Antipin, S. V., Kazarovets, E. V., & Pastukhova, E. N. 2010, Peremennye Zvezdy Prilozhenie, 10, 33
  • Asplund et al. (2009) Asplund, M., Grevesse, N., Sauval, A. J., Scott, P. 2009, ARA&A, 47, 481
  • Baluev (2008) Baluev, R. V. 2008, MNRAS, 385, 1279
  • Bender & Simon (2008) Bender, C. F. & Simon, M. 2008, ApJ, 689, 416
  • Blake et al. (2008) Blake, C. H., Torres, G., Bloom, J. S., Gaudi, B. S. 2008, ApJ, 684, 635
  • Casagrande et al. (2010) Casagrande, L., Ramírez, I., Meléndez, J., Bessell, M., Asplund, M. 2010, A&A, 512, A54
  • Chabrier et al. (2007) Chabrier, G., Gallardo, J., & Baraffe, I. 2007, A&A, 472, L17
  • Claret (2000) Claret, A. 2000, A&A, 363, 1081
  • Claret (2004) Claret, A. 2004, A&A, 428, 1001
  • Clausen et al. (2003) Clausen, J. V., Storm, J., Larsen, S. S., Giménez, A. 2003, A&A, 402, 509
  • Collier Cameron et al. (2007a) Collier Cameron, A., et al. 2007, MNRAS, 375, 951
  • Collier Cameron et al. (2007b) Collier Cameron, A., et al. 2007, MNRAS, 380, 1230
  • Collier Cameron et al. (2010) Collier Cameron, A., et al. 2010, MNRAS, 407, 507
  • Coughlin et al. (2010) Coughlin, J. L., Lopez-Morales, M., Harrison, T. E., Ule, N., Hoffman, D. I. 2010, arXiv:1007.4295v2
  • Cumming (2004) Cumming, A. 2004, MNRAS, 354, 1165
  • Devor et al. (2008) Devor, J. 2008, ApJ, 687, 1253
  • Díaz et al. (2007) Díaz, R. F., Cincunegui, C., & Mauas, P. J. D. 2007, MNRAS, 378, 1007
  • Erskine (2002) Erskine, D. J. 2002, PASP, 115, 255
  • Erskine et al. (2003) Erskine, D. J., Edelstein, J., Feuerstein, W. M., Welsh, B. 2003, ApJ, 592, 103
  • Fernandez et al. (2009) Fernandez, J. M., et al. 2009, ApJ, 701, 764
  • Fleming et al. (2010) Fleming, S. W., et al. 2010, ApJ, 718, 1186
  • Fleming et al. (2011) Fleming, S. W., et al. 2011, in prep
  • Ge (2002) Ge, J. 2002, ApJ, 571L, 165
  • Ge et al. (2002) Ge, J., Erskine, D. J., Rushford, M. 2002, PASP, 114, 1016
  • Ge et al. (2006a) Ge, J., et al. 2006a, Proc. SPIE, 6269, 75
  • Ge et al. (2009) Ge, J. et al. 2009, Proc. SPIE, 7440, 18
  • Ghezzi et al. (2010) Ghezzi, L., Cunha, K., Smith, V. V., de Araújo, F. X., Schuler, S. C., & de la Reza, R. 2010, ApJ, 720, 1290
  • Gunn et al. (2006) Gunn, J. E., et al. 2006, AJ, 131, 2332
  • Hauschildt et al. (1999) Hauschildt, P. H., Allard, F. & Baron, E. 1999, ApJ, 512, 377
  • Hebb et al. (2010) Hebb, L., et al. 2010, ApJ, 708, 224
  • Høg et al. (2000) Høg, E., et al. 2000, A&A, 355L, 27
  • Imbert (2002) Imbert, M. 2002, A&A, 387, 850
  • Kane et al. (2007) Kane, S. R., Schneider, D. P., Ge, J. 2007, MNRAS, 377, 1610
  • Kane et al. (2009) Kane, S. R., Mahadevan, S., von Braun, K., Laughlin, G., Ciardi, D. R. 2009, PASP, 121, 1386
  • Koch et al. (2010) Koch, D. G., et al. 2010, ApJ, 713, 79
  • Landolt & Uomoto (2007) Landolt, A. U. & Uomoto, A. K. 2007, AJ, 133, 768
  • Landolt (2009) Landolt, A. U. 2009, AJ, 137, 4186
  • López-Morales (2007) López-Morales, M. 2007, ApJ, 660, 732
  • Lyra & Porto de Mello (2005) Lyra, W. & Porto de Mello, G. F. 2005, A&A, 431, 329
  • Marigo et al. (2008) Marigo, P., Girardi, L., Bressan, A., Groenewegen, M. A. T., Silva, L., Granato, G. L. 2008, A&A, 482, 883
  • Martin et al. (2005) Martin, D. C., et al. 2005, ApJ, 619, 1
  • Mazeh (2008) Mazeh, T. 2008, in Groupil, M.-J., Zahn, J.-P., eds, EAS Publ. Ser. Vol. 29, Observational Evidence for Tidal Interactions in Close Binary Systems. EDP Sciences, Paris, p. 1
  • Meibom et al. (2009) Meibom, S., et al. 2009, AJ, 137, 5086
  • Montes et al. (1997) Montes, D., Martin, E. L., Fernandez-Figueroa, M. J., Cornide, M., de Castro, E. 1997, A&AS, 123, 473
  • Montes & Martin (1998) Montes, D. & Martin, E. L. 1998, A&AS, 128, 485
  • Morales et al. (2010) Morales, J. C., Gallardo, J., Ribas, I., Jordi, C., Baraffe, I., Chabrier, G. 2010, ApJ, 718, 502
  • Nelson & Davis (1972) Nelson, B. & Davis, W. D. 1972, ApJ, 174, 617
  • Norton et al. (2007) Norton, A. J., et al. 2007, A&A, 467, 785
  • Pollacco et al. (2006) Pollacco, D. L. et al. 2006, PASP, 118, 1407
  • Popper & Etzel (1981) Popper, D. M. & Etzel, P. B. 1981, AJ, 86, 102
  • Popper (1997) Popper, D. M. 1997, AJ, 114, 1195
  • Porto de Mello et al. (2008) Porto de Mello, G. F., Lyra, W., Keller, G. R. 2008, A&A, 488, 653
  • Prsǎ & Zwitter (2005) Prsǎ, A. & Zwitter, T. 2005, ApJ, 628, 426
  • Prsǎ et al. (2010) Prsǎ, A., et al. 2010, arXiv:1006.2815v1
  • Queloz et al. (2010) Queloz, D., et al. 2010, A&A, 517, L1
  • Santos et al. (2000) Santos, N. C., Mayor, M., Naef, D., Pepe, F., Queloz, D., Udry, S., Blecha, A. 2000, A&A, 361, 265
  • Schlegel et al. (1998) Schlegel, D. J., Finkbeiner, D. P. & Davis, M. 1998, ApJ, 500, 525
  • Skrutskie et al. (2006) Skrutskie, M. F. et al. 2006, AJ, 131, 1163
  • Smith et al. (2002) Smith, J .A., et al. 2002, AJ, 123, 2121
  • Southworth et al. (2004) Southworth, J., Maxted, P. F. L. & Smalley, B. 2004, MNRAS, 349, 547
  • Southworth et al. (2011) Southworth, J., Pavlovski, K., Tamajo, E., Smalley, B., West, R. G., Anderson, D. R. 2011, arXiv:1103.1519
  • Torres et al. (2008) Torres, G., Vaz, L. P. R., & Sandberg Lacy, C. H. 2008, AJ, 136, 2158
  • Torres et al. (2010) Torres, G., Andersen, J. & Giménez, A. 2010, A&A Rev., 18, 67
  • Tripicchio et al. (1997) Tripicchio, A., Severino, G., Covino, E., Terranegra, L., & Garcia Lopez, R. J. 1997, A&A, 327, 681
  • Valenti & Piskunov (1996) Valenti, J. A. & Piskunov, N. 1996, 118, 595
  • Valenti & Fischer (2005) Valenti, J. A. & Fischer, D. A. 2005, ApJS, 159, 141
  • van Hamme (1993) van Hamme, W. 1993, AJ, 106, 2096
  • Wang et al. (2003) Wang, S. et al. 2003, Proc. SPIE, 4841, 1145
  • Wils (2010) Wils, P. 2010, Information Bulletin on Variable Stars, 5955, 1
  • Zucker & Mazeh (1994) Zucker, S. & Mazeh, T. 1994, ApJ, 420, 806
  • Zucker et al. (2003) Zucker, S., Mazeh, T., Santos, N. C., Udry, S., Mayor, M. 2003, A&A, 404, 775
HJD RV (formal) (scaled)
(km s) (km s) (km s)
2454101.98701 22.075 0.086 0.297
2454106.01738 -13.138 0.082 0.282
2454158.78668 -3.698 0.100 0.343
2454186.75584 -19.438 0.081 0.278
2454188.74844 70.056 0.088 0.302
2454189.76185 87.249 0.090 0.310
2454194.79179 80.652 0.122 0.420
2454195.81803 79.963 0.119 0.409
2454215.66891 -16.141 0.086 0.297
2454221.71415 -6.716 0.098 0.339
2454224.70505 69.756 0.106 0.364
2454249.66233 -19.993 0.100 0.343
2454251.67471 65.708 0.101 0.349
2454255.67428 -18.045 0.145 0.500
Table 1: MPP RV Observations - TYC 272
Parameter Value 1- Uncertainty
Period (days) 5.7282 0.0003
K () 54.584 0.294
5.548 0.006
2454102.9539 0.0074
Sys. vel. (km s) 34.430 0.231
2.62 -
Table 2: Orbital Properties - TYC 272
Parameter Value 1- Uncertainty
(J2000)aaTycho-2 Catalog (Høg et al., 2000) 176.54255245  (deg) 11:46:10.21  (HH:MM:SS)
(J2000)aaTycho-2 Catalog (Høg et al., 2000) 1.68753726  (deg) +01:41:15.13  (DD:MM:SS)
bbGALEX (Martin et al., 2005) 21.443 0.113
bbGALEX (Martin et al., 2005) 15.398 0.0042
10.945 0.035
10.266 0.018
9.849 0.022
9.460 0.035
12.055 0.100
10.547 0.010
10.045 0.025
9.880 0.020
9.800 0.026
cc2MASS (Skrutskie et al., 2006) Point Source Catalog 8.893 0.021
cc2MASS (Skrutskie et al., 2006) Point Source Catalog 8.520 0.036
cc2MASS (Skrutskie et al., 2006) Point Source Catalog 8.443 0.023
aaTycho-2 Catalog (Høg et al., 2000) 146.9 1.6
aaTycho-2 Catalog (Høg et al., 2000) -138.3 1.6
5.417 -
Table 3: Stellar Properties - TYC 272 A
Parameter Value 1- Uncertainty
(Surface brightness ratio at disk center) 0.14 0.01
0.096 0.004
0.60 0.12
(degrees) 86.6 0.3
0.0004 0.0009
-0.002 0.001
Period (days) 5.72840 0.00001
(HJD) 2454534.7020 0.0007
0.060 0.002
0.036 0.006
(400-700 nm band) 0.048 0.026
0.002 0.0008
(degrees) 284 101
Table 4: Light Curve Properties - TYC 272
HJD RV (formal) (scaled)
(km s) (km s) (km s)
2454101.86376 10.848 0.090 0.276
2454102.03097 -12.295 0.128 0.393
2454105.97331 67.537 0.109 0.335
2454128.81685 49.394 0.099 0.304
2454136.77639 62.745 0.092 0.282
2454136.81223 53.464 0.093 0.286
2454164.75574 -11.161 0.106 0.325
2454165.75527 33.872 0.090 0.276
2454165.79120 33.716 0.091 0.279
2454186.65694 52.251 0.091 0.279
2454191.72814 53.452 0.099 0.304
2454194.74076 13.140 0.090 0.276
2454195.73280 -11.254 0.090 0.276
2454221.62566 34.133 0.094 0.289
2454224.61948 33.477 0.094 0.289
2454254.63340 59.181 0.129 0.396
Table 5: MPP RV Observations - TYC 1422aaThese data are dominated by flux contamination and are not used in the final RV analysis.
Parameter Value 1- Uncertainty
(J2000)aaTycho-2 Catalog (Høg et al., 2000) 152.41204974  (deg) 10:09:38.89  (HH:MM:SS)
(J2000)aaTycho-2 Catalog (Høg et al., 2000) 17.58113689  (deg) +17:34:52.09  (DD:MM:SS)
10.115 0.010
9.797 0.008
9.726 0.013
bb2MASS (Skrutskie et al., 2006) Point Source Catalog 8.875 0.026
bb2MASS (Skrutskie et al., 2006) Point Source Catalog 8.641 0.028
bb2MASS (Skrutskie et al., 2006) Point Source Catalog 8.563 0.023
aaTycho-2 Catalog (Høg et al., 2000) -30.6 0.9
aaTycho-2 Catalog (Høg et al., 2000) -5.6 1.0
Table 6: Stellar Properties - TYC 1422 System
Component 1 Component 2 Component 3
(km s) (km s) (km s) (km s) (km s) (km s)
2455236.67830 43.107 0.298 -18.420 0.602 132.754 1.650
2455241.64797 41.238 0.450 67.401 0.601 9.110 2.507
2455254.64277 31.593 10.625 31.593 10.625 72.102 1.688
2455270.64841 40.024 0.488 90.128 0.490 -12.069 1.648
2455278.59449 38.709 0.901 78.638 0.655 0.463 1.490
2455284.60109 39.142 0.304 87.278 0.607 -8.814 1.602
2455293.53225 38.567 0.381 -10.746 1.071 108.512 2.509
Table 7: SMARTS RV Observations - TYC 1422
Component 2 Component 3
Parameter Value 1- Uncertainty Value 1- Uncertainty
Period (days) 6.2005 0.002 6.2005 0.002
K () 65.366 0.516 84.046 1.297
-0.1494 0.0123 0.1494 0.0123
0.0576 0.0170 -0.0576 0.0170
2455236.8133 0.0554 2455236.8133 0.0554
Sys. vel. (km s) 45.026 0.356 45.026 0.356
Linear Trend () -0.048 0.018 -0.048 0.018
4.135 - 4.135 -
Table 8: Orbital Properties - TYC 1422
Parameter Value 1- Uncertainty
(Surface brightness ratio at disk center) 0.88 0.01
0.162 0.002
0.43 0.03
(degrees) 87.7 0.5
-0.1283 0.0003
0.092 0.011
(Contamination from Comp.) 0.25 0.08
Period (days) 6.199450 0.000007
(HJD) 2454194.3663 0.0008
0.114 0.003
0.049 0.002
(400-700 nm band) 0.16 0.02
0.158 0.006
(degrees) 144 3
Table 9: Light Curve Properties - TYC 1422
Figure 1: Best-fit orbital solution for TYC 272 from the MPP data. The RV uncertainties have been scaled by a factor of 3.44 to account for additional systematic errors.
Figure 2: Top: Monte Carlo calculation of effective temperature using the infrared flux method. We estimate a K from this analysis. Bottom: RPM-J reduced proper motion diagram. The dividing line between giants and dwarfs is the solid curve. Approximate regions where dwarfs, subgiants and giants inhabit are labeled. The location of TYC 272 is consistent with a dwarf or subgiant.
Figure 3: SED fits using NextGen models to TYC 272. Top: Model that only includes flux from the primary, assuming stellar parameters from the combined spectroscopic analysis. Bottom: Model that includes flux contribution from a K5 dwarf secondary, demonstrating the near-IR excess can be explained by a K-type companion.
Figure 4: Radius (top) and mass (bottom) as a function of effective temperature for different Padova isochrones and metalicity . The solid line is for an assumed young age (100 Myr) while the dashed line is for an assumed old age (9 Gyr). The horizontal, solid lines are the 1- confidence intervals for the various stellar parameter techniques that were applied. Only stars with are shown, since there are no indications the primary has a lower . We use these isochrones to constrain the mass of the primary () in order to derive , and .
Figure 5: Phase-folded SuperWASP lightcurve and best-fit model for the primary (top) and secondary (bottom) eclipses of TYC 272.
Figure 6: Comparison between the R 12,000 FOE spectrum of the G2.5V standard 9 Cet (top of each plot) with the R 30,000 ARCES spectrum of TYC 272 (bottom of each plot). The two wavelength ranges contain absorption lines present in stars mid-F though mid-M. No secondary set of spectra can be seen.
Figure 7: Radii and masses for the two components of TYC 272. Both stars are consistent (1-) with model predictions. Only a metalicity of is shown, since it has only a minor effect on the masses and radii of dwarfs and subgiants.
Figure 8: ARCES spectrum of TYC 1422 centered on the Na D doublets. The three sets of doublets can be clearly seen, and components are marked with the same names assigned to them while measuring their RVs with the SMARTS data. The components were identified in the SMARTS CCF’s by examining the width and shape of their CCF peaks.
Figure 9: Photometry and RVs for the three components of TYC 1422 phase-folded to the best-fit period and epoch of central transit in Table 9. The connecting lines between the measured RVs for each component serve as a visual aid. Components 2 and 3 are identified as the EB pair. Their measured RVs are antiphased as one would expect for an orbiting pair. They also phase appropriately such that the RV zero-crossing occurs near primary and secondary eclipses. Note the connecting lines do not cross at the precise primary and secondary eclipse phases because they do not represent the actual Keplerian orbit.
Figure 10: RVs for Component 1 in TYC 1422. The large error bar for the third epoch is caused by blending. The data are well-fit with a gradual slope of , which could be residual instrument drift given that a similar slope is found when fitting Components 2 and 3.
Figure 11: Best-fit orbital solution for Components 2 and 3 of TYC 1422 from the SMARTS observations. The black data points and curve are the RVs for Component 2 and the best-fit orbital solution, while the grey points are for the less massive Component 3. The jagged nature of the best-fit curve is a result of phase-folding the small linear-trend included in the fitting. The bottom panel shows the residuals to the fit.
Figure 12: Phase-folded SuperWASP lightcurve and best-fit model for the primary (top) and secondary (bottom) eclipses of TYC 1422.
Figure 13: Mass-Radius diagram for TYC 1422, Components 2 and 3. Padova isochrones are overplotted using a . This age was selected as the best-fit by visual inspection of tracks with differing ages. We find excellent agreement using a solar-metalicity of . For comparison, we also show tracks with (left track, dotted line) and (right track, dotted line).
Figure 14: Suggested analysis flow diagram that demonstrates the steps one would take to obtain masses and radii for an EB with a small flux ratio, starting with exoplanet survey data and catalog information. A dashed line represents a “no” decision. The analysis depicted here only works for EBs that undergo both primary and secondary eclipses.
Comments 0
Request Comment
You are adding the first comment!
How to quickly get a good reply:
  • Give credit where it’s due by listing out the positive aspects of a paper before getting into which changes should be made.
  • Be specific in your critique, and provide supporting evidence with appropriate references to substantiate general statements.
  • Your comment should inspire ideas to flow and help the author improves the paper.

The better we are at sharing our knowledge with each other, the faster we move forward.
The feedback must be of minimum 40 characters and the title a minimum of 5 characters
Add comment
Loading ...
This is a comment super asjknd jkasnjk adsnkj
The feedback must be of minumum 40 characters
The feedback must be of minumum 40 characters

You are asking your first question!
How to quickly get a good answer:
  • Keep your question short and to the point
  • Check for grammar or spelling errors.
  • Phrase it like a question
Test description