planets that transit their stars offer us the opportunity to study the physics of planetary atmospheres and interiors , which may help constrain theories of planet formation . from the photometric light curve
, we can measure the planetary radius and also the orbital inclination , which when combined with radial velocity ( rv ) observations , allows us to measure the mass and density of the planet .
infrared observations of planets during secondary eclipse can be used to measure the planet s thermal spectrum @xcite , and spectroscopic observations around the time of transit can constrain the composition of the planet s atmosphere through transmission spectroscopy @xcite .
additionally , optical observations of the secondary eclipse can probe a planet s albedo @xcite .
all of the known transiting planets are hot jupiters or hot neptunes , which orbit so close to their parent stars that the stellar flux plays a major role in heating these planets .
in contrast , the detection of transiting planets with longer periods ( @xmath1 ) and consequently lower equilibrium temperatures , would allow us to probe a completely different regime of stellar insolation , one more like that of jupiter and saturn whose energy budgets are dominated by their internal heat .    currently , it is not feasible to find long - period , transiting planets from the ground .
while rv surveys have found @xmath2 planets with periods @xmath3 , it is unlikely that any of these transit their parent star , given that for a solar - type host , the transit probability is very small , @xmath4 .
thus the sample of long - period planets detected by rv will need to at least double before a transiting planet is expected .
given the long lead time necessary to detect long - period planets with rv observations , this is unlikely to happen in the near future .
furthermore , the expected transit times for long - period systems are relatively uncertain , which makes the coordination and execution of photometric follow - up observations difficult . putting all this together ,
radial velocity surveys are clearly an inefficient means to search for long - period , transiting planets .
ground - based photometric transit surveys are similarly problematic . a long - period system must be monitored for a long time because typically , at least two transits must be observed in order to measure the period of the system .
furthermore , if there are a couple of days of bad weather at the wrong time , the transit event will be missed entirely , rendering years of data worthless .
thus transit surveys for long - period planets require many years of nearly uninterrupted observations .
additionally , many thousands of suitable stars must be monitored to find a single favorably - inclined system .
because of its long mission lifetime ( @xmath5 ) , nearly continuous observations , and large number of target stars ( @xmath6 ) , the _ kepler _ satellite @xcite has a unique opportunity to discover long - period transiting systems
. not only will _ kepler _ observe multiple transits of planets with periods up to the mission lifetime , but it is also likely to observe single transits of planets with periods longer than the mission lifetime . as the period of a system increases beyond @xmath7 , the probability of observing more than one transit decreases until , for periods longer @xmath8 , only one transit will ever be observed . for periods longer than @xmath8 , the probability of seeing a single transit diminishes as @xmath9 .
even with only a single transit observation from _ kepler _ , these long - period planets are invaluable .
we show that planets with periods longer than the mission lifetime will likely be detected and can be characterized using the _ kepler _ photometry and precise radial velocity observations .
furthermore , this technique can be applied to planets that will transit more than once during the _ kepler _ mission , so that targeted increases in the time - sampling rate can be made at the times of subsequent transits .    in @xmath10[sec : np ] we calculate the number of one- and two - transit systems _ kepler _ is expected to observe .
we give a general overview of how the _ kepler _ photometry and radial velocity follow - up observations can be combined to characterize a planet in @xmath10[sec : basic ] .
we discuss in @xmath10[sec : lc ] the expected uncertainties in the light curve observables , and we relate these to uncertainties in physical quantities , such as the period , that can be derived from the transit light curve .
@xmath10[sec : rv ] describes the expected uncertainties associated with the rv curve and how they influence the uncertainty in the mass of the planet .
@xmath10[sec : e ] discusses the potential impact of eccentricity on the ability to characterize long - period planets .
we summarize our conclusions in @xmath10[sec : conclude ] .
details of the derivations are reserved for an appendix ( @xmath10[sec : append ] ) .
the number of single - transit events _ kepler _ can expect to find is the integral of the probability that a given system will be favorably inclined to produce a transit , multiplied by the probability that a transit will occur during the mission , convolved with a distribution of semimajor axes , and normalized by the expected frequency of planets and the number of stars that _ kepler _ will observe ,    @xmath11    here @xmath12 is the total number of transiting systems expected , @xmath13 is the number of stars being monitored , @xmath14 is the probability that the system is favorably inclined to produce a transit ( the transit probability ) , @xmath15 is the probability that a transit will occur during the mission , and @xmath16 is some distribution of the planet semimajor axis @xmath17 , normalized to the expected frequency of planets .    assuming a circular orbit ,
the transit probability is simply , @xmath18 where @xmath19 and @xmath20 are the radius and mass of the parent star , and @xmath21 is the period of the system .
we consider both the probability of observing exactly one transit , @xmath22 , and the probability of observing exactly two transits , @xmath23 .
these probabilities are @xmath24 @xmath25 thus there will be a range in periods from @xmath7 to @xmath8 where it is possible to get either one or two transits during the mission . if we assume a mission lifetime of @xmath26 years and a solar - type star ( @xmath27 ) , we can find the total probability ( @xmath28 ) of observing a planet that transits _ and _ exhibits exactly one or exactly two transits as a function of period or semimajor axis .
these probability distributions as functions of @xmath17 are shown in the middle panel of fig .
[ fig : np ] .    for @xmath29 ,
the total probability , @xmath30 , of observing a single transit is @xmath31 thus the probability that the planet will transit _ and _ that a single transit will occur during the mission is generally small for periods at the mission lifetime , and beyond this is a strongly decreasing function of the period , @xmath32 or @xmath33 .    the observed distribution of semimajor axes for planets with minimum masses @xmath34 is shown in the top panel of figure [ fig : np ] , where @xmath35 is the mass of the planet
, @xmath36 is the orbital inclination , and the data are taken from the `` catalog of nearby exoplanets '' of @xcite .
note that we have included multiple - planet systems .
we normalized the distribution by assuming that 8.5% of all stars have planets with @xmath34 and @xmath37 @xcite .
we can extrapolate this distribution to larger semimajor axes by observing that it is approximately constant as a function of @xmath38 and adopting the value found by @xcite , @xmath39 . using this extrapolation
, we predict @xmath40 and @xmath41 of stars have a @xmath34 planet within @xmath42 and @xmath43 , respectively , in agreement with @xcite .
we find the number of planets expected for _ kepler _ by convolving this normalized distribution with the probability distribution and multiplying by @xmath44 .
we convolve the distribution of semimajor axes with the total transit probability distribution by weighting the probability of observing a transit for each member of the bin to determine an overall probability for the bin . the bottom panel of fig .
[ fig : np ] shows the resulting distribution of one- and two - transit systems . integrating this distribution over all semimajor axes out to the limit of current rv observations
predicts that _ kepler _ will observe a total of 4.0 single - transit systems and 5.6 two - transit systems during the mission lifetime of @xmath45 .
if we include the extrapolation of the observed distribution to larger semimajor axes , the integrated number of single - transit systems increases to 5.7 ; the effect is modest because the probability of observing a transit declines rapidly with increasing semimajor axis ( @xmath46 ) .
thus , _ kepler _ is likely to detect at least handful of single - transit events .
it is important to bear in mind that this is a lower limit .
first , our estimate does not include planets with @xmath47 , primarily because rv surveys are substantially incomplete for low - mass , long - period planets .
as we will show , _ kepler _ will be able to constrain the periods of planets with radii as small as that of neptune with the detection of a single transit .
indeed , microlensing surveys indicate that cool , neptune - mass planets are common @xcite .
second , @xcite and @xcite showed that a distribution in eccentricity increases the probability that a planet will transit its parent star .
finally , if the _ kepler _ survey lifetime is extended , the number of detections of long - period transiting planets will also increase .
one might assume that a single - transit event must simply be discarded from the sample because the event is not confirmed by a second transit , and the period can not be constrained by the time between successive transits .
although we do not expect to see many of them , long - period transiting planets are precious , as the information we could potentially gain by observing these systems when they transit in the future would allow us to greatly enhance our understanding of the physical properties of outer giant planets , and in particular would allow us to compare them directly to our own solar system giants .
thus , these single transit events should be saved , if at all possible .    .
the shaded region indicates @xmath48 , where we expect the sample to suffer from incompleteness .
the dashed line indicates our extrapolation of the distribution , found by assuming the constant distribution in @xmath38 from @xcite .
( _ middle _ ) the total transit probability .
the solid line shows the probability that a planet transits _ and _ exhibits a single transit during the _ kepler _ mission lifetime .
the dotted line shows the same probability for two transits .
( _ bottom _ ) the expected number of one and two - transit events .
the probability distribution has been convolved with the observed fraction of planets and multiplied by the number of stars kepler will observe ( 100,000 ) to give the number of single transit systems ( solid line ) and two transit systems ( dotted line ) that can be expected .
the dashed line shows the convolution of the single transit probability distribution with the extrapolation to larger semimajor axes .
[ fig : np ] ]
with a few basic assumptions , _ kepler _ photometry of a transiting planet can provide a period for the system that can be combined with precise radial velocity measurements to estimate a mass for the planet . here
we briefly sketch the basics of how these properties can be estimated and what their uncertainties are , and we provide details in the following sections .
we assume the planet is much less massive than the star , @xmath49 , and circular orbits .
we consider the impact of eccentricity on our conclusions in  [ sec : e ] .    assuming no limb - darkening and assuming the out - of - transit flux is known perfectly
, a transit can be characterized by four observables , namely the fractional depth of the transit @xmath50 , the full - width half - maximum of the transit @xmath51 , the ingress / egress duration @xmath52 , and the time of the center of transit @xmath53 .
these can be combined to estimate the instantaneous velocity of the planet at the time of transit , @xmath54 for a circular orbit , @xmath55 can be used to estimate the period of the planet @xmath56 where we have employed kepler s third law assuming that the planet s mass is much smaller than the star s mass .
thus , if the stellar density @xmath57 is known , then the planet s period can be estimated from photometry of a single transit @xcite .
note that there is a degeneracy between @xmath57 and @xmath21 , so with a single transit @xmath57 must be determined by other means in order to derive an estimate of @xmath21 .
@xmath57 can be estimated via spectroscopy combined with theoretical isochrones , or asteroseismology .
the total uncertainty in the period is the quadrature sum of the contribution from the estimate of @xmath57 and the contribution from the _ kepler _ photometry ( i.e. , the uncertainty at fixed @xmath57 ) , @xmath58 the contribution from _ kepler _ is dominated by the uncertainty in the ingress / egress time @xmath52 , @xmath59 where we have assumed @xmath60 , and that the number of points taken out of transit is much larger than the number taken during the transit ( and thus the out - of - transit flux is known essentially perfectly ) . here
@xmath61 is approximately equal to the total signal - to - noise ratio of the transit , @xmath62 where @xmath63 is the photon collection rate . for _ kepler
_ we assume , @xmath64 and thus _ kepler _ will detect a single transit with a signal - to - noise ratio of @xmath65 where @xmath66 is the radius of the planet . for a neptune - size planet with these fiducial parameters , @xmath67 .
_ kepler s _ contribution to the uncertainty in the period is , @xmath68 where we have assumed the impact parameter @xmath69 .
due to the strong scaling with @xmath66 and relatively weak scaling with @xmath21 , there are essentially two regimes . for @xmath70 , the uncertainty in @xmath21
is dominated by uncertainties in the estimate of @xmath57 , which is expected to be @xmath71 , whereas for @xmath72 , it is very difficult to estimate the period due to the uncertainty in @xmath52 from the _ kepler _ photometry .
once the period is known from the photometry , the mass of the planet , @xmath35 , can be measured with radial velocity observations . for a circular orbit , and for observations spread out over a time that is short compared to @xmath21
, the radial velocity of the star can be expanded about the time of transit , and so approximated by the velocity at the time of transit @xmath73 , plus a constant acceleration @xmath74 , @xmath75 where @xmath76 , and @xmath77 is the stellar radial velocity semi - amplitude .
thus , @xmath78 where @xmath79 is the surface gravity of the star , and we have assumed @xmath80 .
the uncertainty in the planet mass has contributions from three distinct sources : the uncertainty in @xmath79 , which can be estimated from spectroscopy , the uncertainty in @xmath74 , which is derived from rv observations after the transit , and the uncertainties in @xmath81 and @xmath50 , which are derived from _ kepler _ photometry .
in fact , the uncertainty in @xmath52 dominates over the uncertainties in @xmath51 and @xmath50 .
therefore , we may write , @xmath82 for reasonable assumptions and long - period planets ( @xmath83 ) , we find that the uncertainty in @xmath74 dominates over the uncertainty in @xmath52 .    assuming that @xmath84 equally - spaced rv measurements with precision @xmath85 are taken over a time period @xmath86 after the transit , the uncertainty in @xmath74 is , @xmath87 where we have assumed @xmath88 .
thus , radial velocity observations combined with _
kepler _ photometry can confirm the planetary nature of a jupiter - sized planet in a relatively short time span .
an accurate ( @xmath89 ) measurement of the mass for a jupiter - mass planet , or even a rough characterization of the mass for a neptune - mass planet , will require either more measurements , or measurements with substantially higher rv precision .
of course , additional radial velocity observations over a time span comparable to @xmath21 will further constrain the mass and period .    in @xmath10[sec : lc ] and @xmath10[sec : rv ] , we derive the above expressions for the uncertainty in the mass and the period of the planet using a fisher information analysis .
in the absence of significant limb - darkening , a transit light curve can be approximated by a trapezoid that is described by the five parameters @xmath53 , @xmath51 , @xmath52 , @xmath50 , and @xmath90 , where @xmath90 is the out - of - transit flux .
figure [ fig : lc ] shows this simple trapezoidal model and labels the relevant parameters .
mathematically , the flux @xmath91 as a function of time @xmath92 is given by , @xmath93)/\tau , &   t_1 \le t \le t_2 , \\
f_0-\delta , &   t_2 < t < t_3 , \\     f_0-\delta(1-[(t -[t_c - t/2-\tau/2])/\tau ] ) , & t_3 \le t \le t_4 , \\    f_0 , & t_4 < t ,    \end{array}\right.\ ] ] where @xmath94
 @xmath95 are the points of contact .
