we had two distinct objectives in undertaking this work , dominantly to investigate the merit and reliability of k giant stars as absolute fir and millimeter calibration sources , but secondarily to determine just how shallow in depth simple radiative equilibrium cool star model atmospheres might be trusted .
we have chosen @xmath0  tau and @xmath0  boo for this initial study because these stars , being both bright and nearby , are easily the two most completely and thoroughly studied red giants .
this makes them the ideal test cases for our purposes .
it should be noted that they both lie on the cool , non - coronal side of the  linsky - haisch dividing line " @xcite which separates coronal from non - coronal giants , although @xcite have found evidence for buried high temperature emission in the case of @xmath0  boo .
cohen and colleagues have presented a self - consistent context for absolute calibration in the ir , ideal for use by ground - based , airborne , and spaceborne spectrometers and radiometers @xcite .
their approach is based upon a pair of absolutely calibrated , ir - customized models of vega and sirius calculated by kurucz .
these efforts have furnished the absolute stellar spectrum of sirius that underpins cobe / dirbe bands @xmath3 @xcite , the spectra of k0-m0iiis for the on - orbit calibration of the near- and mid - infrared spectrometers on the joint nasa / isas infrared telescope in space @xcite , likewise for esa s infrared space observatory ( iso ) @xcite instruments ( e.g. , the short wavelength spectrometer ) @xcite , for the spatial infrared imaging telescope ( spirit - iii ) on the us midcourse space experiment ( msx ) @xcite , and for 2mass @xcite . an extension of this work to substantially fainter stars has provided calibrators for the _ spitzer _ space telescope @xcite .
most recently , @xcite have absolutely validated the two calibrated model spectra , and the empirical spectra of the brighter k giants such as @xmath0 tau and @xmath0 boo , in six bands in the mid - infrared ( mir ) , using precision measurements from the msx .
the absolute calibration of isophot @xcite was supported by planets , asteroids , and stars , ranked in order of decreasing flux density @xcite .
the stellar calibrators @xcite necessitated long - wavelength extrapolations , as far as 300  @xmath2 m , of the observed ( 1.235  @xmath2 m ) absolute stellar spectra by model atmosphere spectra ( computed by one of us [ dfc ] : @xcite ) .
we note that the bright , secondary reference stars that have been selected for ir calibration purposes are precisely those  quiet " giants in which it was suggested that radiatively - cooled regions very largely dominate the stellar surfaces so that single component atmospheric models were most likely to be valid .
examples of these are @xmath0  boo , @xmath0  hya , @xmath0  tau , and @xmath4  dra @xcite .
consequently , one of our objectives is to establish whether these stars radiate as predicted , thereby validating the scheme used for iso at the longest wavelengths , or possess long - wavelength chromospheres .      the iso long wavelength spectrometer ( lws )
@xcite used uranus as a calibration standard but also observed neptune , and mars when it became available .
although uranus has remained the lws primary calibrator , the mars spectrum has been shown to be fully consistent with this uranus framework @xcite . in the case of neptune ,
a new model based on the lws data has been adopted @xcite .
lws was able to observe some stellar calibrators @xcite , although these observations were difficult to make due to poor signal - to - noise in the relatively narrow lws bandwidths .
however , the stellar spectra did demonstrate consistency with the uranus - mars framework , in the lws 45170  @xmath2 m region . in this paper
we present the totality of results from the iso lws spectra obtained on @xmath0  tau and @xmath0  boo .
the long - wavelength stellar spectral extrapolations are based upon the assumption that these calibrator stars are not attended by extensive chromospheres .
radiometric testing with isophot @xcite suggests that these extrapolations were meaningful , i.e. consistent with calibrations based on asteroids and planets .
it would be of great interest , both scientific and pragmatic , to determine independently whether these same stellar calibrators radiate as predicted in the millimeter ( mm ) domain .
the basis for mm - wave absolute flux calibration is mars .
therefore , we observed the stellar calibrators against mars , or against secondary microwave standards that were ultimately traceable to a comparison with mars , uranus , or neptune . by this
means we hoped to provide a link between the reference flux density scales of planets and stars and to unify flux calibration across the ir - mm regime .
early efforts in this direction at berkeley were described by @xcite .
@xcite and @xcite have continued this radio absolute calibration work .
our secondary goal is to probe the structure of the outer atmospheres in k and , eventually , m giants .
@xcite demonstrated in a classic paper how observations of the solar fir continua , atomic lines , and uv continua could be brought together to construct a single - component , time - independent mean solar temperature structure .
fir continuum fluxes were of crucial importance in helping to establish this mean structure from the upper photosphere , through the temperature minimum region , and into the chromosphere .
using the fir continuum to probe the shallow depths in solar and cooler stars is possible because the primary ir continuous opacity , h@xmath5 free - free , increases as @xmath6 , pushing the depth of continuum formation into the outermost stellar layers . until now
these layers in cool stars have been accessible only through non - lte analyses of the cores of strong atomic and molecular lines ( e.g. , @xcite , @xcite , @xcite , @xcite ) .
while chromospheric models have become more complex in the years since the @xcite study , the fir and millimeter fluxes still have the potential to be valuable diagnostics of the shallow stellar layers . in accurately determining the fir and millimeter fluxes of selected red giants , we hope to provide useful information for testing models of the cool star stellar photospheres and chromospheres .
we are particularly interested in determining the location of the transition from observed fluxes , reproducible by purely radiative equilibrium models , to fluxes that require invoking a non - radiatively heated chromosphere .
in addition , we wish to trace the observed run of decreasing fluxes with increasing wavelength from the upper photosphere and determine how closely current stellar model atmosphere calculations can match this run .
this information on the location of the temperature reversal and the quality of fit to the fluxes originating in the upper photosphere are vital in defining the wavelength range for which simple radiative models may be used for fir calibrations .
table  [ obs ] summarizes our various sets of mm observations , listing array configurations and frequencies used , integration times on the stars , phase and flux calibrators ( chosen to be as close to each star s direction in the sky as possible ) , and the primary or secondary mm - calibrator observed with each set of stellar observations .
a few measurements were taken in the b - array , with which planetary calibrators like mars were partially resolved , necessitating an extra step in the calibration process to calculate the flux expected on the shortest baseline pairs ( see below ) . consequently , we preferred measurements in c- or d - array configurations .
angular diameters of these stars in the optical to mir regions are of order 20  mas , rendering them unresolved targets in all these configurations , as were the local calibrators chosen . on the basis of estimates of the bima performance at 1 and 3  mm , our observing scripts called for typical tracks of several hours duration , of which between roughly two and four hours were spent on the stars themselves .
this was a continuum experiment , and we selected observing frequencies from considerations of the expected stellar spectra , receiver performance , and atmospheric transmission .
whenever possible we included measurements of mars with each set of stellar and local calibration source observations , in an attempt to unify the calibration .
for those few data sets taken in an array configuration that resolved mars , we calculated visibilities as a function of baseline ( in kilo@xmath7 ) and applied these partial - resolution factors to the mars model @xcite flux densities using our shortest baseline pairs .
sometimes no planet was visible close in time to our observing tracks .
then we used mwc  349 when it was available .
hat creek has a long series of observations of this star , it is relatively bright , has shown no statistically significant variability in our archives , and has been modeled in detail by @xcite so that flux extrapolation and interpolation to our frequencies is viable .
further , @xcite have obtained its mm flux densities to a precision of 1% absolute . for our final observations , in december 2000
, we used w3oh as reference for our @xmath0  tau track , again because of its relative proximity to this star , long - term temporal stability of signal , and the available multifrequency hat creek archives which enable interpolation between known flux densities at other frequencies to our own .
we treated w3oh as having a canonical hii region spectrum and included the frequency dependence of the gaunt factor .
other secondary references such as 1415 + 133 and 0530 + 135 were likewise observed when no other obvious flux reference presented itself .
independent observations of these were extracted from the hat creek archives , preferably with respect to a primary or acceptable secondary standard . in such circumstances , table  [ obs ] carries the primary standard to which the flux pedigree of such a secondary is traceable .
upper and lower sidebands were processed independently , maps made of the stars at each frequency , and peak flux densities extracted ( see the final column of table  [ obs ] ) . sometimes , weather conditions were adequate in the lower sideband but led to unacceptably large noise in the upper sideband .
in such cases we eliminated the noisy datasets which were readily recognized by the presence of artifacts in the maps , such as global striping , or large negative - going blobs . to estimate the uncertainties in the tabulated peak flux densities derived from any map
, we used the standard deviation calculated from substantial regions of the map that excluded the target star itself . to combine the simultaneous upper and lower sideband maps at one epoch
, we used inverse - variance weighting of the two maps based upon these noise estimators , relying upon the miriad task maths .
likewise , we combined results across epochs .
phase calibration of interferometer data relies on transferring phases accurately from a point source ( usually a qso ) to the target source .
the accuracy of the calibration depends on the distance between source and calibrator , baseline errors , atmospheric phase stability , and the signal - to - noise ratio ( snr ) in the phase measurement of the calibrator .
all sources of phase error result in decorrelation of the calibrated image ; i.e. , a lower amplitude in the final map compared with perfectly calibrated data . self - calibration offers a way to eliminate these phase errors , but is limited by the timescale on which phase can be accurately measured on the target source itself .    in the case of @xmath0
boo , the measured peak flux density in the normally calibrated image is about 50  mjy .
given the system temperature ( 350  k ) , bandwidth ( 670  mhz , reduced from 800  mhz due to edge - channel flagging in the wideband average ) , and a correlator efficiency of 0.88 ( 2-bit ) for these observations , the expected noise in a single - interferometer measurement is about 70  mjy in 10  min . using self - calibration with a 9-element array improves the snr for antenna - based phases by about a factor of 3 , due to baseline averaging .
the expected r.m.s .
noise in the antenna - based phases using self - calibration is , therefore , about 25  mjy in 10  min . assuming the true flux of @xmath0  boo to be 80  mjy
, we would expect an snr of about 3 using 10-min averaging on the source .
the expected phase error for an snr of 3 is @xmath8 for each antenna ( arcsin(1/3 ) ) .
the decorrelation expected for an r.m.s .
phase error of 20@xmath9 is about 6% ( reduced amplitude = exp(@xmath10 ) , where @xmath11=the r.m.s .
phase error in radians ) .
it would , therefore , appear safe to use an interval of 10 min or longer to self - calibrate the 1-mm @xmath0  boo data in phase .
of course , any additional phase error ( e.g. , atmospheric ) incurred during the self - calibration interval will increase the decorrelation in the final image .
we can estimate the amount of atmospheric decorrelation expected .
the average r.m.s .
delay path in 10  min for a 100  m baseline on the day of the @xmath0  boo 1-mm observations ( 01 jun 2000 ) was 210  @xmath2 m , or @xmath121 rad at 1360  @xmath2 m wavelength . in d - array ,
the maximum baseline is about 30  m. assuming that the atmospheric r.m.s .
phase scales as ( baseline)@xmath13 , the r.m.s .
phase becomes 0.37  rad on the longest baseline [ ( 30  m/100  m)@xmath13  *  1  rad ] , and about 0.2  rad for the median baseline of 15  m. the expected atmospheric decorrelation in 10  min for the 30  m baseline ( @xmath11=0.37  rad ) is 0.85 , and for the 15  m baseline ( @xmath11=0.2  rad ) is 0.98 .
thus , atmospheric decorrelation on 10-min timescales is not expected to be very important , and phase variations on the longer timescales will be removed by the self - calibration .
we validated these ideas on the june 2000 observations of @xmath0  boo by experimenting with self - calibration timescales in the reduction of the lower side - band data .
the total cleaned flux in the resulting images of this star showed a spread of only 3.5% , while the more relevant quantity , the stellar peak flux densities , were 82.1 , 81.5 , and 81.1 mjy for a 15- , 20- , and 30-min timescale .
we settled on a 20-min self - calibration timescale to reduce the data for @xmath0  boo , and 30  min for @xmath0  tau .
the final rows in table  [ obs ] present the combined peak flux density information based on the totality of the available , good quality , stellar data for @xmath0  tau and @xmath0  boo .
the corresponding combined images from which these final flux densities were determined appear as figure  [ bima13 ] .
observations of @xmath0  tau with full wavelength coverage were made in three revolutions : 818 ( tdt 81803302 ) , 848 ( tdt 84801701 ) , and 861 ( tdt 86101101 ) .
the observations of @xmath0  tau each consisted of 28 scans , 14 in each direction with a spectral oversampling of 4 .
each observation had about 4630  s integration .
lws made two aot observations of @xmath0  boo in revolutions 448 ( tdt 44800303 ) and 608 ( tdt 60800904 ) , each with about 2300  s integration .
these observations consisted of 12 scans , 6 in each direction , taken with an oversampling of 4 .    each lws spectrum was processed with offline processing version 10 ( olp10 ) and all the data were reduced using the infrared spectral analysis package ( isap ) .
the reduction steps consisted of removing outliers and averaging over individual detectors at maximum resolution , clipping at 2.5@xmath14 .
no attempt was made to scale the detectors to make the sub - spectra overlap in the reductions .
we accomplished this outside isap using the splicing techniques described by @xcite for constructing absolutely calibrated spectra .
care was taken to use the same number of scans in each direction , so the effects due to transients in the data are minimal , and no transient correction was applied . the fundamental uranus model for calibrating
the lws is a synthesis of globally averaged voyager infrared interferometer spectrometer ( iris ) data up to 50  @xmath2 m and james clerk maxwell telescope data covering the 0.352  mm wavelength range @xcite with a radiative transfer model similar to the one described by @xcite .    a formal off - source , full - grating
, lws integration is available only for @xmath0  boo ( tdt 60801005 ) , and two fixed - grating ( sampling only ten wavelengths with no detector scans ) observations exist ( tdts 21701113 and 27202302 ) .
however , it is vital to recognize that corrections for off - source emission in lws are inappropriate for faint objects like normal stars in low - cirrus regions .
all observed lws spectra have had the instrumental dark currents subtracted from them during the reductions . in our olp10 reductions , we used the  fixed dark currents " , as determined in dedicated calibration observations .
these lws dark current measurements were essentially dark background measurements .
the fabry - perot was placed in the beam and offset . therefore , there was flux on the detectors , radiated by the 3.5-k blank , but this was negligible .
the signals measured remained constant throughout the mission , and were subtracted from the measured photocurrents on - source .
when the dark signal , measured from the blank , was compared with the signal measured on the dark patch of sky , no difference was found .
therefore , for dark sky , no additional  sky " subtraction is necessary . for regions with high cirrus
, there is a sky signal in the lws beam and , if a point source is being observed , an  off " position is needed to determine the local sky background . in dark regions ,
the sky flux predicted by cobe can form a substantial fraction of the lws dark signal . however , this is deemed undetected , because no additional signal is detected over the signal from the blank .
consequently , we must assess the measured sky background around the two k  giants at those times when the lws spectra were taken . if those backgrounds were below the fixed dark values ,
then there is no need for further subtraction of off - source lws spectra . indeed , that subtraction would produce meaningless stellar spectra by removing the dark currents twice .
both lws and isophot , off - source , measure the sum of zodiacal ( hereafter  zodi " ) , cirrus , and extragalactic emission , as did dirbe .
 irsky " shows that extragalactic and cirrus confusion noises are far below zodi , so the dominant background for the two k  giants is zodi .
we have used two methods to evaluate the fir background near @xmath0  tau and @xmath0  boo : actual measurements in off - source pixels when these stars were observed in the fir arrays of isophot ; and , as a coarse check , dirbe predictions for the lws observations based on their solar elongation angles .
the computation of isophot background sky measurements near @xmath0  tau and @xmath0  boo is summarized in tables  [ ataupht ] and  [ aboopht ] .
these present the relevant isophot observation by filter , nominal wavelength , tdt and measurement number , and the background brightness and its 1-@xmath14 uncertainty .
our conversion from point source flux to surface brightness is different from the one given in the isophot handbook @xcite .
the supporting analysis and justification are given by schulz .
when multiple data exist for a filter , inverse - variance weighted combinations of those data are also given .
conversion from mjy  sr@xmath15 to f@xmath16 ( in w  @xmath17  @xmath2m@xmath18 ) was accomplished using the published values for lws detector beam solid angles , effective apertures , and bandwidths .
table  [ bgnd ] presents these backgrounds , and compares them with the fixed dark currents and dirbe predictions ( interpolated between 140 and 240  @xmath2 m to assess the 160-@xmath2 m values ) , matching each isophot band to the nearest lws detector .
the dirbe predictions are generally in accord with isophot measurements
. clearly , the measured sky backgrounds around the two stars are all below the fixed dark currents , confirming that only low level cirrus appears in their vicinities ( with the sole exception of @xmath0 boo in lw2 , where isophot suggests 50% more sky brightness than the lw2 dark current ) . the dirbe predictions are generally in accord with the isophot measurements .
_ subtraction of off - source lws spectra would be inappropriate for these k  giants_.    figures  [ ataulws ] and  [ aboolws ] illustrate our reduced , spliced lws spectrum for each star , presenting the mean spectrum over all tdts , with the @xmath19 bounds for @xmath0  tau ( which has more lws data , longer integrations , and higher snr ) and @xmath20 bounds for @xmath0  boo ( with its noisier data sets ) .
the sizeable triangular excursions in both spectra near 53 and 106  @xmath2 m represent the difficulty , for nir - bright normal stars , of mitigating the substantial leaks of 1.6-@xmath2 m radiation to which lws is prone .
we have cleaned the observed spectra by 30-point boxcar smoothing of the data ( equivalent to about a 2-@xmath2 m interval ) . @xmath0
tau s combined spectrum contained 1181 wavelength points ; @xmath0  boo s 706 points .
each figure incorporates the absolutely calibrated model continuum spectra that we have computed ( see below ) as the almost flat long - dashed lines that cross the figures , and the mean@xmath21 bounds on our older calibrated energy distributions ( dash - triple - dotted lines ) delivered in january 1996 for the calibration of isophot . to convert the model surface fluxes to predicted fluxes at the earth
, we have adopted the angular diameter , 20.88@xmath220.10 mas , of @xcite ( using lunar occultation ) for @xmath0  tau and the angular diameter , 21.0@xmath220.2 mas , of @xcite for @xmath0  boo ( measured by optical intensity interferometry ) .    our new , more detailed , calculations of stellar continuum spectra fall well within the uncertainties associated with the calibrated spectra originally delivered to iso .
we wish to quantify the degree of accord between measured and modeled spectra in the fir for each star .
the grating simultaneously illuminated the ten lws detectors , and each detector was separately wavelength and flux calibrated .
therefore , lws consists of ten spectrally independent detectors .
the degree of oversampling of these stellar spectra was 4 ; consequently , adjacent points in figures  [ ataulws ] and  [ aboolws ] are not independent .
therefore , to assess the consistency of observations and predictions , and to handle the dependence of adjacent lws points in the fully - plotted spectra correctly , we include the data binned by detector as a series of bold crosses in figures  [ ataulws ] and  [ aboolws ] .
filled circles represent the inverse - variance weighted average of all the meaningful data in each lws detector ( e.g. after the rejection of data contaminated by @xmath23-band leaks ) .
error bars show the associated 1-@xmath14 uncertainties .
asymmetric horizontal error bars represent the actual wavelength ranges usable for each detector and star .
filled circles are plotted at the average wavelength of each detector s set of usable data for that star .
the binned detector data lie within 1@xmath14 of the newly modeled energy distribution for @xmath0  tau , and within 1@xmath14 for @xmath0  boo .
there is good agreement between the underpinning uranus calibration of lws , the totally independently calibrated stellar spectra made from the empirical spectra by our older model extrapolations based on the available observations between 8 and 23  @xmath2 m , and the newly calculated , absolutely calibrated , emergent stellar energy distributions . averaging the results for the independent detectors for each star with inverse - variance weighting
, we determine that , for @xmath0  tau , this mean scale factor is 1.000@xmath220.011 ( over nine detectors ) while , for @xmath0  boo , this factor is 1.012@xmath220.025 ( over seven detectors ) .
therefore , empirical stellar and planetary calibrations used by iso , and stellar modeling , are absolutely reconciled to 3.3% ( 3@xmath14 ) for the better measured @xmath0  tau , and 7.5% ( 3@xmath14 ) for the more poorly measured @xmath0  boo .
the spectrum calculations were carried out with the source model atmosphere program in essentially the same fashion as described in @xcite .
specifics pertinent to the current calculations are described in this section .    for @xmath0  tau
, we adopted an effective temperature of 3920  k @xcite and a log(g ) of 1.5 @xcite .
for the isotopic abundances of c , n and o , we adopted log  @xmath24(c )  =  8.40 , log  @xmath24(n )  =  8.20 , log  @xmath24(o )  =  8.78 , @xmath25c/@xmath26c  =  10 , @xmath27o/@xmath28o  =  562 , and @xmath27o/@xmath29o  =  475 @xcite .
the remaining metal abundances were chosen to be solar , following the lead of @xcite .
the preceding abundances were derived using the moog spectrum synthesis code @xcite and model atmospheres @xcite constructed with the atlas ( e.g. , @xcite ) model atmosphere code and opacity sampling .
our choices of t@xmath30 , log(g ) , and abundances are entirely consistent with current determinations of the parameters of @xmath0  tau . in particular , the values are well within the error bars of the values recently derived directly by decin and her collaborators @xcite from fits to stellar spectral energy distributions generated by using the marcs @xcite model atmosphere code .
bertrand plez @xcite kindly computed an @xmath0  tau model atmosphere for us using the aforementioned parameters and the programs described in @xcite .
thus , this model was computed with essentially the same model atmosphere program and opacities as employed by decin .    for @xmath0  boo
, we adopted a slightly revised version ( @xcite : hereafter pdk ) of the @xcite model atmosphere which has more depths than the original published model .
this atmosphere has an effective temperature of 4300  k , log(g ) of 1.5 , [ fe / h ]  =  -0.5 , and other elemental abundances as determined by @xcite .
the pdk model parameters derived by @xcite were derived entirely from the observed stellar spectrum . as was the case for the decin et al .
investigation of @xmath0  tau , no assumptions were made regarding stellar temperature or gravity .
we adopted all of the @xcite values for our study .
the parameters t@xmath30 , gravity , microturbulent velocity , and [ fe / h ] derived by @xcite for @xmath0  boo are essentially identical to those derived later and independently by @xcite , while the c , n , and o abundances agree within the error bounds .
the alert reader will observe that we have adopted model atmospheres for @xmath0  tau and @xmath0  boo that were derived using entirely different model atmosphere codes ( atlas and marcs ) .
we shortly will describe results obtained by using these atmospheres as input to a third model atmosphere code , our source .
we do not feel that significant improvement would result from insisting that _ only _ atlas _ or _ marcs models be used in our study .
all three codes have been vetted against the solar case and should give very comparable results in the g - k star domain .
we cite as support for this stance the striking agreement between the results of @xcite and @xcite for @xmath0  boo .
this occurs despite @xcite having used a completely different model atmosphere program and blanketing representation ( marcs with opacity sampling ) to compute their @xmath0  boo model than used by @xcite ( atlas with opacity distribution functions ) .
the agreement is even more impressive when one considers that these two investigations derived their results from two entirely different spectral regions , 0.50.9  @xmath2 m _ versus _
2.3812.0  @xmath2 m .
we regard this as evidence that different well - vetted model atmosphere codes and opacity representations now produce comparable model atmospheres and visible - infrared spectra for k - giants * _ within the limited approximations _ * of static , homogeneous atmospheres , lte , and mixing - length - theory radiative / convective equilibrium .
that does * _ not _ * mean that the resultant models are accurate , or even adequate , representations of a real star .
differences between observed and computed spectra are now much more likely to arise from the presence of physics ( _ e.g. _ , nlte , non - radiative heat deposition , and inhomogeneous dynamical structure ) * _ not treated _ * in the simple models we use here .
we will discuss this point in more detail as we proceed .
the column mass ( g  @xmath17 ) _ vs. _ temperature relations for both the plez @xmath0  tau and pdk @xmath0 boo models were linearly extrapolated to shallower depths as required to obtain optical transparency at the longest wavelengths .
for the plez @xmath0  tau model , the extrapolated layers are shallower than 0.020 g  @xmath17 ; for the pdk @xmath0  boo model , shallower than 0.010 g  @xmath17 .
the resultant model atmospheres are shown in figure  [ atau_aboo_models ] .
it should be stressed that both of the above adopted model atmospheres have purely radiative equilibrium surface temperature structures .
specifically , they do _ not _ have chromospheric temperature rises , a point whose importance will become evident later .
we have included in figure  [ atau_aboo_models ] , for comparison , the @xcite ( kel ) @xmath0  tau chromospheric model b and the @xcite ( al ) @xmath0  boo chromospheric model .
these chromospheric models are semiempirical , derived from non - lte analyses of and line fluxes . while these chromospheric models are now over 20 years old , they are very similar from the upper photosphere , through the temperature minimum region , and into the lower chromosphere to models of @xmath0  tau and @xmath0  boo used by @xcite , and the @xmath0  tau model derived by @xcite .    using the plez @xmath0  tau atmosphere and parameters detailed above ,
we have computed the synthetic spectrum for @xmath0  tau over the entire lws range and beyond , 43 to 665  @xmath2 m .
the spectrum was computed using a full set of molecular line opacities which includes all the principal isotopomers of co , sio , oh , and @xmath31 and using a variable wavelength mesh which locally provides @xmath322.5 points per doppler width . to allow comparison with the lws observations , we convolved the high - resolution result with a gaussian instrumental profile function having a fwhm equal to 0.29  @xmath2 m for @xmath33 90.5  @xmath2 m , and 0.60  @xmath2 m , @xmath34 90.5  @xmath2 m ( e.g. , @xcite ) . the resulting spectrum , normalized to the continuum and on a greatly magnified ordinate scale , is shown in figure  [ atau_spectrum ] .
a detailed description of the displayed spectrum may be found in @xcite .
we do note here that almost all the absorption beyond @xmath32  50  @xmath2 m is due to the isotopomers of sio . at shorter wavelengths
the principal contributor is oh with some h@xmath35o and co as well .
it is clear that , _ at the iso lws and the bima spectral resolutions _ , all differences between the pure continuum and the synthetic spectrum computed with the full set of molecular line opacities are less than 2% ; beyond 130  @xmath2 m , they are less than 1% .
this is _ much less _ than the noise level in the observed @xmath0  tau spectra in this region . as a consequence ,
there is no significant advantage gained by using the full synthetic spectrum . for clarity in our presentation
, we will use only the continuous spectrum in our comparisons with the @xmath0  tau observations .
since @xmath0  boo is both hotter and has lower metallicity , we shall use the continuum alone for comparisons there as well .
for both the plez @xmath0  tau and the pdk @xmath0  boo models , the dominant continuous opacity in the continuum - forming layers for wavelengths redward of the h@xmath5 opacity minimum at 1.64  @xmath2 m arises from h@xmath5 free - free .
this opacity increases in strength as @xmath6 . as a result
, the emergent continuum flux is formed at progressively shallower layers in the stellar atmosphere as one proceeds to longer wavelengths .
figure  [ atau_aboo_models ] illustrates this point for the plez @xmath0  tau and the pdk @xmath0  boo models . in this figure
we have marked , for the radiative equilibrium models , the atmospheric depths for which the continuum optical depth , @xmath36 , reaches unity at the indicated wavelengths .
that is , we have marked the column masses at which @xmath36  =  1.0 ; roughly speaking , in the fir , depths near this value and shallower make the greatest contribution to the emergent continuum flux at wavelength @xmath7 .
note , for example , that the continuum fluxes between 1  mm and 3  mm are sampling very different parts of the atmospheric temperature structure in both radiative equilibrium models than , say , the continuum fluxes at 2.3  @xmath2 m .
indeed , the longest wavelengths are sampling , _ in the continuum _ , atmospheric levels reached only at the cores of lines in bluer portions of the spectrum .
only when one proceeds to the vacuum ultraviolet and its strong photoionization continua of metals is it possible to reach in the continuum the shallow atmospheric levels probed by the continuum at long wavelengths ( e.g. figure 1 of @xcite ) .
thus , the fir continuum is an effective probe of the shallowest stellar layers which can supplement information derived from line cores and uv continua .
as noted above , there is very little difference between the full synthetic spectrum and the spectrum of the continuum .
figure  [ atau_aboo_models ] provides an important clue as to why spectral lines generally become less prominent at longer wavelengths .
it is evident that , as one goes to shallower depths in the models , the slope of the temperature - column mass relation becomes smaller .
the significance of this can be appreciated by considering an absorption line in the @xmath0  tau model with a characteristic , depth - independent , line core / continuum opacity ratio of 1000 .
if the absorption line occurs at a wavelength of 2.3  @xmath2 m , the continuum will be formed at @xmath325200  k and the line core will be formed at @xmath323200  k. the same line at a wavelength of 200  @xmath2 m will have the continuum formed at @xmath323100  k and the line core at @xmath322600  k , a much smaller temperature difference between line core and continuum than at 2.3  @xmath2 m .
couple this with the decreasing temperature sensitivity of the planck function as one goes to longer wavelengths and one has a recipe for progressively less prominent spectral features .    figure  [ mass_bt ] provides another window into the wavelength dependence of continuum formation . in the upper panel of figure  [ mass_bt ]
we show the column masses for which @xmath36  =  1.0 over the spectral range 0.5  @xmath2 m to 200  @xmath2 m . as noted above ,
there is a steady march to progressively shallower depths of continuum formation as one goes to longer wavelengths . since the dominant h@xmath5 bound - free opacity peaks at roughly 0.825  @xmath2 m , continuum fluxes at wavelengths between 0.825  @xmath2 m and 1.64  @xmath2 m probe the same atmospheric layers as continuum fluxes between 1.64  @xmath2 m and 6.5  @xmath2 m .
longward of 6.5  @xmath2 m , the continuum flux probes atmospheric layers accessible only in line cores and in the metal bound - free continua of the vacuum ultraviolet .    in the lower panel of figure  [ mass_bt ]
we show the continuum  brightness temperature " over the same spectral range for both models .
the brightness temperature is defined as the temperature of a blackbody which gives the same flux as the model atmosphere at the indicated wavelength ( a precise definition is given in section  [ allbt ] ) .
three points are important .
first , neither @xmath0  tau nor @xmath0  boo could have its continuum fluxes over any large spectral range matched by a blackbody of a single temperature .
second , the blackbody temperatures that do match the model fluxes at specific wavelengths are generally quite different than the effective temperatures of the models ( or the stars ) themselves .
third , comparing with figure  [ atau_aboo_models ] , one sees that the brightness temperature is systematically smaller than t(@xmath36  =  1.0 ) , the temperature of the atmospheric layer at which optical depth unity is reached in the continuum at wavelength @xmath7 .
( for an atmosphere with a linear source function at wavelength @xmath7 , the brightness temperature at @xmath7  and
t(@xmath36  =  2/3 ) would be identical .
unfortunately , in the particular cases of the pdk and plez models , the fir source functions are not linear in @xmath36 but , rather , quadratic . ) the differences between t(@xmath36  =  1.0 ) and brightness temperature vary systematically with wavelength .
the differences are greatest at 1.64  @xmath2 m and decrease to the red . _ beyond 40  @xmath2 m , the differences are always less than 100  k. _ ( had we chosen for illustration the more correct , but _ still _ approximate , t(@xmath36  =  2/3 ) , the temperature differences would be smaller yet . )
thus , the brightness temperature is a good first guess at t(@xmath36  @xmath32  1.0 ) at the longer wavelengths .
later , when we discuss differences between observed and predicted brightness temperatures , the reader may use this approximation and figure  [ atau_aboo_models ] to roughly estimate the implied errors in the radiative equilibrium model temperature structures .
table  [ compare ] shows the comparison between the observed fluxes and the predictions of the radiative equilibrium models . only for the case of @xmath0  tau at 1.4  mm does there appear to be approximate agreement between the observed flux and
the flux predicted from the radiative equilibrium model .
the remainder of the observed fluxes are all significantly _ larger _ than the fluxes predicted by the models .
the discrepancies are particularly striking in the case of @xmath0  boo .
we do not believe that the flux computation for the radiative equilibrium models could have errors sufficiently large to produce the marked differences between computed and observed millimeter fluxes .
our model atmosphere code has been carefully vetted and the dominant continuous opacities are themselves quite well determined ( e.g. , @xcite ) .
( as an example of evidence for this , we refer the reader forward in the text to figure  11 .
there we show the continuum for the plez @xmath0  tau model computed with an entirely independent model code . ) on the observational side , our derived errors for the observed 1.4-mm and 2.8-mm fluxes preclude dismissing the differences between observed and computed fluxes as simply due to measurement error .
we believe other explanations must be considered .    since the adopted angular diameters for these two stars
appear to be extremely well determined , with standard deviations of @xmath321% or less , they would not seem to be a significant source of error .
moreover , a major change in the angular diameters would compromise flux agreement for the @xmath33 200  @xmath2 m region .
nevertheless , it should be noted that the @xcite diameter is an arithmetic mean of values measured from 0.45  @xmath2 m to 2.2  @xmath2 m and the @xcite diameter is a error - weighted mean biased toward values obtained at 1.6  @xmath2 m and 3.8  @xmath2 m . at any particular monochromatic wavelength
, an observed angular diameter represents the apparent diameter of the @xmath37  @xmath32  1.0 surface of the star .
since , as shown in figures  [ atau_aboo_models ] and  [ mass_bt ] , there is a progressive march in depth of continuum formation to ever shallower atmospheric layers as one moves redward of 1.64  @xmath2 m , the angular diameters themselves should increase with wavelength . if the angular diameters at 1.4  mm and 2.8  mm are substantially larger than those at the wavelengths of the @xcite and @xcite observations , our predicted fluxes in table  [ compare ] could be seriously underestimated .
we can roughly estimate the magnitude of this source of error , at least in the case of the radiative equilibrium models . in table
[ extents ] we show the physical radii of @xmath0  tau and @xmath0  boo based on the indicated angular diameters and hipparcos parallaxes @xcite .
we also show the physical distances in the radiative equilibrium models between @xmath38  =  1.0 and @xmath39  =  1.0 .    given the tabulated numbers ,
it is clear that , for the pdk and plez models , the atmospheric extent in the continuum at 2.8  mm is not sufficient to seriously affect our predicted fluxes .
we find a similar result in the case of the only model chromosphere with a published height scale , mcmurry s chromospheric model for @xmath0  tau @xcite .
the mcmurry atmosphere is plane - parallel and hydrostatic , constructed with log g = 1.25 .
note that this gravity leads to a more distended atmosphere than does the higher gravity we adopted .
the mcmurry atmosphere has an extent of @xmath322.6  @xmath40  10@xmath41  cm between the layers where the 1.64  @xmath2 m continuum is formed and the layer with t  =  7120  k in the chromosphere .
this extent is only 8% of the stellar radius and again seems too small to explain the discrepancy between computed and observed fluxes .
finally , as we will describe in section  [ allbt ] , @xcite determined that @xmath0  boo has essentially the photospheric radius at 2  cm .
given the expected increase in angular diameter with wavelength , this would imply that @xmath0 boo would have the photospheric radius at @xmath42  2  cm as well .
based on the above arguments , we believe that the flux excesses observed at 1.4  mm and 2.8  mm , relative to the predictions of the purely radiative models , are real and significant . furthermore ,
as we shall show , there is good evidence that these excesses are signatures of the chromospheric temperature rise in the two red giants .
figures  [ ataummcm ] and  [ aboommcm ] present our millimeter flux densities ( filled circles ) and the color - corrected iras measurements of each star at 12 , 25 , 60 , and 100  @xmath2 m ( crosses ) .
each figure includes the plez @xmath0  tau or pdk @xmath0  boo spectrum , as appropriate , from 10  @xmath2 m to 3  mm .
there are additional high - frequency radio continuum observations of these two k giants in the literature .
@xcite made single - dish 86-ghz measurements of @xmath0  boo ( 21.4@xmath227.5 mjy ) with the iram 30-m telescope , while @xcite likewise detected both @xmath0  tau and @xmath0  boo at 250 ghz , obtaining 51@xmath226 and 78@xmath228 mjy , respectively .
there are also vla centimeter measurements on both stars at 4.9 , 8.4 , and 14.9 ghz by @xcite , and unpublished data by drake , linsky , & judge from the compilation of radio data by @xcite .
these long wavelength measurements all have substantial uncertainties , in part due to the faintness of the stars , in part to the potential for confusion with extragalactic background objects .
nonetheless , in combination with our own radio data ( figures  [ ataummcm ] ,  [ aboommcm ] ) one sees that these high - frequency data ( open diamonds ) are in good agreement with our aperture synthesis measurements ( filled circles ) .
it is difficult to fully appreciate the significance of the excess fluxes from the logarithmic scales of figures  [ ataummcm ] and  [ aboommcm ] .
an insightful alternative view may be obtained by combining the mm and cm fluxes from other observers with our data and looking at the whole dataset as brightness temperatures over a broad spectral range . in figures
[ atauallbt ] and  [ abooallbt ] we show the brightness temperatures corresponding to all the observed flux data between 43  @xmath2 m and 6  cm used in our previous figures . for comparison
we have also plotted the predictions of the radiative equilibrium models . in deriving brightness temperatures from the observed fluxes
, we used the angular diameters of table  [ extents ] in the following relation : @xmath43 where @xmath44 is the observed flux in jy , @xmath0 is the angular diameter in mas , @xmath7 is the wavelength in microns , and @xmath45 is in kelvins .
the error bars in figures  [ atauallbt ] and  [ abooallbt ] correspond _ only _ to the uncertainties in the observed fluxes . in the case of the @xmath0  tau data at 2  cm and 3.6  cm ,
no error values were reported by @xcite and we used our own rough estimates of the plausible errors . before proceeding , a caveat needs to be noted regarding the longer wavelength points especially .
the observed fluxes presumably come from progressively shallower stellar layers as @xmath7 increases .
all of the brightness temperatures in figures  [ atauallbt ] and  [ abooallbt ] have been deduced assuming the photospheric angular diameters of table  [ extents ] .
atmospheric extension may play a role in the chromospheric layers of @xmath0  tau and @xmath0  boo .
@xcite estimated the extension of the chromospheres of these stars .
they argued that the uv 0.01 lines which they studied were formed in chromospheric layers with temperatures of @xmath3210000  k in @xmath0  tau and @xmath328000  k in @xmath0  boo .
they deduced that these chromospheric layers had radii of @xmath322.0  r@xmath46 and @xmath321.4  r@xmath46 , respectively , where r@xmath46 is the stellar photospheric radius .
although these estimates are , in the authors words , `` relatively crude '' , they do suggest that atmospheric extension may play a role at least in the hotter layers of the chromosphere .
different evidence of extension was found by @xcite in their study of 2  cm and 6  cm fluxes .
they treated all the radio emission as coming from an optically thick ionized wind for which a `` half - power radius '' could be deduced . in the case of @xmath0
boo they found the 2  cm half - power radius was @xmath32r@xmath46 ; at 6  cm , the half - power radius had grown to @xmath321.7  r@xmath46 . taking the long wavelength approximation to the brightness temperature expression above , we see that , for a given observed flux , @xmath47 , where @xmath0 is the angular diameter .
since the shallowest layers may have angular diameters significantly larger than those of the photosphere , the corresponding brightness temperatures in figures  [ atauallbt ] and  [ abooallbt ] should be regarded as upper limits .
we now describe why we believe that the flux excesses at the longer wavelengths arise from the presence of stellar chromospheres in @xmath0  boo and @xmath0  tau .
circumstantial evidence suggests the possibility that stellar chromospheres could explain the flux excesses .
this evidence comes by examining the depths at which @xmath36  =  1.0 . in the case of @xmath0
boo , as shown in figure  [ atau_aboo_models ] , the al chromospheric model departs from the purely radiative equilibrium model at a column mass of @xmath32  0.75  g  @xmath17 @xcite .
this is just in the column mass range of the pdk radiatve equilibrium model where optical depth unity is reached between 200  @xmath2 m and 1  mm . given this , it is not surprising that the pdk model might significantly _ underpredict _ the flux at both 1.4  mm and 2.8  mm in @xmath0  boo .
we believe that the progressively larger differences between the predictions of the pdk model and the observations beyond 120  @xmath2 m are due to the contributions of @xmath0  boo s chromosphere .
the circumstantial evidence in the case of @xmath0  tau is not so clear .
as shown in figure  [ atau_aboo_models ] , the kel @xmath0  tau chromospheric model begins its temperature rise at a column mass of @xmath480.3  g  @xmath17 @xcite .
this is slightly _ shallower _ in the plez @xmath0  tau model than where @xmath49  =  1.0 occurs . nevertheless , the 2.8  mm observed brightness temperature is distinctly higher than that predicted by the plez model in figure  [ atauallbt ] .
in fact , with the exception of our 1.4  mm point , the general run of observed points beyond 200  @xmath2 m suggests that some enhanced ( chromospheric ? ) emission may be present at 1  mm or even shorter wavelengths . whether this is true or not will take additional observations in this spectral region with better snr .      while nlte chromospheric modeling is well beyond the intent of this investigation , one of us ( adm )
has computed the infrared fluxes predicted by his chromospheric model of @xmath0  tau .
details of the model and the nlte computation may be found in @xcite .
this model is essentially _ identical _ to the kel model from the photosphere through the lower and middle chromosphere , differing from the kel model only at column masses shallower than  10@xmath50  g  @xmath17 ( see fig .  5 , @xcite ) .
the nlte radiative transfer computation was carried out with the independent code multi .
we show in figure  [ atau_nlte ] the results of this nlte computation with a chromospheric model compared with the observations , and with the predictions of the plez model .
it is immediately apparent that , again with the exception of our 1.4  mm point , the observed data beyond 200  @xmath2 m is better fit by the nlte chromosphere model than the lte plez radiative equilibrium model .
we believe this is good evidence that the stellar chromosphere of @xmath0  tau has a significant impact on the observed fluxes beyond 200  @xmath2 m .
indeed , it should be noted that the fit is even improved between 100  @xmath2 m and 200  @xmath2 m .      for @xmath0
boo ( figure  [ abooallbt ] ) it appears that the transition from purely radiative occurs beyond 120  @xmath2 m , and certainly well before 1.4  mm . in our opinion
this makes any purely radiative model for @xmath0  boo a poor choice for calibration in the sub - mm spectral region .
this appears to be the case for @xmath0  tau as well , especially considering the predictions of the mcmurry chromospheric model ( figure  [ atau_nlte ] ) .
given the evidence described above , we do not believe that it is prudent to use cool star radiative models for calibration purposes at the longer wavelengths _ unless alternatively calibrated observations already confirm _ the validity of the purely radiative models .
one must be certain that the flux predicted by the purely radiative model at the long wavelength end of the region being calibrated by the theoretical stellar spectrum does indeed correspond to the actual star .
in other words , we believe that when one is approaching the spectral regions where a chromospheric contribution can be significant , one can safely use purely radiative cool star models only as an interpolation tool for calibrations , not as an extrapolation tool .
we would like to emphasize the importance of the spectral region from @xmath3230  @xmath2 m to several millimeters to furthering an understanding of cool stellar atmospheres .
fluxes from this spectral region reflect the mean time- and spatial - averaged emission from the upper photosphere and chromosphere for stars like those we consider here ( see figure  [ atau_aboo_models ] ) .
accurate fluxes measured through this spectral region should allow one to test directly model atmosphere structures , even complex ones that include nlte and dynamics , over the very interesting depth range in which non - radiative processes like shocks begin to alter the structure .
by selectively sampling progressively longer wavelengths , one can follow the important transition from photosphere to chromosphere , locating the depth at which purely radiative models like the pdk and plez models begin to fail .
our original , theoretically extrapolated , continuum spectra for the two k giants , supplied to isophot as absolutely calibrated energy distributions , are in extremely good agreement with our new calibrated synthetic spectra ( figures  [ ataulws ] ,  [ aboolws ] ) , well within the @xmath226% uncertainties assigned to the original calibrators @xcite .
the lws spectra , reduced with respect to uranus as calibrator , are consistent with the new ( and old ) synthetic spectra , thereby uniting fir planetary and stellar calibrations at the few percent level ( figures  [ ataulws ] ,  [ aboolws ] ) .
there are independent data on @xmath0  boo , derived from the validation of isophot s final pipeline ( olp10 ) @xcite , that show it conforms to our 1996 model spectrum .
in addition , @xcite has described the performance through the iso mission of a significant number of isophot calibrators . by removing the contributions made by any particular source to the calibration process , and calibrating that source with respect to all other calibrators
, one can assess the ratio of observed to modeled irradiance .
using the database described by @xcite , we have been able to extend the range over which our new model of @xmath0  boo can be tested , and to validate @xmath0 tau independently of our lws data . for @xmath0
boo , we have determined the observed to modeled irradiance ratios for the isophot c100 , c120 , c135 , c160 , c180 , and c200 filters (  c " denotes the fir camera within the isophot instrument , and each filter is designated by its reference wavelength in  @xmath2 m ) . expressing these observed irradiance ratios with respect to our new modeled spectrum for @xmath0  boo we find : c100 , 1.02@xmath220.02 ; c120 , 1.14@xmath220.08 ; c135 , 1.01@xmath220.02 ; c160 , 1.06@xmath220.02 ; c180 , 1.00@xmath220.03 ; c200 , 1.00@xmath220.26 . for @xmath0  tau ,
measurements are available solely in the c160 filter and yield a ratio of 0.96@xmath220.05 , with respect to our new model .
we conclude that the separate absolute calibrations of lws and isophot are in accord within the uncertainty of the observations , despite their different dependences on planets , asteroids , and stars .
this should not be interpreted to mean that more accurate observations might not reveal discrepancies between fir observations and the predictions of static , lte model atmospheres computed in simple radiative / convective equilibrium .
our primary objective was to investigate whether purely radiative models for two normal k giants could reproduce the observed fir and microwave fluxes . we have found that for both @xmath0  tau and @xmath0  boo the radiative equilibrium temperature structures of the upper photospheres produce fluxes in agreement with the observations , but only up to a certain point . for @xmath0
boo and @xmath0  tau the agreement extends to at least 125  @xmath2 m , corresponding to a column mass of @xmath325  gm  @xmath17 in the pdk atmosphere and @xmath3210  gm  @xmath17 in the plez atmosphere .
both k giants clearly show definite chromospheric emission at 2.8  mm . in @xmath0
boo this is evident by 1.4  mm .
only additional observations can pinpoint precisely where between 125  @xmath2 m and 2.8  mm the chromospheric contribution becomes significant .
for this reason , we can not rely on purely radiative model atmospheres to form a basis for stellar calibration in the sub - mm region .
purely radiative models of the continuum might be used to provide calibrations _ between _ those shorter fir wavelengths where the validity of the model fluxes already has been independently verified by observations .
we wish to thank bertrand plez for computing an @xmath0  tau model for our particular set of input parameters and ruth peterson for the extended version of the @xcite @xmath0  boo model .
their help is greatly appreciated .
it is a pleasure to thank ullie klaas and phil richards of the isophot team for sending us data on their validation of olp10 as it bears on @xmath0  boo s spectrum , and sumita jayaraman for assistance in extracting dirbe data for the elongation angles at the times of lws and isophot observations of these k  giants .
mc acknowledges support for this work through subcontract uf99025 between the university of florida and uc berkeley , supporting the calibration of isophot .
dfc gratefully acknowledges access to the resources of the nasa advanced supercomputing ( nas ) division , nasa / ames research center , moffett field , ca , as well as the support of the nas division chief walter brooks and nas branch chiefs dochan kwak and william w. thigpen .
sgi / cray c90 supercomputer cycles for this research were provided through the consolidated supercomputing management office ( cosmo ) at nasa / ames .
this material is based upon work by dfc and dg supported by the national aeronautics and space administration under nra-99 - 01-adp-073 issued through the astrophysics data program .
we acknowledge support for bima and carma through nsf grants ast-9981308 and ast-0228963 .
altenhoff , w. j. , huchtmeier , w. k. , schmidt , j. , schraml , j. b. , & stumpff , p. 1986
, , 164 , 227 altenhoff , w. j. , thum , c. , & wendker , h. j. 1994 , , 281 , 161 ayres , t.r . , brown , a. , harper , g.m . ,
bennett , p.d . ,
linsky , j.l . ,
carpenter , k.g . , & robinson , r.d .
1997 , , 491 , 876 ayres , t.r . &
linsky , j.l .
1975 , , 200 , 660 blackwell , d.e . , lynas - gray , a.e . , &
petford , a.d .
1991 , , 245 , 567 van der bliek , n.s . ,
gustafsson , b. & eriksson , k. 1996 , , 309 , 849 carbon , d.f . ,
langer , g.e . , butler , d. , kraft , r.p . , trefzger , ch.f . ,
suntzeff , n.b . , kemper , e. , & romanishin , w. 1982 , , 49 , 207 carbon , d.f . , idehara , h. , & goorvitch , d. 2005 , aj , in preparation carpenter , k.g . ,
brown , a. , & stencel , r.e .
1985 , , 289 , 676 clegg , p. e. et al .
1996 , , 315 , l38 cohen , m. 1998 , , 115 , 2092 cohen , m. & davies , j. k. 1995 , mnras , 276 , 715 cohen , m. , megeath , t.g . ,
hammersley , p.l .
, martin - luis , f. , & stauffer , j. 2003 , , 125 , 2645 cohen , m. , wheaton , w.a . , & megeath , t.g .
2003 , , 126 , 1090 cohen , m. , walker , r.g . ,
barlow , m.j . , & deacon , j.r .
, 1992a , , 104 , 1650 cohen , m. , walker , r.g . , carter , b. , hammersley , p. , kidger , m. , & noguchi , k. , 1999 , , 117 , 1864 cohen , m. , walker , r.g . , jayaraman , s. , barker , e. , & price , s.d
. 2001 , , 121 , 1180 cohen , m. , walker , r. g. , & witteborn , f. c. 1992b , , 104 , 2030 cohen , m. , witteborn , bregman , j. d. , wooden , d.h . , salama , a. , and metcalfe , l. 1996 , , 112 , 241 cohen , m. , witteborn , f. c. , carbon , d. f. , augason , g. c. , wooden , d.h . , bregman , j. d. , & goorvitch , d. 1992c , , 104 , 2045 cohen , m. , witteborn , f. c. , carbon , d. f. , davies , j. k. , bregman , j. d. , & wooden , d. h. 1996b , , 112 , 2274 cohen , m. , witteborn , f.c . , walker , r.g . ,
bregman , j.d . , &
wooden , d.h .
1995 , , 110 , 275 decin , l. , waelkens , c. , eriksson , k. , gustafsson , b. , plez , b. , sauval , a.j . , van assche , w. , & vandenbussche , b. 2000 , , 364 , 137 decin , l. , vandenbussche , b. , waelkens , c. , decin , g. , eriksson , k. , gustafsson , b. , plez , b. , & sauval , a.j .
2003 , , 400 , 709 drake , s. a. & linsky , j. l. 1986 , , 91 , 602 dreher , j.w .
& welch , w.j .
1983 , , 88 , 1014 esa 1997 , the hipparcos and tycho catalogues , esa sp-1200 , or http://astro.estec.esa.nl/hipparcos/catalog.html gibson , j. , welch , w. j. , & de pater , i. 2005 , icarus , in press griffin , m. & orton , g. 1993 , icarus , 105,537 gry , c. et al .
2000 , in  the iso handbook , vol .
iii : lws - the long - wavelength spectrometer " ( esa sp-1262 ) gustafsson , b. , bell , r.  a. eriksson , k. , & nordlund , a. 1975 , , 42 , 407 harris , m.j .
& lambert , d.l .
1984 , , 285 , 674 johnson , h.  r. bernat , a.  p. krupp , b.  m. 1980 , , 42 , 501 laureijs , r. j. , klaas , u. , richards , p. j. , schulz , b. , & abraham , p. 2003 ,  the iso handbook vol .
iv : pht - the imaging photo - polarimeter " ( esa sp-1262 ) kelch , w.l . , linsky , j.l . ,
basri , g.s . ,
chiu , h .- y . ,
chang , s .- h .
, maran , s.p . , & furenlid , i. 1978 , , 220 , 962 kessler , m.f .
1996 , , 315 , l27 kurucz , r. 1993 , atlas9 stellar atmosphere programs 2 km / s grid .
kurucz cd - rom
cambridge , mass . : smithsonian astrophysical observatory , 1993 lemke , d. et al .
1996 , , 315 , l64 lim , t. , swinyard , b.m . ,
liu , x .- w . ,
burgdorf , m. , gry , c. , pezzuto , s. , & tommasi , e. 1997 ,  first iso workshop on analytical spectroscopy " , eds . a.m. heras , k. leech , n.r .
trams & m. perry , esa sp-419 , pg .
281 linsky , j.l . , & haisch , b.m .
1979 , , 229 , l27 mcmurry , a.d .
1999 , mnras , 302,37 mill , j. et al .
1994 , aiaa , 31 , 900 mitchell , k. j. et al .
1996 , in  unveiling the cosmic background " , ed .
e. dwek ( aip , cp 348 ) , p.301 muller , t.g . &
lagerros , j.s.v .
1998 , , 338 , 340 murakami , h. et al .
1996 , pas japan , 48 , l41 orton , g. , burgdorf , m. , davis , g. , sidher , s. , griffin , m. , swinyard , b. , & feuchtgruber , h. 2000 , aas , dps meeting # 32 , # 11.06 peterson , r.c .
1995 , private communication peterson , r.c . , dalle ore , c.m . & kurucz , r.l .
1993 , , 404 , 333 .
plez , b. 1999 , private communication .
plez , b. , brett , j.m . &
nordlund , a. 1992 , , 256 , 551 .
price , s.d . , paxson , c. , engelke , c. & murdock , t.l .
2004 , , 128 , 889 quirrenbach , a. , mozurkewich , d. , buscher , d.f . ,
hummel , c.a . , & armstrong , j.t .
1996 , , 312 , 160 richards , p.j .
2003 , in  the calibration legacy of the iso mission " , eds .
l. metcalfe , a. , salama , s.b .
peschke , m.f .
kessler ( esa sp-481 ) , p.279 ridgway , s.t . , jacoby , g.h . , joyce , r.r . , siegel , m.j .
, & wells , d.c .
1982 , , 87 , 1044 schaeidt , s. et al .
1996 , , 315 , l55 schulz , b. 2003 , in  the calibration legacy of the iso mission " , eds .
l. metcalfe , a. , salama , s.b .
peschke , m.f .
kessler ( esa sp-481 ) , p.83 schulz , b. et al . 2002 , , 381 , 1110 sidher , s.d . , swinyard , b m. , griffin , m.j . , orton , g.s . ,
burgdorf , m. , & lim , t.l .
2003 ,  the calibration legacy of the iso mission " , eds .
l. metcalfe , a. , salama , s.b .
peschke , m.f .
kessler ( esa sp-481 ) , p.153 smith , v.v . & lambert , d.l . 1985 , , 294 , 326 smith , v.v . & lambert , d.l .
1990 , , 72 , 387 sneden , c. 1973 , ph.d .
thesis , univ . of texas , austin swinyard , b. et al .
1996 , , 315 , l43 vernazza , j. e. , avrett , e. h. , & loeser , r. 1976 , , 30 , 1 welch , w. j. , & gibson , j. 2005 , in preparation wendker , h. 1995 , vizier catalog ii/199a wiedemann , g. , ayres , t. y. , jennings , d. e. , & saar , s. h. 1994 , , 423 , 806 wright , e. l. 1976 , , 210 , 250 wrixon , g. t. , welch , w. j. , & thornton , d. t. 1971 , , 169 , 171      cccccccccc @xmath0 tau & 14nov97 & c & 106.6845 & 4.32 hr & 0449 + 113 & 3c111 & mars , & & 17.14@xmath222.19 + & & & & & & & uranus , neptune & & + @xmath0 boo & 15nov97 & c & 106.6845 & 4.22 hr & 1415 + 133 & 3c111 & mars , & & 20.27@xmath222.95 + & & & & & 1357 + 193 & & uranus , neptune & & + @xmath0 boo & 15nov97 & c & 110.1282 & 4.22hr & 1415 + 133 & 3c111 & mars , & & 18.47@xmath222.73 + & & & & & 1357 + 193 & & uranus , neptune & & + @xmath0 boo & 15nov97 & c & both & & & & & & 19.73@xmath222.16 + @xmath0 tau & 06jun98 & c & 106.6845 & 2.99hr & 0449 + 113 & 3c111 & mars & & 13.10@xmath223.40 + @xmath0 boo & 11mar99 & b & 106.6845 & & 1357 + 193 & 1415 + 133 & & mars & 24.57@xmath221.20 + @xmath0 boo & 11mar99 & b & 110.1282 & 4.37 hr & 1357 + 193 & 1415 + 133 & & mars & 19.68@xmath221.23 + @xmath0 boo & 11mar99 & b & both & & & & & & 22.03@xmath220.97 + @xmath0 boo & 12mar99 & b & 106.6845 & 3.13 hr & 1357 + 193 & 1415 + 133 & mwc  349 , mars & mars & 21.02@xmath221.74 + @xmath0 boo & 12mar99 & b & 110.1282 & & 1357 + 193 & 1415 + 133 & mwc  349 , mars & mars & 18.90@xmath221.20 + @xmath0 boo & 12mar99 & b & both & & & & & & 19.92@xmath221.02 + @xmath0 tau & 05sep99 & d & 216.0982 & 1.85 hr & 0530 + 135 & 0530 + 135 & & mwc  349/mars & 24.23@xmath226.87 + @xmath0 tau & 05sep99 & d & 219.5418 & & 0530 + 135 & 0530 + 135 & & mwc  349/mars & 27.34@xmath228.12 + @xmath0 tau & 05sep99 & d & both & & & & & & 25.78@xmath225.64 + @xmath0 boo & 01jun00 & d & 216.0982 & 2.91 hr & 3c273 & 3c273 & mwc  349 & mars & 81.55@xmath222.35 + @xmath0 boo & 01jun00 & d & 219.5418 & 2.91 hr & 3c273 & 3c273 & mwc  349 & mars & 85.45@xmath222.44 + @xmath0 boo & 01jun00 & d & both & & & & & & 83.50@xmath221.71 + @xmath0 tau & 09dec00 & c & 106.6845 & 3.19
hr & 3c111 & 3c111 & w3oh , mars , & & 12.55@xmath222.16 + & & & & & & & uranus & & + @xmath0 tau & 09dec00 & c & 110.1282 & 3.19 hr & 3c111 & 3c111 & w3oh , mars , & & 12.77@xmath222.95 + & & & & & & & uranus & & + @xmath0 tau & 09dec00 & c & both & & & & & & 12.40@xmath221.91 + @xmath0 tau & all & all & 108.40 & & & & & & 13.97@xmath221.46 + @xmath0 tau & all & all & 217.82 & & & & & & 25.78@xmath225.64 + @xmath0 boo & all & all & 108.40 & & & & & & 20.09@xmath220.69 + @xmath0 boo & all & all & 217.82 & & & & & & 83.50@xmath221.71 +    lcccc c1 - 100 & 100 & 86401601 & 2 & 42.32@xmath223.57 + c1 - 100 & 100 & 86401602 & 2 & 44.53@xmath222.68 + c2 - 160 & 170 & 86002102 & 1 & 43.06@xmath224.39 + c1 - 60 & 60 & 84300501 & 5 & 34.54@xmath224.68 + c1 - 100 & 100 & 84300501 & 2 & 36.69@xmath224.14 + c1 - 60 & 60 & 84300501 & 5 & 40.89@xmath225.53 + c1 - 60 & 60 & combined & & 37.19@xmath223.57 + c1 - 100 & 100 & combined & & 42.24@xmath221.90 +    lcccc c1 - 50 & 50 & 27501511 & 1 & 18.16@xmath221.42 + c1 - 90 & 90 & 27502117 & 1 & 14.62@xmath220.87 + c1 - 105 & 105 & 27500602 & 1 & 12.54@xmath220.54 + c1 - 105 & 105 & 27500602 & 1 & 11.56@xmath220.94 + c2 - 120 & 120 & 27503008 & 1 & 9.10@xmath221.16 + c2 - 135 & 135 & 27503311 & 1 & 8.73@xmath221.12 + c2 - 160 & 170 & 27503614 & 1 & 9.23@xmath220.96 + c2 - 160 & 170 & 27503614 & 1 & 8.34@xmath220.27 + c2 - 200 & 200 & 27502402 & 1 & 6.61@xmath220.39 + c1 - 105 & 105 & combined & & 12.30@xmath220.47 + c2 - 160 & 170 & combined & & 8.41@xmath220.26 +    lcllll @xmath0  tau & c1 - 60 & sw2 & 4.2e-18 & 4.5e-19 & 3.8e-19 + @xmath0  tau & c1 - 100 & lw1 & 8.0e-19 & 1.3e-19 & 1.1e-19 + @xmath0  tau & c2 - 160 & lw4 & 9.1e-20 & 4.5e-20 & 5.9e-20 + @xmath0  boo & c1 - 50 & sw2 & 4.2e-18 & 2.0e-19 & @xmath511.3e-19 + @xmath0  boo & c1 - 90 & sw5 & 3.2e-18 & 8.7e-20 & + @xmath0  boo & c1 - 105 & lw1 & 8.0e-19 & 5.2e-20 & 2.7e-20 + @xmath0  boo & c2 - 120 & lw2 & 7.7e-21 & 1.2e-20 & 1.2e-20 + @xmath0  boo & c2 - 135 & lw3 & 2.0e-20 & 1.7e-20 & 7.0e-21 + @xmath0  boo & c2 - 160 & lw4 & 9.1e-20 & 1.5e-20 & 1.8e-20 + @xmath0  boo & c2 - 200 & lw5 & 4.4e-19 & 8.6e-21 & +    llcc 217.820 ghz ( 1.376 mm ) & observed flux ( mjy ) & 83.50@xmath221.71 & 25.78@xmath225.64 + & predicted flux ( mjy ) & 33.43 & 32.38 + & @xmath52predicted  observed flux@xmath52/@xmath14(observed ) & 29.3 & 1.17 + & & & + & observed brightness temperature ( k ) & 7040@xmath22140 & 2200@xmath22480 + & predicted brightness temperature ( k ) & 2820 & 2760 + 108.406 ghz ( 2.765 mm ) & observed flux ( mjy ) & 20.09@xmath220.69 & 13.97@xmath221.46 + & predicted flux ( mjy ) & 7.94 & 7.80 + & @xmath52predicted  observed flux@xmath52/@xmath14(observed ) & 17.6 & 4.23 + & & & + & observed brightness temperature ( k ) & 6840@xmath22230 & 4810@xmath22500 + & predicted brightness temperature ( k ) & 2700 & 2690 +    llcc adopted angular diameter ( mas ) & 21.0@xmath220.2 & 20.88@xmath220.10 + adopted parallax ( mas ) & 88.85@xmath220.74 & 50.09@xmath220.95 + deduced stellar radius ( cm ) & 1.77 x 10@xmath25 & 3.12 x 10@xmath25 + 1.64 @xmath2 m - 2.8 mm atmospheric thickness ( cm ) & 4.38 x 10@xmath53 & 4.22 x 10@xmath53 + percentage of stellar radius & 2.5% & 1.4% +