current observations are more consistent with low - density cold dark matter ( cdm ) dominated cosmogonies . for recent discussions
see , e.g. , ratra et al .
( 1999a ) , kravtsov & klypin ( 1999 ) , doroshkevich et al .
( 1999 ) , colley et al .
( 2000 ) , freudling et al .
( 1999 ) , sahni & starobinsky ( 1999 ) , bahcall et al .
( 1999 ) , and donahue & voit ( 1999 ) .
the simplest low - density cdm models have either flat spatial hypersurfaces and a constant or time - variable cosmological  constant " @xmath3 ( see , e.g. , peebles 1984 ; peebles & ratra 1988 , hereafter pr ; efstathiou , sutherland , & maddox 1990 ; stompor , grski , & banday 1995 ; caldwell , dave , & steinhardt 1998 ) , or open spatial hypersurfaces and no @xmath3 ( see , e.g. , gott 1982 , 1997 ; ratra & peebles 1994 , 1995 ; grski et al .
1998 ) .
while these models are consistent with most recent observations , there are notable exceptions .
for instance , recent applications of the apparent magnitude versus redshift test based on type ia supernovae favor a non - zero @xmath3 ( see , e.g. , riess et al . 1998 , hereafter r98 ; perlmutter et al .
1999a , hereafter p99 ) , although not with great statistical significance ( drell , loredo , & wasserman 1999 ) .
on the other hand , the open model is favored by :    * measurements of the hubble parameter @xmath4
(= @xmath5 ) which indicate @xmath6 at 2 @xmath7 ( see , e.g. , suntzeff et al . 1999 ; biggs et al . 1999 ; madore et al .
1999 ) , and measurements of the age of the universe which indicate @xmath8 gyr at 2 @xmath7 ( see , e.g. , reid 1997 ; gratton et al .
1997 ; chaboyer et al .
this is because the resulting central @xmath9 value is consistent with a low - density open model with nonrelativistic - matter density parameter @xmath10 and requires a rather large @xmath11 in the flat - constant-@xmath3 case .
* analyses of the rate of gravitational lensing of quasars and radio sources by foreground galaxies which require a large @xmath12 at 2 @xmath7 in the flat - constant-@xmath3 model ( see , e.g. , falco , kochanek , & muoz 1998 ) .
* analyses of the number of large arcs formed by strong gravitational lensing by clusters ( bartelmann et al .
1998 ; meneghetti et al . 1999 ) . * the need for mild anti - biasing to accommodate the excessive intermediate- and small - scale power predicted in the dmr normalized flat - constant-@xmath3 cdm model with a scale - invariant spectrum ( see , e.g. , cole et al .
1997 ) .
while the time - variable @xmath3 model has not yet been studied to the same extent as the open and flat-@xmath3 cases , it is likely that it can be reconciled with most of these observations ( see , e.g. , pr ; ratra & quillen 1992 , hereafter rq ; frieman & waga 1998 ; perlmutter , turner , & white 1999b ; wang et al .
1999 ; efstathiou 1999 ) .
we emphasize that most of these observational indications are tentative and certainly not definitive .
this is particularly true for constraints derived from a @xmath13 comparison between model predictions and cosmic microwave background ( cmb ) anisotropy measurements ( see , e.g. , ganga , ratra , & sugiyama 1996 ; lineweaver & barbosa 1998 ; baker et al .
1999 ; rocha 1999 ) , see discussions in bond , jaffe , & knox ( 1998 ) and ratra et al .
( 1999b ) .
more reliable constraints follow from models - based maximum likelihood analyses of the cmb anisotropy data ( see , e.g. , grski et al .
1995 ; yamamoto & bunn 1996 ; ganga et al . 1998 ; ratra et al .
1999b ; rocha et al .
1999 ) , which make use of all the information in the cmb anisotropy data and are based on fewer approximations .
however , this technique has not yet been used to analyze enough data sets to provide robust statistical constraints .
kamionkowski & kosowsky ( 1999 ) review what might be expected from the cmb anisotropy in the near future .    in this paper
we examine constraints on a constant and time - variable @xmath3 that follow from recent type ia supernovae apparent magnitude versus redshift data and recent measurements of @xmath4 and @xmath14 .
we focus here on the favored scalar field ( @xmath0 ) model for a time - variable @xmath3 ( pr ; ratra & peebles 1988 , hereafter rp ) , in which the scalar field potential @xmath15 , @xmath2 , at low redshift .
this scalar field could have played the role of the inflaton at much higher redshift , with the potential @xmath16 dropping to a non - zero value at the end of inflation and then decaying more slowly with increasing @xmath0 ( pr ; rp ) .
see peebles & vilenkin ( 1999 ) , perrotta & baccigalupi ( 1999 ) and giovannini ( 1999 ) for a specific model and observational consequences of this scenario .
a potential @xmath17 could arise in a number of high energy particle physics models , see , e.g. , bintruy ( 1999 ) , kim ( 1999 ) , barr ( 1999 ) , choi ( 1999 ) , banks , dine , & nelson ( 1999 ) , brax & martin ( 1999 ) , masiero , pietroni , & rosati ( 1999 ) , and bento & bertolami ( 1999 ) for specific examples .
it is conceivable that such a setting might provide an explanation for the needed form of the potential , as well as for the needed very weak coupling of @xmath0 to other fields ( rp ; carroll 1998 ; kolda & lyth 1999 , but see periwal 1999 and garriga , livio , & vilenkin 1999 for other possible explanations for the needed present value of @xmath3 ) .
a scalar field is mathematically equivalent to a fluid with a time - dependent speed of sound ( ratra 1991 ) .
this equivalence may be used to show that a scalar field with potential @xmath15 , @xmath2 , acts like a fluid with negative pressure and that the @xmath0 energy density behaves like a cosmological constant that decreases with time .
this energy density could come to dominate at low redshift and thus help reconcile low dynamical estimates of the mean mass density with a spatially - flat cosmological model .
alternate mechanisms that also rely on negative pressure to achieve this result have been discussed ( see , e.g. , zer & taha 1986 ; freese et al .
1987 ; turner & white 1997 ; chiba , sugiyama , & nakamura 1997 ; zer 1999 ; waga & miceli 1999 ; overduin 1999 ; bucher & spergel 1999 )
. however , these mechanisms are either inconsistent or do not share a very appealing feature of the scalar field models .
for some of the scalar field potentials mentioned above , the solution of the equations of motion is an attractive time - dependent fixed point ( rp ; pr ; wetterich 1995 ; ferreira & joyce 1998 ; copeland , liddle , & wands 1998 ; zlatev , wang , & steinhardt 1999 ; liddle & scherrer 1999 ; santiago & silbergleit 1998 ; uzan 1999 ) .
this means that for a wide range of initial conditions the scalar field @xmath0 evolves in a manner that ensures that the cosmological constant contributes a reasonable energy density at low redshift .
of course , this does not resolve the ( quantum - mechanical ) cosmological constant problem .
there have been many studies of the constraints placed on a time - variable @xmath3 from type ia supernovae apparent magnitude versus redshift data ( see , e.g. , turner & white 1997 ; frieman & waga 1998 ; garnavich et al .
1998 ; hu et al . 1999 ; cooray 1999 ; waga & miceli 1999 ; p99 ; perlmutter et al . 1999b ; wang et al . 1999 ; efstathiou 1999 ) . however , except for the analysis of frieman & waga ( 1998 ) , who use the earlier perlmutter et al .
( 1997 ) data , these analyses have mostly made use of the time - independent equation of state fluid approximation to the scalar field model for a time - variable @xmath3 ( perlmutter et al .
1999b and efstathiou 1999 do go beyond the time - independent equation of state approximation by also approximating the time dependence of the equation of state , however , they do not analyze the exact scalar field model ) from gravitational lensing statistics ( see , e.g. , bloomfield torres & waga 1996 ; jain et al . 1998 ; cooray 1999 ; waga & miceli 1999 ; wang et al .
1999 ) , with the exception of the analysis of rq . ] .
a major purpose of this paper is to use the newer supernovae data of r98 and p99 to derive constraints on a time - variable @xmath3 without making use of the time - independent equation of state fluid approximation to the scalar field model .    in the analyses here
we use the most recent r98 and p99 data to place constraints on a constant and time - variable @xmath3 .
we note , however , that this is a young field and insufficient understanding of a number of astrophysical processes and effects ( the mechanism responsible for the supernova , evolution , environmental effects , intergalactic dust , etc . )
could render this a premature undertaking ( see , e.g. , hflich , wheeler , & thielemann 1998 ; aguirre 1999 ; drell et al .
1999 ; umeda et al . 1999
; riess et al . 1999 ; wang 1999 ) . furthermore , other cosmological explanations ( time - variable gravitational or fine structure  constants " ) could be consistent with the data ( see , e.g. , amendola , corasaniti , & occhionero 1999 ; barrow & magueijo 1999 ; garca - berro et al .
1999 ) .
in addition to analyzing just the supernovae data , we also incorporate constraints from recent measurements of @xmath4 and @xmath14 and derive a combined likelihood function which we use to constrain both a constant and a time - variable @xmath3 .
we also examine the effect on the model - parameter constraints caused by varying the prior for @xmath18 .
we emphasize that the tests examined in this paper are not sensitive to the spectrum of inhomogeneities in the models considered .
this spectrum ( possibly generated by quantum mechanical zero - point fluctuations during inflation , see , e.g. , fischler , ratra , & susskind 1985 ) is relevant for some of the other tests mentioned above ( cmb anisotropy , anti - biasing , etc . ) .    in @xmath192
we summarize the computations .
results are presented and discussed in @xmath193 .
conclusions are given in @xmath194 .
we examine three sets of supernovae apparent magnitudes .
we use the mlcs magnitudes for the r98 data , both including and excluding the unclassified sn 1997ck at @xmath20 ( with 50 and 49 sne ia , respectively ) .
the third set are the corrected / effective stretch factor magnitudes for the 54 fit c sne ia of p99 . in all three cases
we assume that the measured magnitudes are independent .
we also assume that sne ia at high and low @xmath21 are not significantly different ( drell et al .
1999 , riess et al . 1999 , and wang 1999
consider the possibility and consequences of evolution ) , and that intergalactic dust has a negligible effect ( aguirre 1999 considers the effects of dust ) .
our analysis of the r98 data is similar to that described in their @xmath194 , with a few minor modifications . for the time - independent @xmath3 model
we compute the likelihood function @xmath22 for a range of @xmath18 spanning the interval 0 to 2.5 in steps of 0.1 , for a range of @xmath23 spanning the interval @xmath24 to 3 in steps of 0.1 , and for a range of @xmath4 spanning the interval 50 to 80 km s@xmath25 mpc@xmath25 in steps of 0.5 km s@xmath25 mpc@xmath25 .
our analysis of the p99 data is similar to theirs , with some modifications .
specifically , we work with their corrected / effective magnitudes and so need only compute a three - dimensional likelihood function @xmath26 . here
@xmath27 is their  hubble - constant - free " b band absolute magnitude ( see p99 ) related to @xmath4 by @xmath28 ( determined by us from results in p99 ) . for the time - independent @xmath3 model
we compute this likelihood function for a range of @xmath18 spanning the interval 0 to 3 in steps of 0.1 , for a range of @xmath23 spanning the interval @xmath29 to 3 in steps of 0.1 , and for a range of @xmath27 spanning the interval @xmath30 to @xmath31 in steps of 0.05 ( which corresponds to about the same interval in @xmath4 used in our analyses of the r98 data ) .
note that in our p99 time - independent @xmath3 model plots we do not show the likelihood for the whole @xmath32 region over which we have computed it .
the spatially - flat time - variable @xmath3 model we consider is the scalar field model with potential @xmath15 , @xmath2 .
when @xmath33 the model tends to the flat - constant-@xmath3 case and when @xmath34 it approaches the einstein - de sitter model .
it is discussed in detail in pr , rp , and rq , and we derive the predicted distance moduli for the sne using expressions given in these papers .    as shown in these papers , the scalar field behaves like a fluid with a constant ( but different ) equation of state in each epoch of the model .
for instance , in the cdm and baryon dominated epoch of the model , it obeys the equation of state @xmath35 ( relating the pressure and energy density of the scalar field ) , where @xmath36 see , e.g. , rq eq . ( 2 ) or see zlatev et al .
( 1999 ) for a more recent derivation .
we shall also have need for an average equation of state parameter , used by perlmutter et al .
( 1999b ) and wang et al .
( 1999 ) , @xmath37 where @xmath38 is the scalar field density parameter and @xmath39 is the scale factor ( with @xmath40 being the present value ) .    in this model , for the r98 data we compute the likelihood function @xmath41 , and for the p99 data the likelihood function @xmath42 . in both cases
we evaluate the likelihood function for a range of @xmath18 spanning the interval 0.05 to 0.95 in steps of 0.025 , for a range of @xmath43 spanning the interval 0 to 8 in steps of 0.5 , and for the same range of @xmath4 or @xmath27 as in the time - independent @xmath3 cases discussed above .
we marginalize these three - dimensional likelihood functions by integrating over @xmath4 ( for the r98 data ) or @xmath27 ( for the p99 data ) and derive two - dimensional likelihood functions , @xmath44 for the constant @xmath3 model and @xmath45 for the spatially - flat time - variable @xmath3 case .
these two - dimensional likelihood functions are used to derive highest posterior density limits ( see ganga et al .
1997 and references therein ) in the @xmath46 or @xmath47 planes .
in what follows we consider 1 , 2 , and 3 @xmath7 confidence limits which include 68.27 , 95.45 , and 99.73% of the area under the likelihood function .    when marginalizing over a parameter or deriving a limit from the likelihood functions , we consider a number of different priors .
we first consider a uniform prior in the parameter integrated over , set to zero outside the range considered for the parameter .
we also incorporate constraints from measurements which indicate @xmath48 at 1 @xmath7 ( see , e.g. , biggs et al . 1999 ; madore et al .
1999 ) , by using the prior @xmath49 , \ ] ] and from measurements which indicate @xmath50 gyr at 1 @xmath7 ( i.e. , 0.5 gyr added to the globular cluster age estimate of chaboyer et al .
1998 ) , by using the prior @xmath51 , \ ] ] with @xmath43 replacing @xmath23 in this expression in the spatially - flat time - variable @xmath3 case .
note that figure 2 of chaboyer et al .
( 1998 ) may be used to establish that a gaussian prior of the form of eq .
( 4 ) is a good approximation to the shifted monte carlo globular cluster age distribution .
see ganga et al .
( 1997 ) for a discussion of gaussian priors in a related context .
this method of incorporating constraints from measurements of @xmath4 and @xmath14 is not identical to the classical @xmath9 cosmological test ( see , e.g. , @xmath1913 of peebles 1993 ) .
finally , since @xmath18 is a positive quantity , we also consider the non - informative prior ( berger 1985 ) here ( see @xmath19b.2 of drell et al .
1999 ; @xmath19vii of jaynes 1968 gives a more general discussion ) .
it is however still of interest to determine how the conclusions depend on the choice of prior .
] , @xmath52 since @xmath18 is bounded from below by the observed lower limit on the baryon density parameter ( @xmath53 , depending on data used , see , e.g. , olive , steigman , & walker 1999 ) , this prior does not result in an infinity .
figure 1 shows the posterior probability density distribution function ( pdf ) confidence contours for the time - variable @xmath3 scalar field model , derived using the three data sets discussed above .
figure 1@xmath54 shows that constraints from the three different data sets are quite consistent . at 2 @xmath7 a large region of the parameter space of these spatially - flat models ( with a constant or time - variable @xmath3 )
is consistent with the sne ia data .
this data favors a smaller @xmath18 as @xmath43 is increased from zero .
a fluid with a time - independent equation of state @xmath55 , @xmath56 , has often been used to approximate the scalar field with potential @xmath15 in the time - variable @xmath3 model .
the solid lines in figure 2@xmath39 show the confidence contours for such a fluid model , derived using the p99 fit c data .
these confidence contours are consistent with those shown in figure 1 of perlmutter et al .
( 1999b ) ( but see discussion below ) , figure 10 of wang et al . ( 1999 ) , and figure 4 of efstathiou ( 1999 ) .
the dashed lines in figure 2@xmath39 show the exact scalar field model contours of figure 1@xmath57 , transformed using the relation between @xmath58 and @xmath43 in the cdm and baryon dominated epoch ( eq .
[ 1 ] ) , which comes to an end just before the present .
the two sets of contours agree near @xmath59 , as they must since at @xmath60 this is the flat-@xmath3 model and @xmath3 does behave exactly like a fluid with a time - independent equation of state .
however , the two sets of confidence contours differ significantly at larger @xmath61 .
since the scalar field @xmath58 eventually switches over to @xmath24 ( in the scalar field dominated epoch , rp ) , it is unclear what significance should be ascribed to this difference .
we stress , however , that since the time - independent equation of state fluid model is an approximation to the time - variable @xmath3 scalar field model , constraints on model - parameter values that are based on the fluid model approximation are only approximate ( and possibly indicative ) . in passing ,
we note that perlmutter et al .
( 1999b ) use @xmath62 ( eq . [ 2 ] ) and not @xmath61 to parametrize the fluid model constraints .
figure 2@xmath63 shows contours of constant @xmath62 as a function of @xmath43 and @xmath18 in the time - variable @xmath3 scalar field model . in a large part of model parameter space
@xmath62 is a sensitive function of both @xmath43 and @xmath18 and hence is not the best parameter to use to describe the time - variable @xmath3 scalar field model .
figure 3 shows the effects of incorporating constraints based on @xmath4 and @xmath14 measurements ( eqs .
[ 3 , 4 ] ) .
panel @xmath64 shows that adding the @xmath4 constraint does not significantly alter the contours derived from the sne ia data alone .
this is expected since the value of @xmath4 used here is very close to the value that is indicated by the sne ia data ( see r98 ) .
however , incorporating the @xmath14 constraint , @xmath50 gyr at 1 @xmath7 , does significantly shift the contours , panel @xmath65 .
this is because the sne ia data alone favor a higher @xmath14 , @xmath66 gyr ( r98 ) , or @xmath67 gyr ( p99 ) .
figure 4 shows constraints on the time - variable @xmath3 model , from the three different sne ia data sets used in conjunction with the @xmath4 and @xmath14 measurements ( eqs . [ 3 , 4 ] ) .
the main effect of incorporating the @xmath4 and @xmath14 constraints is to increase the favored values of @xmath18 ; a weaker effect is the disfavoring of larger values of @xmath43 .
even this extended set of data does not tightly constrain model - parameter values .
figure 5 shows the corresponding constraints on the constant @xmath3 model ( from the sne ia , @xmath4 , and @xmath14 measurements ) .
again , the major effect of including the @xmath4 and @xmath14 data is to increase the favored values of @xmath18 .
a weaker effect is that it reduces the odds against reasonable open models ( grski et al .
1998 ) , but not by a large factor .    figures 6 and 7 show the effects of using the non - informative prior , @xmath68 , of eq .
the major effect is a decrease in the favored values of @xmath18 . in the time - variable @xmath3 case
there is also a slight increase in the favored values of @xmath43 ( see figure 6 ) while in the constant @xmath3 case there is a mild reduction in the odds against reasonable open models ( see figure 7 ) .
if the pdf was narrower ( i.e. , if the error bars on the data were smaller ) , changing from the flat to the non - informative prior would not result in as large a change in the confidence contours .
recent sne ia data do favor models with a constant or time - variable @xmath3 over an open model without a @xmath3 . however , this is not at a very high level of statistical significance .
also , the incomplete understanding of a number of astrophysical effects and processes ( evolution , intergalactic dust , etc . ) means that these results are preliminary and not yet definitive .
the constraints on the time - variable @xmath3 model derived here are based on the exact scalar field model equations of motion , not on the widely used time - independent equation of state fluid approximation equations of motion .
we acknowledge the advice and assistance of r. gott , b. kirshner , l. krauss , v. periwal , s. perlmutter , a. riess , e. sidky , m. vogeley , and i. wasserman , and are specially indebted to j. peebles and t. souradeep .
we acknowledge helpful discussions with i. waga who has also computed the pdf confidence contours for the time - variable @xmath3 scalar field model , with results consistent with those found in this paper .
we thank the referee , i. wasserman , for a prompt and detailed report which helped us improve the manuscript .
we acknowledge support from nsf career grant ast-9875031 .
aguirre , a.n .
1999 , , 512 , l19 amendola , l. , corasaniti , s. , & occhionero , f. 1999 , astro - ph/9907222 bahcall , n.a . ,
ostriker , j.p . ,
perlmutter , s. , & steinhardt , p.j .
1999 , science , 284 , 1481 baker , j.c . ,
et al . 1999 ,
mnras , submitted banks , t. , dine , m. , & nelson , a.e . 1999 , j. high ener .
phys . , 9906 , 014 barr , s.m .
1999 , phys .
b , 454 , 92 barrow , j.d . , & magueijo , j. 1999 , astro - ph/9907354 bartelmann , m. , huss , a. , colberg , j.m . , jenkins , a. , &  pearce , f.r .
1998 , a&a , 330 , 1 bento , m.c . , & bertolami , o. 1999 , gr - qc/9905075 berger , j.o . 1985 , statistical decision theory and bayesian analysis ( new york : springer - verlag ) , 82 biggs , a.d . ,
browne , i.w.a .
, helbig , p. , koopmans , l.v.e . ,
wilkinson , p.n .
, & perley , r.a .
1999 , , 304 , 349 bintruy , p. 1999
, , 60 , 063502 bloomfield torres , l.f . , &  waga , i. 1996 , mnras , 279 , 712 bond , j.r . ,
jaffe , a.h .
, &  knox , l. 1998 , phys .
d , 57 , 2117 brax , p. , & martin , j. 1999 , phys .
b , in press bucher , m. , & spergel , d. 1999 , , 60 , 043505 caldwell , r.r . , dave , r. , &  steinhardt , p.j .
1998 , phys .
, 80 , 1582 carroll , s.m .
1998 , phys .
lett . , 81 , 3067 chaboyer , b. , demarque , p. , kernan , p.j .
, &  krauss , l.m .
1998 , apj , 494 , 96 chiba , t. , sugiyama , n. , & nakamura , t. 1997 , , 289 , l5 choi , k. 1999 , hep - ph/9902292 cole , s. , weinberg , d.h . ,
frenk , c.s . , & ratra , b. 1997 , , 289 , 37 colley , w.n . ,
gott , j.r . ,
weinberg , d.h . , park , c. , & berlind , a.a .
2000 , , 529 , in press copeland , e.j . ,
liddle , a.r .
, & wands , d. 1998 , , 57 , 4686 cooray , a.r . 1999 , , 341 , 653 donahue , m. & voit , g.m . 1999 , , 523 , l137 doroshkevich , a.g . , mller , v. , retzlaff , j. , & turchaninov , v. 1999 , , 306 , 575 drell , p.s . ,
loredo , t.j . , & wasserman , i. 1999 , , in press efstathiou , g. 1999 , , submitted efstathiou , g. , sutherland , w.j .
, & maddox , s.j .
1990 , nature , 348 , 705 falco , e.e . ,
kochanek , c.s . , &  muoz , j.a .
1998 , apj , 494 , 47 ferreira , p.g . , &  joyce , m. 1998 , phys .
d , 58 , 023503 fischler , w. , ratra , b. , &  susskind , l. 1985 , nucl .
b , 259 , 730 freese , k. , adams , f.c . ,
frieman , j.a . ,
& mottola , e. 1987 , nucl .
b , 287 , 797 freudling , w. , et al .
1999 , , 523 , 1 frieman , j.a . , hill , c.t . ,
stebbins , a. , & waga , i. 1995 , , 75 , 2077 frieman , j.a .
, &  waga , i. 1998 , phys .
d , 57 , 4642 ganga , k. , ratra , b. , gundersen , j.o .
, &  sugiyama , n. 1997 , , 484 , 7 ganga , k. , ratra , b. , lim , m.a . , sugiyama , n. , &  tanaka , s.t .
1998 , apjs , 114 , 165 ganga , k. , ratra , b. , & sugiyama , n. 1996 , , 461 , l61 garca - berro , e. , gaztaaga , e. , isern , j. , benvenuto , o. , & althaus , l. 1999 , astro - ph/9907440 garnavich , p.m. , et al .
1998 , , 509 , 74 garriga , j. , livio , m. , & vilenkin , a. 1999 , astro - ph/9906210 giovannini , m. 1999 , astro - ph/9903004 grski , k.m .
, ratra , b. , stompor , r. , sugiyama , n. , & banday , a.j .
1998 , apjs , 114 , 1 grski , k.m . , ratra , b. , sugiyama , n. , & banday , a.j .
1995 , , 444 , l65 gott , j.r . 1982 ,
nature , 295 , 304 gott , j.r .
1997 , in critical dialogues in cosmology , ed .
n. turok ( singapore : world scientific ) , 519 gratton , r.g . ,
fusi pecci , f. , carretta , e. , clementini , g. , corsi , c.e .
, &  lattanzi , m. 1997 , apj , 491 , 749 hu , w. , eisenstein , d.j . , tegmark , m. , &  white , m. 1999 , phys .
d , 59 , 023512 hflich , p. , wheeler , j.c . , &
thielemann , f.k .
1998 , , 495 , 617 jain , d. , panchapakesan , n. , mahajan , s. , &  bhatia , v.b .
1998 , apj , submitted jaynes , e.t .
1968 , ieee trans .
ssc-4 , 227 kamionkowski , m. & kosowsky , a. 1999 , ann .
nucl . part .
sci , in press kim , j.e .
1999 , j. high .
phys , 9905 , 022 kolda , c. & lyth , d.h . 1999 ,
b , 458 , 197 kravtsov , a.v . , &
klypin , a.a .
1999 , , 520 , 437 liddle , a.r . , &  scherrer , r.j .
1999 , phys .
d , 59 , 023509 lineweaver , c.h . , & barbosa , d. 1998 , , 496 , 624 lucchin , f. , & matarrese , s. 1985 , , 32 , 1316 madore , b.f .
, et al . 1999 , , 515 , 29 masiero , a. , pietroni , m. , & rosati f. 1999 , , submitted meneghetti , m. , bolzonella , m. , bartelmann , m. , moscardini , l. , & tormen , g. 1999 , , submitted olive , k.a . ,
steigman , g. , & walker , t.p .
1999 , , in press overduin , j.m .
1999 , , 517 , l1 zer , m. 1999 , , 520 , 45 zer , m. , & taha , m.o .
1986 , phys .
b , 171 , 363 peebles , p.j.e .
1984 , , 284 , 439 peebles , p.j.e .
1993 , principles of physical cosmology ( princeton : princeton university press ) peebles , p.j.e . , &  ratra , b. 1988 , apj , 325 , l17 ( pr ) peebles , p.j.e . , &  vilenkin , a. 1999 , , 59 , 063505 periwal , v. 1999 , astro - ph/9906253 perlmutter , s. , et al .
1997 , , 483 , 565 perlmutter , s. , et al .
1999a , , 517 , 565 ( p99 ) perlmutter , s. , turner , m.s .
, & white , m. 1999b , , 83 , 670 perrotta , f. , & baccigalupi , c. 1999 , , 59 , 123508 ratra , b. 1991 , , 43 , 3802 ratra , b. 1992 , , 45 , 1913 ratra , b. , ganga , k. , stompor , r. , sugiyama , n. , de bernardis , p. , &  grski , k.m .
1999a , apj , 510 , 11 ratra , b. , & peebles , p.j.e .
1988 , , 37 , 3406 ( rp ) ratra , b. , & peebles , p.j.e .
1994 , , 432 , l5 ratra , b. , & peebles , p.j.e .
1995 , phys .
d , 52 , 1837 ratra , b. , & quillen , a. 1992 , mnras , 259 , 738 ( rq ) ratra , b. , stompor , r. , ganga , k. , rocha , g. , sugiyama , n. , &  grski , k.m .
1999b , apj , 517 , 549 reid , i.n . 1997 , aj , 114 , 161 riess , a.g .
, et al . 1998 , aj , 116 , 1009 ( r98 ) riess , a.g . ,
filippenko , a.v . , li , w. , & schmidt , b.p .
1999 , aj , in press rocha , g. 1999 , in dark matter in astrophysics and particle physics 1998 , ed . h.v .
klapdor - kleingrothaus & l. baudis ( bristol : institute of physics publishing ) , 238 rocha , g. , stompor , r. , ganga , k. , ratra , b. , platt , s.r . , sugiyama , n. , & grski , k.m .
1999 , apj , 525 , in press , astro - ph/9905127 sahni , v. , & starobinsky , a. 1999 , int
d , in press santiago , d.i . , & silbergleit , a.s .
1998 , phys .
a , submitted stompor , r. , grski , k.m . , & banday , a.j .
1995 , , 277 , 1225 suntzeff , n.b . ,
et al . 1999 , , 117 ,
1175 turner , m.s .
, & white , m. 1997 , phys .
d , 56 , 4439 umeda , h. , nomoto , k. , kobayashi , c. , hachisu , i. , & kato , m. 1999 , , 522 , l43 uzan , j .-
1999 , , 59 , 123510 waga , i. , & miceli , a.p.m.r .
1999 , , 59 , 103507 wang , l. , caldwell , r.r . , ostriker , j.p . , & steinhardt , p.j . 1999 ,
apj , submitted wang , y. 1999 , , submitted wetterich , c. 1995 , , 301 , 321 yamamoto , k. , &  bunn , e.f .
1996 , apj , 464 , 8 zlatev , i. , wang , l. , & steinhardt , p.j .
1999 , , 82 , 896      fig .
1. pdf confidence contours for the spatially - flat time - variable @xmath3 scalar field model , with potential @xmath15 , derived using the three sne ia data sets .
the @xmath69 axis corresponds to the spatially - flat time - independent @xmath3 case .
confidence contours in panels @xmath70 run from @xmath71 to @xmath72 @xmath7 starting from the lower left hand corner of each panel .
panel @xmath64 shows those derived from all the r98 sne , while @xmath65 is for r98 sne excluding the @xmath20 one , and @xmath73 is for the p99 fit c data set .
panel @xmath74 compares the @xmath75 @xmath7 limits from the three data sets : all r98 sne ( long - dashed lines ) ; r98 sne excluding the @xmath20 one ( short - dashed lines ) ; and the p99 fit c sne ( dotted lines ) .    fig .
2. @xmath64 pdf confidence contours derived from the p99 fit c sne , for a spatially - flat fluid model with equation of state @xmath55 ( solid lines show contours from @xmath76 to @xmath72 @xmath7 starting from the lower left hand corner ) and for the spatially - flat time - variable @xmath3 scalar field model with potential @xmath15 which behaves like a fluid model with equation of state @xmath35 in the cdm and baryon dominated epoch ( dotted lines show contours from @xmath71 to @xmath72 @xmath7 starting from the lower left hand corner ) .
the scalar field model contours were derived by using eq .
( 1 ) to transform the vertical axis of figure 1@xmath57 .
the fluid model likelihood function was computed for the full range of @xmath61 and @xmath18 shown .
the dot - dashed lines bound the @xmath77 region that corresponds to the @xmath78 region of figure 1@xmath57 over which the scalar field model likelihood was computed .
@xmath65 contours of constant @xmath62 ( eq . [ 2 ] ) in the @xmath78 plane for the time - variable @xmath3 scalar field model with potential @xmath15 .
starting at the lower left hand corner , the contours correspond to @xmath79 and @xmath80 .
3. pdf confidence contours derived from the p99 fit c sne , for the spatially - flat time - variable @xmath3 scalar field model with potential @xmath15 .
the dotted lines in panels @xmath70 show the @xmath71 to @xmath72 @xmath7 contours of figure 1@xmath57 .
the solid lines in these panels show the corresponding contours derived from the sne data in conjunction with : @xmath64 @xmath4 measurements ( using eq . [ 3 ] ) ; @xmath65 @xmath14 measurements ( using eq .
[ 4 ] ) ; and @xmath73 @xmath4 and @xmath14 measurements ( using eqs . [ 3 , 4 ] ) .
panel @xmath74 compares the @xmath75 @xmath7 confidence contours from the sne data in conjunction with : the @xmath4 constraint ( long - dashed lines ) ; the @xmath14 constraint ( short - dashed lines ) ; and both constraints ( dotted lines ) .
4. pdf confidence contours for the spatially - flat time - variable @xmath3 scalar field model , with potential @xmath15 , derived using the three sne ia data sets .
panel @xmath64 shows those derived from all the r98 sne , while @xmath65 is from the r98 sne excluding the @xmath20 one , and @xmath73 is for the p99 fit c data .
solid ( dotted ) lines in @xmath70 are the 1 , 2 , and 3 @xmath7 contours from the sne data with ( without ) the @xmath4 and @xmath14 constraints ( eqs .
[ 3 , 4 ] ) ; the dotted lines here are the solid lines in figure 1@xmath81 .
panel @xmath74 compares the @xmath75 @xmath7 limits from the @xmath4 and @xmath14 constraints used in conjunction with : all r98 sne ( long - dashed lines ) ; r98 sne excluding the @xmath20 one ( short - dashed lines ) ; and the p99 fit c sne ( dotted lines ) .
5. pdf confidence contours for the time - independent @xmath3 model , derived using the three sne ia data sets .
panel @xmath64 shows those derived from all the r98 sne , while @xmath65 is from the r98 sne excluding the @xmath20 one , and @xmath73 is for the p99 fit c sne .
solid ( dotted ) lines in @xmath70 are the 1 , 2 , and 3 @xmath7 contours from the sne data with ( without ) the @xmath4 and @xmath14 constraints ( eqs .
[ 3 , 4 ] ) .
panel @xmath74 compares the @xmath75 @xmath7 limits from the @xmath4 and @xmath14 constraints used in conjunction with : all r98 sne ( long - dashed lines ) ; r98 sne excluding the @xmath20 one ( short - dashed lines ) ; and the p99 fit c sne ( dotted lines ) . in all panels , models with parameter values in the upper left hand corner region bounded by the diagonal dot - dashed curve
do not have a big bang .
the horizontal dot - dashed line demarcates models with a zero @xmath3 and the diagonal dot - dashed line running from @xmath82 to @xmath83 corresponds to spatially - flat models .
6. pdf confidence contours ( 1 , 2 , and 3 @xmath7 ) derived from the p99 fit c sne , for the spatially - flat time - variable @xmath3 scalar field model .
panel @xmath64(@xmath65 ) ignores ( accounts for ) the @xmath4 and @xmath14 constraints ( eqs .
[ 3 , 4 ] ) .
the solid ( dotted ) lines use the non - informative @xmath68 ( flat ) prior .
the dotted lines in @xmath64(@xmath65 ) are the same as the solid lines in figure 1@xmath57 ( 4@xmath57 ) .
7. pdf confidence contours ( 1 , 2 , and 3 @xmath7 ) derived from the r98 sne excluding the @xmath20 one , for the time - independent @xmath3 model .
panel @xmath64(@xmath65 ) ignores ( accounts for ) the @xmath4 and @xmath14 constraints ( eqs .
[ 3 , 4 ] ) .
the solid ( dotted ) lines use the non - informative @xmath68 ( flat ) prior .
the dotted lines in @xmath64(@xmath65 ) are the same as the dotted ( solid ) lines in figure 5@xmath63 .
the dot - dashed lines are described in the caption of figure 5 . in the non - informative prior cases the likelihood function
is computed down to @xmath84 .


