in relatively close primordial binaries which contain at least one star that can evolve substantially by the current epoch , mass transfer from the primary to the secondary is assured .
the course of the evolution of each system depends on its initial orbital period , the two stellar masses , whether the binary is part of a triple system , spin and tidal interactions , and so forth .
some of the possible evolutionary channels have been explored and even systematically mapped out by rappaport et al .
( 1994 ) , iben & tutukov ( 1995 ) , nelson & eggleton ( 2001 ) , han et al .
( 2002 ; 2003 ) , podsiadlowski et al .  ( 2003a ) , and nelemans ( 2010 ) .
depending largely upon the mass ratio and the initial orbital period , the mass transfer may be stable and mostly conservative , stable and mostly non - conservative , or dynamically unstable .    in the case of a close binary that comprises a primordial primary star with mass @xmath5 and a secondary star that has @xmath6 , the mass transfer is most likely to be stable .
however , depending on the ratio of thermal timescales of the two stars , the mass transfer may range from highly non - conservative for small @xmath7 ratios to mostly conservative for values of @xmath7 close to unity . whether the orbit ultimately evolves to be longer or shorter than the original orbital period depends on how much specific angular momentum is carried away with the matter ejected from the binary
.    in principle
the best way to model such binary evolution is via a stellar evolution code that can simultaneously handle the nuclear and mass loss / gain evolution of _ both _ stars as well as the orbital evolution .
the recently developed _
code ( paxton et al .  2011 ) is capable of doing such binary evolution studies , but , to our knowledge , no systematic studies have yet been carried out of ordinary close primordial binaries .
we know empirically , and from some theoretical studies , that such evolution can lead to systems such as algol binaries ( see , e.g. , nelson & eggleton 2001 ) where the mass transfer is stable and the envelope of the primary has not yet been exhausted ; sdb binaries ( see , e.g. , han et al .  2002 ; 2003 ) which may have evolved from unstable mass transfer in a common envelope ; systems with hot , bloated , low - mass white dwarfs ( see , e.g. , van kerkwijk et al .  2010 ; carter et al .
2011 ) , where the envelope of the primary has been completely removed ; and even double - degenerate binaries ( marsh et al .  2004
; iben , tutukov , & yungelson 1995 ; nelemans 2010 ) .    in low - mass binaries where the mass transfer proceeds stably
, regardless of whether it is or is not conservative , there is a fairly tight relation between the white dwarf mass , @xmath8 , and the final orbital period , @xmath9 ( see , e.g. , joss & rappaport 1983 ; savonije 1983 ; paczyski 1983 ; savonije 1987 ; rappaport et al .  1995 ; fig .  5 of lin et al .
2011 ; tauris & van den heuvel 2014 ) . for orbital periods
@xmath10 years , and white dwarf masses , @xmath8 , up to 1.4 @xmath0 , in systems where the initial masses of the primaries , @xmath11 , were no larger than @xmath3@xmath12 , the following expression gives the final orbital period in terms of the white dwarf mass : @xmath13 for orbital periods up to @xmath320 days , which apply to systems with primaries of initial masses up to @xmath3@xmath12 , a more recently derived expression ( lin et al .
2011 ) is : @xmath14 these expressions have been checked against a number of binary radio pulsars ( fig .  3 of rappaport et al .  1995 ; fig .  2 of tauris & van den heuvel 2014 ) , and against a number of _ kepler _ binaries containing low - mass hot white dwarfs ( van kerkwijk et al .
2010 ; fig .  6 of carter et al .
2011 ) .
prior to this work , only four low - mass white dwarf systems in the _ kepler _ field that can be used to study this relation between orbital period and white dwarf mass were identified .
they are koi 74 , koi 81 ( rowe et al .
2010 ; van kerkwijk et al .
2010 ) , khwd3 ( koi 1375 ; carter et al .
2011 ) , and koi 1224 ( breton et al .
2012 ) . in this work
we report the discovery of new low - mass , highly thermally bloated white dwarfs in the _ kepler _ binaries kic 9164561 and kic 10727668 .
the orbital periods of these binaries are 1.267 and 2.305 days , the shortest periods among the six _ kepler _ binaries containing low - mass hot white dwarfs .
the white dwarf in kic 9164561 is by far the most thermally bloated of any of the six white dwarfs .    in section [ sec : data ]
we discuss the search that led to these two systems among the _ kepler _ binaries , and indicate how the data were prepared for analysis
. folded light curves are presented in section [ sec : lcs ] along with a simple geometric analysis , including a harmonic decomposition of the out - of - eclipse portions of the light curve .
the masses of the primary a stars are inferred from information on their temperatures and mean densities in sect .
[ sec : astars ] . in sect .
[ sec : rv ] we present the radial velocity and spectral studies of kic 9164561 .
the mass ratio of kic 10727668 is examined in the context of the doppler boosting effect in sect .
[ sec : db ] .
inferences about the mass ratios in the two systems are drawn from evaluating the ellipsoidal light variations in sect .
[ sec : elv ] . in sect .
[ sec : wds ] the system parameters for these two binaries are evaluated with emphasis on the white dwarf properties .
evidence for a third star in kic 9164561 is presented in sect .
[ sec : third ] . finally , in sect .
[ sec : discuss ] we discuss kic 9164561 and kic 10727668 in the context of different binary evolution scenarios and white dwarf cooling models .
in the process of searching for third - body companions to the binary stars in the _ kepler _ sample ( see , e.g. , rappaport et al .  2013 ; see also conroy et al .
2014 ) , we produced both an @xmath15 curve for the eclipses ( `` observed minus computed '' eclipse times ) as well as a folded light curve for each of the @xmath32600 cataloged _ kepler _ binaries ( matijevic et al .
2012 ; kirk et al .
we used the _ kepler _ long - cadence time - series photometric data ( 29.4  min samples ) obtained from quarters 0 through 17 .
for all the known binary systems for which at least one quarter of photometry is available , the version 5.0 fits files , which are available for all quarters , were downloaded from the stsci mast website .
we used the data that had been processed with the pdc - map algorithm ( stumpe et al .  2012 ; smith et al .
2012 ) , which is designed to remove many instrumental artifacts from the time series while preserving any astrophysical variability .
we also made use of the information in the file headers to obtain estimates of the combined differential photometric precision ( cdpp ) .
we then examined by eye all 2600 @xmath15 curves for interesting variations .
some of the results have been reported elsewhere ( rappaport et al .
2013 ; tran et al .
the @xmath15 curve for kic 9164561 is discussed in sect .
[ sec : third ] .
we also searched for folded light curves that ( i ) exhibit a flat - bottomed eclipse that is deeper than the other eclipse , and ( ii ) are associated with hot _ kepler _ targets , i.e. , targets that have a reported @xmath16 k. the deep flat - bottomed eclipse - high temperature combination is a signature of a hot white dwarf in orbit with a normal , f , a , or b star .
in addition to the four known _ kepler _ low - mass bloated hot white dwarfs , we found two others , previously unreported , in the binaries kic 9164561 and kic 10727668 .
the nondegenerate components of kic 9164561 and kic 10727668 are , according to the kepler input catalog , hot with effective temperatures of @xmath38000 k and 10,000 k , respectively , and thus are not susceptible to the starspot modulations that are exhibited by cooler stars .
an examination of the fourier transform of the kic 10727668 flux time series shows no evidence for any periodicities , including those from stellar pulsations , other than that of the orbit . in the case of kic 9164561 ,
the fourier transform of the flux time series shows a set of pulsations at frequencies @xmath17 cycles / day , which does not significantly affect our light curve or timing analysis .
we therefore made no corrections to the data except for normalizing each quarter to its median value .
the orbital periods given in the _ kepler _ binary catalog ( matijevic et al .  2012 ; kirk et al .
2015 ) were used to produce initial folds of the data . from these folds
, we established a trial template for the orbital light curve .
this profile was then used to fit each individual eclipse in the time series , thereby determining a set of 893 and 108 eclipse times for kic 9164561 and 10727668 , respectively .
we then fit a linear function to the sequence of eclipse times of the form : @xmath18 , where @xmath19 is the nth eclipse .
the two orbital periods thus determined are in good agreement with those given in the _ kepler _ binary catalog , but have somewhat smaller uncertainties .
the periods and their uncertainties for these two systems are listed in table [ tab : sys ] .
we also fit the eclipse times for each system including a quadratic term to allow for possible changes in the orbital period or light travel time effects that might occur if the binary is part of a triple system .
these are discussed in sect .
[ sec : third ] below .
the data for both kic 9164561 and kic 10727668 were folded at the orbital periods we determined from fitting the individual eclipses , as described above .
the results were averaged into 3-minute bins , and are shown in figures [ fig:916 ] and [ fig:107 ] .    due to the small curvature observed in the eclipse timing residuals for kic 9164561 ( see sect .  [
sec : third ] ) we removed a quadratic term from the eclipse phases before producing the folded light curve in fig .
[ fig:916 ] . the phase corrections amount to less than 0.002 days ( or 0.0016 orbital cycles ) .    in each case , the light curves are normalized to unity at the bottom of the deeper eclipse when the hot white dwarf has been totally eclipsed .
this is then the light level of the primary a star at a minimum of its ellipsoidal variations .
the relatively small separation of the two components of kic 9164561 is manifest by the short orbital period , the large amplitude of the ellipsoidal light variations ( `` elv '' ) , and the duration of the eclipse relative to the orbital period .
the light curves of these two binaries were analyzed with a relatively simple geometrical model which involves 13 free parameters . in brief
, there are 7 free parameters associated with the eclipses : @xmath20 , @xmath21 , @xmath22 , @xmath23 , two limb - darkening coefficients , and @xmath24 , where subscript `` 1 '' refers to the primary a star , @xmath25 is the orbital separation , @xmath22 the orbital inclination angle , and @xmath24 , the time of the middle of a primary eclipse .
the out - of - eclipse phases are fitted with a set of 6 terms , i.e. , a sine and a cosine term for each of the orbital frequency and its next two higher harmonics .
each of these can be associated with different contributions to the elvs , the illumination term , and the doppler boosting effect , as described in some detail in carter et al .
( 2011 ) .
once the best - fitting parameters were found with a levenberg - marquardt fitting algorithm , the allowed region in parameter space was defined with a markov chain monte carlo ( mcmc ) routine , which used the gibbs sampler and the metropolis - hastings algorithm , and a likelihood proportional to @xmath26 ( see , e.g. , holman et al .
2006 ; sanchis - ojeda et al .
2013 ) .
the results for the fitted model are plotted as red curves in figs .
[ fig:916 ] and [ fig:107 ] .
the high quality of the fit in each case is indicated by the white - noise appearance of the residuals .
the results for the fitted parameters and their uncertainties are summarized in table [ tab : sys ] .
the mass and radius of the primary star in each system are essential to quantitative models of the two binaries .
there are two accessible parameters that can be utilized to estimate the primary mass : the effective temperature , @xmath27 , and the mean stellar density , @xmath28 .
the _ kepler _ input catalog ( kic ) gives @xmath298059 k for the primary star of kic 9164561 and 9696 k for the primary of kic 10727668 .
these temperatures are based on 5-color photometry . for kic 9164561
we measured @xmath30 k from the same spectra used for the radial velocity measurements ( see sect .  [
sec : spec ] ) . in the case of kic 10727668 , where we have no additional spectral information
, we estimate the uncertainty in the kic value of @xmath27 at approximately @xmath31 300 k ( huber et al .
2104 ) .
the second readily , and more accurately , determined parameter associated with each primary star is its mean density , @xmath32 .
this mean density is a function only of the light curve - determined parameters @xmath33 and @xmath34 , with a weak dependence on @xmath35 ( the mass ratio @xmath36 ) .
the dependence of @xmath28 on these parameters may be found by dividing kepler s third law by @xmath37 to obtain : @xmath38 this relation was pointed out by seager & malln - ornelas ( 2003 ) in the context of transiting exoplanets . for exoplanets ,
the value of @xmath35 is much smaller than the error in @xmath34 so that @xmath28 can be determined independently of any mass information . while the same is not quite true with our two binaries , the value of @xmath35 is still sufficiently small , i.e. , @xmath39 that a rather accurate value of @xmath28 may be determined nearly independently of @xmath35 . in particular , by using eqn .
( 3 ) with the results from our photometric analysis of @xmath21 and with , as a first approximation , @xmath40 , we find that @xmath41 g @xmath42 for kic 9164561 , and @xmath43 g @xmath42 for kic 10727668 .
these values are refined below using estimates for @xmath35 that are based on estimates of the masses of the two component stars .
figure [ fig : yy ] shows stellar evolution tracks for stars of a range of initial masses in the @xmath44@xmath45 plane .
the locations of the two target primary stars are indicated by the rectangular boxes .
we conclude from this that the mass of the primary star in kic 9164561 must be @xmath46 ( kic 9164561 ) and that of the primary in kic 10727668 must be @xmath47 @xmath0 .
the spectroscopic follow - up of kic9164561 was carried out with the fibre - fed chelle spectrograph ( fies ; * ? ? ?
* ; * ? ? ?
* ) mounted at the 2.56-m nordic optical telescope ( not ) of roque de los muchachos observatory ( la palma , spain ) .
spectra were taken under clear and stable sky conditions at 11 epochs between may and september 2014 .
we used the @xmath48 _ med - res _
fibre , which provides a resolving power of about 47000 in the spectral range 36007400 .
the observations taken at each epoch were split into two or three consecutive sub - exposures of 1200 seconds each to remove cosmic ray hits .
following the observing strategy described in @xcite , we monitored the radial velocity drift of the instrument by acquiring long - exposed ( t@xmath49=15 - 30 sec ) thar spectra in  sandwich mode " , i.e. , right before and after each epoch observation .
the data were reduced using a customized idl software suite , which includes bias subtraction , flat fielding , order tracing and extraction , and wavelength calibration .
radial velocities ( rvs ) were determined by cross - correlating the target spectra with the spectrum of vega , for which we adopted a heliocentric radial velocity of @xmath50@xcite .
table  [ rv - table ] lists the fies rv measurements , the signal - to - noise ( s / n ) ratio per pixel at 4650   , and the full - width at half maximum ( fwhm ) of the fies cross - correlation function ( ccf ) .
full - widths at half maximum were extracted by fitting a gaussian to the ccfs between the observed spectra and the atlas9 model spectrum having the same photospheric parameters as the target star ( see sect.[sec : spec ] ) and smoothed to the fies spectral resolution .
figure [ fig : rvs ] shows the phase - folded rv curve of kic9164561 , based on the phase established from the _ kepler _ photometry , and the keplerian circular - orbit fit to the results .
the k velocity of the primary a star is 21.64 @xmath31 0.35 km s@xmath51 , and the corresponding @xmath52-velocity is -28.43 @xmath31 0.25 km s@xmath51 .
the corresponding mass function is @xmath53 for an orbital inclinational angle of @xmath54 ( see table [ tab : sys ] ) , this implies : @xmath55 the mass for the a star of @xmath56m@xmath57 obtained above may be combined with the mass function to obtain a white dwarf mass of @xmath58 .
the mass ratio for the kic 9164561 system is therefore @xmath59 .    because kic 10727668 is faint , i.e. , 2.6 _ kepler _ magnitudes fainter than kic 9164561
, we were unable to obtain any useful rv measurements for that system .
[ rv - table ]    ccccccc bjd@xmath60 & rv & @xmath61 & ccf@xmath62 fwhm & @xmath63 & s / n / pixel@xmath64 & date of obser.@xmath65 + ( @xmath66 2450000 ) & & & & & @4650 & ( at start ) + 6783.67056 & -9.397 & 0.673 & 103.6 & 1.4 & 23 & 06 may 2014 + 6833.58987 & -38.015 & 0.714 & 106.7 & 1.2 & 19 & 25 jun 2014 + 6833.67106 & -43.987 & 0.768 & 102.9 & 1.2 & 17 & 25 jun 2014 + 6844.66674 & -7.539 & 0.968 & 112.1 & 1.2 & 11 & 06 jul 2014 + 6845.43519 & -38.741 & 0.794 & 109.1 & 1.0 & 15 & 06 jul 2014 + 6845.71215 & -11.787 & 0.690 & 102.8 & 1.2 & 23 & 07 jul 2014 + 6901.37427 & -22.004 & 0.847 & 107.6 & 1.2 & 15 & 31 aug 2014 + 6902.40137 & -43.375 & 0.743 & 104.6 & 1.3 & 17 & 01 sep 2014 + 6915.38485 & -12.215 & 0.784 & 106.5 & 1.3 & 15 & 14 sep 2014 + 6915.55901 & -7.337 & 0.820 & 107.2 & 1.1 & 14 & 15 sep 2014 + 6920.43515 & -14.633 & 0.791 & 105.6 & 1.3 & 13 & 19 sep 2014 +    \(a ) width of the cross correlation function .
( b ) s / n ratio per pixel at 4650 .
( c ) calendar date in ut at the beginning of the target observation .
we started the spectral analysis by coadding all the fies spectra , after shifting them to a zero - velocity reference frame , which yielded a s / n of @xmath67 55 per pixel at 4650 .
the spectral parameters of the primary component in kic 9164561 were then derived by fitting this spectrum to a grid of synthetic spectra calculated with the spectrum code @xcite using atlas9 models @xcite . we assumed a secondary - to - primary flux ratio of 0.0285 , as estimated from the eclipse light curve modeling ( table  [ tab : sys ] ) .
we adopted the calibration equations from @xcite and @xcite for the purpose of estimating the microturbulent @xmath68 and macroturbulent @xmath69 velocities .
these equations were extrapolated to the effective temperature of the star .
the spectral type was determined following the method described in @xcite .
the estimated stellar parameters are listed in table [ star - parm ] .
figure  [ fies_spectrum ] shows the co - added fies spectrum in the spectral region encompassing the h@xmath70 , h@xmath71 , and mgi b lines , along with the best fitting synthetic model .
lc parameter & + effective temperature @xmath72 ( k ) & @xmath73 + spectroscopic surface gravity log@xmath74 ( log@xmath75 cms@xmath76 ) & @xmath77 + metallicity [ m / h ] ( dex ) & @xmath78 + microturbulent velocity @xmath79 ( kms@xmath51 ) & @xmath80 + macroturbulent velocity @xmath81 ( kms@xmath51 ) & @xmath82 + projected rotational velocity @xmath83sin@xmath22 ( kms@xmath51 ) & @xmath84 + spectral type & 7 +      we looked for the spectroscopic signature of a third ( bound ) component to kic 9164561 by cross - correlating the fies spectra with a set of a- through m - type observed templates . the templates were gathered between 2011 and 2014 , as part of the on - going rv follow - up of _ corot _ and _ kepler _ stars carried out with fies@not ( e.g. , gandolfi et al .
2013 ) .
assuming that the third star is on or near the zams , its mass can not exceed @xmath32 @xmath0  corresponding to an a7v spectral type  without significantly violating a number of observational constraints ( see below ) .
the eclipse timing variations detected in the _ kepler _ light curve ( see sect .
[ sec : third ] ) constrain the orbital period of a putative third star in the system to lie between 8 and 19  years , implying that its rv is within @xmath312 km s@xmath51 of the @xmath52 velocity of the binary system .
the fies rv follow - up covers a time span of about 4.5 months , which corresponds to less than 4.6% of the orbital period of the third component , constraining its rv variation to be less than @xmath34 km s@xmath51 within the observing window . because of the relatively high projected rotational velocity of the primary star ( see table  [ star - parm ] )
, the almost `` stationary '' ccf peak of the third source is expected to be blended with the broad peak of the a7v component .
the fies ccfs show neither a peak with constant rv superimposed on top of the broad peak of the a7v star , nor a skewed profile that would suggest the presence of a another rapidly rotating star ) . ] . in order to assess the significance of the non - detection and constrain the parameters of a possible third body , we built up a set of synthetic double - lined spectroscopic binary ( ` sb2 ' ) spectra using atlas9 model atmospheres .
we adopted gray s ( 2005 ) calibration for dwarf stars to convert spectral types to stellar parameters .
we combined the rotationally broadened model spectrum of an a7v star  resembling the primary component  with that of dwarf stars with different spectral types ( from a7v to m5v ) and values of @xmath83sin@xmath22 .
we added random noise  to account for the average s / n ratio of the fies spectra  and the contaminant light of the second source  to account for the white dwarf .
we assumed a zero velocity for the third light and radial velocities in the range @xmath85rv@xmath86  kms@xmath51 for the a7v primary star .
finally , we cross correlated the combined sb2 spectra with narrow - lined synthetic spectra covering different spectra types , in the same fashion as for the observed spectra .    since the contrast of a ccf peak decreases with increasing @xmath83sin@xmath22 , we considered the following two cases .
_ spectral type later than f5v_. we used gyro - chronology to constrain the @xmath83sin@xmath22 of the third star .
given the age of the system ( @xmath3925  myr ; see sect .
[ sec : cool ] ) , the rotation period of the third component is expected to be longer than about 4  days for spectral later than f5v ( barnes et al .
this places an upper limit of 18 km s@xmath51 on @xmath83sin@xmath22 .
although this value likely overestimates the true @xmath83sin@xmath22 for the late - type stars , we conservatively assumed 18 km s@xmath51 to be the projected rotational velocity associated with the third light .
we found that any stars with spectral type between k2v and f5v would be detectable in the fies ccf ( for an example see fig .
[ fig : simccfs ] ) . for later spectral types
the peak of the third light is embedded in the noise of the ccf , owing to the relatively low s / n ratio of the spectra and the low flux of the third star .    sin@xmath87km
/ s and a k0v star with @xmath83sin@xmath88km / s ( see text for details ) .
the dashed red and blue lines mark the rv of the k0v ( 0 km / s ) and the a7v star ( 21.6 km / s ) , respectively .
the black curve , shown for reference , is the same ccf , but without the presence of the k0v star . ]    _ spectral type between f5v and a7v_. for @xmath83sin@xmath89  km s@xmath51 the ccf peak of the third component is clearly visible on top of the peak of the primary star , as described in the previous case . for @xmath83sin@xmath90  km
/ s the ccf would assume a variable skewed profile , whose shape changes with the changing rv of the primary star , which we do not observe .
finally , for @xmath91sin@xmath92  km / s the two ccf peaks are merged together and larger actual rv variations of the primary component would be needed to account for the measured rvs ( see below ) . in the case of comparable amounts of light
contributed by the a7v star and a putative third star , we can say two things .
first , the fwhm of the ccf would vary in a correlated way with the measured rvs .
the extra widening can be estimated analytically , and the broadening factor is given by @xmath93 $ ] , where @xmath94 is the observed ( i.e. , ` diluted ' ) instantaneous value of the a7v star rv , and @xmath95 is the gaussian width of the ccf curve .
since the ccf fwhm remains constant at 106.5 km s@xmath51 to within an rms scatter of @xmath33 km s@xmath51 ( see table [ rv - table ] ) , and a 42 km s@xmath51 shift ( i.e. , @xmath96 km s@xmath51 ) is expected to widen the ccf by 12 km s@xmath51 , which is distinctly not observed , we conclude that there is no other a- to early f - type star of comparable brightness in the system .
second , as mentioned above , a stationary line would ` dilute ' the actual doppler shift by a factor of @xmath3@xmath97 , where @xmath98 and @xmath99 are the fluxes of the a7v star and the third body , respectively . in this case , the properly inferred mass of the white dwarf ( based on the rv amplitude ) would be a factor of approximately @xmath100 higher .
this , in turn , would yield an implausibly high mass for the white dwarf when @xmath99 exceeds @xmath98 by even a factor of 2 .
therefore , the mass of an early - type third body is independently constrained to be @xmath101 .    thus , based on our simulations
, we can exclude nearly any plausible near - ms third star that has a spectral type earlier than k2v . while there is still some small ` wiggle room ' left in parameter space
, these constraints likely imply an upper limit for the mass of any third star to be about 0.8  @xmath0 .
in principle , the doppler boosting amplitude ( loeb & gaudi 2003 ; van kerkwijk et al .  2010 ) can serve as a surrogate in place of radial velocity ( ` rv ' ) measurements for measuring the binary k velocity . as formulated by carter et al .
( 2011 ) the doppler boosting ( hereafter ` db ' ) amplitude in a binary is given by : @xmath102\ ] ] where the @xmath103 s are the fluxes of the two stars in the _ kepler _ band .
the @xmath104 s are doppler boosting pre - factors calculated from : @xmath105 where @xmath106 is the spectral flux density at frequency @xmath107 , and @xmath104 must be evaluated using stellar atmosphere models averaged over the _ kepler _ band . the second term in eqn .
( [ eqn : db ] ) represents the contribution to the db amplitude from the white dwarf in the system ( note the minus sign ) .
this method for evaluating the k velocities has been verified by ehrenreich et al .
( 2011 ) and bloemen et al .
( 2012 ) for the case of koi 74 , a system quite analogous to the ones discussed in this paper .
_ kic 9164561 _ : the measured doppler boosting amplitude in this system , i.e. , the amplitude of the @xmath108 term , where @xmath109 is the orbital frequency , and @xmath110 is time , is -10 @xmath31 26 parts per million ( ` ppm ' ) . given the @xmath111 velocity obtained from the rv measurements ( see sect .  [
sec : rv ] ) , for nominal values of @xmath112 , @xmath113 for this system should be @xmath67160 ppm . one possible explanation for this discrepancy is based on the large illumination semi - amplitude , i.e. , the coefficient of the @xmath114 term , that is measured to be 2385 @xmath31 36 ppm . in this case , if the phase of the illumination is shifted in the proper sense from its nominal alignment with the eclipses by as little as 4@xmath115 , the entire db effect would be cancelled .
we note that in koi 74 the illumination term is 600 times smaller because of the wider orbital separation and the lower luminosity of the white dwarf , and is therefore very unlikely to materially affect the observed doppler boosting amplitude . in order for the phase of the illumination effect in kic 9164561 to be shifted by the amount required to cancel most of the doppler boosting signal
, the a - star primary would have to be rotating faster than the rate that is synchronous with the orbit , such that some of the heated surface material is carried into the night side of the star .
the measured rotation velocity of this star is @xmath67 75 @xmath31 5 km s@xmath51 , a value that is to be compared with the equatorial rotation velocity of a star of radius 2.5 @xmath116 ( see table [ tab : sys ] ) rotating with a period of 1.267 days , i.e. , @xmath67 100 km s@xmath51 .
thus , if anything , the primary is rotating subsynchronously and the db signal is therefore expected to be enhanced relative to its nominal amplitude .
this is a potentially important discrepancy which merits follow - up study .
_ kic 10727668 _ : in this system , the doppler boosting signal is measured to have a significantly non - negative value , 143 @xmath31 20 km s@xmath51 . the illumination amplitude , 557 ppm , is 4 times smaller than in the case of kic 9164561 ; this system may therefore be less susceptible to an illumination - induced perturbation of the db amplitude . if eqn .
( [ eqn : db ] ) is used to evaluate @xmath111 we find : @xmath117 where we have taken @xmath118 , @xmath119 ( van kerkwijk , private communication ) , and @xmath120 ( see fig .
[ fig:107 ] ) .
this expression for @xmath111 may be entered into the mass function , and the cube root of the result taken to obtain the following equation involving the masses of the two stars : @xmath121 where we have taken the orbital inclination angle to be 86.4@xmath115 ( see table [ tab : sys ] ) .
this equation can be solved numerically for @xmath122 given any specific value for @xmath123 . for a primary star mass of 2.2 @xmath0 ( see fig .
[ fig : yy ] ) , the mass of the secondary is @xmath124 where the error bar takes into account the uncertainty in @xmath123 as well as the uncertainty in @xmath111 .
again we caution that this estimate for @xmath122 does not take into account any systematic error due to an `` amplitude leak '' in the db term .
the ellipsoidal light variations ( ` elvs ' ) seen in the out - of - eclipse region of the light curves also contain information on the mass ratios , @xmath125 . the basic expression for the elvs ( kopal 1959 ;
morris & naftilan 1993 ; carter et al .
2010 ) is : @xmath126 and @xmath127 where @xmath128 is the amplitude of the cosine term at twice the orbital frequency , @xmath129 is the linear limb - darkening coefficient , and @xmath130 is the gravity darkening coefficient .
for the hot primary stars in kic 9164561 and kic 10727668 we take @xmath131 and @xmath132 .
this leads to a numerical value for @xmath133 of 1.68 .
we can then utilize the measured values of @xmath128 of @xmath134 ppm and @xmath135 ppm for the two systems , respectively ( see table [ tab : sys ] ) .
similarly , we take the values for @xmath21 and @xmath22 for each system from table [ tab : sys ] .
finally , we solve for @xmath35 to find : @xmath136 where the error bars are statistical only and do not include any systematic uncertainties in the use of eqn .
( [ eqn : elv ] ) to evaluate @xmath35 from the elv amplitude .
the value of @xmath35 found for kic 9164561 from the rvs is @xmath137 ( sect .
[ sec : spec ] ) which is in quite satisfactory agreement with @xmath138 .    on the other hand , the value of @xmath35 for kic 10727668 found from the elv analysis ( eqn .
( [ eqn : qelv ] ) ) , while still in marginal agreement with the value obtained from the doppler boosting analysis ( @xmath139 ) , would imply a white dwarf mass as high as @xmath140 .
this is quite on the high side of what we expect for this system ( see sect .
[ sec : discuss ] ) .
in order to determine the sizes of the binary systems we need the total masses . in the case of kic 9164561 , where a direct rv measurement has been made ,
the total system mass is 2.22 @xmath31 0.07 @xmath0 . with the application of kepler s third law , this leads immediately to an orbital separation of @xmath141 . for kic 10727668
the doppler boosting analysis yielded a mass ratio @xmath142 .
if we tentatively accept this mass ratio , and apply kepler s third law , we find the orbital separation is @xmath143 .
the reason that the uncertainty is quite small , in spite of the larger uncertainty in @xmath35 for this system , is that @xmath25 depends only on the cube root of the total system mass , and @xmath35 is a small number .
these values for the orbital separation are summarized in table [ tab : sys ] .
the radii of the primary and secondary stars follow immediately from the @xmath21 and @xmath144 values derived from the light curve analyses and the values of the orbital separations .
these radii are summarized in table [ tab : sys ] .
of particular note are the large radii of the two white dwarfs , 0.277 @xmath116 and 0.151 @xmath116 for kic 9164561 and kic 10727668 , respectively , which indicate that these stars are quite thermally bloated .
the degenerate radii for cold white dwarfs of either he or hybrid heco compositions , are given by eqn .
( 4 ) in nelson & rappaport ( 2003 ) .
the degenerate radii for the white dwarfs in kic 9164561 and kic 10727668 are 0.020 @xmath116 and 0.018 @xmath116 , respectively , implying thermal bloating factors of 14 and 8 .
we carried out an eclipse timing analysis for both binaries discussed in this paper .
we used a code to produce the observed minus calculated curve ( ` @xmath15 ' curve ) that was tested in our search for triple star systems ( rappaport et al .
only the @xmath15 curve for kic 9164561 showed an interesting non - linear behavior .
the results are shown in fig .
[ fig : omc ] with black points for the average of the @xmath15 curves for the primary and secondary eclipses .
the @xmath15 curve shows two notable features : ( i ) a clear curvature , and ( ii ) small - amplitude oscillations with a period of @xmath3170 days .
the small amplitude oscillations are due to an artifact that results from a beat between the orbital period and _ kepler _ long - cadence sampling ( see eqns .
( 4 ) and ( 5 ) in rappaport et al .
2013 ) .
the quadratic behavior in the @xmath15 curve is highly significant ( 35 @xmath95 ) , and can be interpreted as either a change in the orbital period , or as a roemer delay due to the motion of the binary induced by a third body .
the first option requires changes in the binary orbital period on a characteristic timescale of @xmath145 years .
since there is no significant mass transfer taking place in the system at the current time , this can not be the cause of the apparent orbital period changes .
one would have to invoke another mechanism , such as the applegate effect ( applegate & patterson 1987 ; applegate 1992 ; kotze & charles 2010 ) to explain the short timescale changes in the orbit .
we focus here on a more interesting and , in our view , more likely scenario , i.e. , that there is a third body in the system .
it is then possible to constrain the orbital radius of such a third body as well as its mass , if we restrict our attention to coplanar and circular orbits .
otherwise , the possibilities are nearly limitless .
the results of such constraints are shown in fig .
[ fig : constraints ] .
we find that the orbital period must lie in the range of @xmath146 years in order to provide at least a somewhat satisfactory fit to the @xmath15 curve ( @xmath147 limits ) .
the blue curve shows the best fit ( projected ) orbital radius as a function of @xmath148 .
the red curve is the corresponding mass function , while the green curve is the mass of the third body under the assumptions that the mass of the inner binary is @xmath149 and that the orbital inclination of the third body orbit is @xmath150 . since the mass of the third body ,
if it is assumed to be on or near the zams , can not exceed @xmath32 @xmath0 without affecting the spectrum , rv analysis , and the third light of the system in ways that are not observed ( see sects .  [
sec : spec ] and [ sec : ccf ] ) , the orbital period is further restricted to the range of @xmath151 years . if we adopt the formal constraint on the mass of a main - sequence third star in the system from sec .
[ sec : ccf ] of @xmath152 we find that : @xmath153 and @xmath154 years .
to this point we have not yet considered the possibility that the putative third body might be a white dwarf . the orbital period
would then be restricted to the range @xmath155 years since the mass would have to be @xmath156 ( see fig .
[ fig : constraints ] ) .
if there had been a prior episode of mass transfer from the progenitor of the outer body white dwarf to the inner binary ( kic 9164561 ; see , e.g. , tauris & van den heuvel 2014 ) then , according to eqn .
( [ eqn : mwd_p ] ) , the mass would lie in the range of @xmath157 .
such white dwarfs would have cooled to the point where they would contribute very little light to the system ( @xmath158 ) .
based on the inferred luminosity of kic 9164561 , and the appropriate bolometric corrections , we infer that its distance from the earth is @xmath159 pc .
typical orbital separations to the putative third body , for periods in the range of @xmath160 years , are expected to be @xmath161 au . at a distance of 390 pc ,
the corresponding angular separations are then @xmath162 mas .
this may be detectable with _
gaia_.    c c c parameter & kic 9164561 & kic 10727668 + [ 0.5ex ] ra ( j2000 ) & 19h 42 m 27.64s & 19h 21 m 45.58s + dec ( j2000 ) & 45@xmath115 30@xmath163 17.0@xmath164 & 48@xmath115 02@xmath163 42.4@xmath164 + @xmath33 [ days ] & 1.2670400 @xmath165 & 2.305897 @xmath166 + @xmath24 [ mjd]@xmath167 & 64.61956 @xmath168 & 1207.3677 @xmath169 + @xmath170 [ mag . ] & 13.71 & 16.32 + @xmath171 & 0.1091 @xmath172 & 0.0893 @xmath173 + @xmath21 & 0.397 @xmath174 & 0.170 @xmath175 + @xmath25 [ @xmath116 ] & 6.40 @xmath177 & 9.95 @xmath31 0.15 + @xmath178@xmath179 & 0.0285 @xmath180 & 0.0245 @xmath180 + inclin .  [ deg ] & 71.59 @xmath181 & 86.45 @xmath182 +   + @xmath123 @xmath183 $ ] & @xmath184 & @xmath185 + @xmath186 [ @xmath116 ] & 2.54 @xmath31 0.03 & 1.69 @xmath31 0.04 + @xmath27 [ k ] & 7870 @xmath31 150 & 9696 @xmath31 300 + @xmath187 [ @xmath188 & 22.3 @xmath31 1.5 & 22.9 @xmath31 3.1 + @xmath189 [ g / cc ] & 0.171 @xmath310.04 & 0.63 @xmath31 0.04 + @xmath190 [ cgs ] & 3.93 @xmath31 0.02 & 4.32 @xmath31 0.04 +   + @xmath122 [ @xmath0 ] & 0.197 @xmath31 0.005 & 0.266 @xmath31 0.035 + @xmath192 [ @xmath116 ] & 0.277 @xmath31 0.005 & 0.151 @xmath31 0.004 + @xmath193 [ k ] & 10410 @xmath31 200 & 14110 @xmath31 440 + @xmath194 [ @xmath188 & 0.82 @xmath31 0.07 & 0.82 @xmath31 0.11 + @xmath195 [ g / cc ] & 12.8 @xmath31 1.2 & 120 @xmath3119 + @xmath196 [ cgs ] & 4.84 @xmath31 0.04 & 5.55 @xmath31 0.06 +   + @xmath197 [ @xmath198@xmath199 & @xmath200 @xmath31 27 & 143 @xmath31 20 + @xmath201 [ @xmath202@xmath203 & 2385 @xmath31 36 & 557 @xmath31 22 + @xmath204 [ @xmath205 & @xmath206 @xmath31 27 & @xmath207 @xmath31 20 + @xmath208 [ @xmath209@xmath210 & @xmath134 @xmath31 35 & @xmath135 @xmath31 21 + @xmath211 [ @xmath212 & @xmath213 @xmath31 28 & @xmath214 @xmath31 19 + @xmath215 [ @xmath216 & @xmath217 @xmath31 36 & @xmath218 @xmath31 23 +    [ tab : sys ]
here we discuss the properties of the kic 9164561 and kic 10727668 systems and some of the evolutionary paths that may have led to their formation .
we can then check whether the expected system parameters match those determined empirically .
there are two important evolutionary pathways that could lead to the formation of low - mass helium white dwarfs in close binaries ( see , e.g. , van kerkwijk et al .
2010 ; carter et al .
2011 ; breton et al .
the first channel includes one episode of common envelope ( ce ) evolution wherein the more massive star of the primordial binary overflows its roche lobe after having become a giant with a degenerate core , and then undergoes dynamically unstable mass transfer onto its binary companion ( see , e.g. , taam et al .
1978 ; webbink 1984 ; taam and bodenheimer 1992 ) .
the rapid rate of mass transfer causes the companion to spiral inward towards the degenerate core of the giant while concomitantly expelling the envelope of the giant into the ism .
the resulting binary is composed of a degenerate white dwarf in a tight circular orbit with a largely unchanged companion star .
the second channel involves stable mass transfer between the components of the progenitor binary .
unfortunately the degree to which mass transfer is conservative and the effectiveness of orbital angular momentum loss mechanisms ( e.g. , magnetic braking ; see , for example , verbunt & zwaan 1981 ) are uncertain . for reasonable physical assumptions in regard to the mass transfer and angular momentum loss processes , this channel can explain the formation and observed properties of these two systems ( see nangreaves et al .
2015 , and references therein ) .    _
( i ) the common envelope channel _
although a physically realistic 3-d hydrodynamic simulation of ce evolution has not yet been carried out , a reasonably good approximation of the effects of the process can presumably be obtained from first principles using energy conservation . as the companion star moves at high velocity through the envelope of the giant
, a fraction @xmath219 of the orbital energy is transferred to the tenuous envelope material .
the companion star then spirals inward and more and more energy and angular momentum are transferred to the giant envelope .
the gravitational potential energy that is released may be sufficient to eject the envelope .
it is assumed that this process occurs on a timescale much shorter than the thermal timescale of the giant star .
one of the more commonly used ce energy constraints was analyzed by dewi & tauris ( 2000 ) .
it is given by @xmath220 \label{eqn : ce}\ ] ] where @xmath221 , @xmath222 , and @xmath223 are the core and envelope masses of the wd progenitor , and the mass of primordial secondary , respectively ( see , e.g. , webbink 1984 ; pfahl et al .
note that @xmath219 is a dimensionless measure of the efficiency of the conversion of orbital energy that goes toward unbinding the envelope of the giant ( the energy required to eject the envelope , @xmath224 , must satisfy @xmath225 ) and is usually taken to be @xmath226 .
the dimensionless factor @xmath227 is intended to take into account the actual binding energy of the giant at the onset of the ce phase .
it largely depends on the ratio of the core mass to the pre - ce - phase mass of the giant .
values of @xmath227 have been tabulated for many different evolutionary states and masses of giants by tauris & dewi ( 2001 ) and podsiadlowski et al .
( 2003b ) . for the types of giants that could produce helium white dwarfs
, @xmath227 is of order unity , normally having a value between 0.5 and 1.0 , which is possible in certain circumstances , this only makes the argument presented here more robust . ] .
also note in equation ( [ eqn : ce ] ) that @xmath221 , the core mass of the giant , is approximately equal to the final white dwarf mass , @xmath8 .
thus the original mass of the primary can be written as @xmath228 and its radius is given by @xmath229 .
finally , @xmath230 is the separation of the white dwarf and secondary star in their mutual circular orbit after the ce phase while @xmath231 is the separation of the components of the progenitor binary . using kepler s third law , equation ( [ eqn : ce ] )
can be rewritten as @xmath232\left(\frac{r_{1,0}}{a_0}\right)=1 \label{eqn : ce2}\ ] ] where @xmath233 is a dimensionless parameter that is @xmath226 . also @xmath234 ) and @xmath235 ) are the mass ratios of the primordial and present day binaries , respectively . at the onset of ( unstable ) mass transfer , the ratio , @xmath236 , is approximately equal to the roche lobe radius of the progenitor primary divided by the primordial separation ; its value is constrained by roche geometry and can be expressed as a simple function of @xmath237 ( see , eggleton 1983 ) . given that the radius of the red giant will have a very strong dependence on @xmath238 and a very weak dependence on @xmath11 , the constraint of roche geometry allows the approximation of @xmath239 by a function of @xmath221 ( see , e.g. , rappaport et al .
1995 ) .    assuming a specific value of @xmath240 , and using the observed orbital period of the binary ( @xmath241 ) and the inferred mass of the present - day companion @xmath123 ( @xmath242 )
, strong limits can be imposed on the possible values of the white dwarf mass ( we take @xmath243 ) , by solving eqn .
( [ eqn : ce2 ] ) implicitly for @xmath8 and applying the constraint that @xmath244 ( i.e. , @xmath245 ) .
it is reasonable to suppose that @xmath240 will lie in the range 0.5 to 1.0 . taking the upper limit of this range for @xmath240 , we find that the only acceptable values of the core mass for kic 9164561 are @xmath246 while for kic 10727668 they are @xmath247 .
assuming the lower limit of this range of @xmath240 , we find that the only acceptable values of the core mass for kic 9164561 are @xmath248 , while for kic 10727668 they are @xmath249 .
these ranges of values , by definition , while the upper limits correspond to @xmath250 to take into account the fact that we are considering only primaries of masses @xmath251 so that only he wds are produced . ] of @xmath221
are clearly incompatible with the estimated values of @xmath8 for both kic 9164561 and kic 10727668 based on the cooling evolution of the white dwarfs ( see the next section ) , as well as with the empirically determined value for kic 9164561 from rvs .
we conclude that for a reasonable range of common envelope efficiencies , the formation of either or both of these systems via the ce channel is rather unlikely .    _
( ii ) the algol channel _
the other channel involves stable mass transfer during the evolution of the progenitor binary . in this scenario ,
the star that is more massive ( the donor ) evolves more quickly and fills its roche lobe on its nuclear timescale before its less massive companion ( the accretor ) can do so . even if the mass ratio is as large as @xmath252 , it is entirely possible for the mass transfer to proceed stably , if the donor has a substantial radiative envelope ( nelson & goliasch 2014 ) . as an example of evolution that could explain the formation of the two observed systems , consider a 1.5 @xmath0 donor that has evolved beyond the tams and lies in the hertzsprung gap when it fills its roche lobe .
this star can stably transfer mass to a 1.2 @xmath0 companion on approximately a nuclear and/or thermal timescale ( case ab evolution ) .
at this mass transfer rate ( of @xmath3@xmath253 yr@xmath51 ) , the accretor is not likely to expand and fill its own roche lobe . during the mass transfer process
, the donor can evolve up the red giant branch and form a substantial degenerate helium core .
since the radius of the giant has a very strong dependence on @xmath221 but a very weak dependence on its total mass , the giant will continuously expand with a concomitant increase in the core mass until the envelope surrounding the core is so depleted ( to a mass of the order of 0.01 @xmath0 ; the remainder having been transferred to the companion star ) that it can not be sustained and ultimately collapses onto the remnant degenerate helium core .
if approximately 1@xmath0 is transferred conservatively to the accretor during this process , the resulting system may then resemble either of the observed ones .
details concerning the possible progenitor evolutions will be discussed in a future paper ( nangreaves et al .
2015 ) .      assuming that stable mass transfer from a red giant led to the formation of the bloated helium wds in both kic 9164561 and kic 10727668
, it is possible to check whether the observed properties of the wds are consistent with cooling models of degenerate helium dwarfs .
moreover , if the properties of either wd are found to be consistent , it is then possible to use it as a chronometer in order to determine the age of the binary after the termination of mass transfer , i.e. , just after the collapse of the giant to form a helium wd ( hewd ) .
the study of the formation and cooling of hewds from progenitor red giants has been carried out by a number of groups ( see , e.g. , podsiadlowski et al .
2003 ; nelson et al .
2004 ; althaus et al .
2005 , 2013 ; istrate et al .
we use the code described by nelson et al .
( 2004 ) to analyze the cooling evolution of very low - mass helium wds .
it is assumed that the giant ( i.e. , the wd progenitor ) had a solar metallicity ( @xmath254 ) , and that its initial mass was between 1.0 and 1.5 @xmath0 .
the initial hydrogen abundance ( by mass ) of the progenitor star while on the zams is taken to be equal to 0.71 and the ratio of the mixing length to the pressure scale - height ratio ( @xmath255 ) is set equal to 1.5 .
the initial orbital separation was carefully selected so that only very low mass helium wds were created .
diffusion ( gravitational settling and radiative levitation ) was not included in the computations because diffusion is a very fragile physical process that can be negated by turbulence in the outer layers of these rapidly rotating ( proto-)white dwarf stars .
moreover , the inclusion of diffusion has little effect on the cooling times for the models of interest ( i.e. , bloated helium wds ) because these systems are very young and have not even evolved to reach their maximum surface temperatures before joining the usual wd cooling track in the hr diagram .
it is only after the wds have passed through their maximum @xmath45 that the compression of the outer hydrogen - rich layers , due to contraction , can cause the residual hydrogen to start to burn , leading to a rapid expansion of the outer layers as a result of a thermonuclear runaway ( tnr ) .    in figure
[ fig : cooling ] , we show the results of the cooling evolution of several models that have masses near or smaller than the tnr lower mass limit , i.e. , @xmath256 .
the evolution of the radius of the wd is shown as a function of its effective temperature , @xmath45 .
when the envelope of the progenitor giant collapses , the radius of the proto - wd continually decreases , as long as no tnrs occur , until it asymptotically reaches its zero - temperature radius .
the surface temperature of the object initially increases ( this increasing temperature defines the proto - wd phase ) until it reaches a maximum value of @xmath45 where it then follows a cooling track and asymptotically approaches zero temperature .    assuming that the white dwarfs in both kic 9164561 and kic 10727668 formed according to the scenario discussed above
, we can infer reasonably precise values of the mass based on the observed radii and @xmath45 values superposed on the cooling tracks ( see fig .
[ fig : cooling ] ) .
we find that the mass of the wd in kic 9164561 is @xmath257 , and that of the wd in kic 10727668 is @xmath258 .
these masses agree with those inferred from the cooling evolutions of althaus et al .
( 2013 ) to within 5 to 10% ( their values of the mass being smaller ) , and to within @xmath259% for the values inferred from the work of istrate et al .
we also carried out this same type of wd cooling calculations with the mesa code ( paxton et al .
2011 ) and found very similar results .
the ages of these objects as measured from the cessation of mass transfer are approximately 925 myr and 950 myr , respectively .
if the companion star in kic 9164561 has a mass of @xmath260 , then the time required to reduce its average density to the inferred value of 0.17 g @xmath42 is in good agreement with the theoretical predictions for the evolution of a solar metallicity star .
the orbital period predicted by our models is approximately 1.7 days , but it should be noted that our assumptions about the degree of mass conservation during the mass - transfer phase and the magnitude of angular momentum losses leading to the formation of the giant are not precisely known , and this uncertainty could have a significant impact on the final value of the orbital period . as for kic 10727668 ,
the predicted orbital period from our models almost matches the observed one .
however , the predicted age and somewhat higher mass of the companion star are very problematic because the observationally inferred average density implies that the companion star has not significantly evolved from its initial position on the zams .
based on the implied cooling age of the system and a mass of 2.2 @xmath0 , we would expect the star to be much closer to the tams with a correspondingly lower density ( contrary to what is implied in figure 3 ) . if the mass of the companion were substantially less ( say , @xmath671.7 @xmath0 ) , then a wd composed primarily of helium and with a mass of @xmath261 would be entirely compatible with the observations .
based on the cooling evolution , we conclude that the theoretical models for the evolution of kic 9164561 lead to properties that agree well with the observationally inferred properties .
the predictions for kic 10727668 are problematic , especially given the implied cooling age of the white dwarf in the system .
for an age in the range of 100200 myr , a white dwarf mass of @xmath262 ( inferred from the doppler boosting effect ; see sect .
[ sec : db ] ) , and @xmath263 , it is very difficult to model the properties with any type of solar metallicity hewd . to understand why this is true , consider the evolutionarily limiting case for which all of the hydrogen has been stripped from the atmosphere and envelope of the hewd .
based on the wd cooling tracks of nelson & davis ( 2000 ) , a 0.27 @xmath0 wd of this type can reach a temperature of 14,000k at either of two ages : ( i ) considerably less than 10 myr ; or , ( ii ) @xmath264 myr .
the luminosity of the wd for the latter case is two orders of magnitude smaller than what is inferred ( see table [ tab : sys ] ) and can be automatically ruled out . for the other solution ( corresponding to the early - wd phase before the maximum @xmath45 is reached ) ,
the luminosity is approximately correct but a cooling age of only several million years after the birth of the wd seems to be unacceptably small based on probability arguments .
if a thick layer of hydrogen were to be added to the envelope , the cooling would be slowed but incompatibilities with the observations would persist .
for example , the luminosity of the hewd while it is still young would be considerably higher and it would likely undergo a tnr . perhaps we are observing the hewd while it is undergoing one of its tnr cycles ; the probability of observing it at this epoch may be unacceptably small .
another possible explanation might be a difference in the actual metallicity from that which we assumed ; however , making the metallicity smaller generally makes the incompatibilities larger .
thus if the mass of the wd in kic 10727668 is actually about 0.27 @xmath0 , then it might have been derived from a higher mass progenitor ( @xmath265 ) or perhaps via a different evolutionary pathway than we have considered here .
we plan to investigate all of these possibilities in a future paper .    for the sake of completeness
, we have included the robust measurements of effective temperatures and radii of eight other eclipsing or tidally distorted , extremely - low mass wds ( found in double degenerate systems ; see hermes et al .
the observations are plotted in fig .
[ fig : cooling ] as open blue circles .
it should be noted that these additional wds have properties that are reasonably similar to the six hot wds found by kepler except that they are typically less bloated .
based on the cooling models , we expect the masses of all of these systems to reside in the range of 0.15 - 0.25 @xmath0 , which is generally consistent with the values reported in hermes et al .
( 2014 ) . finally , in this regard , we note that these white dwarfs have a very different evolutionary scenario from the systems examined in our work in that they are likely the result of two common - envelope episodes .      here
we examine very briefly what the future holds for systems such as koi 74 , koi 81 , koi 1224 , koi 1375 , kic 9164561 , and kic 10727668 .
when these systems evolve to the point where the primary a star fills its roche lobe and commences mass transfer onto the white dwarf , their he core masses will range from 0.19 to 0.27@xmath266 for the shortest ( kic 9164561 ) and longest ( koi 81 ) orbital periods of these systems , respectively ( as estimated using eqn .
( [ eqn : mwd_p2 ] ) ) .
once the mass transfer onto the white dwarf begins , it will be dynamically unstable due to the very large mass ratio ( typically 10:1 ) , and a common envelope will ensue . using eqn .
( [ eqn : ce ] ) , we can compute an approximate expression for the ratio of the final orbital separation to the initial orbital separation , and we find : @xmath267 where @xmath8 is the mass of the hot white dwarfs being studied in these systems , @xmath222 and @xmath221 will be the envelope and core mass of the a star when it evolves to fill its roche lobe , @xmath268 , @xmath240 is the common envelope parameter of order unity discussed above , and @xmath269 is the roche - lobe radius of the a star in units of the orbital separation ( with an expected numerical value of @xmath270 ) .    if we apply eqn .
( [ eqn : aoai ] ) to estimate the final orbital separation that will result from this upcoming common envelope phase , we find that @xmath271 .
the current orbital separations of the six bloated hot white dwarfs in the _ kepler _ field range from 6.4 to 48@xmath272 . from this
we can infer that the final orbital separations will be @xmath273 , for kic 9164561 to koi 81 , respectively .
all of these will be binaries comprised of two roughly equal mass white dwarfs each of @xmath274 , if the core of the a star does not merge with the initial low - mass white dwarf in the system .
the corresponding orbital periods and roche - lobe radii would range from 3/4 to 14 min , and @xmath275 , respectively .
all but the largest of these roche - lobe radii are too small to harbor the white dwarf cores that will be unveiled by the common envelope .
thus , we conclude that the orbital periods of the group of double degenerate systems recently studied by brown et al .
( 2014 ) and hermes et al .
( 2014 ) result from much wider systems .
this is emphasized by the fact that the masses of the primary in the double degenerate binaries in the latter studies are typically much more massive than would result from the six binaries containing bloated hot white dwarfs in the _ kepler _ field .
in this study we report two new low - mass , thermally bloated , hot white dwarfs in binary systems with companion a stars .
the systems are similar in their basic nature to those in koi 74 , koi 81 , koi 1224 , and koi 1375 .
this brings to six the total number of such systems in the _ kepler _ field . of the six systems , the two new systems have the shortest orbital periods at 1.27 d and 2.3 d , while kic 9164561 has the largest , i.e. , the most thermally bloated , of the six white dwarfs
its white dwarf has a radius that is 14 times larger than its fully degenerate radius  it is more than 1/2 @xmath116 in _ diameter _ !
the mass of the white dwarf in kic 9164561 is 0.197 @xmath31 0.005 @xmath0 and is determined from direct rv measurements of the primary a star as well as some inferences about the mass of the a star ( i.e. , the mean density and @xmath45 lead to an estimate for @xmath123 ) . by contrast , the mass of the white dwarf in kic 10727668 is more uncertain at 0.27 @xmath31 0.03 @xmath0 . since there are no rv measurements of this fainter system , we relied on measurements of the doppler boosting effect to determine the orbital velocity of the a star .
the accuracy of this latter mass is not quite sufficient for detailed comparisons with binary evolution or white - dwarf cooling models .
we strongly advocate that a larger telescope be used to determine an accurate rv for this system .
binary evolution scenarios involving two low - mass , e.g. , @xmath67 1.1 and 1.5 @xmath0 , primordial stars that undergo stable , conservative mass transfer can nicely explain the current mass and evolutionary state of both the companion star and the hot white dwarf in kic 9164561 .
for this system , we estimate an age of @xmath3925 myr since the envelope of the white dwarf progenitor was transferred to the current a star .
by contrast , it is more difficult to explain the current state of both the a star and its white dwarf companion in kic 10727668 .
if the white dwarf actually has a mass as high as is indicated by the doppler boosting technique ( i.e. , 0.27 @xmath0 ) then it is difficult to understand how the white dwarf can still be so thermally bloated . if , on the other hand , we take the uncertainty seriously , and allow for a mass as low as , say , 0.22 @xmath0 , then there is still a discrepancy between the inferred cooling age of @xmath67 950 myr , and the apparent location of the a star in the hr diagram which is still near the zams .
the growing statistical ensemble of these systems should make it worthwhile to ( 1 ) make renewed attempts to measure the rvs in four of the remaining bloated hot white dwarf systems in order to better determine @xmath8 , and ( 2 ) directly model the binary evolution of such systems with a code that can evolve two stars _ simultaneously _
, e.g. , mesa ( paxton et al .
2011 ) .
_ note added in manuscript _ : after this work was completed we became aware of an impressive collection of 18 similar bloated hot white dwarf systems found in the wasp photometric database ( maxted et al .
2013 ; 2014 ; and references therein ) .
two of these have accurate measurements of the mass , radius , and @xmath45 for the bloated hot white dwarf and show multi - periodic pulsations .
these are now also shown in fig .
[ fig : cooling ] .      we are grateful to the entire _ kepler _ team for making this study possible .
we thank jon swift for taking exploratory spectra of kic 10727668 to see if rv measurements can be acquired in the future , and dave latham for his early contribution to this project .
the rv observations were made with the nordic optical telescope , operated by the nordic optical telescope scientific association at the observatorio del roque de los muchachos , la palma , spain , of the instituto de astrofisica de canarias .
we gratefully acknowledge that two of the rv points were acquired by pilar montaes - rodrguez , and an additional one was obtained during time allocated to hans j. deeg .
we thank antonio frasca of catania astrophysical observatory for some important tips in combining the model spectra .
mattia gandolfi was quite inspirational for this work .
we are grateful for the cat and not service time awarded to this project , and to the staff members at the not for their valuable and unique support during the observations .
this research was supported in part by the mineco esp2013 - 48391-c4 - 2-r grant .
thanks the natural sciences and engineering research council ( nserc ) of canada for financial support , calcul qubec for providing computing facilities , and ernie dubeau and heather nangreaves for technical assistance .
r.s.o .  acknowledges nasa support through the kepler participating scientist program .
this research has made use of the nasa exoplanet archive , which is operated by the california institute of technology , under contract with the national aeronautics and space administration under the exoplanet exploration program .
the data presented in this article were obtained from the mikulski archive for space telescopes ( mast ) .
stsci is operated by the association of universities for research in astronomy , inc .
, under nasa contract nas5 - 26555 .
support for mast for non - hst data is provided by the nasa office of space science via grant nnx09af08 g and by other grants and contracts .
we made use of j - band images that were obtained with the united kingdom infrared telescope ( ukirt ) which is operated by the joint astronomy centre on behalf of the science and technology facilities council of the u.k .
brown , w.r . ,
kilic , m. , allende prieto , c. , gianninas , a. , & kenyon , s. 2013 , apj , 769 , 66 allende prieto c. , barklem p.  s. , lambert d.  l. , cunha k. , 2004 , a&a , 420 , 183 althaus , l.g . ,
serenelli , a.m. , panei , j.a . ,
crsico , a.h .
, garca - berro , e. , & scccola , c.g .
2005 a&a , 435 , 631 althaus , l.g .
, miller bertolami , m.m . , crsico , a.h .
2013 , a&a , 435 , 631 applegate , j.h .
1992 , apj , 385 , 621 applegate , j.h . , & patterson , j. 1987 , apj , 322 , l99 barnes , s.a .
2007 , apj , 669 , 1167 batalha , n. , et al .
2010 , apj , 713 , l109 bloemen , s. , marsh , t.r . ,
degroote , p. , et al .
2012 , mnras , 422 , 2600 breton , r. , rappaport , s. , van kerkwijk , m.h . , & carter , j.a .
2012 , apj , 748 , 115 bruntt , h. , bedding , t.r . ,
quirion , p .- o . , et al .
2010 , , 405 , 1907 buchhave , l.a . ,
bakos , g.a . ,
hartman , j.d .
, et  al . 2010 , apj , 720 , 1118 carter , j.a . , rappaport , s. , & fabrycky , d. 2011 , apj , 728 , 139 castelli , f. & kurucz , r.l .
2004 , eprint arxiv : astro - ph/0405087 conroy , k.  e. , pra , a. , stassun , k.  g. , et al .
2014 , , 147 , 45 eggleton , p.p .
1983 , apj , 268 , 368 ehrenreich , d. , lagrange , a .- m . , bouchy , f. et al .
2011 , a&a , 525 , a85 .
frandsen , s. & lindberg , b. 1999 , in `` astrophysics with the not '' , proceedings eds : karttunen , h. & piirola , v. , anot .
conf , 71 gandolfi , d. , alcal , j.m .
, leccia , s. , et  al .
2008 , , 687 , 1303 gandolfi , d. , parviainen , h. , fridlund , m. , et al .
2013 , a&a , 557 , 74 gray , r.o . &
corbally , c.j .
1994 , aj , 107 , 742 gray , d.e .
2005 , ` the observation and analysis of stellar photospheres ' , ( uk : cambridge university press ) .
han , z. , podsiadlowski , ph . , maxted , p.f.l . ,
marsh , t.r .
, & ivanova , n. 2002 , mnras , 336 , 449 han , z. , podsiadlowski , ph . , maxted , p.f.l . , & marsh , t.r . 2003 , mnras , 341 , 669 hermes , j.j . ,
brown , w.r . ,
kilic , m. , et al .
2014 , apj , 792 , 39 holman , m.j . ,
winn , j.n . ,
latham , d.w .
2006 , apj , 652 , 1715 huber , d. , silva aguirre , v. , matthews , j. m. , et al .
2014 , apjs , 211 , 2 iben , i. , tutukov , a.v . , & yungelson , l.r .
1995 , apjs , 100 , 233 istrate , a.g . ,
tauris , t.m .
, langer , m. , & antoniadis , j. 2014 , a&a , in press joss , p.c . , & rappaport , s. 1983 , nature , 304 , 419 kotze , m.m . , & charles , p.a .
2010 , mnras , 402 , l16 kirk , b. , pra , a. et al .
2015 , in preparation .
lin , j. , rappaport , s. , podsiadlowski , ph .
, nelson , l. , paxton , b. , & todorov , p. 2011
, apj , 732 , 70l loeb , a. , & gaudi , b.s .
2003 , apj , 588 , l117 matijevic , g. , prsa , a. , orosz , j.a . ,
welsh , w. , bloemen , s. , & barclay , t. 2012 , aj , 143 , 123 marsh , t.r . ,
nelemans , g. , & steeghs , d. 2004 , mnras , 350 , 113 maxted , p.f.l . ,
serenelli , a.m. , miglio , a. , et al .
2013 , nature , 498 , 463 .
maxted , p.f.l . ,
serenelli , a.m. , marsh , t.r . ,
cataln , s. , mahtani , d.p . , & dhillon , v.s .
2014 , mnras , 444 , 208 morris , s.l . , & naftilan , s.a .
1993 , apj , 419 , 344 .
nangreaves , h. , nelson , l. et al .
2015 , in preparation nelemans , g. 2010 , ap&ss , 329 , 25 nelson , l. & rappaport , s. 2003 , apj , 598 , 431 nelson , c. & eggleton , p.p .
2001 , apj , 552 , 664 nelson , l.a . , dubeau , e. , & maccannell , k.a .
2004 , apj , 616 , 1124 nelson , l.a .
, & goliasch , j. 2014 , in preparation paxton , b. , bildsten , l. , dotter , a. , herwig , f. , lesaffre , p. , & timmes , f. 2011 , apjs , 192 , 3 paczyski , b. 1983 , nature , 304 , 421 pfahl , e. , rappaport , s. , & podsiadlowski , ph .
2003 , 597 , 1036 podsiadlowski , ph .
, han , z. , & rappaport , s. 2003a , mnras , 340 , 1214 podsiadlowski , ph . , rappaport , s. , & han , z. 2003b , mnras , 341 , 385 rappaport , s. , di stefano , r. , & smith , j.d .
1994 , apj , 426 , 692 rappaport , s. , podsiadlowski , ph .
, joss , p.c .
, de stefano , r. , & han , z. 1995 , mnras , 273 , 731 rappaport , s. , deck , k. , levine , a. , borkovits , t. , carter , j. , el mellah , i. , sanchis - ojeda , r. , & kalomeni , b. 2013 , apj , 768 , 33 rowe , j. , borucki , w. , koch , d. , et al .
2010 , apj , 713 , 150 sanchis - ojeda , r. , rappaport , s. , winn , j.n . ,
levine , a. , kotson , m.c .
, latham , d.w . , & buchhave , l.a .
2013 , apj , 774 , 54 savonije , g.j . 1983 ,
nature , 304 , 422 savonije , g.j .
1987 , nature , 325 , 416 seager , s. , & malln - ornelas , g. 2003 , apj , 585 , 1038 smith , j.c . ,
et al . 2012 , pasp , stumpe , m.c . ,
van cleve , j.e . , et al .
124 , 1000 stumpe , m.c . ,
smith , j.c . ,
van cleve , j.e .
, et al . 2012 , pasp , 124 , 985 taam , r.e . ,
bodenheimer , p. , & ostriker , j.p .
1978 , apj , 222 , 269 taam , r.e . & bodenheimer , p. 1992 , in `` nato advanced research workshop on x - ray binaries and the formation of binary and millisecond radio pulsars '' , p.  281
tauris , t.m . & dewi , j.d.m .
2001 , a&a , 369 , 170 tauris , t.m . &
van den heuvel , e.p.j .
2014 , apj , 781 , 13 telting , j.h . , avila , g. , buchhave , l. , et al .
2014 , an , 335 , 41 tran , k. , levine , a. , rappaport , s. , borkovits , t. , csizmadia , sz . , & kalomeni , b. 2013 , apj , 774 , 81 .
van kerkwijk , m.h . ,
rappaport , s. , breton , r. , justham , s. , podsiadlowski , ph .
, & han , z. 2010 , apj , 715 , 51 verbunt , f. , & zwaan , k. 1981 , a&a , 100 , l7 webbink , r.f . 1984 , apj , 277 , 355 yi , s. , demarque , p. , kim , y .- c . , lee , y .- w . , ree , c. h. , lejeune , t. , & barnes , s. 2001 , apjs , 136 , 417