in this paper we study the effects of rotation on the radiation driven winds of early - type supergiants .
we will focus on the explanation for the occurrence of disks around fast rotating b[e ] supergiants .
b[e ] supergiants , also designated sgb[e ] stars ( lamers et al .
@xcite ) , are b type supergiants that exhibit forbidden emission lines in their optical spectra . the observations of hybrid spectra , i.e. spectra with _ broad _
uv p cygni features and _ narrow _ emission lines and dust emission of b[e ] supergiants in the magellanic clouds led to the proposal of a disk wind model by zickgraf et al.(@xcite ) .
this model postulates a dense disk of outflowing material in a fast line driven wind to explain the observed characteristics of the spectra of these stars ( see zickgraf  @xcite ) .
the stellar wind at the equator is about ten times slower than that at the pole .
also , the wind at the equator is about a hundred times denser than at the pole . additional evidence for a two - component outflow has been obtained from polarimetric measurements , e.g. by zickgraf & schulte - ladbeck  ( @xcite ) .
the precise mechanism behind the formation of these disks is still a mystery ( cassinelli 1998 ) .
it is however clear that the origin of an axisymmetric wind structure such as a disk may well be connected to the fast rotation of a star .
two theories using rotation in a different way have been considered for the formation of these disk winds :     ( 1 ) the wind compressed disk ( wcd ) model of bjorkman & cassinelli  ( @xcite ) , that invokes the kinematics of the winds from rotating stars .
the streamlines of the gas in the wind from both hemispheres of a rapidly rotating star cross in the equatorial plane .
the concentration of the gas and the shock in the equatorial plane produce an outflowing equatorial disk , with a thickness on the order of a few degrees .
owocki & cranmer ( 1994 ) have argued that the motion to the equatorial plane may be counteracted by the radiation force perpendicular from the plane due to lines .
however , this `` wind - inhibition effect '' may not be effective in winds with a strong density gradient in the equatorial direction .
we return to this in sect .
( 2 ) the rotationally induced bi - stability model ( rib ) of lamers & pauldrach   ( @xcite ) invokes an increase in the mass flux from the equator and a decrease in the equatorial wind velocity compared to the poles from the bi - stability jump .
this jump in mass flux and in wind velocity is due to the temperature difference between the pole and the equator of a fast rotating b[e ] supergiant with gravity darkening .
the jump will occur for stars with effective temperatures between 20  000 and 30  000 k.    in reality , both effects , i.e. the wind compression and the rotation induced bi - stability , may be operating together and amplifying one another in the line driven winds of rapidly rotating early - b stars .
( for a detailed explanation of both models , see lamers & cassinelli ( 1999 ) ( hereafter @xmath5 ) , chapter 11 . )
the theory of radiation driven winds for non - rotating stars was developed by castor et al .
( @xcite ) ( hereafter cak ) and predicts the mass loss rate and @xmath6 for spherical winds .
the influence of rotation on line driven winds was investigated by friend & abbott ( @xcite ) .
they found an increase in the mass - loss rate at the equator and a corresponding decrease in @xmath6 , but these authors did not consider the effects of gravity darkening and the oblate shape of the star .
cranmer & owocki  ( @xcite ) described the finite disk correction factor in case of an oblate star .
in this paper we will modify the line driven wind theory , including these effects of oblateness and gravity darkening as well as the rotational terms in the equation of motion .
this will then be applied to rapidly rotating b[e ] supergiants to investigate whether these effects can explain the occurrence of disks .
we will also apply recent bi - stability calculations to study the effect of the rib - model to explain the occurrence of outflowing disks around rapidly rotating b[e ] supergiants .    in sect .
2 we describe the theoretical background of the radiation driven wind theory and the bi - stability effect . in sect .
3 we derive the equations for the winds from oblate rotating stars with a temperature gradient between poles and equator due to gravity darkening . in sect .
4 we describe the method for solving the equations and for calculating the winds from rotating stars .
these calculations will be applied to b[e ] supergiants in sect .  [ s_b[e ] ] .
in sect .
[ s_bs ] we will investigate the effect of bi - stability on rotating stars and sect .
[ s_concl ] concludes with a summary of this work and a discussion on the formation of disks of b[e ] supergiants .
the computation of the dynamics of line driven stellar winds is a complicated problem in radiation hydrodynamics .
it involves the simultaneous solution of the equations of motion , the rate equations and the radiative transfer equations in order to calculate the radiative acceleration .    for hot luminous stars the line forces are the most important driving forces in the wind .
cak have shown that the line acceleration can be parameterized in terms of the optical depth parameter @xmath7 as @xmath8 , where @xmath0 and @xmath1 are parameters that depend on composition and temperature of the wind . in this expression @xmath0
is a measure of the number of lines and @xmath1 is a measure of the distribution of the line strengths with @xmath9 or 1 for a pure mix of optically thin or thick lines respectively .    in this paper
we adopt the cak formalism and simplify the equations of motion by assuming a stationary , radial flow of a viscousless fluid .
the possible influence of magnetic fields will be ignored .
these simplifications are subject to the following restrictions ( see e.g. abbott 1980 ) :    * the cak theory treats the line acceleration in a simplified way in terms of the force multipliers , described above .
detailed calculations of the line driving , including millions of lines , show that the total line acceleration can be well described by this simple representation ( pauldrach et al .
both observations and calculations show that line driven stellar winds are not stationary .
even a wind in a stationary solution will develop shocks ( see e.g. lucy  @xcite ) .
however , the time - averaged structure follows the stationary state quite well ( see owocki et al .
@xcite , and feldmeier 1999 ) .
therefore we will restrict ourselves to stationary models . *
a sectorial model will be adopted .
this means that for every latitude a one dimensional problem will be solved .
wind compression as in the wcd model will thus be neglected , but the main effect of wind compression is expected to _ redistribute _ the mass loss and not to change the _ total _ mass that is lost from the rotating star . * the absence of viscosity is a good approximation at the high density of line driven winds ( see cak ) .
the representation of the line force by a simple power law may seem to be a gross simplification of the underlying physics of myriad line absorption processes , but it can be shown to hold for a homogeneously distributed mix of optically thick and thin lines ( abbott  @xcite , gayley  @xcite ) .
calculations of cak and more recently in nlte by vink et al .
( @xcite ) of realistic model atmospheres , based on calculations of unified atmosphere / wind models with the isa - wind code ( de koter et al .
1997 ) , confirm this to be generally true to good accuracy for the important part of the wind ( from about the sonic point up to a few stellar radii ) .
therefore , the complex physics of ion populations can be ignored in the investigation of the various effects of rotation on the stellar wind .
however , there is one notable exception which is called the bi - stability jump .
this refers to a jump in @xmath6 around the temperature of 21000 k where @xmath10 climbs from 1.3 ( lower temperature ) to 2.6 ( higher temperature ) .
this jump was observed by lamers et al .
( @xcite ) in b supergiants .
it is linked to a shift in ionization states of fe in the lower part of the wind .
the line driving in the lower part of the wind is dominated by iron .
below about @xmath11 k , fe iv recombines to fe iii and since fe iii is a more efficient line driver than fe iv , the wind structure changes dramatically ( vink et al .
@xcite ) .
the bi - stability jump may also occur in the temperature difference between the pole and the equator of a fast rotating star with gravity darkening .
therefore two sets of force multiplier parameters for the wind will be adopted to reflect the sudden change in ionization states ( a high - temperature set for the pole and a low - temperature set for the equator ) .
intuitively it is clear that rotation has an effect on the shape of a star and on the motion of the gas in the wind .
an additional effect that the rotation of a star can have is the darkening of the equatorial regions of a star via the von zeipel effect .
these effects will modify the wind of a star considerably and we will incorporate them in the line driven wind theory for rotating stars .
the shape of a uniformly rotating star with all its mass concentrated in the core is determined by the equipotential surfaces of the potential in a rotating frame ( roche model ) :
@xmath12    where @xmath13 , @xmath14 is the co - latitude ( @xmath14 = 0 at the pole ) and @xmath15 .
@xmath16 is of course independent of the longitude @xmath17 .
note that @xmath18 .
so the maximum oblateness is @xmath19 . the critical ( break - up ) velocity is defined in this paper as : @xmath20 where @xmath21 is the effective mass of the star ,    @xmath22    which is the reduced mass due to radiation pressure by electron scattering ( see below ) .    the surface of a rotating star is implicitly given by eq .
[ eq_firot ] . solving this expression for @xmath23
is equivalent to the solution of a 3rd degree polynomial , which results in    @xmath24    because of the difference in radius    @xmath25    and the difference in rotational velocity between pole and equator , the radial escape velocity at the equatorial region is smaller than at the pole    @xmath26    we can then roughly estimate the effect of rotation on @xmath6 as this value has an almost linear dependence on @xmath27 ( see later in  [ s_ssol ] ) .
the von zeipel theorem  ( @xcite ) for distorted stars states that the radiative flux from a point on the star is proportional to the local effective gravity : @xmath28 with @xmath29 the boltzmann constant . as @xmath30 , we can write down the flux of a rotating star as a function of the co - latitude @xmath14 . this ( lengthy ) expression for the flux as a function of @xmath14
is plotted in fig .
[ fg_f ] for a few different values of @xmath31 .
the rotation of the star can induce a @xmath14 component of the flow of the matter in the wind . in this study a possible @xmath14 component is neglected as was discussed in  [ s_physics ] , but the motion in the longitudinal @xmath17 direction must still be taken into account .    in the absence of forces that can exert a torque on the wind , the equation of motion for @xmath32
is given by the conservation of angular momentum .
note however that there could be a torque from the line forces themselves , which is dependent on the velocity gradient of the gas .
we assume that there are no external torques , so conservation of angular momentum gives :    @xmath33    the equation of motion for the radial direction becomes :    @xmath34    where @xmath35 is the radial velocity and @xmath36 is the stellar mass .
the conservation of mass for a non - spherical sectorial wind can be written as    @xmath37    where @xmath38 is the  local mass loss rate ,
i.e. the total mass loss rate if the solution for this latitude where valid for a spherical star
. the total mass loss rate from the star is    @xmath39    the equation of motion together with the conservation of mass governs the dynamics of the wind given the equation of state , @xmath40 ( @xmath41 is the isothermal sound speed ) , and the radiative acceleration @xmath42 . using the conservation of mass , the pressure term can be rewritten as : @xmath43    as in the cak - theory the temperature structure is specified a priori as a function of @xmath44 .
( the results depend only very weakly on the chosen temperature structure ; see pauldrach et al .
@xcite ) .
the radiative acceleration consists of two components .
the continuum component due to the electron scattering , @xmath45 , and a second component , @xmath46 , due to line scattering and absorption processes .
the continuum acceleration is given in terms of @xmath47 : @xmath48 where @xmath49 is the electron opacity , @xmath50 is the stellar luminosity and @xmath51 is the radiation flux .
the contributions by other continuum opacities are small .
nevertheless , they were included in the calculation of the force multipliers by vink et al .
( 1999 ) used in our bi - stable wind models ( sect .
for a homogeneous spherical star , @xmath47 is given by as the ratio between the continuum force and gravity and not as the ratio between continuum force and the critical radiation force for a rotating star ] :    @xmath52    in the more general case of a non - homogeneous , non - spherical star the continuum acceleration can be defined as a correction to the classical continuum acceleration of eq .
[ eq_gam ] : @xmath53 where @xmath54 is the continuum correction factor which is given by :    @xmath55    the line acceleration @xmath56 has a more complicated form , since it is a function of distance @xmath44 as well as the velocity @xmath35 and velocity gradient @xmath57 . combining the equation of motion ( eq .  [ eq_bewvgl ] ) with the rewritten pressure and continuum acceleration terms and multiplying by @xmath58 , one finds :    @xmath59    this equation is valid for each co - latitude @xmath14 .
following cak , the critical point of this equation is found by imposing the singularity condition ( where the subscript @xmath60 indicates values at the critical point ) : @xmath61 and the regularity condition : @xmath62 note that the critical point is not the sonic point @xmath63 ( @xmath64 ) due to the fact that there is an additional dependence on @xmath65 in the line acceleration @xmath46 ( abbott  1980 ; @xmath5 chapt 3.3 ) .
the line forces are described within the framework of cak in the sobolev approximation .
this means that the intristic line absorption profiles are considered to be infinitely sharp .
then the line force becomes a function of local properties of the wind only .
the acceleration due to an ensemble of lines is given by the summation of the contributions of all individual lines with rest frequencies @xmath66 .
@xmath46 is given by ( castor  @xcite ) :
@xmath67    where @xmath68 is the absorption coefficient per gram for the @xmath69th line , @xmath70 , @xmath71 the intensity at the rest frequency @xmath72 and @xmath73 is the sobolev optical depth of @xmath69th line , defined as @xmath74 .
the integration is performed over the complete visible disk of the star .
following cak , the line acceleration can be rewritten in terms of the force multiplier as a function of the optical depth parameter @xmath75 . where @xmath75 is defined as : @xmath76 with @xmath77 being the thermal velocity of the protons and @xmath78 is some reference value for the electron scattering @xmath79 @xmath80  g@xmath81 ( see _ isw _ , chapt .
the line acceleration in terms of the continuum acceleration is give by :    @xmath82    if the star is assumed to be a point source , @xmath83 can be approximated by a simple power law parameterization ( cak , abbott 1982 ) : @xmath84 where the @xmath85 accounts for the effect of the electron density on the ionization balance in the wind .
here @xmath86 is the electron density in units of @xmath87 @xmath88 and @xmath89 is the geometric dilution factor .
friend & abbott ( 1986 ) and pauldrach et al .
( 1986 ) showed the importance of the finite disk correction on the line force in case of an extended source .
this finite disk correction factor @xmath90 is given by ( cak ) :
@xmath91    where    @xmath92    this expression is completely general : the actual shape of the star and its intensity distribution enter through the integration domain and the @xmath93 dependence .
however , here we assume the intensity @xmath94 to be locally isotropic , i.e. we neglect limb darkening .
we have tested the effect of limbdarkening .
the inclusion of limbdarkening results in a small @xmath14-independent correction on @xmath90 .
the line acceleration @xmath46 including all the neccesary correction factors is now given by :    @xmath95    combining eqs .
[ eq_bewvglf ] and  [ eq_gl ] gives the full equation of motion for the sectorial wind of a rotating star    @xmath96
for line driven winds of non - rotating stars solutions of the equation of motion have been retrieved by e.g. cak and pauldrach et al .
( 1986 ) . these will be presented to serve as an illustration of the solution of the equation of motion and as a basis to interpret the more complicated results from the full equation of motion ( [ eq_bewvglrot ] ) .        if the star is considered to be a point soure the integral of the line acceleration in eq .
[ eq_gmeerlijn ] collapses to one point . in this case
@xmath90 disappears from eq .
[ eq_gl ] and the line acceleration becomes a simple function @xmath97 .
the solution of eq .
[ eq_bewvglf ] was found by cak and can easily be explained by neglecting the gas pressure terms ( @xmath98 ) . in this case
[ eq_bewvglf ] reduces to    @xmath99    where @xmath100 is a constant containing the mass loss rate .
the equation is solved by imposing uniqueness of the solution ( kudritzki et al .
1989 ) which then fixes the value of the constant and thus @xmath101 .
the solution for the mass loss and velocity law is given by    @xmath102    @xmath103    with    @xmath104    this simplified solution is equal to the full cak solution , in the limit of small sound speed , @xmath105 ( see also @xmath5 chapt 8) .      in case of a homogeneous spherical star
the finite disk correction factor @xmath90 ( eq .  [ eq_corr ] ) can be calculated analytically ( cak ) :    @xmath106    including the finite disk correction results in an increase in @xmath6 , a decrease in @xmath107 and a modification of the simple scaling laws that were found in the original cak approach ( see friend & abbott  @xcite ; pauldrach et al .
the decrease of the mass loss rate compared to the point source case is due to the decrease in @xmath46 close to the star where @xmath101 is determined .
close to the stellar photosphere the finite disk correction is smaller than one , viz . @xmath108 .
the accompanying increase in @xmath6 of typically a factor of two is due to two effects : ( 1 ) a reduction of @xmath107 results in a smaller amount of material to be accelerated and ( 2 ) far from the photosphere , the correction factor @xmath90 becomes larger than 1
. the resulting dependence between @xmath6 and @xmath27 is approximated by friend & abbott ( 1986 ) :    @xmath109    they also found a relation for the modified mass loss rate :    @xmath110      the solution of eq .
[ eq_bewvglrot ] gives the mass loss rate , or rather the local mass loss rate @xmath38 , and velocity structure of the wind .
the solution is complicated however by the presence of @xmath90 .
this is an integral of a velocity - dependent function times the surface intensity over the visible section of the star .
the appearance of the star varies throughout the wind in intensity distribution as well as in shape , as the star is no longer spherical and @xmath111 is a function of latitude .
since the analytic solution would be cumbersome , a numerical approach of the solution is chosen .
this is done with the analytic solutions of simple , i.e. the non - rotating , models in mind for comparison .
the numerical solution of eq .
[ eq_bewvglrot ] is relatively straightforward in case the function @xmath90 is a given function of @xmath44 .
therefore we solve eq .
[ eq_bewvglrot ] with @xmath112 in an iterative way as follows :    1 .
a @xmath113 law for @xmath114 is assumed , viz .
@xmath112 is calculated using this velocity law .
note that @xmath90 is not dependent on the actual value of @xmath114 but on the velocity gradient .
3 .   eq .
[ eq_bewvglrot ] is solved using this correction factor @xmath112 , obtaining a new velocity law @xmath114 .
step 2 - 3 are repeated until convergence is reached .
typically three iterations are sufficient , since further iterations changed the obtained values by less than 0.5 % .
we checked that the solutions converged to the same value , independent of the starting value .
the calculation of the correction factor @xmath90 from the velocity law ( step 1 . in the scheme described above ) is performed by numerical evaluation of eq .
[ eq_corr ] .
this is a non - trivial task , since both the integrand ( eq .  [ eq_zeipel ] ) and
the shape of the visible `` disk '' ( eq .
[ eq_rtheta ] ) are dependent on @xmath44 .
the main parameters needed to obtain @xmath90 for a fixed co - latitude @xmath116 are : the stellar surface @xmath117 and the surface temperature @xmath118 , both given as a function of the co - latitude @xmath14 , as well as the velocity law @xmath114 .
it is convenient to change the integration variables in eq .
[ eq_corr ] from the wind centered coordinates @xmath14 and @xmath17 to the star centered coordinates @xmath119 and @xmath17 .
for the definition of these coordinates , see fig .
[ fg_ofr ] .
note that here @xmath17 is not the stellar longitude anymore , but the angle of rotation around the line from a point @xmath120 in the wind to the stellar center .
@xmath121    where @xmath94 is the frequency - integrated intensity    @xmath122    the angles are related to one another via    @xmath123    with the stellar surface given by @xmath124 , and @xmath125 equation  [ eq_coint ] can be integrated , considering that the limb angle @xmath126 can be determined numerically , yielding the finite disk correction for the line acceleration of an arbitrary shaped star .    the correction factor @xmath54 for the electron acceleration ( defined in eq .
[ eq_dc ] ) similarly becomes :    @xmath127    this correction is also included in the solution of the equation of motion .
[ tb_res1 ]    [ cols="^,^,^,^,^,^,^,^,^,^,^,^,^,^,^,^,^ " , ]     ( 1 ) : the adopted stellar parameters are : @xmath128 ; @xmath129 ; @xmath130 k.      to solve the equation of motion it is neccesary to determine the conditions at the critical point : @xmath131 , @xmath132 and @xmath133 , or , equivalently , @xmath131 , @xmath132 and @xmath134 , as well as the local mass loss rate @xmath135 .
apart from eqs .
[ eq_sing ] and  [ eq_regu ] and the equation of motion ( [ eq_bewvglrot ] ) itself the additional constraint for the continuum optical depth @xmath136 at the stellar radius is neccesary ( pauldrach et al .
1986 ) to uniquely determine these quantities :    @xmath137    if the condition for @xmath138 is used one may adopt any small value , say @xmath139 km / s , since even a large error in @xmath138 will only yield a relatively small error in @xmath6 and @xmath135 .
this is due to the fact that at the surface of the star @xmath114 and thus also @xmath140 vary exponentially with scale height @xmath141 .
this means that the radius where eq .
[ eq_tau2/3 ] is fulfilled is very close to @xmath142 .
this has been checked a posteriori . in our calculations
we have adopted this boundary condition for @xmath138 in the solution of the full momentum equation .    a first guess for the local mass loss rate obtained from the cak solution
is used to integrate the equation of motion from the stellar surface outward .
this gives an approximate value for the critical point where @xmath143 if @xmath144 .
this approximate value for the critical point can then be used to consecutively improve the guessed value of @xmath135 using eqs .  [ eq_sing ] and  [ eq_regu ] .
this iterative procedure is terminated when the local mass loss rate converges . the calculation is then completed by integrating past the critical point on the solution which extends to infinity ( see cak or abbott 1980 for a description of the topology of the solutions ) . to deal with a numerical singularity at the sonic point the equations
are solved for @xmath145 rather than for @xmath114 .
we have calculated models for rotating b[e ] stars using the above described method . for the first models we adopted the following stellar parameters , which are typical for b[e ] stars : @xmath146 , @xmath147 , @xmath148 , @xmath149 and solar abundances .
the effective temperature is defined as @xmath150 where @xmath151 is the total surface of the distorted star .
the following force multiplier parameters from pauldrach et al .
( @xcite ) were adopted to describe the line acceleration : @xmath152 , @xmath153 , @xmath154 .
an overview of the results for this generic b[e ] star rotating at @xmath31 = 0 , 0.3 , 0.4 , 0.5 and 0.6 times the critical rotation speed is given in table 1 .
the calculation of models with higher values of @xmath31 runs into numerical problems , because in the iteration , the approximate critical point can be below the sonic point .
the table shows for different co - latitudes @xmath14 and different rotation velocities @xmath31 the values of @xmath6 , the local mass loss rate @xmath38 and the value of the velocity law parameter @xmath113 , obtained by fitting the calculated velocity structure from @xmath155 to @xmath156 to a @xmath113-law .
we see that the value of @xmath113 does not change much .
[ fg_vr ] shows the velocity structure of the wind and the location of the critical point for @xmath31 = 0 , 0.3 and 0.6 .
note that the critical point is very close the star .
this is true for all finite disk corrected solutions .
figure  [ fg_vmp ] shows the behaviour of @xmath6 versus @xmath135 for various rotational speeds and latitudes . as expected , @xmath6 at the equator is smaller than at the pole because @xmath6 scales roughly with @xmath27 and due to the smaller value of @xmath27 at the equator .
figure [ fg_mp ] shows the latitude dependence of the local mass - loss rate and of the terminal velocity from equator to pole , for different rotational speeds , indicated by @xmath31 . at the pole ( @xmath157 ) @xmath6 is about constant for various rotational speeds
the mass - loss rate at the equator decreases slightly with increasing rotation .
the increase of the mass loss rate at the equator as calculated by pauldrach et al .
( @xcite ) and friend & abbott ( @xcite ) is offset by a decrease in mass loss rate due to the smaller radiative flux at the equator due to the von zeipel effect .
the mass loss rate at the pole increases strongly with increasing rotation rate due to the increase in luminosity at the pole .
this is because the total luminosity of the star is fixed and the smaller luminosity at the equator must therefore be compensated by a higher luminosity at the pole .
the overall effect of rotation is to increase the mass loss rate of a star .
this can be seen in fig .
[ fg_mp ] where the overall mass - loss rate ( i.e. @xmath135 integrated over the stellar surface ) is plotted vs. the rotational velocity .
the mass - loss rate of the models with @xmath158 vary with @xmath31 as @xmath159 .
figure  [ fg_ldens ] shows the resulting density contrast ( far from the star ) between pole and equator ( @xmath160 ) .
the densities at pole and equator are essentially the same : the smaller @xmath6 at the equator is offset by the larger mass loss rate @xmath135 at the pole .
[ fg_dc ]    we have also calculated models for different luminosities ( see table 2 ) , and different stellar masses ( see table 3 ) . the mass - loss rate for the different luminosities is also plotted in fig .  [ fg_mp ] . the density contrast for various luminosities is plotted in fig .
[ fg_ldens ] .
we see that the density at the pole compared to the equator increases for increasing luminosity and the effect on the mass - loss rate is only visible for the highest luminosity .
this is an effect of the continuum acceleration @xmath54 .
7 illustrates the effect of the continuum acceleration on the forces in the wind for a star with high luminosity ( @xmath162 ) .
the net effect of the continuum radiation pressure and the gravity , i.e. @xmath163 , is an inward force in the equatorial region and an outward force in the polar region where the flux is higher .
for small luminosity the variation in @xmath164 through @xmath54 is not important because @xmath47 itself is small . for high luminosity @xmath165 increases and thus the effect of @xmath166 on @xmath164
becomes noticeable ( this is with all other properties of the star fixed ) . in eq .
[ eq_vinfsimple ] we see that a larger value of @xmath166 at the pole causes a relatively larger decrease ( through @xmath164 ) in @xmath27 and thus in @xmath6 , whereas the mass loss rate at the equator will increase less because of the smaller @xmath54 .
this is only a qualitative explanation and the actual dependence of @xmath167 does not follow the relation suggested by eqs . 29 and 31 .
we have found that the actual values of @xmath167 could be approximated quite well ( within 10 percent ) with the following relation :    @xmath168    where the @xmath169 is the radius where @xmath54 reaches its minimum value .
so we see that although rotation alone modifies the wind structure considerably , it can not be responsible for the density contrast between the equatorial and polar wind observed in b[e ] stars , in case only radial effects are considered .
we conclude that at least one other physical effect must be responsible for the formation of disks of rotating b[e ] stars .
the two most likely additional effects are the bi - stability and the flow of the wind material towards the equator .
in the section above we have shown that rotation alone can not explain the observed large density contrast between the poles and the equator of b[e ] stars .
one of the possible additional mechanisms to enhance this contrast is the bi - stability jump .    to investigate whether the bi - stability jump observed in normal b - supergiants can explain the formation of disks around b[e ] stars , we have calculated bi - stable wind models for a fast rotating typical b[e ] supergiant .
spherical , nlte wind models for normal b supergiants have been calculated by vink et al .
( 1999 ) . from these models
cak - like force multiplier parameters have been derived using a monte - carlo method to simulate photon - gas interactions .
these models show the occurence of a bi - stability jump around @xmath11 k.    for two models on either side of the bi - stability jump ( @xmath170 = 17 500 and 30 000 k ) , the cak parameters were determined by fitting a power law from about the sonic point to about @xmath171 ( from @xmath172 to @xmath173 in optical depth , @xmath75 is defined in eq .
[ eq_t ] ) to the calculated force multiplier values .
the resulting values for @xmath0 , @xmath1 and @xmath2 are listed in table  [ tb_ka ] .
note that the @xmath2 parameter is taken equal to zero as it was not possible to extract it explicitly from the models .
its effect is hidden in the constant @xmath0 . as expected , the values for the force multiplier parameters are quite different for the two cases , because the ionization has changed dramatically over the bi - stability jump .
cccc & @xmath0 & @xmath1 & @xmath2 + 17500 k & 0.57 & 0.45 & 0.0 + 30000 k & 0.06 & 0.65 & 0.0 +    [ tb_ka ]    the predicted bi - stability jump occurs around @xmath174 kk .
therefore we calculated a model for a rotating b[e ] star with the following properties : @xmath175 k , @xmath176 , @xmath177 and solar abundances .
this implies that the pole of the rapidly rotating star will be hotter than 25 kk and the equator will be cooler than 25 kk .    for the pole , the force multipliers @xmath0 and @xmath1 from the hot monte carlo model of 30 000 k were used , whereas for the equator , @xmath0 and @xmath1 from the cool model of 17
500 k were used .
the resulting @xmath178 and mass loss rates @xmath38 are plotted in figs .
[ fg_vtheta ] and [ fg_mdtheta ] for @xmath31 = 0.6 .
clearly visible is the drastic decrease in @xmath6 and the drastic increase of @xmath135 towards the equator .
this occurs around a co - latitude @xmath14 ( angle from the pole ) of about 60 degrees , where @xmath179 k. ( the precise location of the jump depends on the effective temperature of the star ) .
since both @xmath135 and @xmath6 have a equator / pole contrast of about a factor of three , the resulting density contrast is about a factor of 10 in a disklike region with a half opening angle of 30 degrees .
our calculations show that the density contrast in the wind between pole and equator due to the bi - stability jump is significantly larger than without the bi - stability jump . yet , the calculated value of about a factor 10 is not sufficient to explain the observed density contrast of a factor @xmath180 .
so another mechanism is needed to enhance the density contrast even further .
this is most likely the wind compression mechanism ( bjorkman &   cassinelli 1993 ; bjorkman 1998 ) .
we have modified the line driven wind theory for rotating stars by including the effects of oblateness and gravity darkening as well as rotational terms in the equation of motion .
we considered a sectorial wind , i.e. we neglected the effect of the motion of the wind towards the equator or towards the pole .
this assumption is justified close to the stellar surface below the critical point of the momentum equation .
this implies that our method will predict about the correct distibution of mass - loss rates from the star as a function of stellar latitude , but it may not be accurate enough to predict the velocity and density distribution further away from the star , if motions in the @xmath14 direction become important .
the equation of motion was solved using an iterative numerical scheme , that includes the conservation of angular momentum and the correction factors to the radiative acceleration by lines and by the continuum , due to the non - spherical shape of the star and due to the latitude dependence of the radiative flux .
this method was applied to study its possible effects on the formation of disks around fast rotating b[e ] supergiants .
the models with constant force multiplier parameters @xmath0 , @xmath1 and @xmath2 show a decrease of _ both _ the mass loss rate and the terminal velocity from pole to equator .
this is mainly due to two effects : the reduction of the escape velocity from pole to equator , resulting in a higher terminal velocity at the pole , and the reduction of the radiative flux from the pole to the equator due to gravity darkening , which results in a decrease in mass loss rate at the equator .    for a star with a fixed luminosity and a fixed polar radius ,
the temperature at the pole increases and the temperature at the equator decreases with increasing rotation rate .
the terminal velocity of the wind from the poles is almost independent of @xmath181 but at the equator @xmath6 decreases with increasing rotation rate . the mass - loss rate at the pole increases with increasing @xmath181 , due to the increase in @xmath170 at the poles , but the mass loss rate from the equator is almost independent of the rotation .
the combination of these effects alone produce a density contrast between the polar and the equatorial wind of a factor @xmath182 , except for fast rotating stars with high luminosity ( in which case @xmath183 ) .
our results agree quantitatively with those obtained by maeder ( @xcite ) for the latitudinal dependence of @xmath135 .
our predictions of the global mass loss rates are also in general agreement with those of maeder for models with luminosities far from the eddington limit . for high luminosities
we find a strong polar outflow , whereas maeder s models predict an enhanced equatorial outflow on the basis of the original cak scaling laws .
this difference is due to the fact that we have included an improved description of the radiation pressure in a rotating star .
the difference between the polar wind and the equatorial wind is strongly enhanced when the bi - stability of radiation driven winds between @xmath149 and 30  000 k is taken into account . in this case
the force multipliers @xmath0 and @xmath1 of a rapidly rotating star change drastically with stellar latitude if the pole is hotter than 25  000 k and the equator is cooler . in   [ s_bs ] we have applied the newly calculated force multipliers from vink et al .
( 1999 ) above and below the bi - stability jump to the models of rotating b[e ] supergiants .
one might argue that it is not allowed to apply the force multipliers of ( bi - stable ) spherical wind models to aspherical winds of rotating stars .
however , we have shown that a disk formed through bi - stability will be sufficiently thick ( typically about 30 degrees for a star with an effective temperature equal to about 25  000 ) to make these spherical models reasonable approximations for the conditions in the wind close to the star , where the mass loss rate is determined .
we find that rotationally induced bi - stability models of b[e ] stars reach a density contrast of about a factor 10 between the dense equatorial wind and the less dense polar wind .
this is less than the factor @xmath184 that is observed ( see e.g. bjorkman  @xcite ) .
the extra density increase is most likely due to the wind compression .
the gas that leaves the photosphere from a rotating wind , will follow an orbit in a tilted plane defined by the local rotation vector and the center of the star .
if the rotational velocity is large or the wind velocity is small , this orbit will cross the equatorial plane where the streamlines of the wind from different stellar latitudes cross .
the resulting shock will compress the gas into a thin outflowing disk with an opening angle of only a few degrees ( bjorkman & cassinelli 1993 ; @xmath5 chapt 11 ) .
we have not considered this wind compression in our model .
however , it is likely to occur in rotationally induced bi - stable winds , because the wind velocity in the equatorial plane is about a factor three smaller than from the poles , thus facillitating the wind compression .
owocki & cranmer ( 1994 ) have argued that the flow of the wind towards the equatorial plane , predicted in the wcd - theory , may be offset by a @xmath14-component of the line acceleration towards the polar regions .
this `` disk inhibition '' mechanism operates in wind models with constant force multipliers @xmath0 and @xmath1 .
however , it is not clear that this mechanism is sufficiently strong to overcome the wind compression in the rib - model .
this is because @xmath0 and @xmath1 change with latitude in the rib - model and there is a strong density gradient in the wind from the equator to the poles ( owocki et al .
1998 , puls et al . 1999 ) .
the combination of the rib and the wcd mechanism may offer the best possibility for explaining the disks of b[e]-supergiants . whether the compression is strong enough to explain the observed high density contrast between the polar and the equatorial wind , remains to be seen .
it can be calculated by combining the solutions of the wind momentum equation of rotating oblate stars with gravity darkening ( derived in this paper ) , and the calculation of the resulting trajectories of the wind .
the combination of the rotation induced bi - stability model and the wind compressed disk model is promising for explaining the disks of b[e ] stars because the rib - mechanism explains the increased mass loss and the small velocity from the equatorial regions and the wcd - mechanism explains the strong compression of the disk .
abbott d.c . , 1980 ,
apj 242 , 1183 abbott d.c . , 1982 ,
apj 259 , 282 bjorkman j.e . , 1998 , in _
b[e ] stars _ , eds a.m. hubert , c. jaschek kluwer : dordrecht , assl 233 , p 189 bjorkman j.e .
, cassinelli j.p . , 1993 , apj 409 , 429 cassinelli 1998 in _ b[e ] stars _ , eds a.m. hubert , c. jaschek kluwer : dordrecht , assl 233 , p 177 castor j.i .
, 1974 , mnras 169 , 279 castor j.i .
, abbott d.c . ,
klein r.i . , 1975 ,
apj 195 , 157 cranmer s.r . ,
owocki s.p . , 1994 , apj 440 , 308 feldmeier a. , 1999 in _ variable and non - spherical stellar winds in luminous hot stars _ , lecture notes in physics , eds .
b. wolf et al .
springer : heidelberg , p 285 .
friend d.b . ,
abbott d.c . , 1986 ,
apj 311 , 701 gayley k.g . , 1995 ,
apj 454 , 410 kudritzki r.p . , pauldrach a. , puls j. , abbott d.c . , 1989 ,
a&a 219 , 205 lamers h.j.g.l.m . , cassinelli j.p . , 1999 , introduction to stellar winds , cambridge univ . press .
( @xmath5 ) lamers h.j.g.l.m .
, pauldrach a. , 1991 , a&a 244 , l5 lamers h.j.g.l.m . ,
, lindholm d.m . , 1995 ,
apj 455 , 269 lamers h.j.g.l.m . ,
zickgraf f .- j .
, de winter d. , houziaux l. , zorec j. , 1998 , a&a 340 , 117 lucy l.b . , 1982 ,
apj 255 , 286 maeder a. , 1999 , a&a 347 , 185 owocki s.p .
, castor j.i . ,
rybicki g.b . , 1988 , apj 335 , 914 owocki s.p . , cranmer r.c . , 1994 ,
apj 424 , 887 owocki s.p . ,
cranmer r.c .
, gayley k.g . , 1998 , ap&ss 260 , 149 pauldrach a. , puls j. , kudritzki r.p . , 1986 , a&a 164 , 86 puls j. , petrenz p. , owocki s.p . , 1999 in _ variable and non - spherical stellar winds in luminous hot stars _ , lecture notes in physics , eds .
b. wolf et al .
springer : heidelberg , p 131 .
, de koter a. , lamers h.j.g.l.m . , 1999 ,
a&a 350 , 181 von zeipel h. , 1924 , mnras 84 , 665 zickgraf f .- j . , 1992 , in _
astronomical ccd observing and reduction techniques _
howell s.b . , asp conference series vol . 22 , 75 zickgraf f .- j . , schulte - ladbeck r.e . , 1989 , a&a 214 ,
274 zickgraf f .- j .
, wolf b. , stahl o. , leitherer c. , klare g. , 1985 , a&a 143 , 421


the 40 transiting planet systems confirmed as of august 2008 show a remarkable diversity of properties , which is indicative of the complexity of planet formation and evolution processes .
many different follow - up studies enabled by the special orientation of these systems ( e.g. , * ? ? ?
* ; * ? ? ?
* ; * ? ? ?
* ; * ? ? ?
* ; * ? ? ?
* see charbonneau et al .
2007 for a review ) have brought about rapid improvements in evolutionary models of planet interiors and atmospheres @xcite .
the increasing predictive power of these models is beginning to drive even more challenging observations . on the other hand , the precision and accuracy with which the most basic planet properties such as the mass and radius can be determined is currently limited by our knowledge of the properties of the host stars .
significant uncertainties remain in the stellar mass and radius determinations of many systems . in some cases
this is due to poorly determined photospheric properties ( mainly temperature and metallicity ) , and in others to a lack of an accurate luminosity estimate .
additionally , the variety of methodologies used for these determinations and the different approaches towards systematic errors have resulted in a rather inhomogeneous set of planet properties , as discussed by @xcite .
this complicates the interpretation of patterns and correlations that are being proposed ( e.g. , * ? ? ?
* ; * ? ? ?
* ; * ? ? ?
recent improvements in the analysis techniques have the potential to increase the accuracy of the stellar and planetary parameters significantly , especially for the large majority without a direct distance estimate .
in particular , the application of the constraint on the stellar density that comes directly from the transit light curves has been shown to be superior to the use of other indicators of luminosity such as the surface gravity ( @xmath15 ) determined spectroscopically @xcite . @xcite
have recently reanalyzed a large subset of the known transiting planets , incorporating these improvements and applying a uniform methodology to all systems .    in the present work we focus on two of the recently discovered transiting systems , @xcite and @xcite , which lack accurate estimates for the photospheric properties of the parent stars and as a result have more uncertain stellar and planetary parameters . to improve upon these properties , we present new radial - velocity and photometric observations of
with which we refine both the light curve solution and the spectroscopic orbit .
we also carry out a reanalysis of existing photometry utilizing a technique which treats stellar limb darkening using adjustable parameters .
we perform the first detailed spectroscopic determination of the photospheric properties of both stars , and we make use of the constraint on the stellar density mentioned above to infer more accurate values for the stellar and planetary masses and radii . our paper is organized as follows . in  [ sec : obs ] we summarize the observations . in  [ sec : atmpar ] we present our effective temperature and metallicity determinations , along with several consistency checks aimed at establishing the accuracy of the temperatures . in  [ sec : orb ] we report an updated spectroscopic orbital solution for , and the light curve solutions for both systems are discussed in  [ lcanalysis ] .
section[sec : physics ] then describes our determination of the stellar masses and radii , which in turn lead to refined values for the planetary parameters over those reported in the discovery papers .
we conclude in  [ sec : summ ] by providing a summary of our results and by discussing whether the properties of the host stars give any useful clues on the origin of the strongly contrasting densities of their close - in gas giant planets , particularly in comparison with the other known transiting planet systems .
high - resolution , high - snr ( signal - to - noise ratio ) spectroscopic observations of were obtained in july and september 2007 with the hires spectrograph on the keck  i telescope @xcite , using essentially the same setup as in the discovery paper . the spectra cover the effective wavelength range @xmath1632008800 , and an iodine ( i@xmath17 ) gas absorption cell placed in front of the spectrograph slit
was used to superimpose a dense set of narrow molecular lines between @xmath165000  and @xmath166000 .
this iodine spectrum provides a stable wavelength reference and a means of monitoring the instrumental profile that are crucial for achieving precise radial velocity determinations .
the observations were reduced and extracted using the makee software written by t.  barlow .
we obtained 4 spectra with typical exposure times of 15 min resulting in snrs of @xmath16100 pixel@xmath18 , which we combined with the discovery observations reported by @xcite for a total of 11 iodine spectra .
precise radial velocities were measured following a procedure based on the methodology developed for the afoe spectrograph  @xcite , and adapted for the processing of hires spectra .
conceptually , the method models the observed star and i@xmath17 spectrum by using templates for the stellar and i@xmath17 spectra .
this model includes physically motivated parameters that describe the spatial and temporal variations of the instrument , including the instrumental profile , as well as the sought - after relative doppler shift with respect to the stellar template .
the model parameters are adjusted to minimize the difference between the model and the observations in a least - squares sense , down to the poisson noise .
the stellar template is estimated by deconvolving an observation taken without superimposing i@xmath17 , whereas the template for the i@xmath17 is based on a high - resolution and high snr scan of the gas cell using the fourier transform spectrometer on the mcmath solar telescope at the kitt peak observatory  ( see * ? ? ?
each spectral order is modeled as a whole but analyzed independently .
the final radial velocities are estimated from the mean of the doppler shifts computed for each order and their uncertainties from the standard deviation of that mean .
the resulting velocities , expressed in the solar system barycentric frame , along with their associated formal errors are reported in table  [ rvtres3 ] .
these measurements include and supersede the velocities presented in the discovery paper , and they also correct a minor error in the previously published dates of observation .
the revised orbital solution is presented and discussed in  [ sec : orb ] .
in addition to the photometric measurements presented in the discovery paper @xcite , we have collected six other high - cadence , high - precision transit light - curves of .
one was obtained in the @xmath19 band on ut 2007 april 23 using the ccd camera of the iac80 telescope at the observatorio del teide , tenerife , spain , and the other five were gathered using keplercam ( see , e.g. , * ? ? ?
* ) at the flwo 1.2-m telescope : two in the sloan @xmath20 and @xmath21 bands on ut 2007 april 25 and 2008 march 27 , and the remaining three in the @xmath22 band on ut 2008 march 9 , april 12 , and may 8 . on all nights conditions
were quite good , except for the presence of high cirrus and highly variable seeing on ut 2008 april 12 .
the transits were observed at airmasses ranging between @xmath161.00 and @xmath161.80 .
all datasets were reduced using standard calibration techniques ( overscan correction , trimming , bias subtraction , flat fielding ) .
we then performed aperture photometry of and between 10 and 30 comparison stars , depending on the filter and exposure time .
we experimented with different choices for the aperture size , comparison star ensemble , and weighting of the comparison stars , aiming for the smallest out - of - transit ( oot ) rms . in practice ,
the best aperture size was approximately twice the fwhm of the stellar image on each night , and the best results were obtained using a straight average of the normalized light curves of the comparison stars . as for the formal errors on each photometric data point , we used the product of the oot rms of each light - curve and the factor @xmath23 , used to account for departures from gaussian ( `` white '' ) uncorrelated noise ( see , e.g. , * ? ? ? * ) .    for consistency ,
we also re - processed in the same way the two discovery light - curves of presented in odonovan et al .
( 2007 ) , as well as the three photometric datasets utilized by mandushev et al .
( 2007 ) in their discovery announcement of .
the main characteristics of the and light - curves are summarized in table  [ tbl : lctres3 ] .
the final set of photometric time - series of in all filters ( including the discovery data ) is available in machine - readable form in the electronic version of table  [ tbl : phottres3 ] , and is plotted in figure  [ tres3phot ] .
table  [ tbl : phottres4 ] reports all photometric data for , which are shown graphically in figure  [ tres4phot ] .
a detailed abundance analysis was carried out using the keck / hires template spectra of and . the stellar atmosphere parameters
@xmath24 , @xmath15 , and [ fe / h ] were determined using standard methodology ( e.g. , * ? ? ?
* ; * ? ? ?
* ; * ? ? ?
* ) , which we summarize as follows .
equivalent widths ( ews ) for a set of relatively weak and lines were measured manually in the keck spectra using the splot task in iraf  ( see , e.g. , * ? ? ? * and references therein , for details on the specific choice of lines ) .
these were then used together with a grid of kurucz atlas plane - parallel stellar model atmospheres @xcite as inputs to the 2002 version of the moog spectral synthesis code @xcite .
atmospheric parameters were derived under the assumption of local thermodynamic equilibrium ( lte ) , imposing excitation and ionization balance .
formal uncertainties on @xmath24 , @xmath15 , and microturbulent velocity @xmath25 were derived using the approach described in @xcite and @xcite , while the nominal uncertainty for [ fe / h ] corresponds to the scatter obtained from the lines rather than the formal error of the mean .
the resulting set of parameters for is @xmath26  k , @xmath27 , @xmath28 @xmath29 , and [ fe / h ] = @xmath30 , and for we obtain @xmath31  k , @xmath32 , @xmath33 @xmath29 , and [ fe / h ] = @xmath34 . these values are collected in table  [ star_tres3 ] and table  [ star_tres4 ] , respectively .    as a check ,
the manual ew measurements used above were compared against the results from the automated software aressousasag / ares .
] , made available to the community by @xcite .
figure  [ aresiraf ] shows the comparison between the ews measured for using ares and those from the manual approach with iraf .
the top panel indicates excellent agreement between the two ; the linear fit has a slope of 1.011 .
there are no significant correlations with the line strength ( see middle panel ) : the fractional difference between ews measured by the two methods has a probability of no correlation of 0.1 .
finally , the histogram in the bottom panel shows there is no appreciable systematic difference between the methods .
the mean difference and scatter are only 0.67 m and 2.6 m , respectively , which are smaller than found by @xcite from a similar comparison between feros , harps and uves spectra .
consequently , the results of the abundance analysis with moog using ews measured with ares are virtually identical to those presented earlier .    following @xcite we also synthesized a number of unblended lines in the template spectra of both and , and determined projected rotational velocities of @xmath35 @xmath29  and @xmath36 @xmath29 , respectively .    our new determinations of the atmospheric parameters for the parent stars of and are generally in good agreement with the values reported by @xcite and @xcite , the main difference being that those authors assumed [ fe / h ] = 0.0 in their studies , whereas we find significant departures from solar metallicity in both stars .
we describe here additional estimates of the effective temperature for and that serve to test the accuracy of our determinations above .
both planet host stars were observed spectroscopically with the center for astrophysics ( cfa ) digital speedometer @xcite as part of the regular follow - up after the discovery of photometric signals suggesting transits .
these observations cover 45    in a single echelle order centered at 5187   , and have @xmath38 .
while this resolving power is moderately high , the snrs of the spectra are low , ranging from 7 to 13 per resolution element of 8.5  @xmath29  for the 13 exposures of , and snrs of 1113 for the 7 exposures of .
nevertheless , useful information on the stellar properties can be extracted from these spectra as described by @xcite .
briefly , the observed spectra are cross - correlated against a library of synthetic spectra based on kurucz model atmospheres ( see * ? ? ?
* ; * ? ? ?
* ) , calculated over a wide range of values of @xmath37 , @xmath15 , [ fe / h ] , and rotational velocity @xmath39 . the combination of parameters yielding the highest correlation averaged over all exposures is adopted as a representation of the properties of the star . due to the narrow wavelength range and limited snrs ,
the first three of the above properties are typically strongly correlated and are difficult to determine simultaneously . given these constraints , in the discovery papers we initially held the metallicity fixed at the solar value for both stars .
for , which our keck spectroscopy indicates is slightly metal - poor , we repeated the above determination using a fixed metallicity of [ fe / h ] @xmath40 , and interpolated @xmath37 and @xmath15 to the precise composition from table  [ star_tres3 ] .
we obtained @xmath41  k and @xmath42 , where the uncertainties include contributions both from the error in the adopted [ fe / h ] and the internal errors .
the temperature as well as the surface gravity are consistent with the keck determinations , within the uncertainties .
similarly for we repeated the determinations for [ fe / h ] @xmath43 , given the metal - rich composition indicated by the keck spectroscopy , and interpolated to the intermediate value from table  [ star_tres4 ] .
the results are @xmath44  k and @xmath45 , which are once again in good agreement with the more reliable estimate from keck .
while the temperature determinations discussed so far are implicitly based on the _ strength _ of the spectral lines , @xcite have demonstrated that highly precise information can also be extracted using the _ ratio _ of the depths of two spectral lines having different sensitivity to temperature .
however , rather than yielding absolute temperatures , this technique in its original formulation only measures changes in temperature , albeit with extremely high precision often reaching a few kelvin ( see , e.g. , * ? ? ?
* ; * ? ? ?
* ; * ? ? ?
absolute temperatures can still be obtained with recourse to external color - temperature calibrations , since line - depth ratios ( ldrs ) are usually strongly correlated with the color index of the star .
the accuracy of such @xmath37 determinations is then limited by the calibrations themselves .
@xcite have presented ldr-@xmath37 calibrations based on 26 carefully selected lines of fe , v , sc , si , and ni in the spectral interval 61906280   , grouped into 16 line pairs .
these calibrations are valid for stars with temperatures between @xmath163800  k and @xmath166000  k , so unfortunately the technique is not applicable to . for we measured all 26 lines in our high - snr keck template spectrum , and adopted the @xcite calibrations appropriate for stars with rotational velocities of 10  @xmath29  ( very close to the @xmath46 value reported by @xcite and obtained also in this work ) .
the temperature scale of these relations relies on a transformation between @xmath47 and @xmath37 by @xcite , which makes use of a mixture of dwarf and giant temperatures obtained by many different methods and does not account for differences in metallicity .
for the present work we have preferred to use more sophisticated color - temperature relations such as those by @xcite and @xcite , which account not only for luminosity class but include also metallicity terms , and are based on effective temperatures derived homogeneously by the infrared flux method .
the conversion from the ldr - based @xmath37 inferred from the @xcite relations back to an average color for the star using the @xcite prescription gives @xmath48 .
our two preferred color - temperature relations mentioned above then yield a weighted average temperature of @xmath49  k for , in which the uncertainty includes observational errors propagated from the ldr measurements as well as the scatter of the calibrations .
this result is consistent with our more direct estimate in  [ sec : keck ] .
as is well known , the wings of the @xmath51 line ( but not its core , formed higher up in the atmosphere under non - lte conditions ) are very sensitive to @xmath24 variations , but are relatively insensitive to changes in @xmath15 and [ fe / h ] ( see , e.g. , * ? ? ?
* ; * ? ? ?
* and references therein ) .
this allows for a useful consistency check on our @xmath24 estimates above . we compared the observed h@xmath50 line profiles in our keck template spectra of and against synthetic profiles for solar - metallicity dwarfs ( [ fe / h ] = 0.0 , @xmath52 ) from the kurucz database .
the results of this exercise are displayed in figure  [ halpha ] , in which 10    regions centered on h@xmath50 are shown for each star together with four calculated profiles for different values of @xmath24 . in both cases the temperatures one would infer from these comparisons agree well with the estimates reported in table  [ star_tres3 ] and table  [ star_tres4 ] .
an additional check on the effective temperatures is available from the multi - color photometry for and .
measurements in johnson @xmath53 , cousins @xmath54 , and 2mass @xmath55 were used to derive seven color indices for , and nine for when considering also the @xmath56 and @xmath57 measurements from the tycho-2 catalog @xcite .
the calibrations by @xcite and @xcite then yielded weighted average temperatures of @xmath58  k and @xmath59  k for and , respectively , ignoring extinction and adopting [ fe / h ] in each case as determined in  [ sec : keck ] .
the uncertainties include photometric errors and metallicity errors , as well as the scatter of the calibrations , but exclude unquantified systematic errors in the calibrations themselves .
these temperatures are 260  k and 65  k cooler than our spectroscopic determinations above .
they can be reconciled with the values in  [ sec : keck ] if we assume the presence of reddening , and correct each of the indices . for the required value of @xmath60 is hardly significant ( @xmath61 mag ) , but for we obtain @xmath62 mag . ,
in which the ldr - based @xmath47 color is significantly bluer than the measured value of @xmath63 @xcite . ] for comparison , the reddening maps of @xcite and @xcite indicate a _
total _ reddening along the line of sight to each star of only @xmath64 mag , and similar results are obtained from the model of galactic dust distribution by  @xcite .
this is not inconsistent with the small value we infer for , but it is much smaller than our estimate for .
possible explanations include patchy interstellar material combined with the relatively coarse resolution of these maps ( a few arc minutes for * ? ? ?
* ) , or perhaps the presence of circumstellar material in .
the reliability of age indicators for stars older than 12  gyr , such as chromospheric activity and lithium ( li ) abundance , as well as their interagreement , have been the subject of much debate in the literature ( e.g. , * ? ? ?
* ; * ? ? ?
* ; * ? ? ?
* ; * ? ? ?
* ; * ? ? ?
the difficulties are due to many factors , including ( but not limited to ) non - trivial correlations between chromospheric activity , rotation , mass , and age , limited availability of activity estimates averaged over entire stellar activity cycles  that range from @xmath65 to @xmath66 .
these correspond to ages of @xmath67  gyr and @xmath16 2.5  gyr , respectively ( see , e.g. , * ? ? ?
] , insufficient understanding of the temporal evolution of li depletion due to the complex interplay between various processes ( e.g. , convection , mixing , diffusion , mass loss ) , and non - negligible differences in the observed behavior of chromospheric activity and li depletion as a function of mass , age , and chemical abundance between stars in young and old clusters , and in the field .
we have nonetheless attempted to use the activity indicator and the lithium abundance as measured in our hires spectra of and to provide independent constraints for comparison with the formal age estimates determined below from evolutionary models ( see  [ sec : physics ] ) , as well as to search for possible peculiarities of these two planet hosts compared to other samples of stars with and without planets .
figure  [ ca_tres3_4 ] shows a region of the hires template spectra for each star centered on the h line .
clear emission is seen in the core of the line of , but not in . following the procedure outlined in @xcite we measured the chromospheric emission ratio @xmath68 , corrected for the photospheric contribution , from the h and k lines in each of our spectra .
we obtained @xmath69 and @xmath70 for and , respectively .
thus appears moderately active , while is quite the opposite .
the resulting chromospheric age estimates , using the relations summarized in @xcite , are @xmath71 gyr for and @xmath72 gyr for .    in figure  [ licomp ]
we show the results of the spectral synthesis of a 10    region centered on the li  @xmath736707.8 line , using the atmospheric parameters derived from the fe - line analysis and the line list of @xcite .
the two panels display the comparison between the observed spectra and three synthetic spectra , each differing only in the assumed li abundance .
neither star shows a measurable li line , and we can only place upper limits of @xmath74 and @xmath75 for and , respectively . by comparison with average li abundance curves as a function of effective temperature for clusters of different ages @xcite , one would infer a rather old age for of @xmath76 gyr , nominally inconsistent with the estimate above from the activity index . based on its li the object appears decidedly older than hyades stars of the same @xmath24 , and this argument
is corroborated ( e.g. , * ? ? ?
* ) by the small @xmath46 we measure ( see table  [ star_tres3 ] ) . for
, the age inferred from the @xcite relations is @xmath77 gyr , which is in qualitative agreement with the absence of significant chromospheric activity as well as with the inferred value of the surface gravity ( see table  [ star_tres4 ] ) , which indicates the star is somewhat evolved .
the measured li abundance for is in line with values determined by @xcite for a sub - sample of nearby planet hosts with the same @xmath24 , while appears more depleted than other planet hosts of similar @xmath24 , @xmath78 , and [ fe / h ] .
we discuss briefly in  [ sec : summ ] some implications of these findings in the broader context of the existence of chemical peculiarities in planet hosts in comparison with stars without detected planets .    from the measurements in this section we conclude that the above results are consistent with the notion that neither star is very young ( @xmath79 12 gyr ) , that is likely to be of intermediate age ( @xmath163  gyr ) , and that appears older and more evolved .
using the radial velocities presented in  [ sec : rvs ] we have updated the orbital solution given by @xcite , adopting the improved ephemeris described below in  [ ttv ] .
a keplerian circular orbit was adjusted to the data , and the results may be seen in table  [ tab : par_tres3 ] and are shown graphically in figure  [ specorb_tres3_4 ] . as found also by @xcite , the solution gives a scatter that is larger than expected from the internal velocity errors .
this is most likely due to velocity `` jitter '' associated with chromospheric activity . if we model this as excess scatter to be added quadratically to the internal errors , we find that the amount of jitter required to produce a reduced @xmath80 value near unity is @xmath1618.5  @xmath81 .
external estimates of the jitter for can be made on the basis of the spectral type , the measured value of @xmath82  @xmath29 , and the activity index @xmath83 , and vary considerably but generally range from about 5 to 20  @xmath81@xcite .
these estimates are thus consistent with our findings . nevertheless , as a test we also modeled the data with an eccentric orbit but did not obtain much improvement , and the result for the eccentricity was not significantly different from zero ( @xmath84 ) .
there are no indications of long - term variations in the observations at hand .    for we adopt in the following the spectroscopic orbit by @xcite , since the radial - velocity material has not changed . for clarity ,
the orbital solution is reported again in table  [ tab : par_tres4 ] .
the analysis of the differential photometry for and was carried out in essentially the same manner as described in @xcite .
we refer the reader to that paper for the details .
given that we have multiple light curves for each system , we first determined the individual times of transit as described below in order to improve the ephemeris , and we subsequently refined the light - curve parameters .      to determine the center of each measured transit event we initially adopted the light - curve parameters from the discovery papers .
for , we adopted a quadratic limb darkening ( ld ) law with coefficients from @xcite , interpolated to the values @xmath85 k , @xmath86 dex , [ m / h ] @xmath87 dex , and @xmath88 @xmath29 .
all ld coefficients for the various filters are reported in table  [ tbl : ldcoeff ] .
we fitted each light curve individually to solve only for the time of transit center @xmath89 and also a linear function of time ( two parameters ) to describe the oot flux ( the slope in the oot flux accounts for systematic errors , including differential extinction ) .
the measured mid - transit times are presented in table  [ midtrans ] .
we then fitted a straight line to the central times of the eight transits of the form @xmath90 and derived the following new ephemeris ( which we report in table  [ tab : par_tres3 ] ) : @xmath91 ( hjd ) , @xmath92 days . with a number of degrees of freedom @xmath93 ,
the resulting reduced @xmath94 indicates a rather poor fit .
the transit timing residuals for are shown in figure  [ trestiming ] .
there are four outliers at the 23.5@xmath95 level .
this could be seen as evidence suggesting that the period is not constant .
alternatively , the errors might have been underestimated .
the data available are not enough to draw any significant conclusion on the nature of these variations , but clearly additional observations are warranted .
we leave for a future study the evaluation of the dynamical interpretation and significance of the transit times , as was done , for example , by @xcite for tres-1 and by @xcite for ogle - tr-111b .    similarly for , we re - fitted the three discovery light - curves using the @xcite ld coefficients listed in table  [ tbl : ldcoeff ] appropriate for @xmath96 k , @xmath97 dex , [ m / h]@xmath98 dex , and @xmath25 = 2.0 @xmath29 .
the corresponding times of transit center , derived with the same procedure adopted for , are presented in table  [ midtrans ] . given the short time baseline of these observations , we did not derive a new ephemeris for the system , but for the purpose of the analysis presented in the next section we simply adopted the mandushev et al .
( 2007 ) values of @xmath99 and @xmath89 .
next , for both and we locked the transit times and oot baseline functions at the values indicated above and then , under the assumption of a circular orbit , we fitted all light curves simultaneously using the algorithm by @xcite to derive the relevant quantities radius ratio @xmath100 ( where @xmath101 and @xmath102 are the planetary and stellar radius , respectively ) , inclination @xmath22 , projected separation @xmath6 ( where @xmath103 is the semi - major axis ) , and impact parameter @xmath104 . for we found @xmath105 , @xmath106 , @xmath107 , and @xmath108 . for
, we obtained @xmath109 , @xmath110 , @xmath111 , and @xmath112 .
these estimates were derived assuming initial values for the stellar and planetary masses of @xmath113 @xmath114 and @xmath115 @xmath116 for , and @xmath117 @xmath114 and @xmath118 @xmath116 for , although the results are insensitive to these values .
recent studies of transiting exoplanets light curves ( e.g. , * ? ? ?
* and references therein ) have highlighted some issues related to how much effect different treatments of limb darkening can have on the light - curve solutions .
in particular , these do not appear to be significantly affected by the specific choice of the ld law , while fixing the ld coefficients to their theoretical values seems to result in significantly smaller uncertainties in the fitted parameters with respect to the case in which the ld coefficients are also adjusted during the fitting procedure .
the above results were obtained keeping all ld coefficients fixed .
we decided to investigate the effect of fitting for the ld coefficients in the following way .
given the quality of the data available to us , we can only fit for one coefficient ( or one combination of the two coefficients ) .
we chose to fix the quadratic coefficient @xmath119 and solve for the linear coefficient @xmath120 .
we found that allowing complete freedom in the linear coefficient resulted in `` unphysical '' solutions ( e.g. , largely negative coefficients ) .
@xcite found that the fitted ld coefficients are usually within 0.10.2 of the theoretical claret values , and that is consistent with our experience ( see , e.g. , * ? ? ? * ) .
we then decided to use an a priori constraint enforcing an agreement of @xmath121 , modifying the merit function as follows : @xmath122 ^ 2+\left(\frac{u - u_1}{0.2}\right)^2\ , , \ ] ] where @xmath123 is the stellar flux observed at time @xmath124 , @xmath125 its corresponding error , @xmath126 is the model value , @xmath127 is the adjustable linear ld coefficient , and @xmath120 its theoretical value ( appropriate for each band - pass ) .
furthermore , for all three @xmath22-band light curves were required to agree on the ld parameter , and the same requirement was set on the two @xmath128-band light curves for .
the results for the system parameters in this case were @xmath129 , @xmath130 , and @xmath131 . for
, we obtained @xmath132 , @xmath133 , and @xmath134 .    in both cases ,
the agreement between the system parameters derived keeping the ld coefficients fixed and those when the ld coefficients are part of the solution was excellent , with only a slight increase in the estimated uncertainties for the latter case . indeed , as already noted by @xcite , fixing the ld coefficients at their theoretically predicted values does not appear to significantly bias the results .
however , in the interest of providing more conservative error estimates , we believe that , for the purpose of the analysis of high - quality light - curves such as the ones presented in this paper , a procedure that treats ld coefficients as adjustable parameters is preferable . based on the above considerations , for both and we elected to include in table  [ tab : par_tres3 ] and table  [ tab : par_tres4 ] the values of the systems parameters obtained from the light - curve analysis in the case in which the linear ld coefficient was allowed to float .
we note that this approach provides a departure from the @xcite analysis .
the revised spectroscopic determinations of @xmath37 and [ fe / h ] , along with the new spectroscopic orbital solution for and the light curve fits presented above , allow us to refine the determination of the stellar and planetary properties for both systems . to establish the properties of the parent stars we rely on stellar evolution models by @xcite and @xcite .
we make explicit use of the constraint on the stellar density provided by the light - curve quantity @xmath135 , as described by @xcite .
the procedure follows closely that given in the previous citation .
briefly , we seek the best match ( in a @xmath80 sense ) between the measured \{@xmath37 , [ fe / h ] , @xmath135 } and points on a finely interpolated grid of isochrones spanning a wide range of metallicities and ages .
two minor improvements in this procedure , described in more detail by @xcite , have to do with the weighting of each point sampled along the isochrones according to the distance in @xmath37-[fe / h]-@xmath135 space compared to the observed values , and an additional weighting according to the likelihood that the star is in a particular evolutionary state .
the latter effect is accounted for by multiplying the first weight by the expected number density of stars at each location in the h - r diagram , according to an adopted initial mass function ( which in this case is simply a power law with a salpeter index ) .
the impact of these weighting factors is generally minor .
the results for the stellar properties of and are presented in table[star_tres3 ] and table[star_tres4 ] . the distance estimate to accounts for interstellar extinction ( @xmath136 mag ) , as described in
[ sec : phot_teff ] .
for both stars the surface gravity inferred from the models is typically much more accurate than the spectroscopic gravity determination , showing the power of the constraint on the luminosity and size of the star afforded by @xmath135 .
the results for both and supersede those given recently by @xcite because of the new photometry and radial velocities contributed here in the former case , and the new light - curve solutions for both systems .
the evolutionary ages for the two stars are qualitatively consistent with our conclusions from  [ sec : age ] , but are somewhat more uncertain in the case of . finally , using the distance estimates inferred in this work , the ucac2 proper motion components  @xcite reported by  @xcite and  @xcite , and the mean radial velocity values @xmath137
@xmath29 and @xmath138  @xmath29 as measured with the cfa digital speedometers , we obtain galactic space motion vectors [ @xmath139 , @xmath19 , @xmath140 = [ + 27.3 , + 6.7 , + 33.0 ]  @xmath29 and [ @xmath139 , @xmath19 , @xmath140 = [ @xmath14143.9 , @xmath14139.1 , @xmath1416.9 ]  @xmath29 for and , respectively ( where @xmath139 is taken to be positive toward the galactic anticenter ) .
we collect these results along with the other properties derived previously in table  [ star_tres3 ] and table  [ star_tres4 ] .
the planet parameters follow from the stellar properties and the results of the transit light curve and spectroscopic orbits , and are presented in table[tab : par_tres3 ] and table[tab : par_tres4 ] . for , our planetary mass and radius are significantly larger compared with the determinations in the discovery paper and in @xcite .
this is due in part to the increased mass and radius for the parent star , but also to the larger radius ratio based on the new photometry .
our detailed spectroscopic analyses of and have yielded accurate values of the atmospheric properties ( @xmath24 and [ fe / h ] ) , which are critical for establishing the fundamental properties of the hosts . the accuracy of the temperatures is supported by a number of independent checks ( low - resolution spectroscopy , line - depth ratios , h@xmath50 line profiles , color - temperature calibrations ) that gives us confidence that the inferred stellar properties are reliable .
we find that is a main - sequence g dwarf with a metallicity about 1.5 times lower than the sun s , not a very common occurence among exoplanets hosts .
is a somewhat evolved late f star that is nearing the end of its main - sequence phase , and is slightly enhanced in its iron content with respect to the solar abundance .
the agreement we find between age indicators for and based on measurements of the activity levels , the lithium abundances , and rotation , and the evolutionary age inferred from the models is fair , although , as discussed in
[ sec : age ] , the reliability of empirical age estimates for stars that are not young ( @xmath142 gyr ) is somewhat questionable .
the model estimates themselves are not without their problems .
nevertheless , we conclude the stars are 13 gyr old .
neither star stands out as peculiar when compared with other planet hosts with similar physical properties . in the small @xmath46 value and the fact that we can only place an upper limit on the li abundance
are consistent with the notion that planet hosts with @xmath24 similar to the sun appear to rotate more slowly and are more li - depleted than stars without detected planets . as pointed out recently ( * ? ? ?
* and references therein ) , this evidence suggests that a planet - forming disk may induce additional rotational braking , leading to enhanced mixing in the stellar envelope , which in turn accelerates the destruction of lithium .
, on the other hand , does not seem consistent with the claim by @xcite that hotter planet hosts with @xmath143  k have higher li abundances , possibly due to self - enrichment processes .
it is worth keeping in mind , however , that these discrepancies may not be significant given the large spread in the li abundance for field stars with the temperature and mass of ( e.g. , * ? ? ?
new investigations on these issues are clearly needed based on uniform analyses of large samples of planet hosts and statistically significant , well - defined control samples of stars without detected planets .
new radial - velocity measurements for presented here have enabled us to revise the spectroscopic orbit for that system .
we detect no indication of any longer - term variations in the radial velocities that might suggest the presence of another body in the system .
however , the small number of observations and their limited time span of only 6 months emphasize the need for continued doppler monitoring of this and other transiting planet systems to investigate the possibility of additional companions . in stark contrast to the considerable ground - based and space - based efforts invested in studying in great detail the atmospheric properties of many of these objects , which have undoubtedly led to tremendous insights into their structure , formation , and evolution , the amount of radial - velocity data available for transiting planets is meager , and often does not go beyond the handful of observations published in the discovery papers .
interest in the radial velocities seems to be quickly lost .
it should be pointed out that the frequency of close - in giant planets ( @xmath144 day ) with additional massive planets in outer orbits ( up to the detection limit of today s doppler surveys , @xmath145 au ; see , e.g. , * ? ? ?
* and references therein ) is about 12% ( 8 out of 70 systems discovered via rv methods ) , and it would therefore be wise to continue the velocity monitoring of some of the transiting systems .
if additional planets in a transiting system were detected , they might also be found to undergo transits , and such a discovery would allow us to constrain structural models for gas giants akin to jupiter or saturn , and open exciting opportunities for additional investigations with present ( _ spitzer _ ) and upcoming ( jwst ) space - borne observatories .    accurate stellar properties for and have been derived here following the approach described in @xcite and @xcite , comparing the spectroscopically determined @xmath24 and [ fe / h ] values along with the photometrically measured @xmath6 with current stellar evolution models .
these properties have in turn allowed us to refine the determination of the mass and radius of the planets .
in particular , we confirm that is the planet with the largest radius among the currently known transiting hot jupiters .
recently , winn et al . ( 2008a ) derived upper limits on the albedo of tres-3 based on the nondetection of occultations observed at optical wavelengths .
our findings are relevant to this study in two ways .
firstly , they strengthen the case for a circular orbit , which is important because if the orbit is eccentric then it would be possible that winn et al .
( 2008a ) did not observe tres-3 at the actual times of occultations and that their data place no constraint on the albedo .
secondly , our revised light - curve parameters are relevant because the upper limit on the geometric albedo ( @xmath146 ) was inferred from the measured upper limit on the planet - to - star flux ratio ( @xmath147 ) according to      our revised value of @xmath149 therefore leads to revised upper limits on the geometric albedo of tres-3 .
however , this revision turns out to be minor : the new value of @xmath149 is only 1.5% larger than the value used by winn et al .
( 2008a ) .
the upper limits on the geometric albedo become weaker by about 3% . at 99% confidence ,
the revised upper limits are 0.31 , 0.64 , and 1.10 in @xmath150 , @xmath151 , and @xmath152band , respectively .
there is a large spread in the observed radii and densities for transiting planets of comparable mass placed at similar orbital distance from stars of very similar properties ( @xmath24 , @xmath15 , [ fe / h ] , and age ) .
for example , if we consider along with 5 other transiting planet hosts ( excluding ) with similar characteristics ( hd  149026 , hd  209458 , ogle - tr-56 , ogle - tr-132 , and wasp-1 ) , the nominal masses of the attending planets vary by a factor of @xmath153 , but reported densities vary by a factor of @xmath154 ( e.g. , * ? ? ?
many theoretical mechanisms have been proposed to inflate the radius of a strongly irradiated planet , such as additional sources of internal heating due to stellar insolation @xcite , tidal heating due to non - zero eccentricity caused by gravitational interaction with an outer companion @xcite or by rotational obliquity @xcite , elevated interior opacity due to enhanced atmospheric metallicity @xcite , or varying core masses @xcite .
none of these appear capable of explaining the observed spread in density in a natural way ( see also * ? ? ?
* ) .    of the 6 stars
just mentioned , all are more metal - rich than the sun , except for hd  209458 , which has [ fe / h ] close to solar .
some of the above models @xcite predict a positive correlation between the inferred planetary core mass and the host star s metallicity , as in the framework of the core - accretion model of giant planet formation ( e.g. , * ? ? ?
this idea assumes [ fe / h ] closely tracks the metallicity of the protoplanetary disk , so that more metal - rich stars should be orbited by more metal - rich planets , with a larger heavy - element content .
however , among these 6 systems only one planet ( hd  149026b ) has an inferred core mass significantly larger than zero , and four of the other planets have measured radii so large that the cores are likely to be insignificant .
in fact , even in the absence of a core the observed radii can not be reproduced by the models , with being the extreme case .
interestingly , over 40% of the transiting planets reported in table  5 of @xcite do not appear to require any core at all to explain their sizes ( according to the models of * ? ? ? * ) .
taking all this into consideration , the claimed evidence for a core mass  metallicity correlation could indeed be seen as supporting the more widely accepted scenario of formation by core accretion , but from the indications above it may also be that a significant fraction of these objects formed in a different way ( e.g. , * ? ? ?
* and references therein ) .    given the evidence collected so far , we suggest that simply connecting the host star s characteristics to the structural properties of transiting planets may in fact be an over - simplification .
we conclude by stressing the importance of refining our understanding of the complex interplay between the disk environment and a forming giant planet , and its evolutionary history after envelope accretion , which might turn out to be more directly responsible for its final structure and composition than the metal content of the parent star .    as gratefully acknowledges the kepler mission for partial support under nasa cooperative agreement ncc 2 - 1390 .
gt acknowledges partial support for this work from nasa origins grant nng04lg89 g .
dc is supported in part by nasa origins grant nng05gj29 g .
ftod acknowledges partial support for this work provided through the nasa postdoctoral program at the goddard space flight center , administered by oak ridge associated universities through a contract with nasa .
jbl gratefully acknowledges support from nsf grant ast-0307340 .
some of the data presented herein were obtained at the w.  m.  keck observatory , which is operated as a scientific partnership among the california institute of technology , the university of california and the national aeronautics and space administration .
the observatory was made possible by the generous financial support of the w.  m.  keck foundation .
the authors wish to recognize and acknowledge the very significant cultural role and reverence that the summit of mauna kea has always had within the indigenous hawaiian community .
we are most fortunate to have the opportunity to conduct observations from this mountain .
this research has made use of nasa s astrophysics data system abstract service and of the simbad database , operated at cds , strasbourg , france .
durisen , r. h. , boss , a. p. , mayer , l. , nelson , a. f. , quinn , t. , & rice , w. k. m. 2007 , in protostars and planets v , b. reipurth , d. jewitt , and k. keil ( eds . ) , university of arizona press ( tucson ) , 607                                                                                lcc 54187.04136 & 176.0 & 8.2 + 54187.13832 & 99.4 & 10.7 + 54188.02349 & 189.5 & 11.2 + 54188.12067 & 292.4 & 13.3 + 54189.01696 & @xmath141331.6 & 5.5 + 54189.09680 & @xmath141246.4 & 11.4 + 54189.14038 & @xmath141195.8 & 7.3 + 54288.85394 & 241.6 & 7.8 + 54288.98956 & 116.7 & 11.5 + 54289.82337 & @xmath1415.9 & 10.8 + 54372.86102 & @xmath141335.9 & 9.5 +       + 2007 - 03 - 25 & @xmath128 & flwo & 1.73 & 0.0012 & 0.0011 & 1.17 + 2007 - 04 - 08 & @xmath155 & ftn & 1.15 & 0.0011 & 0.0014 & 1.46 + 2007 - 04 - 24 & @xmath19 & ot & 1.34 & 0.0016 & 0.0016 & 1.30 + 2007 - 04 - 25 & @xmath20 & flwo & 0.98 & 0.0014 & 0.0015 & 1.16 + 2008 - 03 - 09 & @xmath22 & flwo & 0.73 & 0.0018 & 0.0016 & 1.00 + 2008 - 03 - 27 & @xmath21 & flwo & 0.73 & 0.0014 & 0.0014 & 1.07 + 2008 - 04 - 12 & @xmath22 & flwo & 1.05 & 0.0015 & 0.0015 & 1.02 + 2008 - 05 - 08 & @xmath22 & flwo & 0.73 & 0.0021 & 0.0020 & 1.16 +   + 2007 - 05 - 03 & @xmath128 & flwo & 0.73 & 0.0015 & 0.0016 & 1.26 + 2007 - 05 - 10 & @xmath128 & flwo & 0.73 & 0.0019 & 0.0017 & 1.23 + 2007 - 05 - 10 & @xmath155 & lowell & 1.54 & 0.0015 & 0.0015 & 1.22 +       + 2454185.850884 & 0.99965 & 0.00115 + 2454185.853303 & 0.99948 & 0.00115 + 2454185.854495 & 0.99852 & 0.00115 +   + 2454198.948010 & 0.99987 & 0.00109 + 2454198.948809 & 1.00054 & 0.00109 + 2454198.949606 & 0.99878 & 0.00109 +   + 2454214.574418 & 1.00159 & 0.00157 + 2454214.575348 & 1.00247 & 0.00157 + 2454214.576280 & 0.99766 & 0.00157 +   + 2454215.850393 & 0.99951 & 0.00150 + 2454215.851076 & 0.99912 & 0.00150 + 2454215.851771 & 0.99842 & 0.00150 +   + 2454535.897759 & 0.99906 & 0.00163 + 2454535.898268 & 0.99996 & 0.00163 + 2454535.898766 & 0.99755 & 0.00163 +   + 2454552.866642 & 0.99987 & 0.00135 + 2454552.867140 & 1.00181 & 0.00135 + 2454552.867638 & 1.00086 & 0.00135 +   + 2454569.875280 & 1.00116 & 0.00149 + 2454569.875778 & 0.99975 & 0.00149 + 2454569.876287 & 1.00090 & 0.00149 +   + 2454594.713269 & 1.00015 & 0.00204 + 2454594.713778 & 0.99628 & 0.00204 + 2454594.714276 & 1.00394 & 0.00204 +       + 2454223.741853 & 0.99446 & 0.00160 + 2454223.745823 & 0.99218 & 0.00160 + 2454223.746321 & 0.99209 & 0.00160 +   + 2454230.705029 & 0.99892 & 0.00170 + 2454230.705550 & 1.00057 & 0.00170 + 2454230.706059 & 0.99974 & 0.00170 +   + 2454230.770620 & 0.99944 & 0.00150 + 2454230.773240 & 0.99791 & 0.00150 + 2454230.774280 & 0.99852 & 0.00150
+       + @xmath155 & 0.6379 & 0.1792 + @xmath19 & 0.4378 & 0.2933 + @xmath20 & 0.5535 & 0.2351 + @xmath21 & 0.3643 & 0.3178 + @xmath22 & 0.2777 & 0.3191 + @xmath128 & 0.2179 & 0.3162 +   + @xmath128 & 0.1483 & 0.3600 + @xmath155 & 0.5377 & 0.2579 +       + 2454185.910430 & 0.000198 & 0 + 2454198.973147 & 0.000223 & 10 + 2454214.646298 & 0.000280 & 22 + 2454215.952080 & 0.000214 & 23 + 2454535.968246 & 0.000166 & 268 + 2454552.948971 & 0.000147 & 281 + 2454569.929089 & 0.000153 & 294 + 2454594.745943 & 0.000253 & 313 +   + 2454223.797215 & 0.000847 & 0 + 2454230.904913 & 0.000656 &
7 + 2454230.905624 & 0.001106 & 7 +    lc @xmath24 ( k ) & 5650  @xmath156  75 + @xmath15 & 4.4  @xmath156  0.1 + @xmath15 & @xmath157 + @xmath158 ( km s@xmath18 ) & 1.5  @xmath156  1.0 + @xmath25 ( km s@xmath18 ) & 0.85  @xmath156  0.1 + @xmath159fe / h@xmath160 $ ] & @xmath1410.19  @xmath156  0.08 + @xmath161 & @xmath162 + @xmath163 & @xmath164 + @xmath165 ( g @xmath166 ) & @xmath167 + @xmath168 ( @xmath169 ) & @xmath170 + @xmath102 ( @xmath171 ) & @xmath172 + age ( gyr ) & @xmath173 + @xmath174 ( @xmath175 ) & @xmath176 + @xmath177 ( mag ) & @xmath178 + distance ( pc ) & @xmath179 + @xmath139 , @xmath19 , @xmath180 ( km s@xmath18 ) & [ + 27.3 , + 6.7 , + 33.0 ] +    lc @xmath24 ( k ) & 6200  @xmath156  75 + @xmath15 & 4.0  @xmath156  0.1 + @xmath15 & @xmath181 + @xmath158 ( km s@xmath18 ) & 8.5  @xmath156  0.5 + @xmath25 ( km s@xmath18 ) & 1.50  @xmath156  0.05 + @xmath159fe / h@xmath160 $ ] & + 0.14  @xmath156  0.09 + @xmath161 & @xmath182 + @xmath163 & @xmath183 + @xmath165 ( g @xmath166 ) & @xmath184 + @xmath168 ( @xmath169 ) & @xmath185 + @xmath102 ( @xmath171 ) & @xmath186 + age ( gyr ) & @xmath187 + @xmath174 ( @xmath175 ) & @xmath188 + @xmath177 ( mag ) & @xmath189 + distance ( pc ) & @xmath190 + @xmath139 , @xmath19 , @xmath180 ( km s@xmath18 ) & [ @xmath14143.9 , @xmath14139.1 , @xmath1416.9 ] +    lc  @xmath99 ( days ) & 1.30618581 ( fixed ) +  @xmath89 ( hjd)&2,454,185.9104 ( fixed ) +  @xmath191 & @xmath192 +  @xmath193 & @xmath194 +  @xmath195 & @xmath196 +  @xmath22 ( deg ) & @xmath197 +  @xmath198 ( @xmath81 ) & @xmath199 +  @xmath200 ( @xmath81 ) & @xmath201 +  @xmath202 & 0(fixed ) +  @xmath203 ( @xmath114 ) & @xmath204 +  rms ( @xmath81 ) & 22.0 +  @xmath205 ( @xmath206 ) & @xmath207 +  @xmath208 ( @xmath209 ) & @xmath210 +  @xmath211 ( @xmath212 ) & @xmath213 +  @xmath103 ( au ) & @xmath214 +  @xmath215 ( cgs ) & @xmath216 +    lc  @xmath99 ( days ) & 3.553945 ( fixed ) +  @xmath89 ( hjd ) & 2,454,230.9053 ( fixed ) +  @xmath191 & @xmath217 +  @xmath193 & @xmath218 +  @xmath195 & @xmath219 +  @xmath22 ( deg ) & @xmath220 +  @xmath198 ( @xmath81 ) & @xmath221 +  @xmath200 ( @xmath81 ) & @xmath222 +  @xmath202 & 0 ( fixed ) +  @xmath203 ( @xmath114 ) & @xmath223 +  rms ( @xmath81 ) & 11.1 +  @xmath205 ( @xmath206 ) & @xmath224 +  @xmath208 ( @xmath209 ) & @xmath225 +  @xmath211 ( @xmath212 ) & @xmath226 +  @xmath103 ( au ) & @xmath227 +  @xmath215 ( cgs ) & @xmath228 +    . each of the light - curves is labeled with the telescope and filter employed .
we have overplotted the simultaneous best - fit solution , adopting the appropriate quadratic limb - darkening parameters for each band pass ( see text for details).[tres3phot ] ]    .
each of the light - curves is labeled with the telescope and filter employed . as in figure
[ tres3phot ] , we have overplotted the simultaneous best - fit model , adopting the appropriate quadratic limb - darkening parameters for each band pass ( see text for details).[tres4phot ] ]