we also define @xmath96 to be the total duration of the observations .
this model is fully differentiable and can be used with the fisher matrix formalism to derive exact expressions for the uncertainties in the parameters @xmath53 , @xmath51 , @xmath52 , @xmath50 , and @xmath90 .
we note that the definition for @xmath50 we use here differs slightly from the definition we adopted in
[ sec : basic ] . here
@xmath50 is the depth of the transit ( e.g. , in units of flux ) , rather than the fractional depth .
these two definitions differ by a factor of @xmath90 such that @xmath97 .
however , when @xmath98 ( as will be the case for _ kepler _ ) , the uncertainty in @xmath90 is negligible , and so one may define @xmath99 , thus making these two parametrizations equivalent .
figure [ fig : models ] compares this simple model with the exact model computed using the formalism of @xcite , in this case for a transit of the sun by neptune with a 3.5-year period and an impact parameter of @xmath100 .
we see that the trapezoidal model provides an excellent approximation to the exact light curve , and as we will show the analytic parameter uncertainties are quite accurate .
we also show the case of significant limb - darkening as expected for a sun - like star observed in the _ r_-band ( similar to the _ kepler _ bandpass ) . in this case , the match is considerably poorer , but nevertheless we will show the analytic parameter uncertainties we derive assuming the simple trapezoidal model still provide useful estimates ( see @xcite for a more thorough discussion ) .
we derive exact expressions for the uncertainties in @xmath53 , @xmath51 , @xmath52 , @xmath50 , and @xmath90 by applying the fisher matrix formalism to the simple light curve model ( a detailed explanation of the fisher matrix is given by @xcite ) .
the full details of the derivation are given in the appendix ( @xmath10[sec : append ] ) .
in the case of the long - period transiting planets that will be observed by _ kepler _ , we can simplify the full expression by making a couple of assumptions .
we assume that the flux out of transit , @xmath90 , is known to infinite precision and that @xmath101 .
the uncertainties in the remaining parameters are , @xmath102    we tested the accuracy of these expressions using monte carlo simulations .
we generate 1000 light curves for a given set of orbital parameters .
we assume that the planet orbits a solar - type star in a circular orbit with an impact parameter of 0.2 .
we repeat the analysis for four of the solar - system planets : jupiter , saturn , neptune , and earth .
we assume the expected photometric precision for _ kepler _ for a stellar apparent magnitude of @xmath103 ( a representative magnitude for _ kepler _ s stellar sample ) , a sampling rate of one per @xmath104 , and a total duration for the observations of @xmath105 hours .
then we fit for the parameters @xmath53 , @xmath51 , @xmath52 , @xmath50 , and @xmath90 using a down - hill simplex method . the uncertainty in each parameter
is taken to be the standard deviation in the distribution of the fits for that parameter .
a comparison of the analytic expressions and the monte carlo simulations is shown in fig .
[ fig : sigobs ] , where we have plotted the uncertainties as a function of period for a planet with radius @xmath106 crossing a star with radius @xmath107 at an impact parameter @xmath100 .
we find that the fisher matrix approximation breaks down as the number of points during ingress / egress becomes small .
we also considered expected uncertainties for the exact uniform source and limb - darkened models for the transit light curve .
@xcite give the full solution for a transit involving two spherical bodies and provide code for calculating the transit light curve .
they also include a limb - darkened solution .
these models may be parametrized by the same observables described in our simplified model . the full and
limb - darkened light curves are plotted in fig .
[ fig : models ] as the dotted and dashed curves , respectively .
we use the limb - darkening coefficients from @xcite closest to observations of the sun ( @xmath108 , @xmath109 , @xmath110 = 0.0 $ ] ) in the r - band .
we apply the fisher matrix method to numerical derivatives of these model light curves to compare the uncertainties for the full and limb - darkened solutions with the analytic uncertainties for the trapezoidal model .
a comparison of the uncertainties from the different models is shown in fig .
[ fig : sigobs ] .
we find that the simplified model is a good approximation for the exact , uniform - source transit model .
while the comparison is less favorable with the limb - darkening model , the uncertainties do not differ by more than a factor of a few .
we also compared how the uncertainties in the transit observables varied with impact parameter @xmath111 for the three models , because varying @xmath111 will affect the relative sizes of @xmath51 and @xmath52 .
these comparisons are shown in fig .
[ fig : impact ] for neptune in a circular orbit around the sun with a 3.5-year period .
as can be seen from the figure , variations in @xmath111 have little impact in the uncertainties in the observables for @xmath112 .    ) , the transit depth ( @xmath50 ) , the ingress / egress time ( @xmath52 ) , the fwhm duration of the transit ( @xmath51 ) , the total duration of observations ( @xmath96 ) , and the time of the transit center ( @xmath53 ) .
[ fig : lc],width=576 ]     and an impact parameter @xmath100 .
typical error bars for kepler assuming @xmath113 and 30-min sampling ( solid ) or 20-min sampling ( dotted ) are indicated .
[ fig : models],width=576 ]    , @xmath51 , @xmath52 , and @xmath50 are plotted versus period for neptune orbiting the sun with an impact parameter of 0.2 .
the different line styles indicate the trapezoidal transit model ( solid line ) , the full solution ( dotted line ) , and a model including the effects of limb - darkening ( dashed line ) .
monte - carlo simulations of the trapezoidal model using 30-minute sampling are shown as crosses .
the gray crosses show simulations for 20-minute sampling .
the dotted line in the bottom panel is not visible because it is nearly identical to the solid line .
the vertical dotted lines mark where the ingress / egress duration @xmath52 is equal to the sampling rate @xmath114 , as well as where @xmath115 for 30-minute sampling.[fig : sigobs],width=528 ]    , @xmath51 , @xmath52 , and @xmath50 are plotted versus impact parameter for neptune orbiting the sun with a period of 3.5 years .
the line styles are the same as for fig .
[ fig : sigobs ] . as in fig .
[ fig : sigobs ] , the dotted line in the bottom panel is not visible because it is nearly identical to the solid line .
[ fig : impact],width=528 ]      assuming the planet is in a circular orbit , and assuming the stellar density , @xmath57 , is estimated from independent information , @xcite demonstrated that the planet period , @xmath21 , and impact parameter , @xmath111 , can derived from a single observed transit , with sufficiently precise photometry .
in fact , provided independent estimates of @xmath20 and @xmath19 are available , then it is also possible to derive the planet radius , @xmath116 , semimajor axis , @xmath17 , and instantaneous velocity at the time of transit , @xmath55 .
the uncertainties in these physical parameters can be calculated through standard error propagation .
equations for these quantities and their uncertainties are given in eqs .
[ eqn : rp][eqn : p ] ( note that we calculate @xmath117 instead of @xmath111 because numerical methods have difficulty when @xmath118 ) . in each case
below , the second approximate uncertainty equation is calculated by substituting the variances and covariances of the observables @xmath119 , and @xmath50 , and simplifying under the assumption that @xmath101 .
@xmath120\nonumber\\ & \simeq & t\delta q^{-2}\left(\frac{6t^2}{\tau^3}\right)\nonumber\\ & & \nonumber\\ a & = & \frac{g_{\star}}{4}\left(\frac{t\tau}{\sqrt{\delta}}\right)\\ \sigma_a^2 & = & a^2\left[\frac{1}{g_{\star}^2}\sigma_{g_{\star}}^2 + \frac{1}{t^2}\sigma_t^2 + \frac{1}{\tau^2}\sigma_{\tau}^2 + \frac{1}{4\delta^2}\sigma_{\delta}^2+\frac{2}{t\tau}\sigma_{t\tau}^2 -\frac{1}{t\delta}\sigma_{t\delta}^2 - \frac{1}{\tau\delta}\sigma_{\tau\delta}^2 \right]\nonumber\\ & \simeq & a^2\left[\frac{1}{g_{\star}^2}\sigma_{g_{\star}}^2+q^{-2}\left(\frac{6t}{\tau}\right)\right]\nonumber\\ & & \nonumber\\ v_{\rm tr , p } & = & 2r_{\star}\left(\frac{\sqrt{\delta}}{t\tau}\right)^{1/2}\\ \sigma_{v_{\rm tr , p}}^2 & = & v_{\rm tr , p}^2\left[\frac{1}{r_{\star}^2}\sigma_{r_{\star}}^2+\frac{1}{4t^2}\sigma_t^2+\frac{1}{4\tau^2}\sigma_{\tau}^2+\frac{1}{16\delta^2}\sigma_{\delta}^2+\frac{1}{2t\tau}\sigma_{t\tau}^2-\frac{1}{4t\delta}\sigma_{t\delta}^2-\frac{1}{4\tau\delta}\sigma_{\tau\delta}^2\right]\nonumber\\ & \simeq & v_{\rm tr , p}^2\left[\frac{1}{r_{\star}^2}\sigma_{r_{\star}}^2 + q^{-2}\left(\frac{3t}{2\tau}\right)\right]\nonumber\\ & & \nonumber\\ \label{eqn : p } p & = & \frac{\pi g}{4}\left(\frac{4}{3}\pi\rho_{\star}\right)\left(\frac{t\tau}{\sqrt{\delta}}\right)^{\frac{3}{2}}\\ \sigma_p^2 & = & p^2\left [ \frac{1}{\rho^2_{\star}}\sigma_{\rho_{\star}}^2 + \frac{9}{4t^2}\sigma_t^2 + \frac{9}{4\tau^2}\sigma_{\tau}^2 + \frac{9}{16\delta^2}\sigma_{\delta}^2+\frac{9}{2t\tau}\sigma_{t\tau}^2 -\frac{9}{4t\delta}\sigma_{t\delta}^2 - \frac{9}{4\tau\delta}\sigma_{\tau\delta}^2\right]\nonumber\\ & \simeq & p^2\left[\frac{1}{\rho^2_{\star}}\sigma_{\rho_{\star}}^2 + q^{-2}\left(\frac{27t}{2\tau}\right ) \right]\nonumber\end{aligned}\ ] ]    of course , the quantities @xmath19 , @xmath79 , and @xmath57 and their uncertainties must be estimated from some other external source of information , such as spectroscopy or theoretical isochrones .
we can generally expect that the fractional uncertainties on these quantities will be of order @xmath121 . for long - period ( @xmath122 ) , jupiter - sized planets , @xmath123 , and
thus the variances and covariances of the transit observables will be small in comparison to the expected uncertainties on @xmath124 , and @xmath79 .
the fractional uncertainties in these derived parameters are plotted in fig .
[ fig : deriv ] .
we compared the analytic uncertainties given in eqs .
[ eqn : rp][eqn : p ] to the monte carlo simulations described in sec .
[ sec : models ] .
the fractional uncertainty in each parameter from the monte carlo simulations is one - half the range of the middle 68% of the data divided by the median .
notice that as @xmath52 decreases from @xmath125 to @xmath114 , the simulations become increasingly disparate from the theoretical expectations .
@xmath126 does not show this behavior because it does not depend on @xmath52 .
in this regime , as @xmath52 decreases , it becomes increasingly probable that only one point will be taken during ingress leaving @xmath52 relatively unconstrained and increasing its uncertainty and the uncertainty of quantities that depend on it .
we now consider in detail the uncertainty in the estimated period . given that @xmath127 , we have that @xmath128 . since @xmath129 and @xmath130 , the contribution to the uncertainty in the period due to the _ kepler _ photometry is a much stronger function of the planet s radius than its period , @xmath131 . as a result , for a fixed stellar radius @xmath19 , the ability to accurately estimate @xmath21 depends almost entirely on @xmath66 . furthermore , because the scaling with @xmath66 is so strong , there are essentially two distinct regimes : for large @xmath66 , the uncertainty in @xmath21 is dominated by the uncertainty in @xmath57 , whereas for small @xmath66 , the period can not be estimated .
the boundary between these two regimes is where the contribution to the uncertainty in the period from @xmath57 is equal to the contribution from the _ kepler _ photometry , i.e.  where the two terms in the brackets in eq .
[ eqn : p ] are equal .
this occurs at a radius of , @xmath132^{1/5},\ ] ] assuming @xmath69 .
note that @xmath133 is a weak or extremely weak function of all of the parameters except @xmath19 .
[ fig : contour ] illustrates this point by showing contours of constant uncertainty in the period as a function of the radius and period of the planet .
the calculations are for systems with @xmath134 , @xmath100 , and @xmath113 .
this figure shows that for planets larger than the radius of neptune , the uncertainty in the period will be dominated by the uncertainty in the inferred stellar density , whereas for planets much smaller than neptune , an accurate estimate of the period from the _ kepler _ light curve will be impossible .    as shown above , the uncertainty in @xmath21 is dominated by the uncertainty in @xmath57 and @xmath52 . our ability to determine @xmath52 , and its uncertainty , depends both on its length and how many points we have during the ingress or egress .
the length of @xmath52 depends on @xmath55 ( which is a proxy for @xmath21 in the circular case ) , @xmath111 , and @xmath66 .
these properties are intrinsic to the system , but may be derived from the observables .
we can explore how the contours shown in fig .
[ fig : contour ] vary with impact parameter .
figure [ fig : impact ] shows that a system with a larger impact parameter will have a smaller fractional uncertainty in @xmath52 .
thus , a system with a larger impact parameter would have a smaller fractional uncertainty in @xmath21 at fixed period and planet radius .    the other effect that influences
the contours is the sampling of @xmath52 .
as the number of points taken during the ingress / egress becomes small , the fractional uncertainty in @xmath52 , and hence in @xmath21 , increases .
in particular , if only one point is taken during the ingress , then the duration is relatively unconstrained .
the probability of taking only one point during the ingress increases linearly from 0 to 1 as @xmath52 decreases from @xmath135 to @xmath114 . in fig .
[ fig : deriv ] , the monte carlo simulations for the fractional uncertainty in @xmath21 diverge from the simple model over this range .
we indicate this region by the shaded portion of fig .
[ fig : contour ] .
where @xmath136 , the fractional uncertainty in @xmath21 can be several times that predicted by the theoretical calculations , but these uncertainties converge as the shading gets lighter towards @xmath115 . below the line where @xmath136 , the uncertainty in @xmath52 increases rapidly , but fortunately , in much of this regime , _ kepler _ will observe multiple transits , and this analysis will be unnecessary . as shown in the figure , a faster sampling rate , such as 1 per 20 min , significantly expands the parameter space over which our theoretical uncertainties are valid .
, @xmath17 , @xmath21 , @xmath55 , and @xmath117 are plotted versus period for neptune orbiting the sun with an impact parameter of 0.2 .
the line styles are the same as for fig .
[ fig : sigobs ] . the vertical dotted lines mark where the ingress / egress duration @xmath52 is equal to the sampling rate @xmath114 , as well as where @xmath115 for 30-minute sampling .
[ fig : deriv],width=528 ]     changes with @xmath137 magnitude ; it shows the 0.20 contour for @xmath138 and @xmath100 .
the diagonal black dotted line represents the boundary below which the assumptions of fisher matrix break down ( the ingress / egress time @xmath52 is roughly equal to the sampling rate of 1/30 min ) .
thus , the contours below this boundary are shown as dash - dotted lines .
the region between @xmath115 and @xmath136 is shaded to indicate the increasingly probability of obtaining only one point during ingress or egress leading to an uncertainty in @xmath21 that is larger than theoretical expectations .
the diagonal gray dotted line indicates the @xmath136 boundary for 20-min sampling .
the vertical dotted line shows the mission lifetime of _ kepler _ ( @xmath139 ) .
solar system planets are indicated .
[ fig : contour],width=576 ]
the mass of the planet comes from sampling the stellar radial velocity curve soon after the transit is observed . near the time of transit we can expand the stellar radial velocity , @xmath140\approx v_0 - a_{\star}(t - t_c ) , \ ] ] where @xmath73 is the systemic velocity
, @xmath141 is the stellar radial velocity semi - amplitude , and @xmath142 is the stellar acceleration . because the planet is known to transit and has a long period , we assume @xmath143 . a fisher matrix analysis of the linear form of eq .
[ eqn : rv ] gives the estimated uncertainty in @xmath74 to be @xmath144 where @xmath85 is the radial velocity precision , @xmath145 is the total time span of the radial velocity observations , and @xmath84 is the number of observations , and we have assumed @xmath146 and that the observations are evenly spaced in time .
the details of this derivation are given in the appendix ( @xmath10[sec : append ] ) .
equations for the mass of the planet and the uncertainty are @xmath147\nonumber\\ & \simeq & m_{\rm p}^2\left[\frac{4}{g_{\star}}\sigma_{g_{\star}}^2+\frac{1}{a_{\star}^2}\sigma_{a_{\star}}^2+q^{-2}\left(\frac{24t}{\tau}\right)\right ] .
\nonumber\end{aligned}\ ] ] the approximate expression for the uncertainty in the mass is taken in the limit @xmath101 .
furthermore , the @xmath148 term can be neglected in most cases because , as we showed for the corresponding term in uncertainty in the period , it will be very small for planets with radii larger than that of neptune .
the scaling for the @xmath149 term is given in eq .
[ eqn : sigbscale ] .
contours of constant fractional uncertainty in the mass of the planet are shown in fig .
[ fig : mpplot ] as a function of @xmath21 and @xmath35 , assuming @xmath150 radial velocity measurements are taken over @xmath151 mos with a precision @xmath152 .
it shows that the planet mass can be estimated to within a factor of two over this time period , thus establishing the planetary nature of the transiting object . by doubling the length of observations
, one can put stronger constraints on the mass of the planet .
there are two points to bear in mind when applying this estimate for the uncertainty in the mass of the planet .
first , for @xmath153 , the @xmath154 term is no longer small .
the second consideration is that as time progresses away from the time of transit , the straight line approximation to the radial velocity curve will break down . in that case
, the period will begin to be constrained by the radial velocity curve itself , and the uncertainties should be calculated from a fisher matrix analysis of the full expression for the radial velocity curve with three parameters : @xmath73 , @xmath141 , and @xmath21 .     as a function of @xmath35 and @xmath21 .
the solid lines show the result for 20 radial velocity measurements with precision of @xmath155 taken over a period of 3 months after the transit .
the dashed line shows the contour for @xmath156 for 40 observations taken over 6 months .
the dotted line indicates the mission lifetime of 3.5 years
. the positions of jupiter and saturn are indicated .
[ fig : mpplot],width=576 ]
the results presented in the previous sections assumed circular orbits . given that the average eccentricity @xmath157 of planets with @xmath158 is @xmath159 , this is not necessarily a good assumption . in this section ,
we assess the effects of non - zero eccentricities on the ability to characterize single - transit events detected by _
kepler_. we note that our discussion has some commonality with the study of @xcite , who discuss the possibility of characterizing the orbital eccentricities of transiting planets with photometric observations .
however , our study addresses this topic from a very different perspective .    in  [ sec :
basic ] we demonstrated that , under the assumption that @xmath160 , the planet period @xmath21 and mass @xmath35 can be inferred from two observables : the velocity of the planet at the time of transit @xmath161 , and the projected acceleration of the star @xmath74 . in the case of a non - zero eccentricity , there are two additional unknown parameters : the eccentricity @xmath157 and the argument of periastron of the planet @xmath162 .
thus , with only two observables ( @xmath161 and @xmath74 ) and four unknowns ( @xmath163 ) , it is not possible to obtain a unique solution .
although a unique solution to the planet parameters does not exist using the observable parameters alone , it may nevertheless be possible to obtain an interesting constraint on @xmath21 and/or @xmath35 by adopting reasonable priors on @xmath157 and @xmath162 .
from the transit observables , one can estimate the velocity of the planet at the time of transit ( eq . [ eqn : vtr ] ) . in the case of @xmath164 ,
this is given by , @xmath165 where here and throughout this section we assume @xmath166 .
solving for @xmath21 , @xmath167 ^ 3.\ ] ] for a fixed value of @xmath157 , the inferred period relative to the assumption of a circular orbit has extremes due to the unknown value of @xmath162 of @xmath168 taking a typical value for the eccentricity of @xmath169 , this gives a range of inferred periods of @xmath170 relative to the assumption of a circular orbit .
in fact , because the probability that a planet with a given @xmath17 transits its parent star depends on @xmath157 and @xmath162 @xcite , a proper bayesian estimate , accounting for these selection effects , might reduce the range of inferred values for @xmath21 significantly .    what can be learned from rv observations immediately after the transit ?
expanding the stellar projected velocity around the time of transit , we can write , @xmath171 the projected acceleration of the star at the time of transit is , @xmath172 this can be combined with @xmath55 to derive the mass of the planet , @xmath173 thus , for a fixed eccentricity , the inferred mass relative to the assumption of a circular orbit has extremes of , @xmath174 which for @xmath169 yields a range of @xmath175 , which is generally _ smaller _ than the contribution to the uncertainty in @xmath35 due to the measurement uncertainty in @xmath74 for our fiducial case of @xmath176 and @xmath177 ( see eq .
[ eqn : sigbscale ] ) .
we may also consider what can be learned if it is possible measure the curvature of the stellar radial velocity variations immediately after transit .
the projected stellar jerk is given by , @xmath178 this can be combined with @xmath161 and @xmath74 to provide an independent constraint on a combination of the eccentricity and argument of periastron , @xmath179 unfortunately , it will be quite difficult to measure @xmath180 . using the same fisher formalism as we used to estimate the uncertainty in @xmath74 (  [ sec : append ] ) , assuming @xmath84 evenly - spaced rv measurements with precision @xmath85 , taken over time span @xmath86 after the transit , the uncertainty in @xmath180 is , @xmath181 where we have assumed @xmath88 . for our fiducial parameters ,
the fractional uncertainty in @xmath180 is , @xmath182 where we have approximated @xmath183 .
we conclude that , in order to obtain interesting constraints on the eccentricity of long - period planets detected by _ kepler _ , higher - precision rv measurements taken over a baseline comparable to @xmath184 will be needed .
the discovery of long - period transiting planets along with subsequent follow - up observations would greatly enhance our understanding of the physics of planetary atmospheres and interiors .
such planets would allow us to gather constraints in a regime of parameter space currently only occupied by our own giant planets , namely planets whose energy budgets are dominated by their residual internal heat , rather than by stellar insolation .
these constraints , in turn , might provide new insights into planet formation . in this paper , we demonstrated that it will be possible to detect and characterize such long - period planets using observations of single transits by the _ kepler _ satellite , combined with precise radial velocity measurements taken immediately after the transit .
indeed , these results can be generalized to any transiting planet survey using the scaling relations we provide , and it may be particularly interesting to apply them to the _ convection , rotation & planetary transits ( corot ) _
mission @xcite .
we calculated that _ kepler _
will see a few long - period , single transit events and showed that , for circular orbits , the period of the system can be derived from the _ kepler _ light curve of a single event .
we derived an expression for the uncertainty in this period and showed that it is dominated by the uncertainty in the stellar density derived from spectroscopy ( which we assume is @xmath185 ) for planets with radii larger than the radius of neptune , rather than being dependent on the properties of the transit itself .
this method can also be applied to planets that _
kepler _ will observe more than once , so that the second time a planet is expected to transit _ kepler _ can make a selective improvement to its time sampling to better characterize the transit .
we have also shown how the mass of the planet can be constrained by acquiring precise ( @xmath186 ) radial velocity measurements beginning shortly after the transit occurs .
we have shown that 20 measurements over 3 months can measure the mass of a jupiter - sized object to within a factor of a few and that extending those observations to 40 measurements over 6 months significantly reduces the uncertainty in the mass . knowing the mass to within a factor of a few in such a short time can distinguish between brown dwarfs and planets rapidly and allow us to maximize our use of radial velocity resources .
we explored the effect of eccentricity on the ability to estimate the planet mass and period .
allowing for a non - zero eccentricity adds two additional parameters , and as a result it is not possible to obtain a unique solution for the planet mass , period , eccentricity , and argument of periastron .
however , by adopting a reasonable prior on the eccentricity of the planet , the period and mass of the planet can still be estimated to within a factor of a few .
detailed characterization of the planet properties will require precision rv measurements obtained over a duration comparable to the period of the planet .
thus , in the interest of `` getting the most for your money , '' we have shown that the sensitivity of _ kepler _ extends to planets with periods beyond its nominal mission lifetime . with the launch of the _ kepler _ satellite ,
we are poised to discover and characterize several long - period transiting systems , provided that we are prepared to look for them .
we are grateful to josh carter , jason eastman , eric ford , andy gould , matt holman , yoram lithwick , josh winn , and andrew youdin for helpful discussions .
we thank the referee ron gilliland for constructive comments .
jcy is supported by a dean s graduate enrichment fellowship from the ohio state university .
19 natexlab#1#1    , a. 2003 , advances in space research , 31 , 345    , j.  w. 2007 , , 119 , 986    , g. , borucki , w.  j. , & koch , d. 2005 , new astronomy review , 49 , 478    , j .- p . , bennett , d.  p. , fouqu , p. , williams , a. , dominik , m. , jrgensen , u.  g. , kubas , d. , cassan , a. , coutures , c. , greenhill , j. , hill , k. , menzies , j. , sackett , p.  d. , albrow , m. , brillant , s. , caldwell , j.  a.  r. , calitz , j.  j. , cook , k.  h. , corrales , e. , desort , m. , dieters , s. , dominis , d. , donatowicz , j. , hoffman , m. , kane , s. , marquette , j .- b . , martin , r. , meintjes , p. , pollard , k. , sahu , k. , vinter , c. , wambsganss , j. , woller , k. , horne , k. , steele , i. , bramich , d.  m. , burgdorf , m. , snodgrass , c. , bode , m. , udalski , a. , szymaski , m.  k. , kubiak , m. , wic eckowski , t. , pietrzyski , g. , soszyski , i. , szewczyk , o. , wyrzykowski ,  . , paczyski , b. , abe , f. , bond , i.  a. , britton , t.  r. , gilmore , a.  c. , hearnshaw , j.  b. , itow , y. , kamiya , k. , kilmartin , p.  m. , korpela , a.  v. , masuda , k. , matsubara , y. , motomura , m. , muraki , y. , nakamura , s. , okada , c. , ohnishi , k. , rattenbury , n.  j. , sako , t. , sato , s. , sasaki , m. , sekiguchi , t. , sullivan , d.  j. , tristram , p.  j. , yock , p.  c.  m. , & yoshioka , t. 2006 , , 439 , 437    , w. , koch , d. , boss , a. , dunham , e. , dupree , a. , geary , j. , gilliland , r. , howell , s. , jenkins , j. , kondo , y. , latham , d. , lissauer , j. , & reitsema , h. 2004 , in esa special publication , vol .
538 , stellar structure and habitable planet finding , ed .
f.  favata , s.  aigrain , & a.  wilson , 177182    , c.  j. 2008 , arxiv:0801.2579 , 801    , r.  p. , wright , j.  t. , marcy , g.  w. , fischer , d.  a. , vogt , s.  s. , tinney , c.  g. , jones , h.  r.  a. , carter , b.  d. , johnson , j.  a. , mccarthy , c. , & penny , a.  j. 2006 , , 646 , 505    , j.  a. , yee , j.  c. , eastman , j. , gaudi , b.  s. , & winn , j.  n. 2008 , arxiv:0805.0238 , 805    , d. , allen , l.  e. , megeath , s.  t. , torres , g. , alonso , r. , brown , t.  m. , gilliland , r.  l. , latham , d.  w. , mandushev , g. , odonovan , f.  t. , & sozzetti , a. 2005 , , 626 , 523    , d. , brown , t.  m. , noyes , r.  w. , & gilliland , r.  l. 2002 , , 568 , 377    , a. 2000 , , 363 , 1081    , a. , butler , r.  p. , marcy , g.  w. , vogt , s.  s. , wright , j.  t. , & fischer , d.  a. 2008 , arxiv:0803.3357 , 803    , d. , seager , s. , richardson , l.  j. , & harrington , j. 2005 , , 434 , 740    , e.  b. , quinn , s.  n. , & veras , d. 2008 , arxiv:0801.2591 , 801    , a. 2003 , arxiv : astro - ph/0310577    , a. , udalski , a. , an , d. , bennett , d.  p. , zhou , a .- y . , dong , s. , rattenbury , n.  j. , gaudi , b.  s. , yock , p.  c.  m. , bond , i.  a. , christie , g.  w. , horne , k. , anderson , j. , stanek , k.  z. , depoy , d.  l. , han , c. , mccormick , j. , park , b .-
, pogge , r.  w. , poindexter , s.  d. , soszyski , i. , szymaski , m.  k. , kubiak , m. , pietrzyski , g. , szewczyk , o. , wyrzykowski ,  . ,
ulaczyk , k. , paczyski , b. , bramich , d.  m. , snodgrass , c. , steele , i.  a. , burgdorf , m.  j. , bode , m.  f. , botzler , c.  s. , mao , s. , & swaving , s.  c. 2006 , , 644 , l37    , k. & agol , e. 2002 , , 580 , l171    , j.  f. , matthews , j.  m. , seager , s. , miller - ricci , e. , sasselov , d. , kuschnig , r. , guenther , d.  b. , moffat , a.  f.  j. , rucinski , s.  m. , walker , g.  a.  h. , & weiss , w.  w. 2007 , arxiv:0711.4111 , 711    , s. & malln - ornelas , g. 2003 , , 585 , 1038
the fisher matrix formalism is a simple way to estimate the uncertainties in the parameters , @xmath187 , of a model @xmath188 , that is being fit to a series of measurements @xmath189 with measurement uncertainties @xmath190 .
the covariance of @xmath191 with @xmath192 is given by the element @xmath193 of the covariance matrix @xmath194 , where @xmath195 and the entries of @xmath111 are given by , @xmath196 where @xmath197 is the number of data points .
in the limit of infinite sampling ( @xmath198 ) and fixed precision , @xmath199 , @xmath200 where the interval of interest is given by @xmath201 $ ] , and in this case @xmath96 is the total duration of observations .
thus , if the partial derivatives of a model with respect to its parameters are known , then the uncertainties in those parameters can be estimated @xcite .    for the simplified transit model described in @xmath10[sec : models ] ,
the observable parameters are @xmath202 $ ] . for a sampling rate @xmath203 ,
the @xmath111 matrix is , @xmath204 the covariance matrix is , @xmath205 where @xmath206 is the duration of the flat part of the eclipse and @xmath207 is the time spent out of eclipse .
thus , the uncertainty on the ingress / egress time , @xmath208 , is given by @xmath209    define @xmath210 .
then , in the limit @xmath60 , the uncertainties in the observable parameters are given by , @xmath211 note that @xmath61 is approximately the total signal - to - noise ratio of the transit . assuming that the photometric uncertainties are limited by photon noise , we have that @xmath212 , where @xmath63 is the photon collection rate .
this recovers the expression for @xmath61 in  [ sec : basic ] .
consider @xmath84 measurements of @xmath216 , each with precision @xmath85 , taken at times @xmath217 .
fitting these measurements to the linear model in eq .
[ eqn : vexpand ] , we can estimate the uncertainties in the parameters @xmath73 and @xmath74 using the fisher matrix formalism , @xmath218 .\ ] ] the covariance matrix is the inverse of @xmath111 , @xmath219-\left[\displaystyle\sum_{k=1}^{n}(t_k - t_c)\right]^2}\left[\begin{array}{cc } \displaystyle\sum_{k=1}^{n}(t_k - t_c)^2 & \displaystyle\sum_{k=1}^{n}(t_k - t_c)\\ \displaystyle\sum_{k=1}^{n}(t_k - t_c ) & n \end{array}\right ] .\ ] ] if the points are evenly spaced by @xmath220 @xmath221 and thus , @xmath222 .\end{aligned}\ ] ] the uncertainty in the projected stellar acceleration is , @xmath223    in the limit as @xmath224 , the covariance matrix reduces to , @xmath225 .\end{aligned}\ ] ] defining the total length of observations as @xmath226 , we can also write , @xmath227,\end{aligned}\ ] ] and so the uncertainty @xmath74 in the limit of @xmath224 becomes , @xmath228