when evaluating justification for new facilities , one may question whether and by how much the new technique is superior to already existing ones .
we describe a general approach which leads to a quantitative answer , expressed as a relative figure of merit of two observational techniques .
the approach is based on the comparison of the corresponding errors in the model parameters determination .
the derived figure of merit can also be used in the closely related problem of optimization of complex observations , thus contributing in maximizing their scientific output .
we illustrate the developed method by comparing the technique of optical long baseline interferometry to the classical spectroscopy , both used to find the parameters of the gaseous envelope around the p cyg star . for this first trial application
, we use simulated data , obtained by solving the radiative transfer problem for grids of physical models of the envelope and then adding the photon noise .
the chosen example of the interferometry vs spectroscopy comparison appears to be timely . indeed , in the recent years
, interferometry has been undergoing a steady progress witnessed by publication of fringe visibilities for circumstellar envelopes of various types sometimes combined with a good spectral resolution ( mourard et al .
@xcite , vakili et al .
@xcite , @xcite ; quirrenbach et al .
@xcite , @xcite , @xcite ; harmanec et al .
@xcite ; stee et al .
@xcite ; see also expected performances of interferometric arrays under construction , vlti , von der lhe et al .
@xcite , petrov et al .
@xcite , and chara , mcalister et al .
@xcite ) . while this effort is fully justified , for the knowledge of the flux distribution on small angular scales can unveil important new features of astrophysical objects , it is also justified to wonder whether this effort is well payed for all classes of objects and for all kinds of physical problems . as a matter of fact , interferometric observations as compared to those at a single telescope are more complex , more expensive , more time consuming , and , when possible to compare , have a larger accumulated measurement error due to their complexity . furthermore , in the foreseeable future the instruments of this type will remain orders of magnitude less numerous and by all given reasons much less available for the community than single telescopes .
thus , it makes sense to evaluate in quantitative terms what are the targets for which the interferometric instrumentation has the highest potential in terms of astrophysical information as compared to single telescopes .
also , the optimization of interferometric observations may be a crucial aspect .
indeed , the existing and forthcoming long baseline interferometers comprise only few individual apertures , covering a small set of spatial frequencies during one observing cycle . on the other hand , the information provided by observations can depend critically on spatial frequencies .
certainly , in such a situation the scientific yield of the interferometric array could be greatly improved if thorough model calculations can indicate the configurations expected to be the most useful for determination of aimed physical parameters of the object .    not surprisingly , the relative figure of merit derived in the present article depends not only on the experimental techniques , but also on the chosen set of model parameters , reflecting the fact that the physical model makes a necessary part of the merit problem . however , in modern astronomy , designing and building instrumentation on the one hand , and modeling objects on the other hand , were pushed to such a level of sophistication that they are as a rule two distinct activities .
the optimum use of interferometry needs a close collaboration of two communities .
after these preliminary remarks , let us specify that in the present first trial application , our discussion will be limited to the comparison of observations with a two - aperture optical long baseline interferometer ( hereafter olbi ) and the classical single telescope spectroscopy ( hereafter shortly spectroscopy ) , both applied to study a circumstellar gaseous envelope formed by the wind of a massive star .
the olbi observations are assumed to have the same spectral resolution as spectroscopy . the discussion will be further limited to spectral profiles and visibilities in the h@xmath1 line of the p  cyg envelope , allowing to use published high quality observations both in olbi and spectroscopic , and thus to provide a clear numerical illustration .    in our earlier attempts to solve the olbi vs spectroscopy merits problem ( burgin & chalabaev @xcite ; bourguine & chalabaev @xcite ) , we used qualitative comparison of observables , visibility and spectral profile , computed for a small number of envelope models .
the results were ambiguous .
analising them , we arrived to the firm conviction that meaningful conclusions can be obtained only if ( 1 ) intercomparison of different types of observations is performed using a clearly defined figure of merit , indicating how much information on the studied object is provided by various techniques , and ( 2 ) a sufficiently large ranges in the space of envelope model parameters are analyzed .    in the case of p
cyg , the required number of computed models ( see sect .  [ results ] ) turned out to be @xmath2 .
an efficient simplified envelope model code , computing the emergent spectral profile and visibilities for the h@xmath1 line , was developed and is described in sect .
[ model ] .
we considered only spherically symmetric outflow models , the choice which is justified in sect .
[ model ] .
we consider astrophysical observations and their interpretation in the framework of the theory of parameter estimation ( e.g. kendall & stuart @xcite ) .
that is , the observed object is assumed to be exactly described by a fixed multiparametric physical model , the errors being solely those of the experimental measurements ( see however sect .  [ systerrs ] ) .
any observable quantity is then a function of the vector of @xmath3 model parameters @xmath4 .
an observation consists in finding the values of @xmath5 observables @xmath6 , where each observable @xmath7 is either measured directly or can be calculated in a model - independent way as a known function of directly measurable quantities .
obviously , observations with instruments of different types provide observables of different number and nature .
the @xmath8-th observable can be represented as @xmath9 , where @xmath10 is the `` theoretical '' value of corresponding observable in the absence of observational errors and @xmath11 is the random vector of observational errors , which is assumed here to be distributed normally .
we consider the case when @xmath12 and the parameters @xmath13 could have been uniquely determined if the vector @xmath14 were known .
we shall call the parameters to be determined the `` target parameters '' and denote a set of @xmath15 target parameters , or target set , as @xmath16 .
it is possible that @xmath17 and the set @xmath16 is a proper subset of the set of all parameters of the model .
the situation of this kind may arise for two reasons .
first , a model parameter may be set to some a priori fixed value because , say , it has been measured earlier by entirely different methods .
this reduces the dimension of the problem and computational difficulties , the general method remaining the same .
let @xmath18 denote the set of @xmath19 fixed parameters , and @xmath20 be the value a priori assigned to a parameter @xmath21 .
secondly , some of the unknown physical parameters may eventually be considered `` not interesting '' if their values are not relevant to the astrophysical problem under study .
following the theory of parameter estimation , we will call them the `` nuisance '' parameters and denote the set of nuisance parameters , or `` nuisance set '' , as @xmath22 , their number being @xmath23 . as a rule
, their influence can be separated in the error analysis only at the final stage of computations , for their values are to be calculated along with the values of target parameters .
interpretation of an observation consists in calculating the vector @xmath24 that provides the best fit to the observed values .
here @xmath25 is the index of the component of vector @xmath26 corresponding to the same model parameter as the @xmath8-th target parameter ( @xmath27 ) .
the precision of the values @xmath28 resulting from interpretation of the observations is characterized by @xmath29 and @xmath30 , where for any random vector @xmath31 expressions @xmath32 and @xmath33 denote its mean and its covariance matrix respectively , and @xmath34    in what follows , we assume that @xmath35 is an unbiased estimate of @xmath14 , i.e. @xmath36 , and that the measurement errors are small enough , so that the error analysis can be performed using the linearized version of the least squares method .
we additionally assume that all a priori fixed parameters are set to their true values ( see however [ systerrs ] ) , that is @xmath37 for any @xmath21 , then @xmath38 and statistical properties of errors in parameter determination are completely characterized by the covariance matrix of errors c = d.      according to kendall & stuart ( ( * ? ? ?
* chap.19 ) ) , the covariance matrix of errors c is related to @xmath14 and @xmath39 by the following expression : @xmath40 where @xmath41 is the @xmath42 matrix with elements @xmath43 for @xmath44 and @xmath45 is the vector of experimental errors defined in sect .
[ general - theory ] .
the natural scalar characteristics of the precision of parameter determination for a given set @xmath16 of target parameters is the following principal subdeterminant of the covariance matrix : @xmath46 geometrically , this subdeterminant is proportional to the hypervolume of the scattering ellipsoid in the space of target parameters .
it depends not only on the physical model and errors of measurements , but also on the analytical form used for description of the model : two physically equivalent but mathematically different ( e.g. interrelated by a reversible substitution of variables ) analytical representations could yield entirely different values of @xmath47 .    however , a pair of observational methods can well be compared if one makes use of the ratio @xmath48 where @xmath49 and @xmath50 are the subdeterminants calculated for the observation of the same object with instruments i and ii . this ratio depends only on the physical model used , and the obtained precision on the model parameters .
we shall call the quantity @xmath51 , introduced by the eq .
( [ rdef ] ) , the `` random error ratio '' , or the `` relative figure of merit '' .    if the instruments are of the same kind , and differ from each other only in precision , the value @xmath51 is merely the ratio of observational errors .
if the instruments are of different kinds , providing observables of different nature and number , and the target set consists of only one model parameter , the value @xmath51 is merely the ratio of resulting random errors in the parameter determination .
however , in the general case of instruments of arbitrary kinds and multiparametric models , no simple ratio of errors exists , and the evaluation of relative merits can be done only using the quantity @xmath52 defined by eq .
( [ rdef ] ) .
it allows to compare various observational techniques applied to objects described by various multiparametric models : equality @xmath53 means that the instrument i is better suited for determination of parameters from the target set than the instrument ii .
the fact that the description of an object by a multiparametric physical model is only an approximation to the reality implies that we have to consider the robustness of the method , that is the stability of the results it yields with respect to deviations of the real situation from the model .
the present framework offers a way to obtain certain quantitative characteristics of the robustness .
indeed , let us consider a multiparametric model with @xmath54 . if @xmath55 for some @xmath21 then , in general , @xmath56 for @xmath57
that is , in addition to random errors of observational origin , the result is biased by systematic errors due to inaccurate interpretation .
the value of that bias characterizes the robustness of the method with respect to deviations of @xmath20 from its true value .    in the linear approximation , @xmath58 , where @xmath59 is the @xmath60 matrix with elements @xmath61 where @xmath62 .
a comprehensive study of systematic errors requires a joint analysis of individual elements of the matrix @xmath59 and uncertainties in model parameters from @xmath63 . in the present paper
, we will develop a simplified approach providing semiquantitative indications concerning the relative robustness of different observational techniques .
let us first define for each observational method under consideration the value @xmath64 when the estimates for target parameters are uncorrelated with estimates for parameters from @xmath63 , equality @xmath65 takes place : the method is robust ( with respect to deviations of specified form from the model ) .
when such a correlation exists , @xmath66 exceeds unity ; the closer the correlation , the larger its value , eventually implying poor robustness of the method .
further , to compare the robustness of methods i and ii with respect to inaccuracies in parameters from @xmath63 , we will define the `` robustness ratio @xmath67 '' as follows : @xmath68 where function @xmath69 is defined in eq .
( [ rdef ] ) . note that once the values of @xmath70 , which are necessary for analysis of random errors , are found , comparison of robustness is straightforward .
in this section we obtain the expressions for the elements of the covariance matrix @xmath71 in the form that can be directly used to compare the olbi and spectroscopy .
first , the measurement errors have to be specified .
for this first trial application , the only considered source of random errors is the photon noise .
other sources of error , in particular those in olbi arising from atmospheric seeing and calibrations of the modulation transfer function ( see roddier & lna @xcite ; mourard et al .
@xcite ) are important and has to be incorporated in a future work .
also , the spectral coverage and the spectral resolution have to be specified .
we assume that they are the same for both kinds of observations ; the covered spectral region is @xmath72 $ ] , the spectral resolution is sufficiently high , i.e. the monochromatic intensity received from each point of the observed object remains nearly constant across each spectral channel . in other words , we assume that @xmath73 , the width of the spectral channels , in terms of velocity is lesser than the thermal velocity of the emitting atoms . for the p  cyg type envelopes
this corresponds to @xmath74 , which is easily satisfied with modern spectrographs .
then one can formally consider the convenient limit case @xmath75 . in this case , the resulting covariance matrices tend to a limit which is independent of the spectral resolution ( cf .
( [ cspec ] ) and ( [ colbi ] ) ) .
the physical quantity provided by spectroscopy is the flux density as a function of wavelength . in terms of observables
, it is given by the vector @xmath76 with the components @xmath77 where @xmath78 is the number of photons recorded in the _
i_-th spectral channel during the exposure , and @xmath79 is the number of channels .
the random values @xmath78 are related to the physical parameters of the object by @xmath80 where @xmath81 @xmath82 is the monochromatic flux from the object at wavelength @xmath83 , @xmath84 is the width of @xmath8-th spectral channel , and @xmath85 is the measurement error due to the photon noise , with the variance given by @xmath86 finally the factor @xmath87 is given by @xmath88 where @xmath89 is the instrument effective aperture and @xmath90 is the exposure time .
hereafter we suppose that @xmath91 , where @xmath92 and @xmath93 pertain correspondingly to the spectrometer and interferometer .
since the values @xmath92 and @xmath93 enter the final results only through the ratio @xmath94 , we can for the sake of simplicity set @xmath95 .
thus , @xmath96 and the elements of the covariance matrix of measurement errors are given by @xmath97 substituting eqs .
( [ ys])([nncovar ] ) into eq .
( [ covmateval ] ) and proceeding to the limit @xmath98 we obtain : @xmath99 for spherically symmetric objects considered in the present paper , @xmath100 where @xmath101 is the monochromatic intensity at the angular distance @xmath102 from the center of the object .
we consider the case of a two - aperture olbi with an adjustable baseline .
the physical quantity eventually provided by olbi measurements is the fringe visibility v as a function of the baseline vector .
the visibility v is equal to the real part of the complex degree of coherence , and therefore is the fourier transform of the brightness distribution in the focal plane by virtue of the van cittert - zernicke theorem ( e.g. mariotti @xcite , textbooks of perina @xcite , and goodman @xcite ) .
however , the value estimated in practice is @xmath103 .
we refer the reader for details to the thorough discussion by mourard et al .
( @xcite ) . in the present analysis
we consider the estimates @xmath104 of the quantity @xmath105 where @xmath8 refers as previously to the @xmath8-th spectral channel , and @xmath106 is the projected baseline length of the interferometer . in the case of
spherically symmetric objects the intensity distribution is circularly symmetric and @xmath107 is given by the normalized hankel transform as follows ( bracewell @xcite ) : @xmath108 where @xmath109 is the bessel function , and @xmath110 .    statistical errors @xmath111 obey the equations @xmath112 and @xmath113 the validity of eqs .
( [ wwcovar ] ) and ( [ nwcovar ] ) for @xmath114 is evident .
the case @xmath115 is treated in [ olbidiag ] .    as a rule ,
the spectrum of the object , i.e. the values @xmath78 , is also recorded .
consequently , the vector of observables for the olbi and the elements of its covariance matrix of errors are given by @xmath116 and @xmath117 substituting this expressions into eq .
( [ covmateval ] ) , we obtain in the limit @xmath98 that @xmath118 where @xmath119 is defined in eq .
( [ cspec ] ) .      as discussed in sect .
[ randerr ] , the relative figure of merit of olbi as compared to spectroscopy is the random error ratio : @xmath120 from eqs .
( [ colbi ] ) and ( [ risdef ] ) , it follows that @xmath121 , which merely reflects the fact that the observables of spectroscopy constitute a subset of the olbi observables .
consequently , in the adopted comparison , the statistical errors in parameter determination by olbi are always lesser than those obtained when only spectroscopic data are used , if however the spectroscopic data are of equal precision .    as explained in sect .
[ spectroscopy ] , the eq .
( [ risdef ] ) is obtained in the assumption that @xmath91 . to compare a pair of instruments for which that equality does not hold , the rhs of eq .
( [ risdef ] ) should be multiplied by @xmath122 .
the star p  cyg ( b1 iape ) is characterized by a high mass - loss rate and exhibits bright optical emission lines , formed in the dense and nearly fully ionized stellar wind .
it belongs to the class of luminous blue variables , a short - living transition phase in the evolution of a massive star at the end of the hydrogen burning ( see e.g. humphreys & davidson @xcite , langer et al .
@xcite , maeder @xcite ) .
there were at least two strong reasons to choose p  cyg as the astrophysical object of the present study .
first , this bright object , @xmath123 , is one of the best studied emission line stars .
although in this first trial we use simulated data , it is important to note that for p  cyg there exists a rich literature providing not only high quality spectroscopic data ( e.g. scuderi et al .
@xcite ) , but also interferometric data ( vakili et al .
@xcite ) , as well as thorough theoretical analysis ( drew @xcite ; pauldrach & puls @xcite ) .
this insures that the present work can be followed up by a practical application .
secondly , it happens that in the growing list of outflows known to be non - spherical ( e.g. wolf , stahl , fullerton @xcite ) , that of p  cyg is an exception exhibiting spherical symmetry to a good degree of accuracy ( nota @xcite ) , becoming clumpy only on short time and small flux scales ( taylor et al .
@xcite ; vakili et al .
@xcite ; nota @xcite ) .
this implies that the assumption of spherical symmetry of the envelope , used in the present work , is realistic .
let us remind that it allows to compute a model in a reasonable amount of time , and thus to explore a large domain of parametric space , which is a requirement for the evaluation results to be meaningful .    for the distance to the star @xmath124 , its radius @xmath125 and the effective temperature @xmath126 we adopted the following values : @xmath127 , @xmath128 , @xmath129 ( lamers et al .
@xcite ; pauldrach & puls @xcite ) .
the spectrum of the star was assumed to be blackbody .
our study is limited to the hydrogen h@xmath1 line , which is the most prominent and the best studied feature in the spectrum of the star .
the mass distribution within a spherically symmetric and stationary outflow is described by its mass loss rate @xmath130 and the outflow velocity @xmath131 , which we assume to depend on radial distance @xmath132 in the following way : @xmath133 where @xmath134 is the terminal velocity of the wind , @xmath135 is velocity of the wind at the base of the photosphere , and @xmath1 is a dimensionless parameter that characterizes the rate at which the velocity approachs its asymptotic value at large @xmath132 .
the radial distance @xmath132 is measured in units of stellar radius @xmath125 , while @xmath130 , @xmath135 , @xmath134 , and @xmath1 are parameters of a model .    the temperature in the envelope is often assumed to be constant .
however , as it was shown by drew ( @xcite ) , across the region of the h@xmath1 line formation @xmath136 the temperature in the envelope can vary by as much as 6000k .
therefore , along with isothermal models , we also computed the emergent emission for non - isothermal models .    to keep a finite number of scalar parameters , the dependence of envelope temperature @xmath137 on radial distance is approximated in the following manner
: it is assumed that @xmath138 for @xmath139 , that @xmath137 is a linear function of @xmath140 at each interval @xmath141 $ ] for @xmath142 , and that @xmath143 for @xmath144 . the values @xmath145 were fixed for each model of the family .
we consider only non - increasing temperature laws @xmath137 , that is @xmath146 for @xmath147 , where @xmath148 .
since the computations are organized in such a way that the model parameters vary independently on each other , it is the values @xmath149 that are used along with @xmath150 as the parameters defining the envelope temperature , and the vector of model parameters is given by @xmath151      the `` supersonic '' sobolev approximation is a well known efficient method of calculating line spectral profiles ( e.g. castor @xcite ) .
however , it would have a low accuracy if applied in the p  cyg case , for a noticeable fraction of h@xmath1 flux is emitted in the inner part of envelope , where the outflow velocity is comparable to that of sound .
we adopted a mixed method computing the source function in the sobolev approximation , and then finding the emergent intensities by exact numerical solution of the transfer equation .
as it was shown by hamann ( @xcite ) , the errors in line profiles computed in the original sobolev approximation come mainly from calculations of the emergent intensities , whereas the source function is accurate for a wide range of physical conditions .
the choice of the model of the hydrogen atom was based on the fact that the regions of the envelope emitting the major fraction of the h@xmath1 flux are nearly completely ionized and opaque to lyman continuum and l@xmath1 , so that direct recombinations to and photoionizations from the ground level cancel out , and l@xmath1 is saturated ( drew @xcite ) .
the populations of levels @xmath152 and @xmath153 are mainly defined by collisional and radiative transitions between these two levels , radiative ionizations due to stellar radiation , and radiative recombinations ( including indirect ) to the level @xmath153 .
we therefore adopted the three - level + continuum model of hydrogen atom .
the balance equations for @xmath154 and @xmath155 , the number densities of hydrogen atoms respectively at levels 2 and 3 , take then the following form : @xmath156 where @xmath157 is the spontaneous emission coefficient , @xmath158 and @xmath159 are collisional excitation and de - excitation coefficients for the transition , @xmath160 and @xmath161 are respectively ionizations coefficients and recombination rates for the levels involved , and @xmath162 is the escape probability .
since the hydrogen is nearly completely ionized , the values @xmath158 , @xmath159 , @xmath163 and @xmath164 can be considered independent on @xmath154 and @xmath155 , the only source of non - linearity in eqs .
( [ balance ] ) being the terms containing @xmath165 .
the escape probability @xmath162 is a rather complex function of the population of the lower level of the transition , so that commonly the eqs .
( [ balance ] ) are solved by iterations , and this and , in particular , the recalculation of the escape probability , is by far the most time consuming operation .
to speed up computations , we developed a new approximate method for solving the equations of statistical equilibrium .
it is based on the fact that , as shown in [ epapprox ] , the function @xmath166 can be approximated with sufficiently good accuracy by a simple analytic expression as follows : @xmath167 where @xmath168 @xmath169 and @xmath170 are statistical weights of the levels , and @xmath83 is the line wavelength .    introducing the dimensionless variables @xmath171 where @xmath172 and @xmath173 is the energy of transition , and substituting the approximation ( [ betapp ] ) into eqs .
( [ balance ] ) we obtain after some algebra that @xmath174 is the positive solution of the quadric equation @xmath175 where @xmath176 here @xmath177 and @xmath178 when @xmath174 is calculated , the value of @xmath179 can be found using the relation @xmath180 the source function is then easily computed from level populations , which can be obtained using eq .
( [ pdef ] ) .
integration of the transfer equation was performed using the code developed by bertout(@xcite ) , who kindly provided it to the authors . in computing the line profile
, the code assumes the envelope to be isothermal , so that its use for a non - isothermal case requires some comments .
the envelope temperature enters the calculations at two points : ( 1 ) in calculating the source function , through coefficients @xmath158 and @xmath159 of eq .
( [ balance ] ) . since the bertout code is applicable to arbitrary source function , variations in temperature does not cause any difficulties here .
( 2 ) in integrating the transfer equation , @xmath137 enters the result through the local doppler line width .
since we consider only the envelopes with relatively low temperature contrast @xmath181 , and dependence of the doppler width on temperature is rather week , the error induced by this isothermal code structure is still negligible .
once the intensity @xmath101 is computed , the observable quantities , i.e. the emergent flux @xmath182 and the visibility squared @xmath183 are obtained using eqs .
( [ fdef])([vaxisym ] ) .
the relative figure of merit @xmath184,@xmath16,@xmath185,@xmath186 and the robustness ratio @xmath187,@xmath188,@xmath63,@xmath186 are calculated in two stages , implemented as separate programs .
first , for the given grids of model parameters @xmath189 and the projected baseline lengths @xmath190 , we calculate and store the matrices @xmath191 and @xmath192 , defined in eqs .
( [ cspec ] ) and ( [ colbi ] ) .
the defining parameters of the grids are entered as input data , the grid of models being constructed as the direct product of uniform grids for individual parameters .
the derivatives entering definitions of the matrices are approximated by finite differences .
all the model parameters that are not constant on the grid are treated as adjustable , that is either target or nuisance .
the distinction between those two types of parameters is irrelevant at this stage .    at the second stage , for given partitions of set @xmath193 of all adjustable parameters on the subsets @xmath16 and @xmath22 , we compute the values of @xmath184,@xmath16,@xmath22,@xmath186 and @xmath187,@xmath16,@xmath63,@xmath186 ( see eqs .
( [ rdef ] ) and ( [ sdef ] ) ) on the grid . in this way
, the dependence of the results on @xmath16 and @xmath22 can be studied without repeating the time consuming physical modeling of the envelope .
in total , the values @xmath194 and @xmath195 were calculated as a function of projected baseline length for various sets of target and nuisance parameters at more than @xmath0 points of the parametric space .
table [ grids ] presents the principal characteristics of the model grids .
the ranges of variations of model parameters on the grids were chosen so as to include the values of parameters from earlier works on p  cyg cited in sect .
[ stpars ] .    lrllllllrrrrrrrr grid & @xmath196&@xmath197&@xmath1 & @xmath130&@xmath134&@xmath135 & @xmath150&@xmath163&@xmath198&@xmath164&@xmath199&@xmath200&@xmath201&@xmath202 & @xmath203 + no . & & & & @xmath204 & & & k & & k & & k & & k & & k + 1 & 10000 & 3&3.5 & 1.010@xmath205 & 200 & 15.0 & 11000 + & & & 4.5 & 2.010@xmath205 & & 55.0 & 15000 + 2 & 10000&3 & 3.5 & 1.010@xmath205 & 200 & 15.0 & 11000 & 2.0 & 1300 & 3.4 & 2100 & 5.0 & 200 & 10.0 & 3000 + & & & 4.5 & 2.010@xmath205 & & 55.0 & 15000 + 3 & 32805&4&3.5 & 1.010@xmath205 & 150 & 15.0 & 11000 + & & & 4.5 & 2.010@xmath205 & 250 & 55.0 & 15000 + 4 & 15625 & 6&3.5 & 1.010@xmath205 & 200 & 15.0&11000 & 2.0 & 0 & 5.0 & 0 + & & & 4.5 & 2.010@xmath205 & & 55.0 & 15000 & & 1300 & & 2300 + 5 & 46656 & 6&3.25 & 1.010@xmath205 & 200 & 15.0&11000 & 2.0 & 0 & 3.4 & 0 + & & & 4.25 & 2.010@xmath205 & & 55.0 & 15000 & & 1300 & & 2100 + 6 & 6561 & 8&3.5 & 1.010@xmath205 & 200 & 15.0 & 11000 & 2.0 & 0 & 3.4 & 0 & 5.0 & 0 & 10.0 & 0 + & & & 4.5 & 2.010@xmath205 & & 55.0 & 15000 & & 1300 & & 2100 & & 200 & & 3000 + 7 & 65536 & 8&3.5 & 1.010@xmath205 & 200 & 15.0 & 11000 & 2.0 & 0 & 3.0 & 0 & 4.0 & 0 & 5.0 & 0 + & & & 4.5 & 2.010@xmath205 & & 55.0 & 15000 & & 1300 & & 2100 & & 200 & & 3000 + 8 & 6561 & 8&4.0 & 1.510@xmath205 & 180 & 20.0 & 14000 & 2.0 & 0 & 3.0 & 0 & 4.0 & 0 & 5.0 & 0 + & & & 4.5 & 2.310@xmath205 & 220 & & 16000 & & 1000 & & 1000 & & 1000 & & 1000 + 9 & 19683 & 9&3.5 & 1.510@xmath205 & 220 & 15.0 & 16000 & 2.0 & 0 & 3.0 & 0 & 4.0 & 0 & 5.0 & 0 + & & & 4.5 & 2.310@xmath205 & 300 & 25.0 & 18000 & & 1300 & & 2100 & & 200 & & 3000 +    the main result of
the present study is that the increase in accuracy provided by olbi strongly depends on the number of free parameters in the model : it varies from almost negligible for models of low dimension to very significant for models of high dimension . since in our case
the boundary between `` low '' and `` high '' is nearly coincident with division onto models with fixed and adjustable thermal structure , we discuss the results pertaining to those two classes of models separately .      among our grids
, the grid  1 most closely corresponds to the models often used in interpreting h@xmath1 observations of p  cyg stars , where @xmath134 is excluded from the set of adjustable parameters of the model ( see e.g. scuderi et al .
@xcite ) and its value is taken from the analysis of absorption lines of metals in the ultraviolet ( casatella et al .
@xcite ; lamers et al @xcite ) .
the typical results for an individual point in the parametric space are illustrated by fig .
[ isothtpc ] , where the random error ratio @xmath206 is plotted as a function of the normalized baseline length @xmath207 for various @xmath16 and @xmath22 . here
@xmath106 is the projected baseline length , @xmath208 , and @xmath209 is the angular radius of the central star . for the adopted values of @xmath210 , @xmath127 , and @xmath211
one gets @xmath212 .
the principal features of dependence of @xmath213 on @xmath106 are independent on details of physical model used and can be easily explained qualitatively . as shown in sect .
[ relfigmer ] , inequality @xmath121 holds for all @xmath214 .
when the projected baseline length is very small , the interferometry does not provide any additional information as compared with the spectroscopy , for the object gets unresolved .
consequently , @xmath215 for @xmath216 . in the opposite case of large @xmath106
the fringe contrast is close to 0 for all physically realistic intensity distributions and its dependence on model parameters is impossible to measure because of observational errors .
that is in that case the interferometry again provides no additional information and @xmath215 for @xmath217 .
it follows from the above stated general properties of @xmath213 that for any @xmath16 and @xmath22 it has a ( possibly non - unique ) minimum at a certain finite value of @xmath106 .
the location of the minimum indicates the projected baseline length at which interferometric observations are most informative , and the corresponding value of @xmath213 characterize to what degree the interferometry can reduce the errors in model parameters determination as compared to spectroscopy .
let us note that it would be insecure to compare the merits of the methods using values of @xmath218 obtained for an isolated point or even for few arbitrarily chosen points of parametric space .
since the model parameters are not known a priori , it is necessary to use integral characteristics describing the behavior of @xmath213 on the whole grid .
we will use the following three values calculated as a function of @xmath219 for various @xmath16 and @xmath22 : @xmath220 @xmath221 where @xmath222 is the number of points on the grid for which function @xmath223 reaches its minimum at @xmath224 .    the function defined in eq .
( [ poptdef ] ) is immediately related to the problem of optimal choice of the baseline length : for given @xmath16 and @xmath22 , the higher the value @xmath225 , the higher the probability that the interferometric observations at the baseline length @xmath106 would yield most accurate model parameters .
the dependence of @xmath226 and @xmath227 on the normalized baseline length and set of target parameters for grid  1 is shown in fig .
[ isothrdt ] . for this grid ,
the value of @xmath228 is close to unity for all @xmath16 and @xmath22 and is not shown in the plot .
it can be easily seen that for this grid of models olbi provides only relatively small reduction of the random error in parameter determination as compared with spectroscopy .
this conclusion remains also valid if we add the terminal velocity @xmath134 to the set of adjustable parameters ( grid 3 ) .
as @xmath229 for all @xmath16 and @xmath22 , the same approximate equality is evidently valid for the robustness ratio @xmath230 , which indicate ( see sect .
[ systerrs ] ) that , for models with a priori fixed thermal structure , interferometric data will neither reduce significantly the level of systematic errors .
the number of possible target sets is an exponentially increasing function of the dimension @xmath197 of the model parameter space , so that when the values @xmath149 are added to the set of unknowns , it becomes impossible to make an exhaustive presentation of even the integral results , i.e. the functions @xmath231 , + @xmath232 , and @xmath225 . in fig.[3lmerr ] , we show the dependence of @xmath226 , @xmath228 , and @xmath227 on the number of nuisance parameters for one of our simplest grids of models with adjustable thermal structure , the target set consisting of one parameter , the mass loss rate @xmath130 . from many points of view , it is the most important parameter of the stellar wind , characterizing both its effect on the evolution of the star and the influence of the outflowing matter on the surrounding interstellar medium .
as it can be seen , the lower limit of the random error ratio @xmath226 rapidly decreases as the number of adjustable temperature parameters increases , so that when the thermal structure of the envelope becomes a nuisance parameter , interferometry can reduce the random error in @xmath130 by as much as an order of magnitude .
note , that the decrease in @xmath226 is caused by the fact that parameters @xmath149 are adjustable , rather then by deviation from isothermicity .
this can be seen from the upper two plots in fig .
[ 3lmerr ] , where @xmath233 for a grid containing high proportion of models with noticeable temperature gradient , and is also supported by the results obtained for grid  2 , which differs from grid  1 only in that its thermal structure approximates that of drew s model b , which corresponds to @xmath234 , @xmath235 , @xmath236 , and @xmath237 in our notations , instead of being isothermal .    for other single - parameter target sets the dependencies of @xmath226 and @xmath228 on @xmath22
are qualitatively the same .
[ 3l1 + 5 ] displays the results for the most important case when all the parameters of the model except for @xmath134 are adjustable .
as one can see , again , in favorable circumstances interferometry can increase accuracy by an order of magnitude .
in contrast to @xmath226 , the value @xmath228 is close to unity for all projected baseline lengths and combinations of sets @xmath238 and @xmath22 that we have studied .
this signifies that for any projected baselength there exists a combination of model parameters for which interferometry practically do not provide reduction in error .
it is interesting to notice that for target sets containing several parameters ( see fig .
[ 3lmultyt ] ) , the influence of individual parameters on @xmath226 and @xmath228 to some extent averages out : @xmath226 is systematically higher , and @xmath228 is systematically lower than the corresponding values for individual parameters . for this reason ,
when more than one model parameter is to be determined , interferometry at a single projected baselength is unlikely to yield overall gain in accuracy in excess of factor of two even at optimal baselines .
we studied further the effect of increasing the number of sublayers @xmath239 , which permits to investigate the influence of fine details of the thermal structure of the envelope .
this appears to be justified for drew ( @xcite ) showed that in certain cases this structure is quite complex , and gradient of @xmath240 varies strongly with radial distance .
consequently , a model aiming to closely approximate the realistic envelope should be parameterized using a large @xmath239 ( see sect .
[ stpars ] ) and hence requires computations on grids of high dimension @xmath197 .
since the number of mesh points of the grid depends on @xmath197 exponentially , for computations based on such complex models ( grids 59 from table  [ grids ] ) the limited computer resources dictate sparser mesh points along each axis of the parameter space and a tuning of parameters controlling the computational process in such a way as to increase the speed of computation at the expense of precision .
although numerically less precise than computations for simpler models , our results for more realistic grids indicate that the conclusions derived for low @xmath239 can be safely extrapolated for higher @xmath239 .
this point is illustrated in fig .
[ 5l1 + 7 ] showing the results for single - parameter target sets for grid  7 .    as it can be seen in figs .
[ 3lmerr][5l1 + 7 ] and is confirmed by our other results not presented here , the distribution of optimal projected baseline lengths in the majority of cases has a maximum in the range 0.40.8 in dimensionless units . in some cases this maximum shifts to higher baseline lengths , but this occurs only for a very special choice of @xmath16 that correspond to measurements of wind temperature gradient at the distances of several stellar radii from the star with all other model parameters considered nuisance ( see e.g. bottom plots in figs .  [ 3l1 + 5 ] and [ 5l1 + 7 ] ) .
thus , the aforementioned range , corresponding to 4590  m in linear scale , can be recommended as the optimal choice for interferometric observations aimed at determination of global characteristics of the p  cyg wind .    the typical dependence of the robustness ratio + @xmath241 on the baseline length and the set of fixed parameters for a model with adjustable thermal structure is shown in fig .
[ rbst ] . at nearly all points of parametric space @xmath242 , and @xmath243 is generally lower when the set of fixed parameters @xmath63 includes the values @xmath149 . for a small fraction of models and for short projected baselines where interferometry can not significantly reduce random error , the value @xmath244 by a negligible amount .
thus , our results indicate that if the baseline length and the model used for interpretation are chosen in such a way as to reduce random errors , olbi appears also able to reduce the systematic errors .
using the theory of model parameter determination , we developed a general method which permits a quantitative comparison of observational techniques and optimization of complex observations through the relative figure of merit , defined as a generalization of the ratio of random errors of model parameters .
the method was applied to compare the outcome of olbi and classical spectroscopic observations , both used to determine the parameters of the outflow of the p  cyg star .
the observable quantities were computed using an efficient radiative transfer code and realistic grids of the envelope model .
the olbi and spectroscopic measurements errors were simulated assuming that the only source of errors is the counting statistics , i.e. we considered a nearly - ideal interferometer .
other main assumptions and simplifications of the present work were as follows : ( 1 ) the p  cyg wind is stationary and spherically symmetric ; ( 2 ) the only olbi observable considered here is the visibility modulus ; ( 3 ) only the hydrogen h@xmath1 line is considered .
the main conclusions of the present work are as follows :    1 .
a meaningful evaluation of the olbi vs spectroscopy relative figure of merit requires exploration of large domain of the model parametric space ( about @xmath0 points computed in the present work ) .
2 .   if the p  cyg outflow can be accurately described by a simple model of a stationary spherically symmetric isothermal envelope , then the use of interferometry does not substantially reduce the random error in the model parameter determination .
3 .   if the p  cyg outflow requires more complex models , e.g. with a priori unknown dependence of the envelope temperature on radial distance , use of a nearly - ideal interferometer can reduce errors in model parameters by an order of magnitude as compared to the spectrometer of the same collecting power and noise level .
the study of the relative robustness , i.e. the stability of the obtained model parameter determination with respect to systematic biases in the model , indicates that when the olbi and the physical model reduce random errors , they also tend to reduce the systematic error .
the optimal projected baseline lengths for observation of p  cyg with a nearly ideal interferometer lie in the range 4590 m. 6 .
if several model parameters are to be determined simultaneously , the olbi observations at a single projected baseline can not provide a substantial error reduction .
in this section we discuss to what extent our results are relevant for real observations , and how the methods developed here could be generalized to more complex cases .
although general formalism presented in sect .
[ statth ] remains valid for any source of errors provided that the linearized least - square method is applicable , eqs .
( [ cspec ] ) and ( [ colbi ] ) need to be modified if the accuracy of measurements is not shot noise limited .
in the simplest case of uncorrelated errors in spectral channels , this modification reduces to multiplication of the integrands in eqs .
( [ cspec ] ) and ( [ colbi ] ) by corresponding ratios ( shot noise error)/(total error ) . if this ratios are nearly constant across the observed spectral range , then after appropriate scaling our results can be directly applied to real observations .
the assumption of uncorrelated errors may be invalid for interferometric observations , the observables @xmath104 being obtained as a result of much more complicated processing of raw observational data ( cf .
mourard et al @xcite ) then that considered in [ olbidiag ] .
this point must be further analysed , and in future works , it may be necessary to use directly eq .
( [ covmateval ] ) instead of eq .
( [ colbi ] ) .    in any case
, a generalization on arbitrary measurement errors could be easily incorporated into existing code with no considerable increase in the required computer resources .      the theory developed in sect .
[ statth ] can be applied to the case of multiple observed lines without any modifications .
formally , it reduces to changing integration domains in eqs .
( [ cspec ] ) and ( [ colbi ] ) .    rapid evaluation of source function , on the contrary , would require a detailed study and development of efficient methods specific to each spectral line involved , because different lines can be emitted in different regions of the envelope and originate from different physical processes .
the method used here can be generalized to visibility measurements done simultaneously at more than one baseline : the only modification required is the replacement of rhses of eq .
( [ colbi ] ) by the sum of corresponding expressions over given set of projected baseline lengths , the directions of the baselines being unimportant because of circular symmetry of intensity distribution .
this can be done by modifying only the code performing the second stage of calculations ( see sect .  [ compfinres ] ) and using the data prepared for the single - baseline case during the most time consuming first stage of computation .
although time - averaged structure of p  cyg wind can be considered spherically symmetric , polarimetric ( taylor et al .
@xcite ) and interferometric ( vakili et al .
@xcite ) observations revealed existence of time - dependent `` clumpy '' structures in the envelope , giving rise to deviation from spherical symmetry .    in this case
, measurements of the fringe phase can yield important information .
this contrasts with the spherically symmetric case , where the fringe phase is irrelevant to physical modeling of the object .    on the other hand
, a unique determination of model parameters using only spatially unresolved spectroscopy appears impossible for the complex physical models resulting in two - dimensional brightness distributions , so that even the formulation of the comparison problem for spectroscopy and interferometry can be hardly done .
the choice of the optimal set of baseline lengths for interferometry still remains the important problem , and the method developed here can be used for such an optimization after certain modifications .
first , the fringe phases ( or observable linear combinations thereof ) should be included into the set of observables defined in eq .
( [ yolbi ] ) , and eq .  ( [ colbi ] ) should be modified accordingly .
this would require comprehensive error analysis depending on particular method used for extracting phase information from measurements .
second , errors in parameter determination for observations performed with interferometer configurations @xmath245 and @xmath246 ( here @xmath245 and @xmath246 denote corresponding sets of simultaneously used projected baseline lengths ) can be compared using the ratio + @xmath247 instead of @xmath248 .
the authors are grateful to a. labeyrie , d. mourard and f.vakili for encouraging discussion at the beginning of this project .
we thank c.bertout who provided the code for computation of the emergent intensities .
msb is grateful to the staff of observatoire de grenoble for the hospitality during his stay in grenoble .
this research made use of simbad astronomical database , created and maintained by the cds , strasbourg , and nasa s astrophysics data system abstract service .
part of this work was supported by cnrs through the grants pics @xmath249  194  france  russie , astronomie  haute rsolution angulaire " , gdr  milieux circumstellaires " , and gdr and pn  haute rsolution angulaire en astronomie " .
in this appendix , we compute the non - trivial elements of the covariance matrix of errors used in sect .
[ interferometry ] .
since we consider here a single spectral channel , the index @xmath8 will be omitted everywhere .
thus , the values @xmath252 , @xmath253 , @xmath254 , @xmath255 , @xmath256 , and @xmath257 in this appendix are identical with the values denoted as @xmath258 , @xmath259 , @xmath85 , @xmath260 , @xmath261 , and @xmath262 respectively in sect .
[ statth ] .    in practice
, the error in @xmath263 depends on too many details of experimental techniques , atmospheric conditions and data processing methods to be analyzed in general form ( mourard et al @xcite ) .
the present analysis is restricted to a highly idealized situation .
namely , as well as in sect .  [ spectroscopy ] , the only source of measurement errors taken into account is the photon shot noise .
we assume that the static fringes are detected by counting the photons in nonintersecting channels that cover in total the interval of length @xmath264 of optical path difference ( opd ) .
if z is the opd , which is measured here in units of @xmath265 , and @xmath266 is the photon counting rate per unit @xmath267 , then @xmath268 where @xmath269 is unknown fringe phase shift and @xmath252 is proportional to the total received flux .    if we define @xmath270 and @xmath271 then assuming that @xmath272 and neglecting the terms that decrease as @xmath273 for @xmath274 we obtain that @xmath275 and @xmath276 an estimate for @xmath277 can be constructed by replacing the values in the rhs of eq .
( [ wvalue ] ) by their estimates that can be easily deduced from eqs .
( [ fcos])([nolbi ] ) .
if the interval of opds covered by measurements is divided in @xmath278 subintervals of length @xmath279 and @xmath280 for @xmath281 , then estimate @xmath263 of @xmath255 based upon the number of photons received in each of these subintervals is given by @xmath282 where @xmath283 , @xmath284 is the number of photons received in the @xmath102-th subinterval , @xmath285 , and @xmath286 .    if we assume further that @xmath263 as a function of @xmath287 can be linearized in the vicinity of zero and take into account that @xmath288 and @xmath289 , then after some algebra we obtain for the statistical characteristics of @xmath263 @xmath290 which proves the validity of eqs .
( [ wwcovar ] ) and ( [ nwcovar ] ) for @xmath115 .
in the sobolev approximation , the escape probability for a line photon in the spherically symmetric accelerating outflow is given by ( see e.g. mihalas @xcite , sect .
14.2 ) @xmath291 where @xmath292 @xmath293 @xmath294 , @xmath131 and @xmath295 are the outflow and scattering atoms thermal velocities respectively , @xmath296 , @xmath297 is the line frequency , @xmath298 is the oscillator strength , and @xmath254 is the number density of atoms on the lower level of the transition .
here and in what follows we consider the fixed point of the envelope at the distance @xmath132 from the center and do not show explicitly the dependence of @xmath299 , @xmath107 , @xmath254 , and other physical parameters of the envelope on @xmath132 .
also , in expression for @xmath300 we neglected the stimulated emission , which do not introduce a noticeable error when the formation of h@xmath1 in the p  cyg envelope is considered .
our purpose is to find an approximation @xmath301 to @xmath299 that permits fast evaluation . in order
that in solving the equations of statistical equilibrium the approximation could be used instead of exact definition , the former should retain principal mathematical properties of the latter reflecting underlying physics .
namely , we require that @xmath302 and @xmath303 should monotonically decrease to zero as @xmath304 .
the simplest function that obey the above stated requirements is @xmath305 where @xmath306 is a parameter . to choose the optimal value of @xmath306
, we consider the behavior of @xmath307 for large @xmath254 .
when @xmath304 , the numerator of the integrand in eq .
( [ epdef ] ) can be replaced by 1 and asymptotic form of @xmath299 can be easily computed .
comparison of that asymptotic form with eq .
( [ epappdef ] ) shows that if we set @xmath308 then @xmath309 up to the terms of order @xmath310 when @xmath304 . here
@xmath311 and @xmath312 are statistical weights of the lower and upper levels of the transition respectively , @xmath83 is the line wavelength , and @xmath313 is the spontaneous emission coefficient .    to study the precision of approximation ( [ epappdef ] ) , we performed numerical computations for a wide range of values @xmath314 and @xmath315 ( it can be easily shown that the error is uniquely determined by the values of these two parameters ) .
our results show that the relative error reaches its maximum of 0.23 for @xmath316 , that is when @xmath317 , and @xmath318 .
the approximation found here is applicable if the escape probability is given by eq .
( [ epdef ] ) , that is if the frequency change in scattering is described by complete frequency redistribution .
of course , it is useful only if the usual conditions of applicability of the sobolev approximation are satisfied .
the expressions ( [ epappdef ] ) and ( [ nas ] ) can be generalized in two ways .
first , the error can be reduced if rational approximations of higher order are used instead of eq .
( [ epappdef ] ) .
second , a similar approximation can be obtained for a general three - dimensional flow as it will be showed in a future work .
p. , morand f. , bonneau d. , jiang y. , yang s. , guinan e.f .
, hall d.s . , mourard d. , hadrava p. , bozic h. , sterken c. , tallon - bosc i. , walker g.a.h . ,
mccook g.p . , vakili f. , stee p. , le contel j.m . , 1996 , a&a 312 , 879