the cosmic microwave background ( cmb ) is a tremendous tool for studying cosmology .
the anisotropies in the cmb provide a wealth of cosmological information ( for a recent review see hu & dodelson 2002 ) and upcoming experiments will provide high sensitivity ( approaching 1 @xmath6 ) and high angular resolution ( approaching 1 ) . such precise measurements will provide interesting constraints on secondary anisotropies in the cmb , imprinted by material along the line of sight at redshifts @xmath7 .
the largest secondary anisotropy is expected to be caused by galaxy clusters , through the sunyaev - zeldovich ( sz ) effect @xcite .
compton scattering of cmb photons by intra - cluster electrons leaves a spatial and spectral imprint in the cmb .
recent measurements of the sz effect have provided detailed maps of the electron distribution in galaxy clusters ( see birkinshaw 1999 and carlstrom 2002 for recent reviews ) but new experiments offer the promise of an order of magnitude increase in sensitivity . with @xmath6 sensitivity it should be possible to probe other cluster properties , such as bulk velocities and electron temperatures .    in this paper
we investigate the feasibility of measuring peculiar velocities and cluster temperatures at cm and mm wavelengths in the presence of contaminating sources , focusing on the effects of dusty high - redshift galaxies .
these issues have been partly addressed previously by @xcite and @xcite but both our understanding of dusty distant galaxies and cmb experimental capabilities have advanced tremendously in recent years .
we forecast errors by calculating the sz - parameter fisher matrix for multi - frequency maps that contain signals from the clusters as well as noise from cmb anisotropy , dusty galaxies and the measurement process .
our work extends previous studies of sz parameter reconstruction @xcite by considering point source and cmb anisotropy contamination and the influence of the atmosphere on the relative sensitivity of different frequency bands .
while @xcite did include the effect of cmb anisotropy contamination , explicitly subtracting it with an optimal filter , point source contamination was ignored . and
while @xcite included both cmb and point source contamination , forecasts were for one specific experiment : planck .
here we vary experimental parameters in order to guide experimental design .    in  2 we outline the sz effect and how it relates to cluster properties .
 3 summarizes the relevant known properties of dusty galaxies observed in the direction of galaxy clusters , while  4 outlines our methods for estimating the effects of contaminating sources on cluster measurements .