measurements of the masses and radii of neutron stars provide unique information on the equation of state of matter at densities above the nuclear saturation density .
extremely accurate mass measurements in binary systems can be obtained through dynamical methods ( see , e.g. , stairs 2006 ; demorest et al .
however , without independent determinations of the radii of these neutron stars , these measurements provide only limits for the equation of state .
in contrast , an accurate measurement of the radius of even a single neutron star can provide strong constrains on the equation of state of neutron star matter ( lattimer & prakash 2001 ; zel & psaltis 2009 ; read et al . 2009 ) .
thermonuclear x - ray bursts observed from low mass x - ray binaries have proven to be one of the most suitable events that can be used to measure both the radii and the masses of neutron stars ( see , e.g. , van paradijs 1978 , 1979 ; damen et al . 1990 ; lewin , van paradijs , & taam 1993 ; zel 2006 ) . using high signal - to - noise x - ray data and
time - resolved spectroscopy allows a measurements of the apparent radius during the cooling segments of the bursts .
in addition , time - resolved spectroscopy during a photospheric radius expansion burst enables a measurement of the eddington luminosity for each source . combining these two measurements leads to an uncorrelated determination of the neutron star mass and radius .
thanks to the rossi x - ray timing explorer ( rxte ) , which provided more than 15 years of observations of bursts ( galloway et al .
2008 ) , it has now become possible to not only make neutron star mass and radius measurements using observations of specific systems ( see , e.g. , zel , gver , & psaltis 2009 ) but also to quantify the systematic uncertainties in these measurements . in paper
i of these series ( gver , psaltis , & zel 2012a ) , we analyzed more than 13,000 x - ray spectra from 12 different sources and showed that the errors in the spectroscopic determination of neutron - star radii due to systematic effects in the cooling tails of x - ray bursts are in the range of @xmath23%@xmath08% . in paper
ii ( gver , zel , & psaltis 2012b ) , we used all the archival photospheric radius expansion bursts and found that for sources from which a large number of bursts are observed , systematic uncertainties in the determination of the eddington flux are at the 5%@xmath010% level .    in this paper , we focus on a third possible source of systematic uncertainty , which is the absolute flux calibration of the proportional counter array ( pca ) on board rxte .
absolute flux uncertainties affect the normalization of the spectra obtained during bursts and therefore the measurement of both the apparent radii and of the eddington flux .
determining the absolute flux calibration of an x - ray detector requires a major effort from any instrument team and all x - ray observatories put large amount of resources into in - flight calibrations of their detectors .
however , these efforts are hindered by the fact that almost all x - ray sources in the sky that may be used as standard candles often show unpredicted variability over both long and short time scales ( see , e.g. , wilson - hodge et al .
2011 , and references therein for the recent activity of the crab nebula ) .
additionally , complications due to significantly different responses of different detectors , both in terms of effective areas and the energy ranges that the detectors are sensitive to , become a major hurdle and restrict the inter - calibration efforts .
furthermore , the often unknown intrinsic spectral shape of the astronomical objects being used for the calibration makes the problem even more complex .
simultaneous observations of sources with well known intrinsic spectra by two or more x - ray detectors that have different energy responses allow us to better understand and quantify the systematic uncertainties in the absolute flux calibrations of individual detectors .    in the early days of high energy astrophysics ,
the crab nebula was seen as the ideal stellar source for the calibration of various instruments . compiling a number of rocket , balloon , and satellite observations performed between 1964 and 1972 ,
toor & seward ( 1974 ) concluded that the diffuse x - ray emission should be steady at least for the time period they analyzed , since any variation in the pulsar emission larger than 10% at a few kev would have been detected . with these limits
, they assumed that the crab is a steady source and provided a `` universal '' power - law amplitude and photon - index with their relevant uncertainties . in the following years
, these values were often cited and used as a reference point to calibrate the performance of various instruments .
norgaard - nielsen et al .
( 1994 ) even proposed a specific mission called _ excalibur _ to continuously monitor the crab nebula using two pinhole cameras one with a ccd and another with a proportional counter to measure the spectrum of the nebula with an accuracy @xmath32% over the energy range 0.8@xmath020  kev .
such a mission could , in principle , be a reliable inter - calibration reference for subsequent detectors .    using the crab as a standard source , kirsch et al .
( 2005 ) compared the xmm - newton / epic - pn measurements of the spectrum of crab with those from all the other modern x - ray satellites .
in particular , they found that although the photon index inferred from xmm - newton / epic agrees within errors with the toor & seward ( 1974 ) value , the normalization derived from epic - pn is significantly lower than that from the rest of the instruments .
they also noted that , with the new generation of instruments , crab is becoming a source that is too bright and too extended to be used as a standard calibration source . in a similar , later study , weisskopf et al .
( 2010 ) investigated possible departures from a single power - law of the intrinsic spectrum of crab for rosat / pspc , rxte / pca , and the xmm - newton / epic - pn , which would affect the instrument calibrations . using the high statistical quality of the rxte / pca data
, these authors concluded that they would have been able to detect significant deviations from a simple power - law spectrum but they did not find any evidence for such a deviation .    with the advancement of imaging detectors and considering the brightness of the crab for the onboard ccds , searches for new calibration sources and techniques started .
in order to coordinate and improve the efforts on inter - calibration , a group of instrument teams formed the international astronomical consortium for high - energy calibration ( iachec , sembay et al . 2010 ) .
one of the recent results of this group demonstrates that in the 2.0@xmath07.0  kev range , epic - pn measures systematically lower fluxes than the chandra acis , by 11.0@xmath10.5% ( nevalainen et al .
2010 ) , using a sample of 11 nearby clusters of galaxies .
the difference is smaller yet still significant when acis is compared with epic mos1 and mos2 , with differences at 3.0@xmath10.5% and 6.0@xmath10.5% , respectively ( nevalainen et al .
2010 ) .
the spectroscopic measurements of neutron star masses and radii have been predominantly performed with rxte .
launched in 1995 and operated until 2012 , rxte / pca has always been calibrated against the crab nebula ( jahoda et al .
2006 ; shaposnikov et al .
2012 ) together with an onboard calibration source used mainly for channel - to - energy conversion . throughout the mission , any possible variations in the crab , as the one reported by wilson - hodge et al .
( 2011 ) , were assumed to be averaged out ( shaposnikov et al .
2012 ) ; this issue will be discussed further in section  [ variable_source ] . because the pca was intrinsically calibrated against the crab and because it is a non - imaging instrument , inter - calibration of the pca has always been a challenge , especially as more and more sources such as clusters of galaxies or similarly soft and relatively dim objects ( i.e. , rx  j1856.5@xmath03754 ) started to be used as standard candles .
time resolved x - ray spectra extracted during thermonuclear x - ray bursts observed from low mass x - ray binaries offer a different opportunity in determining the systematic uncertainties caused by calibrations of different detectors . since the early studies of x - ray bursts , it has been shown that the burst x - ray spectra can statistically be modeled with blackbody functions ( see , e.g. , swank et al .
1977 ; lewin , van paradijs , & taam 1993 ; galloway et al .  2008 ) . in paper
i , we analyzed more than 13,000 x - ray spectra extracted from the cooling tails of x - ray bursts observed from 12 different sources and showed that a planckian function can statistically fit the vast majority of the observations .
the high count rates of rxte and the large number of bursts allowed us to show that any deviation from a planckian spectrum is limited to @xmath4 in the 2.5@xmath025.0  kev range .
this almost unique characteristic of burst spectra make them ideal for testing our ability to convert countrate spectra to bolometric fluxes .
furthermore , because the obtained spectral parameters are directly related to neutron star mass - radius measurements , results obtained from an analysis of x - ray burst data can be directly converted to uncertainties in the measurements of neutron star masses and radii .    using simultaneous rxte , chandra x - ray observatory ( chandra ) , and xmm-_newton _ observations of x - ray bursts from gs  1826@xmath0238
, we quantify the differences in the measurements of the apparent surface area of the thermal emission when using different detectors with independent calibrations .
we describe our method in section  2 , which is based on predicting the observed chandra and xmm-_newton _ countrates using time resolved spectral parameters obtained from rxte / pca .
section  3 details the analysis of the data from different x - ray satellites .
finally , in sections  4 and 5 we present our results and discuss their implications .
the rapid evolution of the effective temperature during the cooling tails of thermonuclear x - ray bursts demands very short exposure time observations for direct comparison of spectral data that are obtained simultaneously .
the short exposures require very high countrates and , hence , detectors with large photon collecting areas .
this has been a unique property of rxte / pca but can not be achieved by the other instruments that we are going to compare its results against .
therefore , we formulated an approach that does not rely on directly comparing the extracted spectra of x - ray bursts simultaneously observed with two ( or more ) x - ray satellites .
instead , we first analyze the rxte data by performing time resolved x - ray spectroscopy .
we extract spectra for short intervals ( 0.25@xmath02.0  s depending on the countrate ) and fit these data with blackbody functions . based on the obtained @xmath5 values
, we then calculate the likelihood for a wide range of spectral parameters that cover the @xmath15@xmath6 range of the best fit values .
we use this range of individual parameter pairs ( blackbody temperature and normalization ) to create simulations of chandra or xmm-_newton _ x - ray spectra and then to calculate the total rxte - predicted countrate in a given energy range .
we calculate the uncertainty in these countrates by taking into account the posterior likelihoods of the rxte / pca parameters for each simulation . finally , we compare this predicted countrate with the countrate actually observed by chandra or xmm-_newton _ in the same energy range , taking into account all the associated uncertainties .    assuming that there are no calibration uncertainties , the chandra and/or xmm-_newton _ count rates that were predicted by analysing the rxte observation should be in perfect agreement with the actual count rates observed with these satellites . a systematic offset on the other hand ,
would imply a discrepancy in the overall calibration of either or both detectors .
it is evident that this method allows for the detection of a systematic uncertainty caused by the overall flux calibration of detectors and specifically in the determination of the normalization in a particular energy range .
however , the method does not allow us to constrain any calibrational uncertainties in the determination of the spectral shape .
we searched the archival rxte observations to find the best source that can be used for this purpose .
an ideal source should be a frequent burster with a number of simultaneous observations and predictable burst properties .
two low mass x - ray binaries in the archive conform to these conditions : the so - called clockwork burster gs  1826@xmath0238  and 4u  1636@xmath0536 .
the problem with the latter source is that its bursts do not last long enough ( only about 20  s on average ) to obtain sufficient amount of data to infer any systematic trends during the cooling tails of x - ray bursts . on the other hand , with bursts that last on average 130  s ,
gs  1826@xmath0238  provides a significant amount of counts to be used for this purpose .
gs  1826@xmath0238   was first discovered with the ginga satellite ( makino et al .
1988 ; tanaka   1989 ) .
first conclusive detections of thermonuclear x - ray bursts from this source came from bepposax observations ( ubertini  et  al .
1997 ) , although x - ray bursts may have been observed from this source earlier by oso  7 ( markert et al .  1977 ) and oso  8 ( becker et al .
1976 ) . using 2.5 years of monitoring observations by bepposax , during which 70 x - ray bursts were detected , ubertini  et  al .
( 1999 ) reported that the x - ray bursts show a quasi - periodicity of 5.76  hr , a remarkable characteristic of this source .
later , using 260 x - ray bursts observed with the wide field camera ( wfc ) on board bepposax over 6 years , cornelisse et al .
( 2003 ) showed that the burst recurrence time is indeed stable but decreases as the persistent flux increases .
a total of 40 x - ray bursts were observed with the rxte between 1997 - 2003 ( galloway et al .
2008 ) . during this time interval ,
the persistent flux steadily increased as the burst recurrence time decreased ( galloway et al .
another remarkable feature of the x - ray bursts observed from this source is their similarity . although the initial decay time - scale and the burst fluence show slight variations with persistent flux , the observed x - ray burst light curves are nearly identical and reproducible ( galloway et al .
2004 ) .
the nearly identical light curves of individual bursts , predictable recurrence times , and , finally , already existing simultaneous observations with other x - ray satellites ( chandra and xmm- _ newton _ ) make gs  1826@xmath0238   an ideal target to understand and quantify the systematic uncertainties in the measurement of neutron - star masses and radii introduced by the uncertainties in the absolute flux calibration of rxte / pca .
in this study , we use data from three different x - ray satellites .
gs  1826@xmath0238  was observed simultaneously by chandra and rxte for the first time in 2002 , where three x - ray bursts were recorded .
we present the analysis of these bursts below . furthermore , rxte observed the system in 2003 simultaneously with xmm-_newton_. during this observation , 7 x - ray bursts were detected with rxte .
we analyze here 5 of these events because , during one of the remaining bursts , xmm-_newton _ was affected by high energy particle background ( burst # 28 ) while the other burst was only partially observed and did not have high time resolution data with rxte / pca ( burst # 29 ) .
we present a log of these bursts in table  [ bursts ] .
cccc 22 & 52484.41775 & 70044 - 01 - 01 - 000 & chandra + 23 & 52484.56684 & 70044 - 01 - 01 - 00 & chandra + 24 & 52485.00672 & 70044 - 01 - 01 - 02 & chandra + 26 & 52736.53279 & 80048 - 01 - 01 - 00 & xmm-_newton _ + 27 & 52736.66488 & 80048 - 01 - 01 - 010 & xmm-_newton _
+ 28 & 52736.80048 & 80048 - 01 - 01 - 010 & xmm-_newton _
+ 29 & 52738.34118 & 80048 - 01 - 01 - 04 & xmm-_newton _ + 30 & 52738.47686 & 80048 - 01 - 01 - 04 & xmm-_newton _
+ 31 & 52738.61013 & 80048 - 01 - 01 - 07 & xmm-_newton _ + 32 & 52738.74636 & 80048 - 01 - 01 - 07 & xmm-_newton _ [ bursts ]      the rxte observations were analyzed following the method detailed in paper  i. time resolved , 2.5@xmath025.0  kev x - ray spectra were extracted from science event mode data ( 125  @xmath7s time resolution and 64 energy channels ) using the ftool _ seextrct _ from all of the rxte / pca layers . for each burst , spectra were integrated over 0.25 , 0.5 , 1 or 2  s time intervals so that the total number of counts in each spectrum remain roughly the same while the countrate decreases during a burst .
a 16  s x - ray spectrum was extracted from a time interval just prior to each burst , as the spectrum of the persistent emission , and subtracted from the x - ray burst spectra as background .
we took the pca deadtime correction into account following the procedure described in paper i.    response matrix files were created for each burst separately using the pcarsp version 11.7 , heasoft 6.7 and heasarc s remote calibration database .
x - ray spectral analysis was performed using the interactive spectral interpretation system ( isis ) , version 1.4.9 - 55 ( houck & denicola 2000 ) and , for each fit , a systematic uncertainty of 0.5% was included as suggested by the rxte calibration team .
we note that we also investigated the possibility of photon pile - up in the rxte / pca data following tomsick & kaaret ( 1998 ) and jahoda et al .
( 2006 ) . during the decay of the x - ray bursts investigated here , the observed count rates are @xmath81600  cts / s / pcu , which is far below the threshold of 10000 cts / s / pcu , where pile - up can become important according to those studies .
we , therefore , conclude that photon pile - up in the pca does not affect our results .
we fit each spectrum with a blackbody function using the _ bbodyrad _ model ( as defined in xspec , arnaud 1996 ) .
pinto et al .
( 2010 ) reported an analysis of high resolution x - ray spectroscopy of the interstellar medium toward gs  1826@xmath0238   and found an equivalent hydrogen column density of n@xmath9 . in order to take into account the effect of the interstellar extinction in the x - ray spectra
, we used this value with the _ phabs _ model in xspec .
we set the elemental abundances to the protosolar abundances of lodders ( 2003 ) , following pinto  et  al .
the resulting histogram of @xmath5/dof values obtained from the cooling tails of all the bursts and the spectral evolution during one of the bursts ( burst # 23 ) are given in figures [ chi_hist ] and [ rxte_burst23 ] .
similar to the results we obtained for a number of sources in paper  i , the blackbody model could adequately fit the data with a systematic uncertainty @xmath10 of only @xmath11 .
chandra observed gs  1826@xmath0238   in 2002 for approximately 70  ks ( obsid 2739 ) .
we analyzed the data using the _ chandra_repro _ tool with ciao version 4.3 and caldb 4.4.5 .
chandra_repro _ script automates the data processing following the standard pipeline steps but applies the latest available calibrations . using the _ dmextract _ tool we created a 2.5@xmath08.0  kev range background - subtracted source lightcurve using only the data from the grating arms ( + /-
1st orders ) .
the zeroth order data for this observation is saturated due to the high countrate of the source and the nominal frame time of 1.7  s.    even for the grating arms , pile - up was significant during the peak of the bursts . in order to correctly carry out the comparison between the rxte and chandra observations , we opt to correct for pile - up the _ predicted _ chandra countrates based on the rxte spectra , as opposed to correcting the observed chandra countrates .
because rxte observations have better statistics , this forward folding of the correction does not unnecessarily inflate the uncertainties .
we applied a pile - up correction procedure that is described in the appendix , which was developed based on the understanding of the detector properties and calibrated against independent data .
the evolution of the applied correction is shown in figure  [ ch_pileup ] for burst # 23 . at the burst peak ,
photon pile - up causes a net loss of 25% of the flux , which can be recovered by the pile - up correction model that was developed independently of the present data ( see the appendix ) .
nevertheless , the rise and the peak of each burst , which approximately correspond to the first 30  s , are most heavily affected by pile - up .
furthermore , the rise and the peak are also the periods where the intrinsic spectrum may show some deviation from a planckian shape ; this is expected from theoretical model spectra and was discussed in paper i using rxte observations .
because of these reasons , we exclude the first 30  s of all chandra bursts from our analysis such that in the time interval that we use , the pile - up fraction decreases to 2% .
we used the _ mkgrmf _ and _ fullgarf _ tools to create observation specific response and auxiliary response files for the first orders of the high energy transmission grating ( hetg ) and used these when creating the simulations based on the spectral parameters determined from time resolved x - ray spectroscopy .
xmm-_newton _ observed gs  1826@xmath0238  in 2003 for approximately 200  ks with two different back - to - back exposures that lasted 107,850 and 91,916  s ( obsids : 0150390101 , 0150390301 ) .
the observation was performed with epic - pn in timing mode and the epic - mos detectors were turned off .
we analyzed the epic - pn data using the sas version 11.0.0 ( xmmsas_20110223_1801 - 11.0.0 ) and the calibration files as of june 2011 .
we calibrated the observation with the _ epproc _ tool .
due to the high energy particle background we could only use the first 67  ks of the first observation and @xmath880  ks of the second .
similar to the chandra observation , we extracted background subtracted source light curves in the 2.5  @xmath0  10.0  kev range using the sas tool _ epiclccorr_. the observatory guide gives a limit of 450  ct / s for telemetry of the epic - pn timing mode , which can be increased if some of the other instruments are turned off as in this case .
in the burst lightcurves that were observed with xmm-_newton _ , two gaps that last approximately 10  s ( see figures [ comp_xmm1 ] and [ comp_xmm2 ] ) are seen .
these are most likely due to this limit . while creating the simulations we calculated observation specific effective area and redistribution matrix files using the sas tools _ rmfgen _ and _
arfgen_.    in order to explore the possibility of pile - up in the epic - pn data , we created two event lists : one including the full psf for the source and another where the central two pixels , which should be the most affected from pile - up , are excised .
if pile - up were significant , the ratio of the countrates in the two event lists would evolve as the count rate changes during the burst because the two event lists would be subject to different amounts of pile - up .
the ratio of the two lightcurves is shown in figure  [ xmm_pileup ] for bursts # 30 , 31 , and 32 . as expected ,
a large fraction of counts are collected in the central two pixels , making this ratio approximately equal to two . however , as the burst eveolves , the ratio remains constant at the few percent level , implying that the pile - up effect will be at most of the same order .
we consider this as strong evidence for the lack of any pile - up effect on the measured fluxes .
nevertheless , in order to avoid any residual effects of photon pile - up and any possible deviations from a planckian shape , as with the chandra bursts , we exclude the first 30  s of each burst when determining the systematic offsets between the instruments .
we show these time intervals in the plots for comparison and completeness .
we followed the method outlined in section  [ method ] for comparing all of the x - ray bursts observed simultaneously with rxte and chandra or with rxte and xmm-_newton_. briefly , our analysis steps are as follows :    1 .
we first obtained the best fit values of the parameters for each burst spectrum observed with rxte / pca .
we then calculated the two - dimensional confidence contours for the resulting color temperature and the blackbody normalization parameters covering a @xmath12 range from the best fit values .
an example confidence contour is shown in figure  [ example_conf ] for the spectrum obtained approximately 9.7  s after burst  # 26 started . note that we show 68% and 95% confidence contours in the figure but use the entire likelihood in the calculations .
2 .   for each pair of spectral parameters , we calculated the predicted countrates for the chandra gratings ( only the @xmath13  1st orders for heg and meg ) in the 2.5@xmath08.0  kev range and for xmm- _ newton_/epic - pn countrates in the 2.5@xmath010.0  kev range ( i.e. , the energy range that matches the rxte / pca sensitivity ) . for the case of the chandra simulations
, we applied the pile - up correction , as described in the previous section .
3 .   in order to take into account the uncertainties in the rxte / pca modeling
, we used the posterior likelihoods over the parameters for each spectral fit to calculate the uncertainty in the total predicted countrates .
4 .   finally , we compared the observed chandra hetg / acis - s or xmm-_newton_/epic - pn countrates to the countrates predicted using the rxte / pca spectral fits . for this comparison , we binned the countrates to match the time resolution of each data set and subtracted the persistent flux level of the source from the burst emission . during the chandra observation , the nominal frame time was 1.7  s. the time resolution of the epic - pn observation was much higher .
we , therefore , binned the epic - pn data to match the varying exposure time of rxte spectra .        in figures
[ comp_ch ] , [ comp_xmm1 ] , and [ comp_xmm2 ] , we show the time evolution of the source countrate during each x - ray burst simultaneously observed with rxte and chandra in the 2.5@xmath08.0  kev range or with rxte and xmm-_newton _ in the 2.5@xmath010.0  kev range . in principle
, all of the calculated ratios should be equal to unity within their formal uncertainties , if there are no systematic calibration differences .
it is already evident from this figure that this is indeed the case between rxte and chandra but that there is a systematic calibration offset between rxte and xmm-_newton_.    to explore in more detail the calibration differences between the three instruments and whether such differences evolve with countrate , we performed two analyses .
first , for the xmm-_newton_/epic - pn data , for which we infer the largest systematic calibration offset , we created histograms of the ratios of the predicted and observed countrates with varying start times compared to the burst start times ( see figure [ phtest2 ] ) .
if there were a significant countrate dependence of the calibration offset , we would have seen a significant change in the distribution of calculated ratios at later stages of these bursts .
however , it is evident from figure  [ phtest2 ] that there is no such time dependent change in the ratios .
second , in order to quantify the level of offset between the three instruments , we followed the bayesian procedure described in papers i and ii and outlined below for completeness .
we calculated the ratio @xmath14 between the observed and predicted countrates for each time bin in the 30@xmath0120  s interval after the start of every burst together with its formal uncertainty .
we assume that the uncertainty in each measurement is described by a gaussian with a centroid @xmath15 and a dispersion @xmath16 i.e. , @xmath17\;,\ ] ] where @xmath18 is a set of normalization constants .
we will also assume that the ensemble of these ratios for each burst is drawn from a gaussian distribution with a centroid at @xmath19 and a dispersion @xmath20 , i.e. , @xmath21\;,\ ] ] where @xmath22 is a proper normalization constant .
this distribution is our model for the underlying systematic flux calibration offset between the instruments .
if there are no calibration differences , the centroid should be equal to one , to within @xmath23 , where @xmath24 is the number of measurements , and the dispersion should tend to zero .
a significant offset in the flux calibration between two instruments will move @xmath25 away from unity , but keep the dispersion consistent with zero . on the other hand ,
a large dispersion can not be produced by an overall flux calibration offset .
instead , it likely points to effects such as an inadequate modeling of the pile - up correction , of the energy response of the detector(s ) , or of the astrophysical properties of the source .
our goal is to calculate the posterior likelihood for the parameters @xmath19 and @xmath20 , @xmath26data@xmath27 , given an ensemble of @xmath24 measurements .    using bayes theorem ,
we can write @xmath28 here @xmath29 is a normalization constant , and @xmath30 and @xmath31 are the priors over the centroid and dispersion of the underlying gaussian distribution .
we take these prior distributions to be flat in the range of parameters shown in figures  [ res_ch ] , [ res_xmm1 ] , and [ res_xmm2 ] , which is large enough to incorporate any range of interest .
the quantity @xmath32 measures the likelihood that we will make a particular set of measurements given the values of the parameters of the underlying gaussian distribution .
because these measurements are uncorrelated , we write @xmath33 inserting equation  ( [ eq : data ] ) into equation  ( [ eq : bayes ] ) we obtain the posterior likelihood @xmath34 where @xmath29 is an overall normalization constant .
we show two - dimensional contours of these likelihoods for each burst in figures  [ res_ch ] , [ res_xmm1 ] , and [ res_xmm2 ] and display the most likely values of the parameters in table  [ res_all ] .
all figures show that the underlying distributions of offsets are extremely narrow and that the most likely values of their dispersions are consistent with being zero .
this indicates that there are no systematic problems in the modeling of the energy response of the detectors and that the corrections described earlier adequately account for any pile - up effects .
the most likely values of the centroids are consistent with unity for the comparison between rxte and chandra but as low as 0.86 for the comparison between rxte and xmm-_newton_.    because the offsets are stable both throughout the cooling tails of individual bursts and from burst to burst , we combined all the measurements from all the bursts observed by each x - ray satellite to quantify the overall offset and its uncertainty .
we computed the posterior likelihood for the offset @xmath25 given the two dimensional posterior likelihood of the centroid value @xmath19 and dispersion @xmath20 of the underlying gaussian distribution as @xmath35    using 5 bursts and 447 individual measurements and 3 three bursts and 170 individual measurements for comparison with xmm-_newton _ and chandra , respectively , we found that the ratio of the observed count - rate with the chandra and xmm-_newton _ compared to the predictions based on the rxte spectral analysis in the 2.5  @xmath0  8.0  kev and 2.5  @xmath0  10  kev range is 1.012@xmath36 and 0.857@xmath10.003 , respectively .
in other words , the fluxes measured with chandra hetg / acis are in statistical agreement with the rxte / pca while the epic - pn measures one that is 14% less . in figure
[ final_overplot ] , we show the combined likelihoods of the resulting ratios of countrates between the different instruments
.                  lcc 22 & 1.01 & hetg - acis / pca + 23 & 1.01 & hetg - acis / pca + 24 & 0.97 & hetg - acis / pca + 26 & 0.86 & epic - pn / pca + 27 & 0.85 & epic - pn / pca + 30 & 0.86 & epic - pn / pca + 31 & 0.86 & epic - pn / pca + 32 & 0.86 & epic - pn / pca [ res_all ]
in the previous section , we addressed the cross calibration between rxte and _ chandra _ as well as between rxte and xmm- _ newton _ , using observations of the burster gs  1826@xmath0238 .
we obtained the most likely fractional offsets between the fluxes measured using the various instruments on simultaneous observations and showed that they are consistent with being constant throughout individual bursts as well as across different bursts .
however , the offsets we inferred correspond to a single point in time for each comparison .
if the efficiencies of the detectors were constant in time , the results of this comparison would be applicable to all the bursts observed thoughout the durations of the missions .
however , the responses of the detectors evolve with time for a number of reasons ( e.g. , gas leakage , charge transfer inefficiency ) . this evolution is typically accounted for using smooth , time - dependent models of the response parameters of the detectors .
these parameters are determined by fitting the models to regular observations of the cosmic sources with known brightness and spectral shape or using internal calibrations sources ( e.g. epic cameras ) . as an example
, rxte / pca had been calibrated against the crab nebula ( shaposnikov et al .
2013 ) throughout the mission .
the smooth models of the evolution of the detector response clearly can not reliably capture the true evolution , which is expected to have a stochastic component .
this will , therefore , introduce an uncertainty in the measurements that will change in time .
[ line - through]*t*his complication is exacerbated if the calibration source shows intrinsic variability , as is the case with the crab nebula , which showed variability at a level of @xmath8  7% ( wilson - hodge et al .
it is , therefore , possible that were we to repeat our cross calibration procedure between rxte and chandra at different epochs , we would obtain different offsets .
even though randomly catching an epoch of identical flux measurements between rxte and chandra might seem unlikely , we can not base our conclusions on such a posterior likelihood .
for this reason , we now explore the possibility that there is a time - dependent systematic uncertainty in the calibration of the instruments and how this uncertainty would affect the scatter in the measured burst fluxes .
if there were single - epoch measurements of the fluxes during the cooling tails of bursts , this potential systematic error would have to be added to the error budget of the measurements .
however , for many of the sources we used for neutron - star radius measurements in papers i and ii , the burst flux measurements were carried over many epochs spanning the lifetime of rxte . in this case , as we will show below , this potential calibration uncertainty should not be treated as a separate source of systematic uncertainty as it is already inherently part of the measured scatter of observed fluxes .    to demonstrate this point , we compute the mean and the variance of the ( burst ) flux measured by an instrument ( e.g. , rxte / pca ) . let the source flux for each burst be @xmath37 , with a mean of @xmath38 and an intrinsic dispersion @xmath39 .
let the instrument have an inferred overall efficiency factor @xmath40 ( e.g. , as reported by a calibration team ) .
if this efficiency factor has an unknown systematic bias , which may also be variable in time , then the observed countrate flux will be @xmath41 where we have denoted the mean bias in the calibration over the time span of the observations by @xmath42 and the potential time - variable part by @xmath43 . by definition ,
@xmath44 .
when an observation is carried out and standard calibration procedures are applied , the observed source flux will be inferred to be equal to @xmath45 which will be related to the true source flux by @xmath46 clearly , the observations reproduce the true source flux if there is no bias in the calibration .
the mean inferred flux of the source is equal to @xmath47 where we have used the fact that @xmath48 because the time evolution of the detector response and the variability of the source are uncorrelated .
the variance in the inferred burst flux is @xmath49 and the mean of the square of the inferred fluxes is given by @xmath50 ^ 2    \right\rangle \nonumber \\    & = & \left(1+\frac{b}{r}\right)^2 \langle    s_{i,\rm true}^2\rangle + \left\langle \left(\frac{v_i}{r}\right)^2    s_{i,\rm true}^2 \right\rangle + 2\left(1+\frac{b}{r}\right ) \left\langle    \frac{v_i}{r }
s_{i,\rm true}^2 \right\rangle.\end{aligned}\ ] ] since @xmath51 , we find that @xmath52      + \frac{1}{r^2}\left\langle v_i^2s_{i,\rm true}^2 \right\rangle \;.\end{aligned}\ ] ] the term in the square bracket is equal to the true variance of the source flux @xmath53 = \sigma_{\rm true}^2 $ ] .
the average in the second term is @xmath54 , where we have used the definition of the variance of the time - variable bias @xmath55 as a result , @xmath56 or in terms of fractional errors , @xmath57\;.\end{aligned}\ ] ] keeping only terms up to second order in the dispersions , the final result is @xmath58 this expression shows that the measured fractional variance in the burst fluxes is equal to the quadrature sum of the intrinsic fractional variance and the fractional variance of the calibration bias .
in other words , the observed variance in the burst fluxes is an upper limit to the intrinsic scatter .
if the statistical properties of the calibration uncertainties could be measured independently , this contribution could be subtracted from the dispersion of measured fluxes .
in this paper , we quantified the systematic uncertainties that arise from the flux calibration of the x - ray instruments used in neutron star mass and radius measurements .
we used simultaneous observations of the burster gs  1826@xmath0238  with rxte / pca and chandra hetg / acis - s as well as with rxte / pca and xmm-_newton_/epic - pn .
we performed a bayesian statistical analysis to measure the level of offset in the flux calibration between these instruments .    in the case of the comparison between chandra and rxte
, we found no evidence for a significant offset between the measured fluxes and concluded that the most likely value of the flux ratio is @xmath59 .
as we describe in the appendix , the effect of pileup is unlikely to change the value of this ratio in a statistically significant way .
in contrast , we found a 14.0@xmath10.3% difference between the fluxes inferred using xmm-_newton_/epic - pn and rxte / pca .
this discrepancy is persistent in all five analyzed bursts and is independent of the source spectrum or countrate .
the flux discrepancy between xmm-_newton_/epic - pn and other detectors has been noted in earlier studies of different sources .
nevalainen et al .
( 2010 ) reported a difference between flux measurements with the epic - pn , chandra acis ( zeroth - order data ) , and epic mos detectors at a level of 5@xmath010% . according to their results ,
acis flux measurements are significantly higher than epic - pn , by 11.0@xmath10.5% .
similarly , the difference in the flux measurements between the epic - pn and mos detectors has been noticed earlier ( see e.g. , mateos et al .
2009 and xmm - newton calibration documents xmm - newton - soc - cal - tn-0018 and xmm - newton - soc - cal - tn-0052 ) ; mos yields fluxes that are approximately 7% higher than epic - pn ( see also read et al . 2014 ) .
most recently , using 64 clusters of galaxies , schellenberger et al .
( 2014 ) showed that the effective area cross - calibration uncertainty between acis and epic - pn is energy dependent and in the 2@xmath07  kev range stabilizes at a 15% level a similar result is also presented by guainazzi et al .
( 2015 ) .
given that the calibration errors of the detectors may be variable in time , it is possible that the excellent agreement between rxte and chandra is a coincidence .
we showed , however , that even in this case , the spread of the measured burst fluxes reported in papers i and ii already accounts for this potential systematic uncertainty .
we , therefore , conclude that , while an overall bias may still exist in the measurements , potential time - dependent flux calibration uncertainties do not increase the error budget of the radius measurements we reported in these earlier papers .
one final consideration is whether the difference in the measured flux between xmm-_newton _ epic - pn and the other two instruments also affects the neutron star radii determined through the observations of quiescent low - mass x - ray binaries .
the flux levels in question are very different in the observations of these dim sources : they are approximately six orders of magnitude below the peak burst fluxes and the crab nebula , and about three orders of magnitude below other typical calibration sources ( e.g. , clusters of galaxies ) . given the signal - to - noise ratio of the spectra that is achieved with even the longest observations of quiescent low - mass x - ray binaries , it is not possible to discern any measured flux differences between instruments , even if this difference persisted at the 14% level .
this supports the lack of evidence for flux differences presented in the earlier studies ( catuneanu et al .
2013 ; heinke et al .
2014 ; guillot et al .
2013 ) , which compared neutron star radii derived from the quiescent sources using xmm-_newton _ and chandra data .
to account for pileup in the chandra hetg / acis - s spectra , we used a correction that was developed independently by one of the co - authors , herman marshall , for mildly piled - up calibration data .
the approach relies on the measurement of counts per pixel around the ir - m and si - k edges , at which the spectral response changes rapidly . denoting the observed counts per detector pixel per frame in wavelength bin @xmath60 as @xmath61
, we define the degradation of the effective area as        the value of the parameter @xmath63 and the uncertainty associated with its determination was found using an observation of mk  421 ( observation i d 4148 ) performed on 26@xmath027 october 2002 with the letgs / acis - s for over 91  ks .
the data were downloaded from _ tgcat _ ( huenemoerder et al . 2011 ) .
the countrate spectra were rebinned by a factor of two so that each spectral bin corresponds closely to one acis pixel .
the source was extremely bright so that there were of order 5000 counts per spectral bin in both the + 1 and -1 orders near 2  kev .
the count spectrum was divided by the arf obtained from tgcat to generate spectra in photons/@xmath64/s /   for both the @xmath65 and @xmath66 orders . when the two orders are compared ( see the top panel of figure  [ final_overplot ] ) without being corrected for pileup , there are clear edges at ir - m ( 6.0   ) and si - k ( 6.74   ) in the effective - area - corrected spectra .
these features arise from differential pileup , where @xmath61 may increase significantly over a few spectral resolution elements at an edge in the effective area of the instrument , causing increased pileup and a larger fractional loss of events with a corresponding decrease in apparent photon flux .
we devised a model - independent method for estimating the parameter @xmath63 in equation  ( [ eqa1 ] ) , which depends on the very smooth and slowly varying photon spectrum of mk  421 . using the ratio of fluxes in the 5.5 - 5.9    range to those in the 6.1 - 6.5   range ( straddling the ir - m edge ) , we compute the value of @xmath63 that matches the ratios based on the average @xmath40 values in the two narrow bands .
we find @xmath67 for @xmath66 and @xmath68 for @xmath65 and , therefore , adopt @xmath69 and a range of acceptable values that are @xmath70 away from this value .
figure  [ final_overplot ] shows the sensitivity of the spectral inversion to @xmath63 ; for @xmath71 , the spectra still show edges , while , for @xmath72 , the data are over - corrected to form an excess over a smooth spectral model .
( a global spectral fit to a double power law indicates that the photon spectrum should drop monotonically by only 1% per   through this wavelength range , so the assumption that the fluxes are equal should be good to 1% . )    using equation  ( [ eqa1 ] ) , which is verified to @xmath73 cts / frame / pixel , a @xmath74% uncertainty in @xmath63 results in at most a @xmath75% uncertainty in @xmath76 , which is smaller than the expected calibration uncertainty in @xmath77 .
this method required the hrma effective area correction implemented in 2005 , which corrected the ir - m edge caused by a thin contaminant layer on the hrma .
note that , the correction given in equation  ( [ eqa1 ] ) is valid for low to modest pileup , i.e. , @xmath78  cts / frame / pixel and the data used to compute the coefficient @xmath63 ranged up to about 0.10 cts / s / pixel .
this yields a maximum correction of about 50% .    in our comparison between rxte and chandra
, we took into account the pileup losses in the predicted countrates rather than correcting the observed rates for pileup . if we denote by @xmath79 the number of photons per detector pixel per frame that should have been recorded in wavelength bin @xmath60 in the absence of pileup , we can then invert the above expression to obtain                                      houck , j.  c. , & denicola , l.  a. 2000 , in astronomical society of the pacific conference series , vol . 216 , astronomical data analysis software and systems ix , ed .
n.  manset , c.  veillet , & d.  crabtree , 591+