detached , double - lined eclipsing binaries provide a largely model - independent means to precisely and accurately measure fundamental stellar properties , particularly masses and radii . in the best - observed systems
the precision of these can be at the @xmath7 per cent level , and thus place stringent constraints on stellar evolution models ( e.g. @xcite ) .    despite this , as far as we are aware , there are at present only four known systems with one or more components between @xmath8 and the hydrogen burning limit on the main sequence ( i.e. old , field stars ) : cm  dra @xcite , cu  cnc  b @xcite , lp  133 - 373 @xcite , and sdss  j031824 - 010018 @xcite .
although jw  380 @xcite , nsvs  02502726b @xcite , and the ngc  1647 system of @xcite also satisfy this mass criterion , these objects are still contracting on the pre  main - sequence .
of the known systems , the only two with parameters determined to better than @xmath9 per cent are cm  dra and cu  cnc . for sdss
j031824 - 010018 , the knowledge of the parameters is limited largely by radial velocity errors , since this system is extremely faint ( sdss @xmath10 ) and has an extremely short orbital period ( @xmath11 ) , meaning long integrations can not be used to obtain better signal - to - noise .
hence there is little possibility for substantial improvement in the parameters in the near future .
it is clear that , in order to better - constrain the stellar mass - radius relation on the main sequence , more bright , low - mass eclipsing binaries are needed to yield extremely precise masses and radii .
furthermore , observations of these systems have indicated significant discrepancies with the stellar models .
this is particularly the case below @xmath12 , and the components of cm  dra ( the lowest - mass system with @xmath13 observational errors ) have radii @xmath14 per cent larger than the theoretical predictions from state - of - the - art stellar evolution models .
it has been suggested ( e.g. @xcite ) that the reason for this discrepancy may be that close binaries are not actually representative of single - stars at this level of precision .
the effect of the close companion and tidal locking is likely to significantly increase activity levels in close binaries , and these authors suggest that it could be this effect that is responsible for the inflated radii of cm  dra .
@xcite examined the available sample of single - star and eclipsing binary measurements to search for correlations of radius with activity levels and metallicity , finding that such a correlation of activity with radius does appear to exist for members of close binary stars .
this result was based on a small number of measurements for late m - dwarfs , and the conclusions would be strengthened by the availability of additional precise estimates of the parameters for binaries with a range of activity levels and orbital periods , to explore the available parameter space .
we present the discovery of a new bright , low - mass eclipsing binary system , gj 3236 .
the components show high activity levels , as evidenced by h@xmath6 emission lines in the spectra of both stars , and show rotational modulations in the light curves that are synchronized with the binary orbital period of @xmath15 .
this system thus has the potential to yield an additional precise test of stellar evolution models for high activity and short orbital periods .
eclipses in gj 3236 were detected during 2008 january in the first two weeks of routine operations of the newly - commissioned mearth observatory , a system designed primarily to search for transiting super - earth exoplanets orbiting around the nearest @xmath16 mid to late m - dwarfs in the northern hemisphere ( @xcite ; @xcite ) .
exposure times on each field observed by mearth are tailored to achieve sensitivity to a particular planet size for the assumed stellar parameters of the target star , and were @xmath17 for gj 3236 .
after the initial detection of the eclipses , we switched to a follow - up mode , observing at the highest possible cadence ( i.e. continuously , resulting in a cadence of approximately @xmath18 including overheads ) during eclipse and for an additional @xmath19 window surrounding the eclipse , and at the normal cadence for the mearth survey of approximately @xmath20 for the remainder of the time . in this way
, we sample as well as possible the rapid flux decrement during eclipse , and the out of eclipse portions of the light curve .
we note that sampling the latter is of crucial importance for obtaining accurate parameters from light curve analyses of binary systems showing out of eclipse modulations , such as the present one .
we obtained data both in the initial season after the eclipses were first detected ( 2008 january to may , inclusive ) , and following the summer monsoon ( 2008 october onwards ) . during this time , there were a number of software improvements ( including the ability to update the telescope pointing after taking each science image , resulting in substantially improved light curves due to minimizing the drift across the detector , and hence the effect of flat fielding errors , fringing , etc . ) and our filter system was changed from cousins @xmath21 to a custom long - pass filter with transmission from @xmath22 , and limited at the long wavelength end by the tail of the ccd quantum efficiency curve ( see @xcite ) .
we have elected to use only the latter data - set for modeling , despite this having the disadvantage of being observed in a non - standard bandpass , because it is approximately contemporaneous with our @xmath23-band follow - up photometry .
we see evidence for evolution of the amplitude and phase of the out of eclipse modulation between the two observing seasons , so this simplifies the analysis by allowing the use of a single set of spot parameters to describe both bandpasses .
a total of @xmath24 observations were obtained between ut 2008 october 2 and ut 2008 december 11 , including 6 primary and 3 secondary eclipses sampled at high cadence .
data were reduced using the standard mearth reduction pipeline , which is at present identical to the monitor project pipeline described in @xcite .
we used a smaller aperture radius of @xmath25 ( @xmath26 ) than the usual value of @xmath27 adopted for bright stars , due to the presence of a nearby star approximately @xmath28 from gj  3236 at the epoch of our observations , and @xmath29 fainter in the mearth bandpass . this star does not share a common proper motion with gj  3236 and is therefore not physically associated . reducing the aperture size in this way
reduces the flux contributed by this star in the light curve of gj  3236 to a negligible level ( we estimate a @xmath30 error in the measured fluxes , and thus also in the eclipse depths ) .
the high proper motion of the gj  3236 system allows us to constrain the contribution of any additional background stars in the photometric aperture that are not co - moving with gj  3236 itself , using previous epochs of imaging .
we show in figure [ third_light_images ] a series of three images centered on the position of the photometric aperture used for gj  3236 in the mearth images , demonstrating that there are no such objects . if any additional third light is indeed present , it must therefore be co - moving with gj  3236 , meaning it is highly likely to be physically associated , given the extremely small probability of such a chance alignment .
( radius ) photometric aperture used to derive our light curves .
data are from the first and second epoch palomar sky surveys as provided by the digitized sky survey ( top and center panels ) , and the mearth stacked master image ( bottom panel ) .
the approximate epochs of the images are @xmath31 ( poss-1 ) , @xmath32 ( poss-2 ) , and @xmath33 ( mearth ) .
all three panels have the same center , scale and alignment on - sky , with north up and east to the left , covering @xmath34 in the horizontal direction.,width=268 ]    the photometry was calibrated using observations of standard star fields from @xcite taken automatically each night by the mearth observation scheduling software . by fitting all the standard star observations from several photometric nights , we derived the following color equation to convert between the mearth bandpass and the standard johnson - cousins system : @xmath35    the full mearth light curve is reproduced in table [ mearth_lc_table ] .
lrrrrrrr    @xmath36 & @xmath37 & @xmath38 & @xmath39 & @xmath40 & @xmath41 & @xmath42 & @xmath43 + @xmath44 & @xmath45 & @xmath46 & @xmath47 & @xmath48 & @xmath49 & @xmath50 & @xmath51 + @xmath52 & @xmath53 & @xmath46 & @xmath54 & @xmath55 & @xmath56 & @xmath57 & @xmath58 + @xmath59 & @xmath60 & @xmath46 & @xmath61 & @xmath62 & @xmath63 & @xmath64 & @xmath65 + @xmath66 & @xmath67 & @xmath46 & @xmath68 & @xmath69 & @xmath70 & @xmath71 & @xmath72 +      observations centered around the primary eclipse of ut 2008 september 25 and the secondary eclipse of ut 2008
september 27 were obtained using the keplercam instrument on the flwo @xmath73 telescope .
we used the standard binning @xmath74 readout mode , since the plate scale of @xmath75 per unbinned pixel significantly oversamples the typical seeing at flwo .
the resulting plate scale was @xmath76 per summed pixel .
we used the @xmath23 filter and an exposure time of @xmath77 .
we also used observations of standard star fields from @xcite taken on ut 2008
september 28 to calibrate the keplercam photometry onto the standard johnson - cousins system .
a total of @xmath78 observations were taken , where possible starting @xmath79 before first contact and finishing @xmath79 after last contact to sample the out - of - eclipse portions of the light curve and thus allow the eclipse depths to be properly measured .
the end of the primary eclipse was missed due to twilight .
these photometric data were reduced using the same pipeline as described in  [ mearth_obs_sect ] .
we used an aperture radius of @xmath80 binned pixels , corresponding to @xmath81 on - sky ( recalling from
[ mearth_obs_sect ] that the size of this aperture is important due to the nearby star ) .
the difference in magnitude between gj  3236 and the fainter star is only @xmath82 in @xmath23-band , which leads to approximately a @xmath83 error in the measured fluxes of gj  3236 , and thus a similar error in the measured eclipse depths .
since the point spread functions on our images are not very well - behaved it will be challenging to reduce this using psf - fitting photometry , and it is still substantially smaller than the other sources of error in the present light curve models , so we have not pursued this issue further .    the full @xmath23-band light curve is reproduced in table [ v_lc_table ] .
lrrrrrrr    @xmath84 & @xmath85 & @xmath86 & @xmath87 & @xmath88 & @xmath89 & @xmath90 & @xmath91 + @xmath92 & @xmath93 & @xmath94 & @xmath95 & @xmath96 & @xmath97 & @xmath98 & @xmath99 + @xmath100 & @xmath101 & @xmath94 & @xmath102 & @xmath103 & @xmath104 & @xmath105 & @xmath106 + @xmath107 & @xmath108 & @xmath94 & @xmath109 & @xmath110 & @xmath111 & @xmath90 & @xmath112 + @xmath113 & @xmath114 & @xmath94 & @xmath115 & @xmath116 & @xmath117 & @xmath118 & @xmath119 +      observations in the near - infrared @xmath120-band were obtained using the robotic peters automated infrared imaging telescope ( pairitel ) from ut 2008 february 17 to ut 2008 march 3 ( inclusive ) .
we scheduled observations around the times of primary and secondary eclipse , obtaining @xmath121 data points spread over @xmath122 .
exposure times were @xmath123 .
the observations were automatically scheduled , collected , and reduced by the fully robotic pairitel observing system @xcite .
we then produced differential photometry using a set of comparison stars chosen from the 2mass catalog .
we estimate a photometric error of approximately @xmath124 per data point from the scatter of the out - of - eclipse measurements . the full @xmath120-band light curve is reproduced in table [ j_lc_table ] .
lrrrrrrr    @xmath125 & @xmath126 & @xmath127 + @xmath128 & @xmath129 & @xmath127 + @xmath130 & @xmath131 & @xmath127 + @xmath132 & @xmath133 & @xmath127 + @xmath134 & @xmath135 & @xmath127 +      spectroscopic observations were obtained using the tres fiber - fed chelle spectrograph on the flwo @xmath136 tillinghast reflector .
we used the medium fiber ( @xmath137 projected diameter ) throughout , yielding a resolving power of @xmath138 .
table [ spec_obs_table ] summarizes these measurements , including the radial velocities of both components of the binary derived therefrom .
the tres instrument is extremely stable , so we acquired thar wavelength calibration exposures before or after the target exposures , rather than using the simultaneous calibration fiber , which can lead to contamination of the target spectrum in the red - most orders from the very strong ar lines in the thar lamp spectrum .
a second fiber was placed on sky , but this was not used due to a problem causing extremely poor throughput in the tres sky fiber , which has since been resolved .
lrrrrrr    @xmath139 & & & @xmath140 & @xmath141 & @xmath142 & @xmath143 + @xmath144 & & & @xmath145 & @xmath146 & @xmath147 & @xmath148 + @xmath149 & & & @xmath150 & @xmath151 & @xmath152 & @xmath153 + @xmath154 & @xmath155 & @xmath156 & @xmath157 & @xmath158 & @xmath159 & @xmath160 + @xmath161 & @xmath162 & @xmath163 & @xmath164 & @xmath165 & @xmath166 & @xmath167 + @xmath168 & @xmath169 & @xmath170 & @xmath171 & @xmath172 & @xmath173 & @xmath174 + @xmath175 & @xmath176 & @xmath177 & @xmath178 & @xmath179 & @xmath180 & @xmath181 + @xmath182 & & & @xmath183 & @xmath184 & @xmath185 & @xmath186 + @xmath187 & & & @xmath188 & @xmath189 & @xmath190 & @xmath191 + @xmath192 & & & @xmath193 & @xmath194 & @xmath195 & @xmath196 + @xmath197 & & & @xmath198 & @xmath199 & @xmath200 & @xmath201 + @xmath202 & @xmath203 & @xmath204 & @xmath205 & @xmath206 & @xmath207 & @xmath208 + @xmath209 & @xmath210 & @xmath211 & @xmath212 & @xmath213 & @xmath214 & @xmath215 + @xmath216 & @xmath217 & @xmath218 & @xmath219 & @xmath220 & @xmath221 & @xmath222 + @xmath223 & @xmath224 & @xmath225 & @xmath226 & @xmath227 & @xmath228 & @xmath229 + @xmath230 & & & @xmath231 & @xmath232 & @xmath233 & @xmath234 + @xmath235 & & & @xmath236 & @xmath237 & @xmath238 & @xmath239 + @xmath240 & & & @xmath241 & @xmath242 & @xmath243 & @xmath244 + @xmath245 & @xmath246 & @xmath247 & @xmath248 & @xmath249 & @xmath250 & @xmath251 + @xmath252 & @xmath253 & @xmath254 & @xmath255 & @xmath256 & @xmath257 & @xmath258 + @xmath259 & @xmath260 & @xmath261 & @xmath262 & @xmath263 & @xmath264 & @xmath265 + @xmath266 & @xmath267 & @xmath268 & @xmath269 & @xmath270 & @xmath271 & @xmath272 + @xmath273 & @xmath274 & @xmath275 & @xmath276 & @xmath277 & @xmath278 & @xmath279 + @xmath280 & @xmath281 & @xmath282 & @xmath283 & @xmath284 & @xmath285 & @xmath286 + @xmath287 & @xmath288 & @xmath289 & @xmath290 & @xmath291 & @xmath292 & @xmath293 + @xmath294 & @xmath295 & @xmath296 & @xmath297 & @xmath298 & @xmath299 & @xmath300 + @xmath301 & @xmath302 & @xmath303 & @xmath304 & @xmath305 & @xmath306 & @xmath307 + @xmath308 & @xmath309 & @xmath310 & @xmath311 & @xmath312 & @xmath313 & @xmath314 + @xmath315 & @xmath316 & @xmath317 & @xmath318 & @xmath319 & @xmath320 & @xmath321 + @xmath322 & @xmath323 & @xmath324 & @xmath325 & @xmath326 & @xmath327 & @xmath328 + @xmath329 & @xmath330 & @xmath331 & @xmath332 & @xmath333 & @xmath334 & @xmath335 + @xmath336 & & & @xmath337 & @xmath338 & @xmath339 & @xmath340 +    spectra were reduced using standard procedures in iraf @xcite from the echelle package to extract the spectra from the target fiber and sky fiber simultaneously .
the tres ccd is read out using two amplifiers , which necessitates combining them before tracing and extraction .
this was done using the mscred package in iraf , applying a multiplicative gain correction to equalize the difference in gain between the two readout electronics chains .
we then divided by a normalized flat field exposure to correct for fringing in the red - most regions of the spectrum , which has a significant effect on the orders we used for the radial velocity analysis .    at the time of our observations , the tres detector suffered from a very high energetic particle hit rate ,
the source of which is under investigation .
the particle hits were removed by median - combining each set of multiple exposures for the flat fields , target and calibrations , and then using the statistics of the median frame to apply an upper envelope clip to the target frames themselves .
this allowed us to avoid having to combine our target exposures to eliminate the particle hits , which would smear out the radial velocity variations slightly as a result of the extremely short orbital period of gj  3236 .
the stability of the tres instrument allowed us to define the aperture trace for each fiber using the high signal - to - noise flat fields , which were then used to extract the target and calibration spectra , employing the `` optimal '' ( in a least squares sense ) weighting scheme of @xcite .
the spectra were wavelength - calibrated using the standard ecidentify procedure , which employs the ( known ) dispersion relation of the chelle to solve for a wavelength solution over all orders of the spectrum simultaneously .
this is necessary due to the relative paucity of sufficiently intense lines from the thar lamp in each order of the spectrum at wavelengths @xmath341 .
table [ photparams ] summarizes the known system properties , from the literature ( principally the proper motion survey of @xcite used to select the target stars for the mearth survey ) and our own @xmath23 and mearth photometry , converted to the standard johnson - cousins system .
lr    @xmath342 & @xmath343 + @xmath344 & @xmath345 + @xmath346 & @xmath347 + @xmath348 & @xmath349 +   + mearth & @xmath350 +   + @xmath351 & @xmath352 + @xmath353 & @xmath354 +   + @xmath355 & @xmath356 + @xmath357 & @xmath358 + @xmath359 & @xmath360 +    the average color of @xmath361 indicates an average spectral type of approximately m4 using the color to spectral type relation of @xcite for young galactic disk stars .
for detached eclipsing binaries with circular orbits , the radial velocity ( rv ) and light curve models are largely independent .
we therefore carried out a preliminary analysis of the available light curves before starting to obtain spectroscopic observations .
the principal purpose of doing this was to determine an extremely precise orbital period to better - target the radial velocity observations , but we can also constrain the system eccentricity by using the phase of the secondary eclipses relative to the primary eclipses ( related to @xmath362 ) .
a simple geometric argument can then be used to obtain a limit on the eccentricity itself , assuming no _ a priori _ knowledge of @xmath363 .
we used jktebop ( @xcite ; @xcite ) , a modified version of ebop ( eclipsing binary orbit program ; @xcite ; @xcite ) , to obtain these parameters by fitting the mearth light curve .
the program was modified to fit simultaneously for the eb model , and a simple form for the synchronized out - of - eclipse modulations , assuming they can be approximated by a truncated fourier series : @xmath364 where @xmath365 , and @xmath366 is the orbital period of the binary .
the normalization term is omitted since this is already taken into account by the standard ebop model .
the revised code yields a fit with very small residuals for the present system , and the use of the form in eq .
( [ ooe_eq ] ) keeps the number of parameters required to describe the modulations to a minimum compared to a full spot model .
doing so will not necessarily reproduce the correct eclipse depths for a given set of physical parameters ( or vice versa ) , but should be adequate for deriving an accurate system ephemeris .
we relax this assumption in  [ wd_sect ] , where we adopt a full spot model for determining the geometric and radiative parameters of the system .
the following parameters were allowed to vary in the fit : the sum of the radii divided by the semimajor axis , @xmath367 ( where @xmath368 , @xmath369 is the stellar radius of star @xmath370 and @xmath371 is the semimajor axis ) , orbital inclination @xmath372 , @xmath362 , central surface brightness ratio @xmath120 ( essentially equal to the ratio of the primary and secondary eclipse depths ) , orbital period @xmath366 , orbital phase zero - point @xmath373 , light curve normalization , and the parameters @xmath374 , @xmath375 , @xmath376 and @xmath377 describing the out of eclipse modulations .
single - band light curves showing grazing eclipses constrain the ratio of the radii only very weakly ( or equivalently , the luminosity ratio ) , so these parameters and the mass ratio ( used for computing the oblateness of the stars and the reflection effect ) were fixed at @xmath378 for the initial analysis .
the results for the orbital parameters were found not to change significantly upon revising these to the measured values from the spectroscopy once they were available .
our light curves are not of sufficient quality to fit for the stellar limb darkening , so we assumed a square root limb - darkening law with coefficients appropriate to the @xmath21 ( cousins ) passband from @xcite using the phoenix model atmospheres .
given the lack of a spectroscopic constraint on the effective temperature , surface gravity and metallicity for either star , we have instead assumed solar metallicity , and derived temperatures and surface gravities of @xmath379 and @xmath380 , and for the secondary , @xmath381 and @xmath382 by iterating the modeling process .
the assumed temperatures were derived using the masses from the combined orbital and light curve solution , interpolating between the compilation of values for field stars by @xcite .
see also  [ kin_sect ] where these temperatures are verified using the final model .
a gravity darkening exponent of @xmath383 , a value appropriate for stars with convective envelopes @xcite , was also assumed , and the option in jktebop to calculate the reflection effect was used rather than fitting for it .
we report only the orbital parameters ( @xmath366 , @xmath373 and @xmath362 ) from the ebop analysis , and have adopted these for the remainder of this work .
the parameter uncertainties were derived using a monte carlo algorithm built - in to jktebop @xcite .
briefly , this algorithm uses the best fit to generate a synthetic light curve , injecting gaussian noise with amplitude determined by the observational errors ( which were scaled such that the reduced @xmath384 of the fit was equal to unity ) to produce a simulated light curve , which is then fit to determine a new set of parameters .
the parameter uncertainties can be estimated using the distribution of the parameters from a large number of realizations ( here , @xmath385 ) of this process .
see also @xcite for more details .
times of minimum light derived for the @xmath386 well - sampled eclipses in our @xmath23 and mearth observations are reported in table [ minlight ] .
these were estimated using the method of @xcite over a @xmath387 region in normalized orbital phase around each eclipse .
we subtracted the truncated fourier series in eq .
( [ ooe_eq ] ) from the light curves before applying this analysis to reduce the effect of any systematic errors arising from the out - of - eclipse modulation .
ll    @xmath366 & @xmath388 + @xmath373 & @xmath389 +   + @xmath362 & @xmath390 + @xmath391 & @xmath392 ( 95% cl ) + & @xmath393 ( 99% cl ) +   + @xmath6 & @xmath394 +   +   +   + @xmath395 & @xmath396 + @xmath397 & @xmath398 + @xmath399 & @xmath400 +   + @xmath401 & @xmath402 + @xmath403 & @xmath404 + @xmath405 & @xmath406 + @xmath407 & @xmath408 +   + @xmath409 & @xmath410 + @xmath411 & @xmath412 +   +   +   + @xmath395 & @xmath413 + @xmath397 & @xmath414 + @xmath399 & @xmath415 +   + @xmath401 & @xmath416 + @xmath403 & @xmath417 + @xmath405 & @xmath418 + @xmath407 & @xmath419 +   + @xmath409 & @xmath420 + @xmath411 & @xmath421 +   +   +   + @xmath395 & @xmath422 + @xmath397 & @xmath423 + @xmath399 & @xmath424 +   + @xmath401 & @xmath425 + @xmath403 & @xmath426 + @xmath405 & @xmath427 + @xmath407 & @xmath428 +   + @xmath409 & @xmath420 + @xmath411 & @xmath429 +    lrrrl    @xmath430 & @xmath431 & @xmath432 & @xmath23 + @xmath433 & @xmath434 & @xmath435 & @xmath23 + @xmath436 & @xmath437 & @xmath438 & mearth + @xmath439 & @xmath440 & @xmath441 & mearth + @xmath442 & @xmath443 & @xmath444 & mearth + @xmath445 & @xmath446 & @xmath447 & mearth + @xmath448 & @xmath449 & @xmath450 & mearth + @xmath451 & @xmath452 & @xmath453 & mearth + @xmath454 & @xmath455 & @xmath456 & mearth + @xmath457 & @xmath458 & @xmath459 & mearth + @xmath460 & @xmath461 & @xmath462 & mearth +
radial velocities were obtained using the two - dimensional cross - correlation algorithm todcor @xcite , which uses templates matched to each component of a spectroscopic binary to simultaneously derive the velocities of both stars , and importantly for the present application , an estimate of their luminosity ratio ( @xmath6 ) in the spectral bandpass .
we performed both a standard absorption line based cross - correlation analysis , and an emission line analysis using the strong h@xmath6 emission observed in both components of gj  3236 .
the signal to noise ratio of the latter is substantially higher in our data , but we caution that the source of the h@xmath6 emission in our target stars is not well - understood
. the emission line analysis results could suffer systematic errors depending on the relative velocities of the regions of the chromosphere emitting the h@xmath6 and the stellar photospheres , and this is exacerbated in the present system because it is highly likely that the stellar spin and binary orbit are synchronized , both from the light curves and expectations from tidal theory @xcite .
we therefore conservatively adopt the results from the absorption line analysis despite the larger errors .    the absorption line analysis used as a template spectrum a single epoch observation of barnard s star ( @xcite ; also known as gj  699 ) taken on ut 2008 october 20 .
we used a wavelength range of 8700 to 8850  in order 50 of the spectrum to derive the velocities , since this region contains a number of reasonably strong metallic lines in m - dwarfs , and is free of telluric absorption lines .
the ca  ii infra - red triplet at 8498 , 8542 , and 8662  was deliberately avoided since the absorption in these lines is saturated , and they are therefore intrinsically very broad , and emission cores are often observed in these lines in active stars , which would severely complicate the use of an inactive field star as a template .
we assumed a barycentric radial velocity of @xmath463 for barnard s star , derived from presently unpublished cfa digital speedometer measurements spanning @xmath464 .
for the emission line analysis , we used as templates spectra of two m4 dwarfs that were found to display h@xmath6 emission , gj  856a and gj  3379 , observed as part of another program using the same instrument .
gj  856a was observed on ut 2008 october 17 and gj  3379 was observed on ut 2008 october 19 .
we used a wavelength range of 6520 to 6645  in order 37 for the emission line cross - correlations , noting that the cross - correlations are dominated by the h@xmath6 lines .
barycentric radial velocities of gj  856a and gj  3379 were determined by cross - correlation with barnard s star , using the 8700 to 8850  region as for our absorption line analysis of gj  3236 .    given that todcor provides a measure of the correlation between the observed and template spectra , it is possible to run cross - correlations against a series of templates and use the correlation value as an indicator of the template match . both the effective temperature and the rotational velocity of the templates affect the solution , so it is vital to explore this parameter space . due to the availability of a limited range of observed templates , the effective temperature could not be varied , but we were able to rotationally broaden the template spectra in order to create primary and secondary templates with @xmath465 .
we used linear limb darkening coefficients @xmath466 and @xmath467 from @xcite for the sdss @xmath468 passband ( a close match to the effective wavelength of the spectroscopy ) , which correspond to our assumed temperatures for the primary and secondary ( see  [ ebop_lc_sect ] ) .
we ran todcor on the grid of @xmath469 template combinations , and sought the maximum of the resulting correlation values to determine the most likely @xmath470 values for the two components of the binary system .
barnard s star is a very slow rotator , so the rotation of the template is negligible given our spectral resolution .
our analysis indicates @xmath471 , which is consistent with the expected value of @xmath472 if the spin and orbit are synchronized .
@xmath473 was poorly - constrained , with the maximum correlation occurring at @xmath474 , whereas we would expect @xmath475 assuming synchronization .
we have simply adopted the latter value , since the light curves and primary velocity indicate that the system is indeed synchronized . using the templates with @xmath471 and @xmath476
, todcor indicates @xmath477 .    for the emission line analysis , we assumed the stellar rotation values from the absorption line analysis , rotationally broadening the observed spectra of gj  856a and gj  3379 .
the intrinsic rotational broadening in both of these templates is again negligible for our purposes .
todcor gave a maximum correlation for @xmath478 , which provides a measure of the relative h@xmath6 emission luminosities of the two stars .
since the orbital period and time of primary eclipse are extremely well - determined from the light curves , we held these values fixed in the final radial velocity orbital solution .
we also assumed a circular orbit .
figure [ rv_fit ] shows the resulting radial velocity curves .
the parameters derived from the radial velocity analysis are given in table [ orbparams ] .
we assumed equal observational errors on each radial velocity point , scaled such that the reduced @xmath384 of the fit was equal to unity .
radial velocities of the components of double - lined spectroscopic binaries are prone to a `` peak pulling '' effect as they approach the @xmath395 velocity ( at times of minimum separation of the two components in the spectrum ) .
the todcor method can largely eliminate this , provided the template spectra match the target sufficiently well . in our experience
, this effect is particularly sensitive to errors in the assumed rotational broadening : if the template contains too much or too little rotation , the velocities would be pushed toward or away from @xmath395 to compensate .
examining the residuals of our solution in figure [ rv_fit ] indicates little or no `` peak pulling '' , which further justifies our assumed rotational velocities and the assumption of negligible rotation in the templates .
in order to properly account for spots on the components of gj  3236 , we use the most recent version ( wd2007 ) of the standard @xcite code to derive the geometric and radiative parameters of the system .
this code incorporates a full physical model including spots and using model atmospheres .
wd normally uses the @xcite model atmospheres to calculate the emergent flux from surface elements on each binary component .
these models span an effective temperature range of @xmath479 , and for temperatures below @xmath480 , wd interpolates between a black body of the appropriate temperature and the @xmath480 model atmosphere .
this is weighted in such a fashion as to produce a pure black body at @xmath481 , and for intermediate temperatures there is a smooth ramp toward a pure model atmosphere at @xmath480 .
molecular sources of opacity increasingly begin to dominate the optical spectra of m - dwarfs in this temperature range , so a black body is clearly a poor approximation for the emergent flux .
since both components of gj  3236 appear to be cooler than @xmath480 , we have modified our copy of wd to substitute phoenix nextgen model atmospheres @xcite in place of these interpolations for effective temperatures cooler than @xmath480 .
our assumptions for the wd light curve model are detailed in table [ wd_pars ] .
we have assumed zero orbital eccentricity , synchronized spin and orbit for both components of the binary , and no third light , which are justified in  [ mearth_obs_sect ] and  [ ebop_lc_sect ] .
llll    & detached +   + @xmath373 & hjd0 & @xmath482 ( hjd ) & phase zero - point + @xmath366 & period & @xmath483 & orbital period + @xmath484 & dpdt & @xmath485 & first derivative of period +   + @xmath403 & sma @xmath486 & @xmath404 & projected semimajor axis + @xmath401 & rm & @xmath402 & mass ratio + @xmath395 & vga & @xmath396 & systemic velocity +   + @xmath391 & e & @xmath485 & orbital eccentricity + @xmath487 & f1 & @xmath378 & primary rotation parameter + @xmath488 & f2 & @xmath378 & secondary rotation parameter +   + & hla & varied & light curve normalization +   + @xmath489 & tavh & @xmath490 & primary effective temperature + @xmath491 & tavc & varied & secondary effective temperature + @xmath492 & phsv & varied & primary surface potential + @xmath493 & pcsv & varied & secondary surface potential +   + @xmath372 & incl & varied & orbital inclination + @xmath494 & el3 & @xmath485 & third light +   + @xmath495 & alb1 & @xmath496 & primary bolometric albedo + @xmath497 & alb2 & @xmath496 & secondary bolometric albedo + @xmath498 & gr1 & @xmath499 & primary gravity darkening exponent + @xmath500 & gr2 & @xmath499 & secondary gravity darkening exponent +   + @xmath501 & xlat1 & @xmath502 & spot 1 latitude + @xmath503 & xlong1 & varied & spot 1 longitude + @xmath504 & radsp1 & varied & spot 1 radius + @xmath505 & temsp1 & 0.96 & spot 1 temperature contrast +   + @xmath506 & xlat2 & @xmath502 & spot 2 latitude + @xmath507 & xlong2 & varied & spot 2 longitude + @xmath508 & radsp2 & varied & spot 2 radius + @xmath509 & temsp2 & 0.96 & spot 2 temperature contrast +   + @xmath6 & & @xmath394 & light ratio +    the limited phase coverage and lack of repeat observations in our @xmath23-band photometry ( leading to an uncertainty in any corrections for color - dependent `` second - order '' atmospheric extinction ) , and poor precision ( large scatter ) in the @xmath120-band photometry means in practice that these are not useful in aiding the determination of the geometric ( and thus physical ) parameters of the system . our experiments with the wd model indicate that these provide minimal additional constraints on any of the parameters when compared to using the mearth data alone .
we have therefore elected to use only the mearth data for fitting for the system parameters , and then evaluated this model with respect to the other passbands as a check .
doing so necessitates fixing one of the component effective temperatures , since the relative eclipse depths measured in a single passband merely determine the relative temperatures of the two components .
we have therefore assumed @xmath379 following  [ ebop_lc_sect ] . in practice
even given multi - band photometry it is usually necessary to fix this parameter : although the color - dependence of the eclipse depths does in principle provide information on the value of @xmath489 , this is normally poorly constrained due to the difficulty of measuring differential colors to high precision .    in order to reproduce the out - of - eclipse variations seen in our optical light curves , at least two spots must be introduced into the model .
the configuration of these spots affects the measured eclipse depths , depending on whether a spotted or unspotted part of the photosphere is eclipsed .
this predominantly affects the derived ratio of surface brightnesses as measured by the ratio of depths of the primary and secondary eclipses ( and thus the ratios of the component radii and effective temperatures inferred from this ) , and a combination of @xmath510 and the orbital inclination inferred from the total eclipse depth .
since our @xmath23-band data do not cover the portions of the light curve out of eclipse , and the scatter in the @xmath120-band data is too large to usefully constrain the model , the spot modeling must be done using single - band light curve information .
spots are introduced in the wd code using a standard four - parameter model , where each spot is parametrized using the spot latitude @xmath511 and longitude @xmath512 , radius @xmath513 and temperature contrast @xmath514 .
the latter quantity gives the ratio of effective temperatures of the spotted and unspotted photosphere .
for a single - band light curve , the spot longitude is well - constrained by the phase of the out of eclipse modulation , but the combination of @xmath511 , @xmath513 and @xmath514 is largely degenerate ( these quantities determine the amplitude and shape of the variation ) . moreover
, we can place each spot on either star .
although the data do rule out some of the possible spot locations , there are still a range of possibilities , which introduces a systematic error into the geometric parameters of the system . in principle , measuring the color - dependence of the out - of - eclipse modulation would allow this degeneracy to be reduced .
given the presently available data , we instead probe this systematic error by considering three characteristic spot configurations that are consistent with the light curve .
these are : ( 1 ) one close to polar spot on each star , with @xmath515 , with both spots located in the same hemisphere as the chord traversed during the eclipse , resulting in the spot on the primary being eclipsed by the secondary during primary eclipse , and vice versa during secondary eclipse ; ( 2 ) likewise , only with the spots located in the opposite hemisphere , resulting in no spots being eclipsed ; and ( 3 ) a case intermediate between these , with a `` symmetric '' spot configuration of two spots on each star , at latitudes of @xmath502 with one located in each hemisphere . in each model , the spot latitudes and temperature ratios were held fixed , allowing only longitude and radius to vary , and in configuration ( 3 ) we enforced symmetry of the spot pattern about the equator , i.e. both spots on one star were forced to have the same longitudes and radii .
this means that all three models have the same number of parameters ( @xmath516 ) .      to derive robust error estimates ,
including the correlations between parameters , we have implemented a markov chain monte carlo ( mcmc ) algorithm around the wd light curve model .
following @xcite , we used the metropolis - hastings algorithm @xcite to estimate the _ a posteriori _ joint probability distribution of the fitted model parameters .
the fitting statistic was the traditional @xmath517 , where the observational error estimates were scaled such that each light curve s reduced @xmath384 was equal to unity .
we briefly summarize the metropolis - hastings algorithm here .
starting from an initial point in parameter space , the algorithm takes the most recent set of parameters and perturbs one or more parameters by a random gaussian deviate .
if the perturbed parameter set has a lower @xmath384 than its progenitor , it is accepted as a new point in the chain .
if it has a larger @xmath384 , it is accepted with a probability @xmath518 . if it is not accepted , the original point is repeated in the chain .
the size of the perturbations were adjusted so that @xmath519 of the proposed points were accepted .
we included the spectroscopic light ratio between the primary and secondary components as an observation in our @xmath384 statistic , which in practice was implemented by adding an extra contribution to @xmath384 of the form : @xmath520^{2}\ ] ] where @xmath521 and @xmath522 are the observed light ratio and its error , taken from table [ wd_pars ] , and @xmath523 is the value calculated from the wd model .
this is necessary for systems with near - circular orbits and grazing eclipses because the light curves only very weakly constrain this quantity ( or equivalently , the ratio of the component radii ) .
uninformative ( uniform ) priors were assumed on all other parameters .
the final parameters and error estimates were determined by taking the @xmath524 confidence interval ( @xmath525 and @xmath526 percentiles , corresponding to @xmath527 standard deviation for a gaussian distribution ) of the samples generated by our mcmc procedure .
the chains were run until they contained at least @xmath528 samples , and we discarded the first @xmath529 of the points in each chain in order to ensure they had converged .
these parameters are reported in table [ params ] , and figure [ lc_model ] shows the light curves with the best - fitting model overplotted .
llll    & @xmath530 & @xmath531 & @xmath532 + @xmath533 & @xmath534 & @xmath535 & @xmath536 + @xmath492 & @xmath537 & @xmath538 & @xmath539 + @xmath493 & @xmath540 & @xmath541 & @xmath542 + @xmath372 & @xmath543 & @xmath544 & @xmath545 +   + @xmath503 & @xmath546 & @xmath547 & @xmath548 + @xmath504 & @xmath549 & @xmath550 & @xmath551 + @xmath507 & @xmath552 & @xmath553 & @xmath554 + @xmath508 & @xmath555 & @xmath556 & @xmath557 +   + @xmath510 & @xmath558 & @xmath559 & @xmath560 + @xmath561 & @xmath562 & @xmath563 & @xmath564 +   + @xmath565 & @xmath566 & @xmath567 & @xmath568 + @xmath569 & @xmath570 & @xmath571 & @xmath572 + @xmath573 & @xmath574 & @xmath575 & @xmath576 + @xmath577 & @xmath578 & @xmath579 & @xmath580 +   + @xmath371 & @xmath581 & @xmath582 & @xmath583 + @xmath584 & @xmath585 & @xmath586 & @xmath587 + @xmath588 & @xmath589 & @xmath590 & @xmath591 + @xmath592 & @xmath593 & @xmath594 & @xmath595 + @xmath596 & @xmath597 & @xmath598 & @xmath599 +   + @xmath600 & @xmath601 & @xmath602 & @xmath603 + @xmath604 & @xmath605 & @xmath605 & @xmath605 + @xmath606 & @xmath607 & @xmath608 & @xmath609 + @xmath610 & @xmath611 & @xmath612 & @xmath611 +     and @xmath120 passbands ( top , middle and bottom panels , respectively ) . in each panel ,
the upper sub - panels show the light curve ( black points ) with the wd model @xmath613 ( see table [ params ] ) overplotted ( solid lines ) , and the lower sub - panels show the residuals after subtracting the model from the data .
the uppermost panel for the mearth data shows a binned light curve ( in 250 bins each spanning @xmath614 in normalized orbital phase ) to allow the features to be more clearly distinguished given the large number of data points .
, title="fig:",width=316 ]   and @xmath120 passbands ( top , middle and bottom panels , respectively ) . in each panel , the upper sub - panels show the light curve ( black points ) with the wd model @xmath613 ( see table [ params ] ) overplotted ( solid lines ) , and the lower sub - panels show the residuals after subtracting the model from the data .
the uppermost panel for the mearth data shows a binned light curve ( in 250 bins each spanning @xmath614 in normalized orbital phase ) to allow the features to be more clearly distinguished given the large number of data points .
, title="fig:",width=316 ]   and @xmath120 passbands ( top , middle and bottom panels , respectively ) . in each panel , the upper sub - panels show the light curve ( black points ) with the wd model @xmath613 ( see table [ params ] ) overplotted ( solid lines ) , and the lower sub - panels show the residuals after subtracting the model from the data .
the uppermost panel for the mearth data shows a binned light curve ( in 250 bins each spanning @xmath614 in normalized orbital phase ) to allow the features to be more clearly distinguished given the large number of data points .
, title="fig:",width=316 ]      the most serious limitations in the present analysis result from the use of only single - band light curves to model the out of eclipse variations .
this is clear from comparing the results for the three solutions reported in table [ params ] , where we find significant ( @xmath615 standard deviation ) differences in the orbital inclination @xmath372 , ratio of effective temperatures @xmath533 , and in @xmath510 , which depend on the assumptions made regarding the locations of the spots on the components of the binary .    given our radial velocity errors , the orbital inclination uncertainty has little effect on the final physical parameters , as evidenced by the component masses and semimajor axis reported in table [ params ] , which are essentially identical to within the observational errors .
the dominant uncertainty in the final solution is therefore in the component radii .    of the three models considered
, the symmetric spot distribution in model ( 3 ) seems the most physically realistic given the known properties of spots on low - mass stars .
these are often found to be polar in light curve models ( e.g. @xcite ) .
moreover , doppler imaging studies ( e.g. @xcite ) indicate that in reality , the surfaces of low - mass stars probably have many small spots distributed over a range of latitudes .
this situation is largely indistinguishable from a small number of large , polar spots in one - dimensional light curve data , although it is important to note that the effects on the measured eclipse depths _ could _ be different .
there is no a - priori reason to expect the spots to be concentrated into one hemisphere given the close to edge - on inclination of the binary orbit .
we therefore favor model ( 3 ) , but report all three solutions to provide a realistic evaluation of the systematic error in our results , and reiterate that the data do not distinguish between the three possibilities , with all of these having comparable @xmath384 values and residuals .    as a check ,
we compare the model fit to the mearth data alone with our @xmath23 and @xmath120-band data in figure [ lc_model ] . in order to do this ,
we have refit the model for these bands allowing only the light curve normalization parameter ( hla ) to vary , with all other parameters fixed to the values determined from the mearth data .
the model is consistent with the @xmath120-band data within the scatter , and is reasonably consistent with the @xmath23-band data , considering that the early parts of the secondary eclipse curve were taken during quite non - photometric conditions , and that we have not accounted for color - dependent or differential atmospheric extinction effects , which are expected to be larger in this band than the redder mearth or @xmath120 bandpasses , and can not readily be determined due to the lack of repeat observations .
the presence of spots and proximity effects in very close binary systems have a small influence on the shape of the radial velocity curve , causing it to depart from a keplerian orbit as we assumed in  [ rv_sect ] .
we have investigated the influence of these effects on our results by comparing radial velocity curves generated using the wd model and the simple keplerian model , finding that the corrections to radial velocities taken out of eclipse are dominated by the spots , the properties of which are largely unknown .
however these are very small ( @xmath616 ) compared to the uncertainties in our radial velocity measurements ( approximately @xmath617 ) , so we are therefore justified in neglecting them for the present analysis .
provided the effective temperature of one component of an eclipsing binary is known , the combined light curve and radial velocity parameters are then sufficient to determine intrinsic , bolometric luminosities of both components of the system .
this follows directly from the definition of the effective temperature in terms of the stefan - boltzmann law : @xmath618    by using bolometric corrections and the measured system magnitudes , we can then infer the distance , provided the reddening can be constrained or assumed to be zero . in the present case ,
the latter is a reasonable assumption , since our target is very nearby ( @xcite give a `` photometric parallax '' of @xmath619 , and @xcite gives @xmath620 ; note that both of these assume the system is a single star , which means they underpredict the distance to a near equal luminosity binary such as the present one by a factor of approximately @xmath621 ) .    in order to determine the effective temperatures
, we assume the intrinsic colors and bolometric corrections for m - dwarfs compiled by @xcite and @xcite to convert the integrated system @xmath622 color presented in table [ photparams ] into @xmath489 , assuming the measured effective temperature ratio and radii for both components from the light curve model .
we used the polynomial fits of @xcite , which cover the required spectral type range , and assume an uncertainty of @xmath623 in these fits as stated in her appendix c. we also assume a @xmath624 systematic uncertainty in the effective temperature scale , which has been included and propagated in the errors we report .
table [ kinparams ] gives our derived parameters for the gj  3236 system .
llll    @xmath489 & @xmath625 & @xmath626 & @xmath625 + @xmath491 & @xmath627 & @xmath628 & @xmath629 + @xmath630 & @xmath631 & @xmath632 & @xmath633 + @xmath634 & @xmath635 & @xmath636 & @xmath637 +   + @xmath638 & @xmath639 & @xmath640 & @xmath641 + @xmath642 & @xmath643 & @xmath644 & @xmath645 +   + @xmath646 & @xmath647 & @xmath648 & @xmath649 + @xmath650 & @xmath651 & @xmath652 & @xmath653 +   + @xmath654 & @xmath655 & @xmath656 & @xmath657 + @xmath23 & @xmath658 & @xmath659 & @xmath660 + @xmath661 & @xmath662 & @xmath663 & @xmath662 +    the refined value of @xmath489 reported here differs by approximately @xmath664 from the assumed value in  [ ebop_lc_sect ] used to determine limb darkening coefficients .
this discrepancy is @xmath665 , and should have a negligible effect on the parameters determined from the light curve analysis , so we elected not to repeat this using the updated value of @xmath489 given the computational expense involved .    given the full system kinematic information ( position , proper motions , @xmath395 velocity and parallax from the eb analysis ) we can also calculate the @xmath666 components of the space velocity relative to the sun .
we use the method of @xcite , but adopt the definition that positive @xmath654 values are away from the galactic center ( a left - handed coordinate system ) for better consistency with the literature .
these are included in table [ kinparams ] . comparing to the velocity ellipsoids derived by @xcite , gj
3236 appears to be a member of the galactic disc , and lies within the bounds of the young galactic disc population as defined by @xcite .
comparing to the @xmath666 space motions of various solar neighborhood moving groups ( e.g. @xcite ) , the velocities are consistent with with those of the hyades moving group , which has @xmath667 @xcite . despite this possible membership in the hyades group ,
we note that the dispersion in stellar parameters for the group members does not provide a very useful constraint on age or metallicity for the present system .
@xcite obtain a metallicity of @xmath668 } = -0.09 \pm 0.17 $ ] for the hyades group .
figures [ mrrel ] and [ mtrel ] show the position of gj  3236 relative to other known objects and theoretical models of the m - dwarf mass - radius and mass - effective temperature relationships .    .
black points with error bars show the two components of the gj  3236 system , for all three solutions reported in table [ params ] .
the gray points are known systems from the literature ( @xcite ; @xcite ; @xcite ; @xcite ) .
the black line shows the theoretical mass - radius relationship from @xcite for an age of @xmath669.,width=316 ]    , but showing the mass - effective temperature relationship .
again , all three solutions are shown for gj  3236 , but the differences are relatively minor in this plane .
we also include several more massive m - dwarf ebs in this figure for comparison : yy  gem @xcite , gu  boo @xcite , and ogle bw3 v38 @xcite .
the black line shows the theoretical relationship from @xcite for an age of @xmath670.,width=316 ]    although many of the best - observed systems are found to exhibit larger radii than the models predict , the slope of the mass - radius relation between the components of each system ( or equivalently , the ratio of the component radii ) is generally found to agree well between the observations and the theoretical models . examining the positions of our solutions for gj  3236 , our models 1 and 3 for the out of eclipse modulations are reasonably consistent with this expectation , whereas the slope defined by the components for solution 2 is substantially different .
although this indicates the assumptions used in model 2 may be unrealistic , the present observational data for gj  3236 do not distinguish between the three possibilities .
using this as a constraint in the modeling would be dangerous because the results would then no longer yield a completely independent test of the theoretical models , so instead we suggest that a better solution to the issue would be obtaining improved multi - band photometry to better - constrain the spot parameters .
therefore , within the present observational uncertainties , the radii of gj  3236 appear to be consistent with the theoretical predictions at approximately @xmath5 , when we account for the systematic uncertainties represented by the three solutions we have presented .
furthermore , the central values appear to lie above the theoretical curve .
this is largely consistent with the results for the well - known systems in the literature , where the radii are typically found to be @xmath671 larger than the theoretical predictions .    in the effective temperature versus mass plane ( figure [ mtrel ] ) ,
the effective temperatures of the components of gj  3236 are found to be cooler than the models predict , by approximately @xmath672 .
this is in good agreement with many of the well - known systems in the literature as shown in the figure , such as cm  dra .
compared to the well - known systems cm  dra and cu  cnc , gj  3236 has a shorter orbital period , and the h@xmath6 emission in the spectroscopy and out - of - eclipse modulations we observe in the light curve are indicative of high activity levels , with both the photometric period and the spectroscopic line broadening apparently consistent with the stellar spin being synchronized to the binary orbit , as expected from tidal theory ( e.g. @xcite ) .
@xcite propose two hypotheses to explain the observed radius discrepancies between theoretical models and eclipsing binaries : ( 1 ) that the inflated radii result from reduced convective efficiency due to high rotation rates and large magnetic fields ; or ( 2 ) that magnetic spot coverage of the surfaces leads to reduced heat flux , and thus larger radii and cooler effective temperatures .
hypothesis ( 2 ) has been found to be consistent with some of the well - known eclipsing binary systems , e.g. yy  gem @xcite , where the radius discrepancy ( and the discrepancy in effective temperature ) can be explained by the presence of starspots covering approximately @xmath673 of the stellar surface ( if the spots are cooler than the photosphere by @xmath674 ; @xcite ) .
gj  3236 shows out of eclipse modulations of comparable amplitude to yy  gem , a compatible radius discrepancy within the present observational errors , and similarly has effective temperatures somewhat cooler than the models predict , so this is an attractive hypothesis . by obtaining precise , multi - band photometry covering the entire orbital phase
, it may be possible to constrain the spot temperatures , and hence observationally test this argument .
the main limitations in the orbital and geometric modeling of the system arise from the lack of out - of - eclipse data in multiple passbands , which leads to large errors ( @xmath675 ) on the radii when we take the degeneracies in the spot configuration into account .
the precision of the mass estimates ( @xmath4 ) is limited by the error in the radial velocity measurements , and an important contribution to the uncertainties in the radii will be from the spectroscopic light ratio if the uncertainty in spot parameters can be resolved .
further high - precision photometric measurements and high - resolution spectroscopy are therefore needed , and combined with careful analysis these should allow the precision of the mass and radius measurements for this system to be improved , potentially to beyond the @xmath124 level as for the well - known systems cm  dra and yy  gem .    the mearth team gratefully acknowledges funding from the david and lucile packard fellowship for science and engineering ( awarded to dc ) .
this material is based upon work supported by the national science foundation under grant number ast-0807690 .
sq and dwl acknowledge support from the nasa kepler mission under cooperative agreement ncc2 - 1390 .
gt acknowledges partial support from the nsf through grant ast-0708229 .
we thank eric mamajek for helpful discussions regarding solar neighborhood kinematics , and the `` astro - comb '' team at cfa : claire cramer , david phillips , alexander glenday , chih - hao li , dimitar sasselov and ronald walsworth , for their generous donation of engineering time to obtain spectra of our target .
the referee is thanked for comments that helped to improve the manuscript .
this research has made extensive 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 nasa and the nsf , nasa s astrophysics data system ( ads ) , and the simbad database , operated at cds , strasbourg , france .
the digitized sky surveys were produced at the space telescope science institute under u.s .
government grant nag w-2166 .
the images of these surveys are based on photographic data obtained using the oschin schmidt telescope on palomar mountain and the uk schmidt telescope .
the plates were processed into the present compressed digital form with the permission of these institutions .
the peters automated infrared imaging telescope ( pairitel ) is operated by the smithsonian astrophysical observatory ( sao ) and was made possible by a grant from the harvard university milton fund , the camera loan from the university of virginia , and the continued support of the sao and u.c .
the pairitel project and jsb are further supported by nasa / swift guest investigator grant nng06gh50 g .
we thank m. skrutskie for his continued support of the pairitel project .
andersen , j. 1991 , , 3 , 91 baraffe , i. , chabrier , g. , allard , f. & hauschildt , p. h. 1998 , , 337 , 403 barnard , e. e. 1916 , , 29 , 181 barnes , j. r. , collier cameron , a. , james , d. j. & donati j .- f .
2000 , , 314 , 162 barnes , j. r. & collier cameron , a. 2001 , , 326 , 950 bessell , m. s. 1991 , , 101 , 662 bessell , m. s. 1995 , in the bottom of the main sequence any beyond , ed . c.  tinney ( springer , berlin , germany ) , p.123 bessell , m. s. & brett , j. m. 1988 , , 100 , 1134 blake , c. h. , et al .
2005 , , 435 , 181 blake , c. h. , torres , g. , bloom , j. s. , & gaudi , b. s. 2008 , , 684 , 635 bloom , j. s. , starr , d. l. , blake , c. h. , skrutskie , m. h. , & falco , e. e. 2006 , in asp conf .
351 , astronomical data analysis software and systems xv , ed . c. gabriel et al .
( san francisco : asp ) , 751 cakirli , o. , ibanoglu , c. & gungor , c. 2009 , new astronomy , in press .
arxiv : 0902.3327 chabrier , g. , gallardo , j. , baraffe , i. 2007 , , 472 , 17 chiba , m. , & beers , t. c. 2000 , , 119 , 2843 claret , a. 2000 , , 363 , 1081 claret , a. 2004 , , 428 , 1001 delfosse , x. , forveille , t. , mayor , m. , burnet , m. , & perrier , c. 1999 , , 341 , 63 eggen , o. j , & sandage , a. 1967 , , 148 , 911 etzel , p. b. 1980 , ebop user s guide , 3rd edn . ,
ucla astronomy & astrophysics ford , e. b. 2005 , , 129 , 1706 gliese , w. & jahreiss , h. 1991 , preliminary version of the third catalogue of nearby stars , astron .
rechen - institut , heidelberg hastings w. k. 1970 , biometrika , 57 , 97 hauschildt , p. h. , allard , f. & baron , e. 1999 , , 512 , 377 hebb , l. , wyse , r. f. g. , gilmore , g. , & holtzman , j. 2006 , , 131 , 555 hewett , p. c. , irwin , m. j. , bunclark , p. , bridgeland , m. t. , kibblewhite , e. j. , he , x. t. , smith , m. g. 1985 , , 213 , 971 hillenbrand , l. a. 1997 , , 113 , 1733 irwin , j. , et al .
2007a , , 380 , 541 irwin , j. , irwin , m. , aigrain , s. , hodgkin , s. , hebb , l. , & moraux , e. 2007b , , 375 , 1449 irwin , j. , charbonneau , d. , nutzman , p. , & falco e. 2008 , to appear in the proceedings of iau symposium 253 : `` transing planets '' , ed . f. pont et al .
( cambridge : cambridge university press ) .
arxiv : 0807.1316 johnson d. r. h. , & soderblom , d. r. 1987 , , 93 , 864 joy a. h. & sanford r. f. 1926 , , 64 , 250 kopal , z. 1959 , semi - detached binaries and stellar evolution , astronomical society of the pacific leaflets , vol . 8 , p.81 kurucz , r. l. 1979 , , 40 , 1 kurucz , r. l. 1993 , vizier on - line data catalog : vi/39 .
kwee , k. k. & van woerden h. 1956 , , 12 , 327 lacy , c. h. 1977 , , 218 , 444 landolt , a. j. 1992 , , 104 , l340 leggett , s. k. 1992 , , 82 , 351 lpine , s. 2005 , , 130 , 1680 lpine , s. , & shara , p. 2005
, , 129 , 1483 lpez - morales , m. & ribas , i. 2005 , , 631 , 1120 lpez - morales , m. 2005 , , 660 , 732 lucy , l. b. 1967 , zeitschrift fr astrophysik , 65 , 89 maceroni , c. & montalbn j. 2004 , , 426 , 577 metcalfe , t. s. , mathieu , r. d. , latham , d. w. , & torres , g. 1996 , , 456 , 356 metropolis , n. , rosenbluth a. w. , rosenbluth m. n. , teller a. h. , teller e. 1967 , , 21 , 1087 morales , j. c. , et al .
2009 , , 691 , 1400 nutzman , p. & charbonneau , d. 2008 , , 120 , 317 popper , d. m. & etzel , p. b. 1981 , , 86 , 102 pra , a. & zwitter , t. 2005 , , 628 , 426 ribas , i. 2003 , , 398 , 239 ribas , i. , morales , j. c. , jordi , c. , baraffe , i. , chabrier , g. & gallardo , j. 2008 , , 79 , 562 rodon , m. 1986 , , 1986 , 165 , 135 soderblom , d. r. & mayor , m. 1993 , , 105 , 226 southworth , j. , maxted , p. f. l. , & smalley , b. 2004a , , 351 , 1277 southworth , j. , zucker , s. , maxted , p. f. l. , & smalley , b. 2004b , , 355 , 986 southworth , j. , smalley , b. , maxted , p. f. l. , claret , a. , & etzel , p. b. 2005 , , 363 , 529 tody , d. 1993 , in asp conf .
52 , astronomical data analysis software and systems ii , ed .
hanisch , r. j. v. brissenden , j. barnes ( san francisco : asp ) , 173 torres , g. & ribas , i. 2002 , 567 , 1140 vaccaro , t. r. , rudkin , m. , kawka , a. , vennes , s. , oswalt , t. d. , silver , i. , wood , m. , smith , j. a. 2007 , , 661 , 1112 van gent , h. 1926 , , 3 , 121 wilson , r. e. & devinney e. j. 1971 , , 166 , 605 zahn , j .- p .
1977 , , 57 , 383 zhao , j. , zhao , g. , & chen , y. 2009 , , 692 , 113 zucker , s. & mazeh , t. 1994 , , 420 , 806