we discuss complications that will arise in analyses of real data in  5 and describe the experiments we consider in  6 .
our primary results are in  7 , where we forecast uncertainties for a four - channel reference experiment as a function of cluster gas temperature and core radius .
we explore the impact of varying experimental parameters such as angular resolution , observing time and number and placement of frequency channels .
we follow that with a discussion of possible radio point source contamination and close with a discussion of our results and implications for future instrumentation .
the sz effect is the change in energy of cmb photons from compton scattering with electrons , primarily in the intra - cluster medium . for recent reviews see @xcite and @xcite .
the high temperature ( several kev ) of the electrons relative to the cmb photons leads to a net increase of energy of the scattered photons , while a bulk motion of the cluster electrons along the line of sight leads to either a net redshift ( moving away from observer ) or blueshift ( moving toward observer ) of the scattered cmb photons . the first effect is the thermal sz effect and it has a distinctly non - thermal spectrum ; the scattering preserves photon number and mainly just shifts the energy of each photon , distorting the original blackbody spectrum . the effect of the bulk motion is the kinetic sz effect , and the emergent spectrum is that of a blackbody with a slightly different temperature .    for the thermal sz effect , the main physical parameters are the optical depth to thomson scattering , @xmath3 and the fractional energy gain per scattering @xmath8
. the combination @xmath9 sets the amplitude of the spectral distortion . for the kinetic sz effect ,
the redshift or blueshift is set by @xmath10 , so the relevant combination is @xmath11 .
typical electron temperatures are on the order of 0.01 @xmath12 , making relativistic effects modestly important @xcite .
the relativistic corrections are sensitive to the electron temperature through the relativistic correction to the thomson cross - section and through the relativistic form of the thermal maxwellian velocity distribution .
these corrections can be on the order of @xmath13 or higher at many observing frequencies .
thus the temperature difference due to the kinetic sz effect is given by @xmath14 while the temperature difference due to the thermal sz effect as a function of frequency @xmath15 ( expressed in dimensionless units @xmath16 ) is given by @xmath17 \quad , \ ] ] where the relativistic corrections have been represented by @xmath18 .
typical values for massive clusters of @xmath19 , @xmath20 and @xmath21 lead to typical amplitudes ( @xmath22 ) for the thermal effect , kinetic effect and relativistic corrections of @xmath23 , @xmath24 and @xmath24 respectively .    in this work ,
we take our fiducial galaxy cluster to have @xmath25 , @xmath26 kev and @xmath27 .
we vary the thermal sz properties using the scaling relations of @xcite .
the central @xmath0 parameter was taken to scale as @xmath28 , which leads to @xmath29 .
self - similar evolution would result in @xmath30 ; we use the slightly steeper relation to be consistent with evidence for excess entropy .
we ignore evolution of intrinsic cluster properties with redshift , such as the expectation that clusters ( of a fixed mass ) become more compact ( and therefore hotter ) at higher redshift .
it is expected that the effect of this evolution is to enhance the thermal sz effect , but the details of this redshift evolution are sensitive to poorly understood gas processes @xcite .
star - forming galaxies can be very luminous in the submillimeter wavelength range , as recent measurements by the sub - millimeter common user bolometer array ( scuba ) and other experiments have discovered ( see blain 2002 for a recent review ) . in particular , galaxy clusters
are often found to coincide with relatively bright submillimeter emitters , probably because of gravitational lensing effects @xcite .
the emission is mainly from dust at a temperature of several tens of kelvin , putting the peak of the radiation spectrum at submillimeter wavelengths , but there can still be significant emission at millimeter wavelengths .
many galaxy clusters have been observed to coincide with submillimeter sources with fluxes of several mjy at @xmath15=350 ghz . assuming a typical spectrum for these sources would lead to expectations of fluxes at 150 ghz on the order of 1 mjy . in a beam of 1x1
this would correspond to nearly 30 @xmath6 of contamination if not correctly taken into account .
this makes dusty galaxies a non - negligible contaminant for studies of galaxy clusters at mm wavelengths @xcite .    throughout
we often characterize a given intensity as the equivalent departure from the mean cmb temperature , @xmath31 .
the conversion factor , for @xmath32 , is ^-1 & = & ( ) ^2 + & = & ( )      the experiments we consider are not sensitive to the absolute flux , but to spatial variations in the flux .
these variations contribute to the variance of temperature fluctuations on the sky , observed with a gaussian beam profile with full width at half maximum of @xmath33 , ( t)^2 = _
l 2l+1 4 c_l where c_l = [ b_t]^-2 s^3 ds .
this variance decreases with beam size as @xmath34 . for the convenience of using a quantity that is a property of the sky only , and not the angular resolution of the telescope , we define @xmath35 where @xmath36 is the solid angle of the beam .
the shot noise at 350 ghz can be inferred from scuba observations .
@xcite used these observations to fit a double power - law form for @xmath37 : dn(>s)/ds = ^-1 with @xmath38 mjy , @xmath39 , @xmath40 and @xmath41 the result is shot noise with an rms of @xmath42k - arcmin where the subscript `` 350 '' is used to denote the frequency .
if the spectral behavior of these sources is known , then it would be straightforward to measure the flux at a higher frequency , where the emission is stronger and the sz effects are smaller , and subtract the appropriate levels from the measurements at other frequencies .
however , the spectral behavior is not perfectly known , so there will be uncertainty due to imperfect subtraction , as well as uncertainty due to noise in the measurement at higher frequencies .    to investigate the homogeneity of the spectral behavior , we used submillimeter observations of local galaxies @xcite .
local galaxies were modeled as having a spectrum of the form @xmath43 , where @xmath44 is the blackbody function and @xmath45 here refers to the emissivity spectral index ( not the bulk velocity of the cluster ! ) .
dunne provide best - fit values of @xmath45 and dust temperatures .
we used these best - fit spectra for the local galaxies , placed them at a range of redshifts and fit them over the range 150 to 350 ghz assuming a power - law @xmath46 . as a function of redshift
the mean spectral index was well - fit by @xmath47 , with a fairly constant scatter at each redshift of @xmath48 .
the moderate flattening with increasing redshift arises because the radiation that is observed at 350 ghz originates at a higher frequency , closer to the peak of the dust emission .    without source redshifts
, it is not clear what to take as the mean spectral index . assuming a uniform distribution in redshift between 0 and 5 , the mean spectral index is 2.6 and the @xmath49
scatter around this mean is 0.4 .
the sources are unlikely to be at very low redshift , simply due to volume considerations , but the redshift distribution of the observed submillimeter - emitting galaxies is very uncertain .
we set our mean value for @xmath50 as @xmath51 and set @xmath52 , unless specified otherwise .
we set @xmath53 since the redshift distribution will be more peaked than the uniform one that leads to 0.4 .
the quantity @xmath54 is the spectral index rms for _ individual _ galaxies . for @xmath55 galaxies of similar brightness the composite intensity will have a scatter in spectral index that is smaller by @xmath56 .
we define an effective number density of galaxies by weighting them according to their contribution to shot - noise variance . thus , n_eff & = & s^3 d / s^3 d + & = & 0.64/(sq .
arcmin ) where the last equality follows from the @xcite @xmath37 . for simplicity we set @xmath57(sq .
finally we define @xmath58 and take our fiducial value to be @xmath59-arcmin .
the @xcite @xmath37 is for an unlensed population of sources .
lensing changes both the brightnesses of galaxies and their number densities .
indeed , the magnification provided by galaxy clusters has been exploited to study the population to fainter limiting magnitudes @xcite .
the magnitude , and even the sign , of the effect on shot noise depends on @xmath37 .
magnification increases @xmath60 which increases shot noise , but also decreases @xmath37 ( at fixed unlensed @xmath60 ) which decreases shot noise . for a power - law , @xmath61 ,
the net effect is an increase in shot noise if @xmath62 . for the relatively bright population at 350 ghz of interest for this work ,
the shape of @xmath63 is such that lensing leads to an enhancement of shot noise .
@xcite pointed out that lensing would significantly enhance the level of contamination by ir luminous galaxies of sz observations .
the key quantity for lensing is the einstein radius , @xmath64 .
a point source with source plane location coincident with the cluster center ( @xmath65 ) will form a ring of infinite magnification appearing in the image plane at distance @xmath64 . in the image plane ,
magnification rises from the center towards @xmath64 and then drops beyond @xmath64 .
thus , as @xcite showed , the enhancement of confusion noise for an observation towards the cluster center increases with increasing beam size until peaking at @xmath66 as more weight is placed on galaxies near @xmath64 , and then drops with further increase in @xmath33 as the effect is diluted by adding in more galaxies from beyond the einstein radius .
for a rough model of the effect of lensing we introduce the confusion noise enhancement factor , @xmath67 which is the ratio of lensed confusion noise to unlensed confusion noise at distance @xmath68 from the cluster center and for observations with angular resolution @xmath33 . for the dependence of @xmath69 on @xmath33 we use the results of @xcite for his galaxy evolution model i1 for a cluster with velocity dispersion of @xmath70 km / s and @xmath71 when it is placed at @xmath72 . for this case
@xmath73 1.1 , 1.5 , 2.5 , 2.1 , 1.8 and 1.5 for @xmath74 10 , 20 , 40 , 60 , 100 and 200 arcsec respectively .
we ignore the redshift - dependence of @xmath69 since for @xmath75 @xcite shows that @xmath69 is highly independent of the cluster redshift .
we also ignore the @xmath76 dependence of @xmath69 since @xcite shows this dependence to be quite slow , covering a range of 30% as @xmath76 varies from 800 km / s to 2000 km / s .
for @xmath77 km / sec one expects a gas temperature of about @xmath78 kev @xcite .
for dependence on @xmath68 we simply set @xmath79 for @xmath80 and @xmath81 otherwise .
for all cases that we consider here the beam size will be larger than a typical einstein radius for a cluster .
we model the data at sky location @xmath82 and frequency @xmath83 as having contributions from both signal and noise : ( _ a,_i ) = s(_a,_i ) + n(_a,_i ) .
we consider each component in detail below .
the signal is due to three components , s(_a , x_i ) = _
m=1 ^ 3 p_m h(x_i)f_m(_a ) , [ eqn : signals ] where p_1 & & y = dl _ t n_e = | _ + p_2 & = & dl _
t n_e v / c = | _  and + p_3 & = & dl_t n_e ^2 = |_^2 = y | _ y. here @xmath84 is the electron number density and @xmath85 is the dimensionless gas temperature .
the subscripts on the average quantities indicate the weighting along the line of sight . for @xmath3
the weighting is with the electron - number density and for @xmath0 the weighting is with the pressure , @xmath86 .
the optical depth , @xmath3 is given by @xmath87 .
the integrals defining @xmath88 and @xmath89 are along the line of sight through the cluster center .
we take the angular profile of the cluster to be h ( ) = ( 1+(||/_c)^2)^-1/2 the beam - convolved cluster profile is denoted by @xmath90 , and we assume gaussian beams in all that follows .
note that the assumption of the same cluster profile for all three signal components in equation [ eqn : signals ] is only valid for an isothermal intracluster medium .
the frequency dependences are given by @xcite f_1 ( ) & & x -4 + f_2 ( ) & & 1 + f_3 ( ) & & -10 + x - x^2 + x^3 + ( -+ x)s^2 where x & & x coth(x/2 ) + s & & x / sinh(x/2 ) .
the first and third components are due to the scattering of cmb photons off of the hot electrons in the cluster and represent the thermal sz effect , while the second component is the kinetic sz effect .
we have neglected contributions to the sz effects that are @xmath3 times higher order products of @xmath91 ( @xmath92 ) and @xmath68 ( @xmath93 ) . for clusters with temperatures of 10 kev or lower the higher order corrections are less than five percent @xcite at all frequencies except near the null of the thermal sz spectrum ( where the thermal sz signal is sub - dominant ) .
because they are the amplitudes of the different frequency shapes , @xmath94 , @xmath95 and @xmath89 are the theoretical parameters most directly related to the data . from them one can estimate the physically interesting quantities : _ y & = & p_3/p_1 =    & = & p_1 ^ 2/p_3 = = + & = & p_2/ = | thus the temperature determinations from sz measurements are measurements of the pressure - weighted average gas temperature .
velocity measurements are of the mass - weighted average velocity ( assuming electron density traces mass ) times a correction factor dependent on differences in different temperature averages .
the optical depth determinations are of the optical depth times another ratio of different average temperatures . throughout ,
we assume isothermality in order to avoid these complications ; the distinctions must be kept in mind when analyzing real data , particularly for comparison of x - ray and sz gas temperatures .    in a recent preprint @xcite it was noted that peculiar velocity measurements of non - isothermal clusters could be biased .
the source of this bias can be clearly seen in the above expression for @xmath96 , where the true @xmath97 is multiplied by the ratio of two different weightings of the electron temperature . in general the optical depth weighted temperature
will not be the same as the pressure weighted temperature , resulting in a biased velocity estimate .
this is a direct consequence of the biased estimate of the optical depth .
the noise has a contribution from measurement error , @xmath98 , from primary cmb fluctuations , @xmath99 and from emission from galaxies in the direction of the cluster @xmath100 .
we assume the composite dusty galaxy spectrum has a power - law form in intensity , @xmath101 where @xmath50 is spatially varying with @xmath102 and @xmath103 as discussed earlier . linearizing the dependence of @xmath104 on @xmath105 we can write @xmath106 as n_gal(_a,_i ) = ( _ a,_i ) e(_i,_b ) where @xmath107 is the shot noise enhancement factor discussed earlier .
the frequency dependences are given by f_4 ( ) & = & ^|-2 x^-2 e^-x(e^x-1)^2/f_4(350  ghz ) + f_5 ( ) & = & f_4 ( )    the cmb noise contribution is independent of frequency , and the instrument noise contributions as a function of frequency depend on the experiment being considered .
the total noise in the map is then @xmath108 .
the noise in the map therefore has the two - point function , considering points in the map @xmath82 at frequency @xmath83 and @xmath109 at frequency @xmath110 : @xmath111 ^ 2 f_5(\nu_a)f_5(\nu_b ) \\
\nonumber   & & + \sigma^2(\nu_a)/\omega_{\rm pix}\delta_{ab}\delta_{i j}\end{aligned}\ ] ]    in the above , @xmath112 is a shot - noise covariance matrix ( with off - diagonal correlations only due to beam - smoothing ) normalized to unit variance , and the last term gives the contribution of instrument noise , which will be considered in detail in  6 .
@xmath113 is the covariance matrix of the cmb fluctuations , given by c^cmb_i j = _
l c_l p_l(_i j ) where @xmath114 is the angular separation between pixels @xmath115 and @xmath116 .
this analysis implicitly assumes that we know the statistical properties of all sources of noise , including the dusty galaxies .
this is an excellent approximation for the cmb fluctuations , but the dusty galaxy noise power spectrum will be affected by both clustering and gravitational lensing .
on arcminute scales the clustering effects should be relatively small , but the effects of gravitational lensing will modify the noise properties .
we forecast how well the three parameters of our model can be measured by calculating a fisher matrix . combining the spatial and spectral indicators @xmath115 and @xmath117 into a combined index @xmath4 ( and @xmath116 and @xmath118 into @xmath15 ) we can write the fisher matrix as f_pp = ( c^noise)^-1 _ to this fisher matrix we usually add a prior fisher matrix which carries the information we have that is not contained in our measurements in the 30 to 350 ghz range . the error covariance matrix for our three parameters is then given by c_pp = _ pp. note that since our data depend linearly on our parameters @xmath119 , this fisher matrix calculation of the expected error covariance matrix is exact , given our model of the data .
we include prior information due to x - ray measurements of the gas temperature with error @xmath120 so that f_prior(1,1 ) & = & ( _ ) ^-2 + f_prior(3,3 ) & = & ( _ y)^-2  and + f_prior(1,3 ) = f_prior(3,1 ) & = & - .
including the prior information in this manner is approximate since it results from taylor expanding the dependence of @xmath68 on @xmath94 and @xmath89 .    to calculate errors on a new set of parameters , @xmath121 , that are functions of our original parameters @xmath119 , we taylor expand the dependence of @xmath121 on @xmath119 to first order about the fiducial value .
we then calculate the fisher matrix in these new coordinates by : f_ii = _
pp r_ipf_ppr^t_pi where the transformation matrix @xmath122 .
for example , such a variable transformation is necessary to get errors on @xmath123 . dropping the higher - order terms in the taylor expansion makes this procedure approximate as well .
calculating these fisher matrices can be computationally expensive as matrices of size @xmath124 have to be inverted and multiplied where @xmath125 is the number of map pixels and @xmath126 is the number of frequencies .
we therefore wish to pixelize no more finely than necessary , and keep the physical area of the map as small as possible .
if the map is too small compared to the beam - convolved cluster profile we will lose the large - scale information necessary for subtracting off the cmb contamination .
if the pixel size is too large we will lose small - scale information .
we have found the following prescriptions to ensure we are neither losing information , nor using much larger @xmath125 than necessary : _ m & = & 4(2_c + _
b ) + _ pix & = & ( _ c+_b)/2.5 where the map and pixels are squares of length @xmath127 and @xmath128 respectively .
our model is fairly detailed and sufficient in many ways , but there are several respects in which real clusters could present some challenges , such as temperature structure in the intra - cluster medium , incomplete knowledge of the lensing properties of galaxy clusters , and internal bulk flows in the cluster
.    temperature gradients in the intra - cluster medium would lead to a mis - match between the signal maps of equation [ eqn : signals ] . observed clusters show significant departures from isothermality @xcite at relatively large radii and at small radii , but the largest scales are already obscured by the primary cmb anisotropies and the small scales contribute relatively little to the sz effect .
it is therefore not expected to be a significant effect but will complicate data analysis .
similarly , we have assumed that we know the cluster template . in the cases where we assume complementary x - ray spectroscopy
it is quite reasonable to assume that the x - ray image gives an estimate of the electron spatial distribution . without x - ray information the best estimate of the spatial template will come from the thermal sz map itself , where a large mismatch between the assumed spatial template and the observed emission should be evident . a parametrized model for the cluster could be introduced into the fit without significantly affecting the constraints on peculiar velocities .
we have assumed that the magnification due to lensing is a simple step function in radius .
the main effect of lensing is to amplify the poisson noise , and the amount of amplification depends on the details of the mass profile . incomplete knowledge of the mass distribution will therefore lead to less efficient component separation .
deep optical images will be useful for the purpose of studying the strong lensing properties of galaxy clusters .    on a related note
, we have assumed that we know the statistical properties of the galaxy contamination .
lensing modifies the noise properties of the background galaxies , and there is uncertainty in the statistics of dusty galaxies .
multi - frequency observations of many fields , both with and without clusters , will lead to a good understanding of the statistics of dusty galaxies .
we have also assumed a single peculiar velocity for the cluster whereas real clusters show evidence for internal flows , sometimes as large as 3 - 4000 km s@xmath129 @xcite .
some striking visualizations of the kinetic sz fluctuations induced by these internal flows are given in @xcite .
it has been shown @xcite that the average peculiar velocity provides an unbiased estimate of the true bulk velocity but with an added dispersion of roughly 100 kms@xmath129 . for our purposes
this can be considered as an extra source of noise which is small compared to the uncertainties due to astrophysical confusion .
we have neglected calibration uncertainty .
since the total signal , anywhere other than near the thermal sz null , is about ten times as large as the kinetic sz signal , accurate calibration of one band relative to another is important .
calibration errors will have to be controlled to better than about 10% for peculiar velocity measurements to be better than about one @xmath130 .
the same argument applies to gas temperatures for clusters with temperatures near 6 kev , since the relativistic correction is also about 10% of the total signal .
but gas temperatures , since they do not suffer from cmb anisotropy confusion , can be measured to much better than 1 @xmath130 .
achieving @xmath131 kev , as can be done if we neglect calibration uncertainty , may require calibration uncertainties as small as 2% .
measurements closer to the thermal sz null will reduce these sensitivities to calibration error .
as our reference experiment we have modeled the sensitivity achievable for a generic ground - based bolometric experiment with three frequency bands selected to lie in the atmospheric windows available at a good millimeter site . we have also assumed that a 30 ghz measurement is available , with sensitivity similar to that expected from the sza .
we assume that the sensitivity in each channel of the bolometric experiment is limited only by fluctuations in the photon background itself and that other potential sources of noise , such as phonon fluctuations in the bolometers or electronic readout noise , are negligible ( this can be achieved with current detector technology , see lange 2002 ) .
we assume that @xmath132 noise introduced by fluctuations in atmospheric water vapor emission can be adequately subtracted from all of the data .
this can be accomplished by some kind of spatial chopping or removal of common - mode signals from array data , or alternatively the spectral properties of the atmospheric noise can be used .
the latter technique has been demonstrated by the sunyaev - zeldovich infrared experiment ( suzie ) which operates at 150 , 220 and 350 ghz and which has been used to set limits to the peculiar velocities of 11 galaxy clusters @xcite .
for the purposes of this paper we do not include in our calculations any of the overhead ( increase in real observation time ) that this process will produce , since this will depend on the individual experiments . for example , since the atmosphere is in the near field of most large telescopes , many detectors view the same column of atmosphere , yet observe different parts of the cluster , or none of the cluster at all . by using detectors that are a large distance from the cluster to subtract common - mode fluctuations , it may be possible to subtract atmospheric noise with little or no degradation in sensitivity .
this is a commonly - made assumption that we also include here .    with the above assumptions the performance at each frequency
is characterized by the noise equivalent power ( nep ; quoted in whz@xmath133 ) , given by : @xmath134 where @xmath135 is the total background loading on the detector in watts , @xmath136 is the detection bandwidth in ghz and @xmath137 is the number of waveguide modes that are detected ( @xmath138 for a diffraction - limited system ) .
the first term in equation  [ model : e1 ] is caused by shot noise due to poisson statistics of the incoming photons while the second term accounts for the effect of photon ( boson ) correlation . for more detail , see @xcite .
sources of background loading include the atmosphere , warm emission from the telescope and surroundings , and of course the cmb itself .
the total background loading is then : @xmath139 where the sum is over all sources of power on the detector . for each source of loading
, the power at the detector can be calculated from : @xmath140 where @xmath141 is the equivalent rayleigh - jeans ( rj ) temperature of the @xmath115th background source , @xmath142 is the fraction of photons from the background source that reach the detector , @xmath143 is the effective telescope area and @xmath144 is the solid angle response of each detector on the sky .
the integral is over the spectral band to which the detector responds .
we assume for simplicity , a square band of width @xmath136 .
we now consider each source of background power in turn .    for the telescope
the equivalent rj temperature @xmath145 is : @xmath146 where @xmath147 is the physical temperature of the telescope , @xmath148 is a frequency - dependent emissivity and @xmath149 .    for the atmosphere : @xmath150\right\ }    \label{model :
e5}\ ] ] where @xmath151 is the physical temperature of the telescope , @xmath152 is the frequency - dependent zenith optical depth and @xmath153 is the zenith angle of the observation .
for the cmb : @xmath154 where @xmath155k and @xmath156 .
the optical efficiency is assumed to be : @xmath157   & \mbox{for $ i=3 $ ( cmb ) } \end{array } \right .
\label{model : e7}\ ] ] which are values typical of those measured for bolometric systems @xcite .
once the nep is determined , the sensitivity to cmb fluctuations , including the sz effect , is obtained by calculating the noise equivalent temperature ( net ) as follows : @xmath158 where @xmath159 is given above .
we further assume that the focal plane is fully sampled with detectors spaced at half the beam size .
each detector will have an angular response that is diffraction - limited at each frequency , in which case the solid angle response of each detector is just @xmath160 . with a suitable scanning strategy , an oversampled map can always be produced with a pixel solid angle @xmath161 .
the rms sensitivity per map pixel is then : @xmath162 where @xmath163 is the integration time per detector and the net is per focal plane element .
the factor @xmath164 accounts for the increase in sensitivity that a fully sampled focal plane can achieve over a focal plane with feed horns that are exactly diffraction - limited .
we assume @xmath165 which is the ideal limit .
the real value of @xmath164 depends upon factors such as detector noise and instrument photon background @xcite .
we include a map pixel size here to allow us to compare experiments with different angular resolutions .
if the array is large enough the scanning can be arranged so that there are always focal plane elements viewing the cluster .
no integration time is lost in this case .
other chopping schemes could degrade the sensitivity by a factor that depends on the details of the particular scheme being used .
thus our final numbers may be taken as best case values .    in order to use equations  [ model : e1 ] through [ model : e9 ] , we now define the details of the experiments . for our reference experiment
we assume a ground - based off - axis 8-m telescope located at a site with an altitude of 17,000@xmath166 and precipitable water vapor of 0.5 mm .
the telescope and instrument are assumed to have an equivalent rayleigh - jeans temperature of @xmath167k that is frequency independent .
figure  [ model : f1 ] shows the dependence of @xmath168 on frequency , for a map pixel size of 1 sq arcmin .
two models are shown , one with a bandwidth of 20% at each frequency and one with a bandwidth of 10% .
there are broad atmospheric windows visible at 135 - 165ghz and at 215 - 290ghz .
the window at 345ghz is sufficiently narrow that only a 10% wide frequency band can be utilized but it is a very useful atmosphere and/or point source monitor .
based on this figure we have placed three of the four channels of the reference experiment at 150 , 220 and 280ghz .
the sensitivities are shown in table  [ tab : experiments ] .
.assumed experimental specifications and sensitivities per square arcminute pixel . [
cols="^,^,^,^ " , ]      we begin with our fiducial cluster , studying the error forecasts as a function of confusion noise and spectral index rms in fig .  [
fig : dusty ] . for the lower @xmath169 values
one can infer from the flattening of the error contours with increasing shot noise amplitude that the data themselves are capable of fitting for the dusty galaxy component with residuals at about the 100 to 200 @xmath4 k level .
if instrument noise were to increase , this transition from dependence on @xmath170 to independence would occur at higher values .
the trend of error level with @xmath169 is as expected : flat at low @xmath171 and steeper at higher values .
peculiar velocities are more affected , particularly at non - zero @xmath169 , than the other parameters .
we now fix the point source contamination parameters and study our results as a function of galaxy cluster core radius , @xmath172 , and gas temperature , @xmath173 .
note that as we vary @xmath172 the map size we assume increases in order to assure complete coverage of the cluster and good cmb subtraction . for
large - format bolometer arrays , expected to have more than 1000 elements , the larger map size does not require longer integration times because the field of view will be much larger than the galaxy cluster .
therefore , despite the added coverage , one hour of observation is still sufficient to obtain the reference experiment sensitivities in table 1 . on the left of fig .
[ fig : reference ] we set the point source contamination parameters to zero and on the right to their fiducial values @xmath174k - arcmin and @xmath175k - arcmin .    because of the linear response of the signal to @xmath0 , @xmath176 is independent of @xmath0 and therefore @xmath173 as seen in the left side of fig .
[ fig : reference ] .
in contrast , at fixed @xmath0 and @xmath177 the signal is proportional to @xmath178 so one expects @xmath179 .
confusion with @xmath0 , more important at high @xmath173 , leads to the error decreasing more slowly than @xmath180 .
the @xmath181 has a similar dependence on @xmath173 since the negligible error on @xmath0 means @xmath182 and our scaling assumes @xmath29 .    for fixed three - dimensional distribution of gas pressure
, central @xmath0 values do not vary as a function of redshift .
core radii do as @xmath183 where @xmath184 is the angular - diameter distance to redshift @xmath185 . thus to study dependence on @xmath185 we vary @xmath172 with fixed @xmath0 .
note that @xmath186 0.5 , 0.7 , 1.3 , 1.4 for @xmath187 , 0.3 , 1 , 2 respectively , assuming @xmath188 and @xmath189 .
note that the number of pixels with sz signal in any given range scales as @xmath190 , so the errors on @xmath3 , @xmath0 and @xmath173 decrease as @xmath191 .
the peculiar velocities on the other hand suffer from contamination from the cmb .
as @xmath172 increases it becomes more difficult to distinguish the signal proportional to @xmath177 from the very red cmb power spectrum , hence @xmath192 goes up . at large @xmath172 , where the error in @xmath97 is dominated by the cmb - induced error in @xmath177 ,
the error in @xmath3 is negligible . in this case , @xmath193 so @xmath194 , which is the scaling we see . at smaller @xmath172 uncertainty in @xmath3
is no longer negligible and the scaling with @xmath173 is more complicated .
including the effect of the fiducial point source contamination ( right side of fig .
[ fig : reference ] ) does not alter the error forecasts dramatically . for @xmath1 the largest impact is at the smaller @xmath172  which is unfortunate since this is where @xmath1 is measured best . at small @xmath172
the cmb noise is smaller and ( less importantly ) the lensing enhancement of confusion noise is largest . at @xmath195
the velocity errors are almost doubled .
figure  [ fig : reference_time ] shows how the error forecasts depend on observing time . on the left side
one can see @xmath196 indicating that for the reference experiment the dominant source of uncertainty on @xmath0 is instrument noise , not dusty galaxies .
similar scalings are seen for @xmath181 and @xmath197 .
the @xmath198 in contrast , quickly saturate with increased observing time producing very little improvement since the dominant source of error is the cmb contamination .
the saturation occurs even sooner for larger @xmath172 .    as can be seen in fig .
[ fig : reference_fwhm ] parameters other than @xmath1 are highly independent of the angular resolution @xmath33 .
this is because none of the signal is lost with increasing @xmath33 : the signal is simply spread out over greater area .
errors would go up with increasing @xmath33 for single - pixel observations as opposed to the map - making observations we consider .
the peculiar velocity dependence on @xmath33 is steep for @xmath199 , since spreading out the signal degrades the ability to subtract off the cmb contamination .
note that to explore the effect of angular resolution independent from sensitivity , we have not scaled the sensitivity with beam size as one would expect from eq .
[ model : e9 ] , but rather fix the weight per solid angle .
the apparent decrease in @xmath198 with increasing @xmath33 at fixed @xmath200 is an artifact of our sparse sampling of the @xmath172 , @xmath33 space .
we now consider the importance of each of the channels of our reference experiment .
on the left side of fig .
[ fig : reference_minus ] we show error forecasts for the reference experiment with no 30 ghz channel , ` reference-30 ' . without this channel it is harder to distinguish a change in @xmath0 from a change in @xmath173 , thus the errors in both of these increase considerably .
their ratio , @xmath201 , suffers similarly .
the impact on peculiar velocities is less dramatic , especially at high @xmath172 where the cmb contamination remains dominant .
changes are significant at smaller @xmath172 where , e.g. , @xmath198 goes from 250 to 400 km / sec at @xmath202 kev .
the strong degeneracy between @xmath0 and @xmath173 for reference-30 can be seen on the right side of fig .
[ fig : reference_minus ] .
the correlation coefficient , @xmath203 increases from 0.93 for reference to 0.998 for reference-30 .
this near - maximal correlation coefficient means that the error in @xmath173 ( @xmath0 ) is increased by a factor of 14 compared to the error in @xmath173 ( @xmath0 ) if @xmath0 ( @xmath173 ) were held fixed . including a reasonable prior on the velocity of 300 km / sec decreases the @xmath173 and @xmath0 errors by about 20% .
the importance of a low frequency measurement for determining @xmath173 has been emphasized by @xcite and @xcite .    after the 30 ghz channel ,
the most critical channels in decreasing order are 150 ghz , 220 and 280 .
removal of the 150 and 220 ghz channels can increase @xmath197 by up to a factor of 5 , although only about a factor of 2 at our small fiducial core radius of @xmath195 assumed in the right side of fig .
[ fig : reference_minus ] . these increases , at the assumed level of instrument noise , are not sufficient to degrade the peculiar velocity errors .
the dashed line on the right side of fig .  [ fig : reference_minus ] is @xmath28 , through our fiducial values of @xmath0 and @xmath173 .
we see that the degeneracy direction is such that assuming the scaling relation and a normalization would lead to much more precise determinations of both @xmath0 and @xmath173 . on the other hand
, the orientation of the degeneracy makes testing the scaling relation more difficult than if it lay along the line .
= .67    we saw that for the reference experiment the loss of the 30 ghz channel increased @xmath197 , as well as @xmath198 for small @xmath172 .
the lack of a low frequency channel can therefore be compensated by an x - ray determination of the temperature .
thus we show forecasts for the act survey and suzie - iii observations both complemented by x - ray temperature determinations of @xmath131 kev .
the results for the two different experiments are very similar .
we assume one hour of integration with suzie - iii for every 48 sq .
arcmin covered .
recall that in order to assure good subtraction of the cmb , we assume a map of linear extent @xmath204 .
we have made no attempt to optimize this observing strategy .
the right - most figure in figure  [ fig : planned ] shows that adding a 30ghz measurement to the millimeter data can break this degeneracy , as we expect from our previous forecasts for the reference experiment in fig .
[ fig : reference ] . in 1 hour of observation with suzie iii ( per 48 sq .
arcmin of map ) and a 12-hour sza observation , the peculiar velocity of a cluster can be measured to 250 km s@xmath129 _ even without an x - ray measurement of temperature_. a low frequency channel will thus be crucial for obtaining peculiar velocity and gas temperatures at high redshift where @xmath205-ray brightness is very low .
we also see in figure  [ fig : planned ] a temperature error for our fiducial cluster of @xmath206 kev ( and even smaller for larger @xmath172 ) . at this accuracy level ,
one expects to begin to see discrepancies between x - ray temperature measurements which are more sensitive to the hotter , inner regions and sz temperature measurements which are more sensitive to the cooler , outer regions @xcite .
deeper integrations can slowly improve these measurements as seen from the study of dependence of forecasts for the reference experiment on observation time .
we have assumed that ( i ) the frequency spectrum of cluster radio sources will be sufficiently steep at high frequencies that measurements at 150 ghz and higher will not be contaminated and ( ii ) that measurements at 30 ghz will be sufficiently high resolution that point sources can be subtracted at that frequency ( this is possible with interferometric experiments such as the sza ) .
consequently we have concentrated only on the contamination introduced by dusty galaxies known to be strong submillimeter - wavelength emitters and have excluded radio sources from our analysis .
we will now justify this omission .
clusters often contain radio point sources at mjy levels , as seen from 30 ghz sz measurements @xcite .
it is important to note that radio point sources and submillimeter point sources are _ not _ the same objects .
the excess of radio sources toward galaxy clusters is most likely due to emission from cluster galaxies themselves @xcite , whereas the dusty galaxies are typically background sources that are not associated with cluster members .
very little is known about the emission from radio point sources at millimeter wavelengths , but their spectra are expected to steepen due to the energy losses of the most energetic electrons to synchotron emission .
this is confirmed to some extent by spectral indices measured between 10 - 90 ghz for bright sources @xcite , although there is substantial scatter in the measured index .
@xcite find a mean index of @xmath2070.5 between 10 and 90 ghz , while @xcite find a mean index at 40 ghz of @xmath208 .
measurements by @xcite of the spectral indices of the 208 point sources detected by wmap at 40 ghz @xcite , suggest a mean of @xmath209 at 20ghz with large scatter , but also show that the spectra of radio sources are often not well - approximated by a simple power law .
sources of only a few mjy at 30 ghz could easily still have mjy fluxes at 150 ghz if the spectrum does not steepen significantly above 90 ghz .
we have used the point source counts from wmap , dasi , vsa and cbi summarized in @xcite to determine a best fit model of : @xmath210 where @xmath211 per sq .
deg . and @xmath212 mjy .
our expression gives number counts for mjy sources that are slightly lower than @xcite , but that are consistent with those determined by @xcite based on sz measurements of clusters at 30ghz .
consequently we expect that the number of sources above 0.1 mjy at 30 ghz is less than one per 4 sq arcmin ( roughly the size of a moderate cluster ) .
repeating our fisher analysis with radio sources at this level , we find that for all the experimental situations we have considered , the peculiar velocity uncertainty for our canonical source is unaffected , even if a spectral index of @xmath213 is assumed , and the number counts are 10 times higher than we have assumed .
the spectra of mjy radio sources will soon be much better understood when the sza and the array for microwave background anisotropy ( amiba ) start operation at 90 ghz .
however , a complete understanding of faint radio sources at mm wavelengths will likely require the sensitivity and resolution of the atacama large millimeter array ( alma ) .
at roughly the level of 1 kev the distinction between x - ray emission - weighted temperatures and electron - weighted temperatures becomes important @xcite , so a comparison of the two `` gas temperatures '' could be valuable for studies of the physics of the intra - cluster medium .    furthermore , the pressure - weighted temperature should be an excellent indicator of cluster mass .
this could allow independent calibration of the mass scale for galaxy cluster surveys that use x - ray , optical or sz selection criteria .
the measurement of peculiar velocities using the kinetic sz effect has been a goal of sz measurements for many years .
peculiar velocities at large distances will allow measurements of cosmological parameters and reconstruction of the large scale gravitational potential .
the two - point function of radial peculiar velocities is sensitive mainly to @xmath214 and largely insensitive to the dark energy equation - of - state parameter , @xmath215 or other cosmological parameters . a linear theory calculation of
how well the two - point function can constrain @xmath214 gives @xcite [ eqn : peelknox ] _ m = 0.04 . the suzie iii and act collaborations plan to measure optical redshifts for 400 clusters with masses greater than @xmath216 . since typical errors for peculiar velocities on these
may be 300 km / sec and the velocity rms is about 400 km / sec , these data may be able to achieve @xmath217 .
however , equation  [ eqn : peelknox ] assumes that the variance in our peculiar velocity measurement errors is perfectly well known . in reality
, the statistical properties of the errors will be difficult to know well because of their dependence on the subtraction of the contamination of dusty galaxies .
if the assumed error variance is different from the actual error variance by @xmath218 then there will be a systematic error in @xmath214 of _ m = f()^-1 0.02 ( f/0.1 ) independent of the number of clusters measured . to control this systematic error
, we will need to learn more about the spectral properties and luminosity functions of sub - mm galaxies . to avoid this source of systematic error
one could use a linear statistic , the average difference between the radial velocities of two clusters separated by distance @xmath219 @xcite .
how this statistic depends on cosmological parameters is under investigation .
another application of peculiar velocities is gravitational potential reconstruction @xcite .
as @xcite point out , due to the low number density of clusters this is only possible on very large scales .
comparison with reconstruction from galaxy number counts can determine galaxy biasing properties .
a good understanding of the large scale potential would be valuable for investigating environmental effects on galaxy formation and would be a starting point for constrained realizations of simulations of the large  scale structure of our universe .
contamination by emission from dusty galaxies will make accurate measurements of peculiar velocities difficult . with four observing frequencies , few @xmath6 sensitivity , and arcminute resolution it will be possible to measure peculiar velocities at the level of roughly 200 km / s for massive galaxy clusters .    for comparison , previous estimates of uncertainties in peculiar velocities by holder ( 2003 ) assumed @xmath220 sensitivity and perfect point source removal , which on figure [ fig : reference ] roughly correspond to @xmath221 and an observing time of several tens of hours .
in that case it was found that uncertainties well below 100 km / s were possible , as could be extrapolated from figure [ fig : reference ] .
such exposure times may be prohibitively long and removing point source contributions to this level will likely require alma .
alma will provide both a better understanding of the physical nature of dusty galaxies in general as well as measurements of the contaminating galaxy fluxes at exactly the frequencies of interest .
the reference experiment , as well as combinations of planned observations , can provide measurements of peculiar velocities at the level of 150 - 200 km / s , even in the presence of contamination by dusty galaxies .
this is sufficient accuracy to constrain @xmath214 to better than 10% . at the same time galaxy cluster temperature measurements
will be possible at the kev level .
this will allow interesting comparisons with x - ray spectroscopic measures of the x - ray emission - weighted temperature .
peculiar velocities measured using the kinetic sz effect , a long - standing challenge to cmb experimentalists , may soon be a reality .
the final limiting obstacle appears to be contamination from dusty galaxies , but sensitive multi - frequency measurements should clear this hurdle .
+          , c.  l. , hill , r.  s. , hinshaw , g. , nolta , m.  r. , odegard , n. , page , l. , spergel , d.  n. , weiland , j.  l. , wright , e.  l. , halpern , m. , jarosik , n. , kogut , a. , limon , m. , meyer , s.  s. , tucker , g.  s. , & wollack , e. 2003 , , 148 , 97                                                                                , m.  c. , ade , p.  a.  r. , bhatia , r.  s. , bock , j.  j. , daub , m.  d. , goldstein , j.  h. , haynes , c.  v. , holzapfel , w.  l. , kuo , c.  l. , lange , a.  e. , leong , j. , lueker , m. , newcomb , m. , peterson , j.  b. , ruhl , j. , sirbi , g. , torbet , e. , tucker , c. , turner , a.  d. , & woolsey , d. 2003 , , submitted , astroph/0303515


the copernican principle implies that the earth , and , by extension , the solar system , do not hold a central or specifically favored position .
this viewpoint is related to the so - called _ mediocrity principle _ @xcite , which notes that an item drawn at random is more likely to come from a heavily populated category than one which is sparsely populated .
these principles , however , have not had particularly apparent success when applied in the context of extrasolar planets .
@xcite used their high precision doppler survey data to deduce that of order 50% ( or more ) of the chromospherically quiet main - sequence dwarf stars in the solar neighborhood are accompanied by a planet ( and in many cases , by multiple planets ) with @xmath3 , and orbital periods of @xmath4 .
taken strictly at face value , this result implies that our own solar system , which contains nothing interior to mercury s @xmath5 orbit , did not participate in the galaxy s dominant mode of planet formation . yet
the eight planets of the solar system have provided , and continue to provide , the _ de - facto _ template for most discussions of planet formation .    indeed ,
where extrasolar planets are concerned , m - dwarfs and mediocrity appear to be effectively synonymous .
recent observational results suggest that low - mass planets orbiting low - mass primaries are by no means rare .
numerous examples of planets with @xmath6 and m - dwarf primaries have been reported by the doppler surveys ( e.g. @xcite , and many others ) , and the _ kepler _ mission has indicated that small planets are frequent companions to low mass stars .
for example , @xcite report that among dwarf stars with @xmath7 , the occurrence rate of @xmath8 planets with @xmath9 is @xmath10 planets per star . improved statistics , however ,
are required for a definitive statement that is couched in planetary masses as well as in planetary radii .
figure [ fig : mdwarfcensus ] shows the current distribution of reported planets and planetary candidates orbiting primaries with @xmath11 , which we adopt as the functional border between `` m - type '' stars and `` k - type stars '' .
the census of low - mass planets orbiting low - mass primaries can be accessed using a variety of techniques . for objects near the bottom of the main sequence
, it appears that transit photometry from either ground @xcite or space @xcite offer the best prospects for planetary discovery and characterization . for early to mid m - type dwarfs
, there is a large enough population of sufficiently bright primaries that precise doppler detection ( see , e.g. @xcite ) can play a lead role . for the past decade
, we have had a sample of @xmath12160 nearby , photometrically quiet m - type stars under precision radial velocity surveillance with the keck telescope and its hires spectrometer . in recent months , this survey has been supplemented by data from the automated planet finder telescope @xcite . here , we present 16.6 years of doppler velocity measurements for the nearby m3 dwarf gj  687 ( including 122 velocity measurements from keck , 20 velocity measurements from the apf , and 5 velocity measurements made with the hobby - eberly telescope ) and we report the detection of the exoplanet that they imply .
we use this discovery of what is a highly archetypal representative of a planet in the milky way  in terms of its parent star , its planetary mass , and its orbital period  to motivate a larger discussion of the frequency of occurrence , physical properties , and detectability of low - mass planets orbiting m - type stars .
the plan for this paper is as follows . in  2 ,
we describe the physical and spectroscopic properties of the red dwarf host star gliese 687 . in  3 , we describe our radial velocity observations of this star . in
4 we describe our keplerian model for these observations , along with an analysis that assesses our confidence in the detection . in  5
, we describe our photometric time series data for the star , which aids in the validation of the planet by ruling out spot - modulated interpretations of the doppler variations . in  6 , we discuss the ongoing refinement of the planet - metallicity correlation for low - mass primaries , in  7 , we discuss the overall statistics that have emerged from more than 15 years of precision doppler observations of m - dwarf stars with the keck telescope , and in  8 we conclude with an overview that evaluates the important future role of the apf telescope in precision velocimetry of nearby , low - mass stars .
gliese 687 ( lhs 450 , bd+68@xmath13946 ) lies at a distance , @xmath14 pc , is the 39th - nearest known stellar system , and is the closest star north of @xmath15 declination .
figure [ fig : hrdiag ] indicates gj  687 s position in the color - magnitude diagram for stars in the lick - carnegie survey s database of keck observations . due to its proximity and its brightness ( v=9.15 , @xmath16=4.548 ) , gliese 687 has been heavily studied , and in particular , the chara array has recently been used to obtain direct interferometric angular diameter measurements for the star .
@xcite find @xmath17 , and derive @xmath18 , @xmath19 , and use the mass - radius relation of @xcite to obtain @xmath20 .
as illustrated in figure [ fig : sindex ] , gliese 687 s mean mt .
wilson s value and the dispersion of its s - index measurements from our spectra indicate that it has a moderate degree of chromospheric activity .
this conclusion is in concordance with our long - term photometric monitoring program , which also indicates that the star is somewhat active .
ccc spectral type & m3 v & @xcite + mass ( @xmath21 ) & @xmath22 & @xcite + radius ( @xmath23 ) & @xmath24 & @xcite + luminosity ( @xmath25 ) & @xmath26 & @xcite + distance ( pc ) & @xmath27 & @xcite + @xmath28 & 1.5 & simbad + @xmath29 mag .
& 9.15 & @xcite + @xmath30 mag .
& 5.335 & @xcite + @xmath31 mag .
& 4.77 & @xcite + @xmath32 mag .
& 4.548 & @xcite + avg .
s - index & 0.811 & this work + @xmath33 & 0.096 & this work + @xmath34 ( days ) & @xmath35 & this work + @xmath36 ( k ) & @xmath37 & @xcite +
doppler shifts from both the keck ( 122 observations ) and apf ( 20 observations ) platforms were measured , in each case , by placing an iodine absorption cell just ahead of the spectrometer slit in the converging beam of stellar light from the telescope @xcite .
the forest of iodine lines superimposed on the stellar spectra generates a wavelength calibration and enables measurement of each spectrometer s point spread function .
the radial velocities from keck were obtained by operating hires at a spectral resolving power r@xmath1270,000 over the wavelength range of 3700 - 8000 @xmath38 , though only the region 5000 - 6200 @xmath38 containing a significant density of iodine lines was used in the present doppler analysis @xcite .
the apf measurements were obtained over a similar spectral range , but at a higher spectral resolving power , r@xmath12108,000 .
for each spectrum that was obtained , the region containing the iodine lines was divided into @xmath12700 chunks , each of @xmath39 width .
each chunk produces an independent measure of the wavelength , psf , and doppler shift .
the final measured velocity is the weighted mean of the velocities of the individual chunks .
all radial velocities ( rvs ) have been corrected to the solar system barycenter , but are not tied to any absolute velocity system . as such , they are `` relative '' velocities , with a zero point that can float as a free parameter within an overall system model .    the internal uncertainties quoted for all the radial velocity measurements in this paper reflect only one term in the overall error budget , and result from a host of errors that stem from the characterization and determination of the point spread function , detector imperfections , optical aberrations , consequences of undersampling the iodine lines , and other effects .
two additional major sources of error are photon statistics and stellar `` jitter '' .
the latter varies widely from star to star , and can be mitigated to some degree by selecting magnetically inactive older stars and by time - averaging over the star s unresolved low - degree surface @xmath40-modes .
all observations in this paper have been binned on 2-hour timescales .
in addition to the radial velocities that we have obtained at keck and apf , we also use five doppler measurements obtained by @xcite at the hobby - eberly telescope located at mcdonald observatory .
these radial velocity observations are presented in the appendix .
the combined radial velocity data sets show a root - mean - square ( rms ) scatter of 7.58 @xmath41 about the mean velocity . this scatter is measured after we have applied best - fit telescope offsets of 0.64 @xmath42 for keck , -1.71 @xmath42 for apf , and 1.27 @xmath42 for het .
a lomb - scargle periodogram of the 149 velocity measurements of gj  687 is shown in figure [ fig : periodogram ] .
false alarm probabilities are calculated with the bootstrap method , as described in @xcite , iterating 100,000 times for a minimum probability of @xmath43 as easily met by the tallest @xmath44=38.14 day peak in figure [ fig : periodogram ] .
this signal in the data is modeled as a @xmath45 planet with an orbital eccentricity , @xmath46 .
lllr period ( d ) & 38.14&(0.015 ) + mass ( @xmath47 ) & 0.058&(0.007 ) + mass ( @xmath48 ) & 18.394&(2.167 ) + mean anomaly ( deg ) & 234.62&(87.962 ) + eccentricity & 0.04&(0.076 ) + longitude of periastron ( deg ) & 359.43&(120.543 ) + semi - major axis ( au ) & 0.16353&(0.000043 ) + & & & + time of periastron ( jd ) & 2450579.11&(9.32 ) + & & & + rv half amplitude ( @xmath42 ) & 6.43&(0.769 ) + & & + first observation epoch ( jd ) & 2450603.97 & + velocity offsets & & +  keck / hires & 0.64 @xmath49 & ( 0.63 ) +  apf / levy & -1.71 @xmath49 & ( 1.68 ) +  het & 1.27 @xmath49 & ( 0.98 ) + @xmath50 & 18.55 & + rms & & +  keck / hires & 6.62 @xmath49 & +  apf / levy & 3.95 @xmath49 & +  het & 2.44 @xmath49 & + jitter & 5.93 @xmath49 & +    using levenberg - marquardt optimization , we obtained a best - fit keplerian model for the system .
this fit , which assumes @xmath51 and @xmath52 for the planet , is listed in table [ tab : fit ] . the phased rv curve for the planet in table [ tab : fit ]
is shown in figure [ fig : phasedcurves ] .
a power spectrum of the residuals to our one - planet fit is shown in figure [ fig : residperiodogram ] and indicates no significant periodicities .
also shown in this figure is a periodogram of our mt .
wilson s - index measurements from the spectra , which are a proxy for the degree of spot activity on the star at a given moment .
none of the peaks in the periodogram of s - index values coincide with the peak that we suspect to be a planet .
the reduced chi - squared statistic for our fit is @xmath53 18.55 and results in a fit with a combined rms of 6.16 @xmath42 and estimated excess variance of @xmath54 = 5.93 @xmath42 ( the estimate of the jitter that is required to bring the reduced chi - squared statistic of the fit down to unity ) .
this value accounts for variance in both the stellar signal and from the telescope itself , though for a moderately active star such as gj  687 a stellar jitter of order 6.0 @xmath42 is reasonable , and could account for the majority of the observed variance .    in order to compute parameter uncertainties for our orbital fit ,
we implement a markov chain monte carlo algorithm @xcite .
the mcmc algorithm returns a chain of state vectors , @xmath55 ( a set of coupled orbital elements , e.g. period , mass , etc . and the three velocity offset parameters ) .
the goal of the markov chain calculation is to generate an equilibrium distribution proportional to @xmath56 $ ] .
we adopt non - informative priors on all parameters ( and uniform in the log for masses and periods ) .
the resulting error correlations are shown in figure [ fig : errorcorr ] , and a set of 100 states drawn randomly from the converged chain are shown in the orbital diagram in figure [ fig : orbit ] .
the error correlation diagram indicates that all parameters are well determined , save the usual degeneracies between mean anomaly and @xmath57 for the low - eccentricity orbit .
the distribution of the residuals relative to the best - fit model shows no evident pathologies .
indeed , a quantile - quantile plot ( shown in figure [ fig : fitplot ] ) indicates that the distribution of residuals is well described by a normal distribution .
we note that the smaller scatter of points obtained with the apf telescope could be a consequence of the fact that they were all taken within a @xmath58 = 140d period , and thus sample only one segment of the stellar activity cycle .    a potentially significant challenge to correctly identifying the orbital period of a proposed exoplanet arises from the discrete and uneven sampling inherent in radial velocity surveys .
the spacing of observations leads to increased noise and the presence of aliases within the star s periodogram which can be mistaken for a true orbital signature . for a real signal occurring at a frequency @xmath59
we expect alias signatures at @xmath60 where n is an integer . in order to aid confirmation
that the periodic signal we observe is actually a planetary signature , we must be able to calculate where aliases due to our observing cadence will occur , and then verify that they are not the source of the signal .
the aliases are determined using a spectral window function as defined by @xcite    @xmath61    where n is the total number of observations and t is the date on which they were taken . plotting this function will result in peaks that are due solely to the sampling cadence of the data .
because our observations are constrained by when the star is visible in the night sky , and because keck telescope time is allocated to doppler surveys primarily when the moon is up , we expect aliases at periods of 1 solar day , 1 sidereal day , 1 synodic month and 1 sidereal year .
examining the window function in figure [ fig : windowfunc ] we do see peaks resulting at these periods , but careful analysis of the periodogram for our radial velocity observations shows no evidence of strong signals occurring at the locations necessary for our @xmath1 day signal to be a potential alias instead of a true keplerian signature
.    with an apparent @xmath16-band magnitude of 4.54 , gliese 687 is brighter ( in the near infrared ) than _ all _ known hosts of transiting extrasolar planets other than 55 cancri . as a consequence , transits by gliese 687 s planetary companion ( which has an equilibrium temperature , @xmath62 ) ,
were they to occur , would be of substantial scientific value .
in particular , transmission spectroscopy with jwst would give insights into what is likely a dynamic and chemically rich planetary atmosphere .
the _ a - priori _ geometric transit probability for gliese 687  b , however , is a scant @xmath63% , and as we describe below , there is no evidence that transits occur . with @xmath64 ,
the currently observed mass - radius range for exoplanets indicates that the planetary radius , @xmath65 could credibly range from @xmath66 to @xmath67 , implying potential transit depths in the @xmath68% to @xmath69% range .
during the 20092013 observing seasons , we acquired a total of 866 photometric observations of gj  687 on 519 nights with the tennessee state university ( tsu ) t12 0.80  m automatic photoelectric telescope ( apt ) at fairborn observatory in arizona .
the t12 apt is one of several tsu automatic telescopes operated at fairborn @xcite .
it is equipped with a two - channel precision photometer that employs a dichroic filter and standard strmgren @xmath70 and @xmath71 filters to separate the two passbands and two emi 9124qb bi - alkali photomultiplier tubes to measure the @xmath70 and @xmath71 count rates simultaneously .
we observed gj  687 , designated our program star ( p ) , differentially with respect to three neighboring comparison stars : c1 ( hd  156295 , @xmath72 , @xmath73 , f0  iv ) , c2 ( hd  160198 , @xmath74 , @xmath75 , f2  v ) , and c3 ( hd  161538 , @xmath76 , @xmath77 , f2  v ) .
a detailed description of the observing sequence and the data reduction and calibration procedures are given in @xcite .
we computed all pairwise differential magnitudes @xmath78 , @xmath79 , @xmath80 , @xmath81 , @xmath82 and @xmath83 in both the @xmath70 and @xmath71 passbands , corrected them for atmospheric extinction , and transformed them to the standard strmgren photometric system .
observations with internal standard deviations greater than 0.01 mag were discarded to remove data taken in non - photometric conditions .
intercomparison of the six sets of differential magnitudes demonstrated that hd  156295 ( c1 ) is a low - amplitude variable while both hd  160198 ( c2 ) and hd  161538 ( c3 ) are constant to the expected measurement precision . to improve our precision , we combined the separate differential @xmath70 and @xmath71 observations into a single @xmath84  passband . "
we also computed the differential magnitudes of gj  687 with respect to the mean brightness of the two good comparison stars : @xmath85 .
the standard deviation of the @xmath81 comparison star differential magnitudes is 0.0020 mag , which we take to be the precision of a single measurement .
a total of 606 nightly measurements in the five observing seasons survived the cloud - filtering process .
these data are plotted as @xmath85 differential magnitudes in the top panel of figure [ fig : phot1 ] .
the five individual observing seasons are plotted in the remaining panels .
the standard deviations for the yearly light curves are given in each panel .
these range from 0.0049 to 0.0092 mag , compared to the measurement precision of 0.0020 mag .
gaps of 1012 weeks in the yearly light curves for 2009 through 2012 are due to southern arizona s july
september rainy season when good photometry is not possible .
low - amplitude variability is seen in gj  687 during each observing season , resembling light curves typical of modestly active stars with spot filling factors of a few percent ( see , e.g. , * ? ? ?
the 2010 light curve has the largest amplitude variability ( @xmath86  mag ) and reveals cyclic variation with a time scale of @xmath87  days .
the other light curves have lower amplitudes and include cyclic variations of @xmath87 and also @xmath88 days .
these year - to - year and cycle - to - cycle variations are also typical of modestly active stars .
we interpret the 60-day variability as the signature of the star s rotation period and the 30-day variability as a sign of spot activity on opposite hemispheres of the star .
frequency spectra of the complete 2009 - 2013 data set and of the 2010 data alone are shown in the top and bottom panels of figure [ fig : phot2 ] respectively .
the rotational modulation signal is seen most clearly in the 2010 data , which matches up with the most coherent light curve in figure [ fig : phot1 ] .
therefore , we take the @xmath89d signal as our best measurement of the star s rotation period
. inspection of the 2010 photometric segment of figure [ fig : phot1 ] clearly shows the overall 60-day modulation that generates the periodogram peak .
departures from perfect periodicity are presumably caused by the evolution of the spot activity on the surface of the star .    finally , we search for transits of gj  687  b by first removing the spot variability from each of the yearly light curves .
we do this by successively subtracting multiple frequencies from each yearly light curve using the method described in @xcite .
we removed three to six frequencies from each light curve until each set of residuals approached the precision of a single observation .
the residuals from all five observing seasons are plotted in the top panel of figure [ fig : phot3 ] , phased with the 38.14-day best - fit planetary orbital period and a time of mid transit computed from the orbital parameters .
the vertical bar represents the 0.0022 mag standard deviation of the residuals from their mean , very close to the measurement precision given above
. a sine fit to the phased data gives a formal semi - amplitude of just @xmath90  mag .
since none of the frequencies removed from the yearly light curves were similar to the orbital frequency or its harmonics , this result limits any periodic brightness variability of the star on the observed radial velocity period to a very small fraction of one milli - magnitude ( mmag ) .
this rules out the possibility that the 38.14-day radial velocity variations in gj  687 are induced by stellar activity , as has been documented in somewhat more active stars , for instance , by @xcite , @xcite , and @xcite .
instead , this lack of photometric variability confirms that the radial velocity variations in gj  687 result from true planetary reflex motion .
the photometric observations within @xmath91p of mid - transit are replotted with an expanded scale in the bottom panel of figure [ fig : phot3 ] .
the solid curve shows the predicted phase , depth ( assuming neptune - like density ) , and duration of a central transit , computed from the stellar radius in table [ tab : stellarparams ] and the orbital elements in table [ tab : fit ] .
the horizontal error bar under the predicted transit time gives the @xmath92 uncertainty in the timing of the transit .
the photometric observations when filtered using the @xcite procedure described above , and when folded at the @xmath1 day best - fit period for the planet , give no indication that transits occur .
we note , however , that the markov chain models generate a five - day window for possible transits , and so a more conservative approach is also warranted . in figure
[ fig : transit ] , we plot the unfiltered photometric data , indicating the range of photometric points that potentially could have been affected by transits were they to occur . because of uncertainties in the orbit , the potential transit duration , the potential size of the planet , and the error in the photometric filtering , we recommend that continued photometric monitoring be carried out to confirm that transits do not occur .
gliese 687 appears to have a slightly sub - solar metallicity .
@xcite use na i , ca i , and h@xmath93o - k@xmath94 calibrations to estimate [ fe / h]=-0.09 for gliese 687 , whereas the m - dwarf metallicity calibration of schlaufman & laughlin 2010 yields a value [ fe / h]=-0.02 .
the connection between the detectable presence of a giant extrasolar planet and the metallicity of the host star was noticed soon after the first extrasolar planets were detected @xcite , and has been studied in many previous works , see , e.g. @xcite
. for m dwarfs , recent work , such as that by @xcite , suggests that the giant planet stellar metallicity correlation holds robustly for m - dwarf primaries , but that for planets with mass , @xmath95 , no correlation is found with host star metallicity , and indeed , @xcite report a hint of _ anti_-correlation between the presence of a low - mass planet and host star @xmath96 $ ] .
our detection of a neptune - mass companion to gliese 687 , and our lick - carnegie database of doppler measurements of m dwarf stars provides an opportunity to revisit this topic .
our database of radial velocity observations taken at the keck telescope contains 142 m - type stars with the necessary spectral information to assess metallicity , 17 of which are known to host planets published in the peer - reviewed literature .
we break the planet - hosting stars into two subgroups based on their masses - stars with @xmath97 planets less than 30 @xmath98 are described as neptune hosting while stars with @xmath97 planets greater than 30 @xmath98 are listed as jupiter hosting .
we replicate the procedure of @xcite and examine how horizontal distance from a field m dwarf main sequence in a @xmath99 vs.  @xmath100 color - magnitude diagram ( cmd ) correlates with metallicity , as noted e.g. , by @xcite .
the top panel of figure [ fig : metallicity ] displays all of the lick - carnegie survey m dwarf stars plotted in @xmath99 vs.  @xmath100 , with grey dots denoting survey stars without known planets , red dots denoting survey stars that host `` neptune - mass '' planets and blue dots representing the survey stars that host `` jupiter - mass '' planets .
it can be seen that most planet hosting stars fall to the right of the field m dwarf main sequence presented in @xcite ( black line ) , which is taken to be a [ fe / h]=0.017 isometallicity contour in this cmd . in order to quantify the likelihood that a star s horizontal distance from the isometallicity contour is related to its propensity to host planets , we compare the distances for our actual planet - hosting stars with randomly drawn samples from the collection of m dwarfs in the survey .
we characterize the position of each m dwarf by obtaining v - band and @xmath101 photometry and then using them to calculate the distance statistic @xmath102 :    @xmath103    to determine if the @xmath102 of our known planet hosting subgroups is significant or , alternatively , if it could be produced by chance , we make use of a monte carlo simulation that calculates the cumulative sample distance of survey m dwarfs from the field m dwarf main sequence presented by johnson and apps 2009 . for the simulation , we randomly select a subset of m dwarfs from the lick - carnegie field star list , setting the sample size equal to the number of m dwarfs known to host either jovian or neptune mass planets .
then we compute the cumulative horizontal distance of those stars from the field m dwarf ms , where stars to the right of the ms add their distance to the sum and stars to the left of the ms subtract their distance .
we repeat this process 10,000 times to determine the distribution of cumulative horizontal distances from the ms given no correlation between whether the star hosts an exoplanet and its location in the @xmath104 cmd .
our results show that for hosts of jupiter mass planets , @xmath105 , which corresponds to a probability of @xmath106 that the stars cumulative distance from the isometallicity contour occurred by chance . for the neptune hosts
, we find that @xmath107 leading to @xmath108 , and for the combination of all planet hosts , we obtain @xmath109 or @xmath110 .
the distributions resulting from the monte carlo simulation and the locations of the actual planet hosting stars @xmath102 values are displayed in the bottom panel of figure [ fig : metallicity ] .
the points plotted on top of each curve in figure [ fig : metallicity ] represent the actual cumulative distance of our planet hosting star samples from the field m dwarf ms .
our results thus indicate that the planet - metallicity correlation is robust for m - dwarf hosts of planets with @xmath111 , but that at smaller masses there is , at present , no evidence a correlation exists .
the lick - carnegie exoplanet survey and its predecessors have carried out a long - term monitoring program of the brightest m - dwarf stars in the sky .
our database of observations contains 159 stars that have more than 10 observations apiece , and which , additionally , have median internal uncertainty @xmath112 . within this group ,
there is a subset with extensive data sets .
for example , 11 stars have @xmath113 observations and median internal uncertainties @xmath114 .
a question of substantial interest , therefore , is the degree to which the observations taken to date have probed the true aggregate of planetary companions to the m - dwarf stars in our survey .
the effort required to obtain the existing data has been substantial . among the m - type stars alone ,
our database contains a total of 5,468 velocity measurements from keck i , totaling 2,579,862 seconds ( 29.86 days ) of on - sky integration .
overheads , including the acquisition of high s / n spectra , ccd readout time , and weather losses , add materially to this time investment .
furthermore , the distribution of total observing time allotted to the stars on the list has been highly uneven .
targets such as gliese 436 and gliese 876 , which harbor planetary systems of particular interest , have received much more attention than the typical red dwarf in the survey .
for example , gliese 436 has 148 observations and gliese 876 has 204 observations obtained with the keck telescope .
the stars themselves also exhibit a range of chromospheric activity levels . the resulting star - to - star dispersion in
 stellar jitter " ( tantamount to a measurement uncertainty , @xmath115 ) complicates the evaluation of threshold levels for @xmath97 as a function of orbital period to which planetary companions can be excluded .
there are a variety of approaches to the measurement of false alarm probabilities ( fap ) in the context of spectral analysis of unevenly sampled data .
see , e.g. @xcite for a recent discussion .
a very simple approach is described by @xcite . for a gaussian random variable . ] , the probability distribution for obtaining a peak at frequency @xmath57 of lomb - normalized power @xcite , @xmath116 , is exponential with unit mean .
if a data set drawn from measurements of a white noise ( gaussian ) distribution supports measurement of @xmath117 independent frequencies , the probability that no peak exceeds power @xmath118 ( the fap ) is @xmath119 .
we adopt a fap of @xmath120 , calculated with the above method ( and using monte - carlo simulations to determine @xmath117 ) as the generic threshold for attributing a given planetary signal to a given dataset . with this detectability threshold ,
we use the systemic console 2.0 software package @xcite to determine the number of readily detectable planets in our m - dwarf data set .
a `` readily detectable '' planet generates a signal that can be isolated algorithmically ( and automatically ) by straightforward periodogram analysis and levenberg - marquardt minimization .
the results of this exercise are shown in figure [ fig : mdwarftable ] , which locates signals corresponding to 19 previously published planets orbiting 14 separate m - dwarf primaries .
other than gliese 667c , there are no stars on our 159-star list for which a planet has been published by another group , and for which the automated algorithm finds no planets .
regarding gj  667c , 40 observations have been made at keck , and these were used in a characterization of the gj  667c system @xcite , however the peak planetary signal for this set fell below our fap threshold when utilizing only the keck data .
the bright planet - hosting red dwarfs gliese 832 , 3634 , and 3470 all have declinations that are too far south to be observed from mauna kea , and hip 79431 ( ra 16h 12 m 41.77s dec -18@xmath13 52@xmath121 31.8@xmath122 ) is not on the list of m - dwarfs being monitored at keck .
the kepler mission s photometric data have been used to infer that small planets orbiting m - dwarfs are very common . for example , @xcite find an occurrence rate of 0.9 planets per star in the range @xmath8 with @xmath123 days .
given the existence of this large number of small - radius planets , it is of interest to make a quantitative analysis of how deep into the expected population of super - earth type planets suggested by the kepler mission the keck radial velocity survey has probed . to answer this question ,
we have created synthetic radial velocity data sets that contain test planets , and which conform with the timestamps , the internal measurement uncertainties , and the stellar properties ( namely mass ) for all 104 m - dwarfs under surveillance at keck with at least 20 radial velocity observations . to address the error source arising from stellar jitter , @xmath124
, we use the median value provided in @xcite of 3.9 @xmath42 as the expected level of @xmath124 for our m - dwarf stars .
this value is then added in quadrature with the internal uncertainty and applied to the synthetic data set to create a more accurate representation of the system .
for each of the 104 stars in the keck survey , we have created 400 synthetic data sets .
each set contains a single planet .
the planets are evenly spaced in log period from 2 to 100 days , and evenly spaced in log mass from 1 @xmath48 to 1 @xmath125 .
we assign a circular orbit to these test planets and assume @xmath126 and @xmath52 in each case .
we then calculate the radial velocity each of these planets would induce on a parent star .
a gaussian distribution with @xmath127 is used to perturb the predicted radial velocity value .
each of the @xmath128 synthetic systems is passed to the planet search algorithm .
figure [ fig : massperpanel ] shows examples of the returned planets for 4 of the 104 stars in this experiment including the star of main interest here , gj  687 . in figure [ fig : massperpanel ] the black lines from bottom to top represent constant @xmath129 values of 1 , 10 , and 100 respectively .
as expected , the detectablity thresholds lie roughly along lines of constant @xmath32 . to determine the smallest @xmath32 value we could reliably detect for each star , we find the smallest value of @xmath32 for which a planet was found for at least 50% of the chosen periods .
this median value generates the green lines seen in the figure .
the top panel of figure [ fig : mink ] shows these minimum @xmath32 values for each of the 104 stars that we tested . for clarity ,
the stars in this figure have been ordered by increasing minimum @xmath32 , and are colored by the number of observations we have for each .
if a test planet lies in its star s habitable zone , ( defined as the semi - major axis at which the flux received by the planet is the solar constant received at earth ) we can ask how large the planet needs to be to be detectable by our radial velocity survey .
figure [ fig : mink ] shows these threshold masses for each of the 104 m - dwarf stars which we analyzed .
these stars maintain the ordering from the top panel , but now have been colored by the mass of the parent star .
we see that while the keck survey has probed substantially into the regime occupied by neptune - mass planets , it has not made significant inroads into the super - earth regime for periods that are of astrobiological interest .
gj  687 is the second planetary system to be detected using data from the apf telescope , with the first being hd 141399 b , c , d , and e @xcite .
apf has successfully navigated its commissioning stage , and , since q2 2013 , it has routinely acquired science - quality data that presents sub - m / s precision on known radial velocity standard stars @xcite . in recent months
, the degree of automation for apf has increased substantially .
the facility currently works autonomously through an entire night s operations , calibration , and observing program .
the apf and its accompanying high - resolution levy spectrograph together form a dedicated , cost - effective , ground - based precision radial velocity facility that is capable of detecting terrestrial - mass planets at distances from their parent stars at which surface liquid water could potentially be present .    unlike other highly successful rv facilities
, the apf uses neither image scrambling nor image slicing . with a peak efficiency of 15% and typical spectral resolutions of @xmath130 ,
the apf represents a critical new resource in the global quest to detect extrasolar planets .
initial speed comparisons indicate that in order to match the signal - to - noise acquired using the keck telescope / hires spectrograph combination , the apf needs only a factor of 6 increase in observing time . since the amortized cost of a night on keck is @xmath131 more expensive than a night on the apf , and because 80% of the apf s nights are reserved for exoplanetary work , the apf ( with its sub - m / s precision and dedicated nightly cadence abilities ) will likely provide key contributions to exoplanet detection and characterization in the coming years .
gliese 687  b s radial velocity half - amplitude , @xmath132 , is substantially greater than the current state - of - the - art detection threshold for low mass planets .
the lowest measured value for @xmath32 in the catalog of doppler - detected extrasolar planets stands at @xmath133 @xcite . on the other hand , gliese 687 s status as one of the nearest stars to the sun imbues it with a great deal of intrinsic interest . in our view , the relatively recent date for gliese 687  b s detection can be attributed both to the substantial amount of stellar - generated radial velocity noise ( as evidenced by figures [ fig : sindex ] and [ fig : phasedcurves ] ) , but also to its location in draco , high in the northern sky , where apf , along with harps north , are the only facilities that can routinely observe at sub-@xmath134 precision .
( as evidenced by the data in this paper , keck can observe at these high declinations , but at significantly higher expense in comparison to stars lying closer to the celestial equator . )    indeed , gliese 687 s stellar coordinates ( r.a .
17h , 36 m , dec @xmath135 ) place it very close to the north ecliptic pole , located at @xmath136h , dec@xmath137 .
this location flags it as a star of potentially great importance for the forthcoming nasa tess mission . as currently envisioned ( and as currently funded ) , tess is a two - year , all - sky photometric survey to be carried out by a spacecraft in a 27-day p/2 lunary resonant orbit .
tess will photometrically monitor @xmath138 bright stars with a @xmath139 1-hour systematic error floor .
( for reference , a central transit of the sun by the earth produces an 86 ppm transit depth . )
the northern ecliptic hemisphere will be mapped during the first year of the mission via a sequence of 13 sectors with 27 days of continuous observation per sector .
these sectors overlap at the north ecliptic pole , and create an area of @xmath121,000 square degrees (  1/50th of the sky ) for which photometric baselines will approach 365 days .
gliese 687 lies at the center of this tess `` overlap zone '' ( which also coincides with jwst s continuous viewing zone ) .
because much longer time series are produced in the overlap zone , the highest - value transiting planets found by the mission will emerge from this part of the sky ( along with the sister segment covering the south ecliptic pole ) .
as mentioned above , however , the tess overlap zone has received relatively little attention from the highest - precision doppler surveys .
about 10,000 target stars from the tess dwarf star catalog , all with v@xmath14012 , are present in the overlap zone .
this , of course , is far too many stars to survey with doppler rv , but there appears to be substantial value inherent in monitoring the _ brightest _ , _ nearest _ , and _
quietest _ members of the cohort of tess overlap stars .
the latest estimates @xcite suggest that @xmath1250% of main sequence stars in the solar vicinity harbor @xmath141 planets with @xmath4 .
assuming a uniform distribution in period between 5 and 100 days , the average transit probability for these planets is @xmath142 , suggesting that of order @xmath143 low mass transiting planets ( and systems of transiting planets ) will be detected by tess within the overlap zone . of these , a small handful , of order 5 systems total ( and perhaps , with probability @xmath144% , including gliese 687  b ) will garner by far the most attention from follow - up platforms such as jwst , due to their having optimally bright parent stars .
our detection of gliese 687  b suggests that by starting now , with a systematic program of doppler observations of a target list of @xmath12200 carefully vetted g , k & m dwarf stars with v@xmath127.5 to v@xmath1210.5 in the 1000-square degree tess overlap zone , apf can ensure that a precise multi - year doppler velocity time series will exist for the most important tess planet host stars at the moment their transiting planets are discovered .
gl acknowledges support from the nasa astrobiology institute through a cooperative agreement between nasa ames research center and the university of california at santa cruz .
ssv gratefully acknowledges support from nsf grants ast-0307493 and ast-0908870 .
rpb gratefully acknowledges support from nasa oss grant nnx07ar40 g , the nasa keck pi program , and from the carnegie institution of washington . s.m .
acknowledges support from the w. j. mcdonald postdoctoral fellowship .
gwh acknowledges support from nasa , nsf , tennessee state university , and the state of tennessee through its centers of excellence program .
the work herein is based on observations obtained at the w. m. keck observatory , which is operated jointly by the university of california and the california institute of technology , and we thank the uc - keck , uh , and nasa - keck time assignment committees for their support .
this research has made use of the keck observatory archive ( koa ) , which is operated by the w. m. keck observatory and the nasa exoplanet science institute ( nexsci ) , under contract with the national aeronautics and space administration .
we thank those who collected the data in the koa including : d. fischer , g. marcy , j. winn , w. boruki , g. bakos , a. howard , and t. bida .
we also wish to extend our special thanks to those of hawaiian ancestry on whose sacred mountain of mauna kea we are privileged to be guests . without their generous hospitality , the keck observations presented herein would not have been possible .
this research has made use of the simbad database , operated at cds , strasbourg , france .
this paper was produced using @xmath145 .    , , .
ccc 2450603.965 & -15.930 & 1.770 + 2450956.065 & -4.740 & 2.040 + 2450982.977 & -12.880 & 1.560 + 2451013.879 & -15.550 & 1.610 + 2451312.036 & -6.870 & 1.910 + 2451368.798 & -17.330 & 1.790 + 2451439.745 & -25.500 & 2.170 + 2451704.957 & -14.920 & 1.810 + 2452007.029 & -3.280 & 1.970 + 2452009.077 & -7.440 & 2.020 + 2452061.913 & 7.450 & 2.040 + 2452062.935 & 4.490 & 1.940 + 2452094.853 & 2.450 & 1.550 + 2452096.894 & 4.560 & 1.850 + 2452097.982 & 5.180 & 1.610 + 2452127.919 & 1.940 & 2.200 + 2452133.737 & 2.830 & 2.010 + 2452160.872 & -6.120 & 2.170 + 2452161.821 & -6.840 & 2.370 + 2452162.787 & -8.300 & 2.370 + 2452445.988 & 17.880 & 2.100 + 2452537.743 & -4.900 & 2.340 + 2452713.113 & 8.860 & 1.590 + 2452806.027 & -4.580 & 1.990 + 2452850.901 & 6.090 & 2.040 + 2453179.985 & 12.360 & 1.730 + 2453479.066 & 8.000 & 1.210 + 2453549.857 & -2.830 & 0.880 + 2453604.881 & -2.530 & 1.220 + 2453838.109 & 3.550 & 1.220 + 2453932.905 & 3.890 & 1.380 + 2453960.873 & 3.880 & 1.630 + 2453961.820 & -0.510 & 1.800 + 2453981.789 & -14.390 & 1.190 + 2453982.904 & -1.880 & 1.900 + 2453983.831 & -8.010 & 1.220 + 2453984.887 & -0.360 & 1.290 + 2454248.023 & 7.800 & 2.090 + 2454248.990 & 8.410 & 2.030 + 2454249.945 & -3.960 & 1.730 + 2454252.032 & -12.260 & 1.870 + 2454255.924 & -12.790 & 1.280 + 2454277.851 & 8.200 & 1.500 + 2454278.896 & 11.300 & 1.700 + 2454279.933 & 11.850 & 1.580 + 2454294.903 & -10.850 & 1.260 + 2454304.888 & -4.110 & 1.570 + 2454306.037 & 2.450 & 1.600 + 2454307.015 & -2.050 & 1.500 + 2454308.072 & 5.540 & 1.600 + 2454309.050 & 7.370 & 1.810 + 2454310.042 & 8.730 & 1.720 + 2454311.024 & 1.600 & 1.480 + 2454312.017 & 0.860 & 1.590 + 2454312.875 & 2.350 & 1.450 + 2454313.875 & -1.170 & 1.320 + 2454314.908 & 6.220 & 1.630 + 2454318.921 & 5.360 & 1.560 + 2454335.826 & -13.560 & 1.940 + 2454338.888 & -16.490 & 1.790 + 2454339.883 & -9.590 & 1.780 + 2454343.791 & -8.170 & 1.950 + 2454396.717 & 3.080 & 2.390 + 2454397.729 & -4.380 & 2.130 + 2454548.052 & 12.600 & 1.330 + 2454549.092 & 3.010 & 1.280 + 2454633.963 & -6.780 & 1.570 + 2454634.903 & -5.330 & 1.650 + 2454635.946 & -2.500 & 1.790 + 2454636.901 & 2.170 & 1.650 + 2454637.939 & -5.830 & 1.610 + 2454638.838 & -5.010 & 1.450 + 2454639.998 & 3.970 & 1.850 + 2454641.954 & 2.760 & 1.540 + 2454674.860 & -0.820 & 1.440 + 2454688.892 & -4.450 & 1.710 + 2454689.917 & 6.340 & 1.890 + 2454717.853 & 4.220 & 2.000 + 2454718.914 & -2.870 & 2.070 + 2454719.866 & 1.210 & 2.040 + 2454720.871 & 0.570 & 2.140 + 2454721.879 & 4.620 & 2.040 + 2454722.797 & 5.810 & 2.210 + 2454723.814 & 2.610 & 2.120 + 2454724.848 & 1.350 & 1.990 + 2454968.017 & 5.210 & 0.970 + 2455016.035 & 0.730 & 1.370 + 2455022.084 & 1.680 & 1.670 + 2455023.044 & 6.360 & 1.070 + 2455024.823 & -11.730 & 0.920 + 2455049.825 & 0.220 & 1.090 + 2455050.858 & 3.990 & 1.650 + 2455051.772 & 6.360 & 1.790 + 2455052.780 & 8.540 & 1.930 + 2455053.925 & -12.540 & 1.510 + 2455143.692 & -3.080 & 2.620 + 2455167.704 & 3.210 & 1.710 + 2455259.119 & 3.760 & 1.610 + 2455260.144 & 0.000 & 1.470 + 2455371.991 & -5.030 & 1.070 + 2455407.955 & -1.480 & 1.060 + 2455463.718 & -8.590 & 1.520 + 2455516.696 & -1.120 & 1.740 + 2455518.722 & -3.780 & 2.320 + 2455609.135 & 11.250 & 1.310 + 2455638.121 & -3.930 & 0.910 + 2455639.143 & -2.250 & 0.910 + 2455665.046 & -5.130 & 1.190 + 2455670.082 & -8.090 & 1.290 + 2455721.057 & 6.840 & 1.850 + 2455825.832 & 3.030 & 2.280 + 2455840.748 & 3.820 & 1.390 + 2456027.098 & 11.390 & 1.940 + 2456116.938 & 16.010 & 1.130 + 2456117.906 & 14.270 & 1.050 + 2456329.167 & 5.420 & 1.970 + 2456432.892 & -0.700 & 1.300 + 2456433.935 & 1.610 & 1.340 + 2456548.855 & -2.340 & 1.490 + 2456549.793 & -7.780 & 1.430 + 2456550.849 & -3.840 & 2.150 + 2456551.724 & -2.520 & 1.570 +      ccc 2456484.738 & 11.650 & 0.910 + 2456486.839 & 9.500 & 0.950 + 2456488.845 & 11.640 & 1.120 + 2456490.726 & 8.300 & 1.260 + 2456492.816 & 7.020 & 0.900 + 2456494.776 & 5.260 & 0.750 + 2456508.764 & -8.650 & 0.780 + 2456510.800 & -9.020 & 0.800 + 2456513.772 & -8.920 & 0.860 + 2456521.823 & -2.120 & 0.730 + 2456522.790 & -2.920 & 0.930 + 2456533.797 & 2.910 & 0.820 + 2456551.743 & -8.410 & 0.810 + 2456552.717 & -6.660 & 0.970 + 2456553.727 & -7.610 & 0.850 + 2456554.720 & -4.640 & 0.890 + 2456564.762 & 1.530 & 0.940 + 2456568.692 & 0.460 & 1.140 + 2456571.657 & 0.580 & 0.940 + 2456622.032 & -9.750 & 1.280 + 2456677.921 & 2.390 & 0.980 + 2456681.011 & 4.660 & 1.360 +                                                                                , s.  s. , allen , s.  l. , bigelow , b.  c. , et  al .
1994 , in society of photo - optical instrumentation engineers ( spie ) conference series , vol . 2198 , society of photo - optical instrumentation engineers ( spie ) conference series , ed .
d.  l. crawford & e.  r. craine , 362+