mass - loss from massive stars in the red supergiant phase is fundamental to the chemical enrichment of the interstellar medium @xcite .
red supergiants , with initial masses of 10@xmath250@xmath3 , are known to lose mass at rates of up to 10@xmath4/yr @xcite .
mass - loss rates are commonly estimated indirectly , using observations of infrared continuum excesses ( e.g. , * ? ? ? * ) or measurements of co rotational line intensities ( e.g. , * ? ? ?
outflowing material is also directly surveyed through observations of systematic spectral line displacements @xcite and interferometric studies over multiple epochs @xcite .
the physical processes involved in red supergiant mass - loss are however not well understood . while mass - loss from less massive stars , such as agb stars , is understood in terms of significant radial pulsations and winds driven by radiation pressure on dust grains , such mechanisms are not applicable to red supergiants @xcite , except those in a superwind phase @xcite .
rotation in stars that are spun - up as they condense from protostellar clouds could play a role in driving mass - loss @xcite , as could alfvn and acoustic waves ( e.g. , * ? ? ?
in the case of betelgeuse ( @xmath5 orionis , spectral type m2iab ) , the elevation of cool photospheric gas @xcite and gas at chromospheric temperatures @xcite has been attributed to large - scale surface convective activity @xcite .
betelgeuse is perhaps the archetypal red supergiant .
its large angular diameter ( 0.047 , * ? ? ? * ) and high apparent magnitude ( @xmath6 ) make it an excellent subject for high spatial resolution studies .
the continuum angular diameter of betelgeuse has been measured using interferometric techniques at radio @xcite , mid - ir ( 8@xmath212@xmath0 m , * ? ? ?
* ; * ? ? ?
near - ir ( e.g. , * ? ? ?
* ; * ? ? ?
* ) , optical ( e.g. , * ? ? ?
* ; * ? ? ?
* ; * ? ? ?
* ) and ultraviolet @xcite wavelengths , and significant variation of the diameter with wavelength is observed .
the apparent mid - ir diameter is also observed to vary significantly on timescales of years @xcite .
numerous interferometric campaigns , from optical to mid - ir wavelengths , have further revealed up to three hotspots on the stellar surface ( e.g. , * ? ? ?
* ; * ? ? ?
* ; * ? ? ?
* ; * ? ? ?
these spots range in intensity from a few to a few tens of percent of the intensity of the stellar disk .
the detection of these spots is generally interpreted as direct evidence for the presence of a few giant convection cells on the stellar surface ( e.g. , * ? ? ?
a profusion of circumstellar components have been identified surrounding betelgeuse .
furthest from the star , an analysis of mid - ir spectra and interferometric data by @xcite showed the presence of two thin shells of dust at approximately 40 and 80 stellar radii ( @xmath7 ) , with temperatures of @xmath1400k and @xmath1200k respectively .
later work by @xcite revealed a third , more compact dust shell at 4@xmath7 that was not apparent in previous data .
the dust condensation zone and acceleration mechanism are not identified .
radio continuum observations by @xcite at wavelengths between 0.7 cm and 6 cm , probing spatial scales of @xmath8 , revealed the presence of neutral gas at temperatures between 3500k and 1300k .
these temperatures are consistent with the temperature profile of the betelgeuse atmosphere in these regions , as derived from mean temperature and density models @xcite . at @xmath9
, the gas distribution is asymmetric , indicating an anisotropic mechanism for elevating the gas above the photosphere .
ultraviolet @xcite and h@xmath5 @xcite imaging observations have shown that this cool gas co - exists with asymmetrically - distributed hot chromospheric material .
observations of [ fe ii ] lines by @xcite that correspond to the dominant cool atmospheric component showed , surprisingly , that this material was at rest with respect to the star .
the circumstellar envelope as viewed across the _ jhk _ wavelength bands was imaged using a ` lucky - imaging ' approach combined with adaptive optics by @xcite , and was found to include a plume extending to 6@xmath7 .
closer to the star , strong evidence exists for the presence of a cool molecular layer or shell directly above the photosphere within @xmath10 ( e.g. * ? ? ?
* ; * ? ? ?
this layer is optically thick at mid - ir wavelengths but is optically thin in the near - ir .
spectral signatures of water from betelgeuse and the red supergiant @xmath0 cephei , found using data from the infrared space observatory ( iso ) , were interpreted by @xcite as originating from a 1500k layer surrounding these stars .
modeling of near- and mid - ir interferometric data by @xcite confirmed the presence of such a layer , with a mid - ir optical depth of @xmath12 and a temperature of @xmath12000k .
further work by @xcite combining both interferometric and spectroscopic data placed the layer surrounding betelgeuse at beyond @xmath11 , with a temperature of 2250k .
the layer was shown not to be solely composed of water by @xcite , who suggested the presence of amorphous alumina ( al@xmath12o@xmath13 ) to explain the high mid - ir opacity .
the spectrum of water lines from betelgeuse is indeed more complex than revealed by the iso data ; signatures of water at photospheric temperatures were detected by @xcite and @xcite .
most recently , wide - band , spectrally dispersed mid - ir interferometric observations by @xcite were fitted by a layer containing sio , amorphous alumina and water .
high resolution measurements of the co overtone band near 2.3@xmath0 m @xcite and one - dimensional imaging @xcite were modeled with a large patch of gas near @xmath14 with motions of 10@xmath215kms@xmath15 with respect to the star .
this emerging model of a photosphere and overlying layer for a selection of red supergiant stars allows a better understanding of dust formation and red supergiant wind acceleration ( e.g. * ? ? ?
* ) , leading towards a convergent picture of red supergiant mass loss .
studies of the extended betelgeuse atmosphere form a basis for a better understanding of the atmospheric dynamics of red supergiant stars . the existence of a cool layer with high opacity at mid - ir wavelengths directly above the photosphere is now established for betelgeuse .
the size , temperature and composition of the layer have been investigated , assuming spherical symmetry .
little , however , is known about the spatial structure and time - variability of this layer .
consistent instrumentation and modeling procedures have also not been applied over multiple observation epochs . in this paper , we present observations of the mid - ir continuum of betelgeuse using the berkeley infrared spatial interferometer ( isi , * ? ? ?
interferometric data using three telescopes were gathered during each year between 2006 and 2010 , with angular resolutions sufficient to model the stellar shape and potential large - scale features .
we describe the observations and image - modeling processes in  2 , and perform fits to a photosphere - and - layer model for our observations in  3 . in  4 , we discuss the nature of the observed non - uniformities in our model images .
we examine the implications of our results for models of the betelgeuse atmosphere and red supergiant mass - loss in  5 , and present our conclusions in  6 .
the isi consists of three 1.65 m aperture telescopes with a heterodyne detection system at each telescope .
the double - sideband detection systems are sensitive to @xmath16@xmath0 m radiation within single sideband widths of @xmath12.7ghz .
the @xmath17co@xmath12 laser local oscillators can operate at numerous wavelengths corresponding to vibrational molecular transitions between wavelengths of 9@xmath0 m and 12@xmath0 m .
for the observations reported here , the local oscillators at each telescope were tuned to a wavelength of 11.15@xmath0 m , chosen so as to exclude strong stellar and telluric spectral lines from the detection band @xcite .
ccc 2006 & 8 , 9 , 10 nov ; 7 dec & 4.2@xmath180.2 + 2007 & 14 , 15 , 16 nov & 4.1@xmath180.3 + 2008 & 22 , 23 , 24 sep & 3.8@xmath180.5 + 2009 & 4 , 5 , 7 , 9 , 18 , 20 , 22 nov & 4.4@xmath180.1 + 2010 & 5 , 17 , 18 nov ; 9 , 13 dec & 4.9@xmath180.3           the van cittert - zernicke theorem equates the complex cross - correlation of the electric fields detected at two spatially - separated telescopes with the two - dimensional fourier transform of the source brightness distribution within the telescope fields of view . instead of measuring complex cross - correlations , the isi , like all optical / ir interferometers , measures the normalized cross - correlation amplitude , known as the michelson fringe visibility .
the michelson fringe visibility , @xmath19 , is given by @xmath20 where the fringe power , @xmath21 is the cross - correlation amplitude between signals measured at two telescopes , @xmath22 and @xmath23 , and @xmath24 and @xmath25 are the source powers measured at telescopes @xmath22 and @xmath23 respectively
. visibility amplitudes are recorded for each isi telescope baseline pair .
the visibility phase in the mid - ir for a given baseline is , however , strongly affected by randomly varying atmospheric pathlengths above each telescope , and can not be recovered from the data .
the isi instead measures the closure phase , or the sum of the measured visibility phases for each baseline .
the closure phase is independent of optical pathlength differences between pairs of telescopes @xcite and it provides information about the point - inversion symmetry of the source . the spatial frequency sampled by the projection of a baseline of length @xmath26 , perpendicular to the direction of a source , at a wavelength @xmath27 , is proportional to @xmath28 . by obtaining measurements at different source hour - angles ,
multiple spatial frequencies at different position angles are sampled .
visibility and closure phase data were recorded for betelgeuse over timespans of five weeks or less each year between 2006 and 2010 ( see table 1 ) . during this time period ,
the three isi telescopes were configured in an approximately equilateral triangle with baseline lengths between 34 m and 40 m . for betelgeuse , this configuration allows spatial frequencies between 20 and 37sfu ( spatial frequency unit ; 1sfu@xmath29cyclesradian@xmath15 ) to be sampled . for all observations ,
tip - tilt corrections were applied to center the source on _ k_-band guider cameras operating at frame rates of 27hz . a standard observing sequence consisted of 4.5-minute dataset intervals of staring at the source and position - switching between the source and blank fields on either side of the source .
the position - switching measurements moved the telescopes between the source and the sky every 15s .
the source powers were measured using additional 140hz chopping between the source , or the sky , and cold loads with temperatures near that of the sky . a single visibility measurement consisted of either : ( 1 ) a fringe power measured during a chopping / position - switching dataset , normalized using the source powers during that dataset ; or ( 2 ) a fringe power measured during a staring dataset , normalized by weighted averages of source powers measured during the nearest chopping / position - switching datasets .
closure phases were calculated for each triplet of visibility measurements .
visibility magnitudes and closure phase zero - point offsets were calibrated using observations of aldebaran ( @xmath5 taurus ) conducted on the same night .
aldebaran was assumed to have a diameter of 20mas and zero closure phase .
the visibility calibrations for the two kinds of visibility measurements were different because of the different times spent on - source during chopping and staring datasets .
data that were affected by systematic errors or poor atmospheric conditions were discarded .
the total flux densities of betelgeuse during each observing epoch were also estimated using aldebaran , assumed to have a flux density at 11.15@xmath0 m of 615jy @xcite , as a calibrator .
these total flux density values are given in table 1 , and represent the total power of betelgeuse within the 5@xmath305  fields of view of the isi detectors .
image models were fitted to the calibrated visibility and closure phase datasets for each epoch .
the data , analytic fits , and coverage of the @xmath31 spatial frequency plane , are shown in figure 1 .
effective model - independent image reconstruction was not possible because of the sparsity of the @xmath31 coverage in all epochs .
our model images include a centered uniform disk ( ud ) , and up to two point sources offset from the image centers .
the uniform disks represent the apparent stellar surface , and the point sources represent localised intensity fluctuations in this surface .
emission from the dust surrounding the star at large angular scales , while contributing to the total detected intensity , was assumed not to contribute to the visibility at the sampled spatial frequencies .
the intensity of each image component was fitted to a fraction of the total intensity .
other free parameters included the ud radius and the positions of the point sources .
the model complex visibility at a point @xmath32 in the @xmath31 plane was given by @xmath33 where @xmath22 is the fraction of the total intensity contributed by the ud , @xmath34 is the uniform disk angular radius , @xmath35 denotes a bessel function of the first kind and of order unity , @xmath36 represents the fractions of the total intensity of the point sources , and @xmath37 and @xmath38 represent angular offsets to the west and north for the point sources , with @xmath39 .
ccccccc 2006 & 0.51@xmath180.02 & 24.5@xmath180.8 & 0.04@xmath180.01 & -2.4@xmath180.1 & -23.7@xmath180.2 & 2.683 + 2007 & 0.54@xmath180.01 & 24.8@xmath180.3 & - & - & - & 2.053 + 2008 & 0.48@xmath180.02 & 24.7@xmath180.8 & - & - & - & 2.416 + 2009 & 0.557@xmath180.008 & 26.4@xmath180.2 & 0.009@xmath180.001 & -25.0@xmath180.3 & 10.5@xmath180.2 & 4.565 + & & & 0.007@xmath180.001 & 18.3@xmath180.3 & -20.2@xmath180.2 & + 2010 & 0.61@xmath180.03 & 26.4@xmath180.3 & 0.025@xmath180.009 & -29.7@xmath180.3 & -9.3@xmath180.1 & 3.65 + & & & 0.009@xmath180.004 & 8.3@xmath180.3 & -25.9@xmath180.1 &    we used a weighted least - squares fitting technique for each epoch , with each baseline and the closure phase data given equal weights .
the differing quality of data from different epochs necessitated using models with different numbers of free parameters to fit different epochs .
data from 2006 were fitted with a single point source in addition to a ud , data from 2007 and 2008 were only fitted with uds , and data from 2009 and 2010 were each fitted with a ud and two point sources . an additional point source was added to the fit for a given epoch only if the reduced @xmath40 values calculated for the full dataset , and for the closure phase data alone , were both decreased upon the addition of the point source .
we also attempted to fit other image models .
these included uniform ellipses , uniform ellipses with point sources , and uds with offset uniform disks or offset two - dimensional gaussian profiles in place of the point sources .
none of these trials , however , converged to a satisfactory image , and the angular sizes of the fitted point sources are not determined .
the fractions of the total intensity contributed by the point sources were not constrained to be positive , as our data are sensitive to both bright and dark features on the stellar disk .
simulations were conducted to assess possible degeneracies in fitted image model parameters given the sparsity of the @xmath31 coverage of our datasets , in order to ensure that our data were not being over - fitted .
no significant degeneracies were found for any epoch .
for example , the 2006 , 2009 , and 2010 datasets constrain both the locations and brightnesses of the point sources .
closure phase measurements at multiple array position angles constrain the position angles of the point sources , and both the closure phase measurements and the visibility measurements constrain the angular offsets of the point sources from the image centers .
the results of the fits , as well as the reduced @xmath40 for each epoch , are given in table 2 .
all the fits have reduced @xmath40 values that are greater than unity .
this implies that the data contain minor features besides gaussian noise that are not being fitted by the free parameters of the image models .
isi data for betelgeuse obtained during 2006 were also analyzed by @xcite . despite our choice of different observing dates , a different method of calibration and slightly different image fitting procedures ,
our results match those of @xcite within the margins of error .
we however do not adopt tatebe et al.s results for the 2006 epoch in order to maintain the same analysis procedure across all epochs .          using isi observations over the period 1993@xmath22009 ,
@xcite reported a systematic decrease of 15% in the mid - ir continuum ud diameter of betelgeuse . here , we augment this finding with new measurements of the ud diameter of betelgeuse . figure 2 shows the data plotted in figure 1 of @xcite , with the new results also included . though the present work re - analyzes some previous observations examined by @xcite , the diameter measurements are obtained from different data sets and analyses from those reported by @xcite .
the measurement of the ud diameter of betelgeuse at 11.15@xmath0 m by @xcite , obtained using the very large telescope interferometer in 2003 , is also plotted .
recently , @xcite showed that k - band diameter measurements , made over a similar time period , vary less than the mid - ir diameters plotted in figure 2 .
while our results from 2006@xmath22008 are consistent to within the margins of error with the trend towards a smaller size @xcite , the ud diameters of betelgeuse in 2009 and 2010 are significantly larger than those during 2006@xmath22008 .
our analysis method for the 2007 and 2008 datasets essentially matched that of @xcite , in that only uds were fitted to the data .
our technique however differed somewhat in the methods of calibration and fitting .
we consider the variability in the apparent size of betelgeuse in the mid - ir further in  5.2 .
effective temperatures of stellar surfaces are among the fundamental parameters of stars @xcite .
the isi measurements of the flux density of betelgeuse , and of the ud sizes and fractions of the total intensity , allow the effective temperature of the 11.15@xmath0 m apparent stellar surface to be determined at each observing epoch .
we use a planck function for the power emitted by a blackbody to express the effective temperature , @xmath41 , as @xmath42 where @xmath43 is planck s constant , @xmath44 is the observing frequency , @xmath45 is boltzmann s constant
, @xmath46 is the solid angle subtended by the stellar disk , @xmath47 is the vacuum speed of light , and @xmath48 is the measured power from the stellar disk in units of wm@xmath49hz@xmath15 . we do not use the widely - applied technique of calculating effective temperatures using estimates for the bolometric flux from the star ( e.g. , * ? ? ?
* ; * ? ? ?
* ) for two reasons : radiation at mid - ir wavelengths is strongly affected by thermal emission from circumstellar dust shells and bolometric correction factors at these wavelengths are unique to individual stars ; and also because the correspondance between the 11.15@xmath0 m scattering surface and the photosphere of betelgeuse is not certain .
the calculated effective temperatures for the observed stellar disk of betelgeuse during each observing epoch are plotted in figure 3 .
we do not account for attenuation by the interstellar medium , estimated as modifying the observed flux density by a factor of 0.975 @xcite .
we also do not account for absorption of the stellar radiation by the dust that surrounds betelgeuse at angular scales of 1  or larger , or for the emission from this dust in the line of sight to the stellar face .
this is because the dust shells are extremely optically thin at mid - ir wavelengths when viewed radially , with a total line - of - sight optical depth of approximately 0.0065 @xcite .
our measured effective temperatures of the mid - ir surface of betelgeuse are consistently lower than the value of 3641k derived by @xcite for the photospheric component of a photosphere - and - layer model .
furthermore , significant variability is evident between the effective temperatures measured at different epochs .
we therefore consider sources of mid - ir continuum opacity in the atmosphere of betelgeuse in an attempt to elucidate the scattering surface probed by our observations and the cause of the observed variability .
stellar atmospheres have opacities that vary with wavelength .
stars appear smaller at wavelengths where the atmospheric opacity is smaller .
the surface of a star , when observed at a given wavelength , can be approximated as the location where the optical depth at that wavelength , measured from the observer , is approximately unity .
the photosphere of a star , however , is the apparent stellar surface observed at a wavelength where the opacity matches the rosseland mean opacity @xcite .
the rosseland mean opacity , @xmath50 , is given by @xmath51 where @xmath52 is the intensity of radiation from a point in a stellar atmosphere , @xmath53 is the temperature at that point , and @xmath54 is the opacity at a frequency @xmath44 .
regions of a stellar spectrum with large radiation intensity but low opacity contribute most to determining the rosseland mean opacity .
@xcite concluded that the mid - ir continuum opacity of a gas with solar atomic abundances under surface conditions present on red giant stars closely matches the rosseland mean opacity of such a gas .
molecular layers surrounding red giant stars have , however , been shown to significantly change the mid - ir opacity from that of a solar abundance atomic gas , while remaining relatively transparent in shorter wavelength ir bands @xcite . the @xcite model for betelgeuse , combining @xmath55-band visibility data and a 11.15@xmath0 m isi visibility curve @xcite , includes a shell with a @xmath55-band optical depth of 0.06 and a 11.15@xmath0 m optical depth of 2.33 .
@xcite showed that this shell was composed of water , amorphous alumina dust and sio gas , and that the 11.15@xmath0 m opacity was dominated by alumina .
the presence of alumina in a layer surrounding betelgeuse , within 1.5@xmath7 , was suggested by @xcite as the only dust species that could condense at the temperatures present in this region .
alumina dust grains forming close to stars possibly play a role as nucleation sites for dust species that condense at lower temperatures ( e.g. * ? ? ?
bremsstrahlung interactions between electrons and neutral hydrogen atoms also contribute to the mid - ir continuum opacity in red giant stars @xcite .
this emission mechanism was adduced by @xcite to explain the variation with wavelength of the apparent sizes of the radio photospheres of long - period variables , and ascribed to betelgeuse by @xcite .
we discuss this opacity model further in  3.4 .
prior work , the dramatic variations in apparent mid - ir diameters of betelgeuse , and the low effective temperatures all provide clear motivation to model the new isi observational results in terms of a photosphere and a molecular layer .
our method is similar to that of @xcite in that we modeled a spherically symmetric layer surrounding an opaque stellar photosphere with an assumed diameter of 0.04371  and an assumed surface temperature of 3641k @xcite .
visibility and closure phase models for any observed point sources were subtracted from the data during each observation epoch in order to allow the model to be spherically symmetric .
we did not attempt to fit an asymmetric layer because of the large number of free parameters involved .
the free parameters of our model were : the temperature , @xmath53 , of the layer ; the optical depth , @xmath56 , of the layer along the line of sight to the center of the star ; and , the outer radius , @xmath57 , of the layer .
the layer was assumed to have uniform temperature and opacity , and was also assumed to have an inner radius that corresponds to the photospheric radius .
radiation within the isi band was assumed to arise only from blackbody emission by the layer .
the modeling procedure involved fitting visibilities derived from a model image , created using a trial set of free parameters , to the measured visibilities .
the measured visibilities for each epoch , and their errors , were scaled by the measurements of the measured power emitted by the ud surface of betelgeuse .
the flux density of a pixel in the model image grid was calculated as : + for @xmath58 , @xmath59 and , for @xmath60 , @xmath61 here , @xmath62 describes the optical path length through the layer , @xmath63 is the absorption coefficient , @xmath64 is the radial distance of a pixel from the image center , @xmath46 is the solid angle subtended by a single pixel , @xmath65 is the specific intensity radiated by a blackbody of temperature @xmath53 at frequency @xmath44 , and @xmath66 is the assumed temperature of the photosphere . @xmath62 and @xmath63 are given by @xmath67 the pixel size of the model images was @xmath68 milli - arcseconds . for each iteration of the fit ,
a model image was generated and integrated along one dimension , and a one - dimensional fourier transform was performed in order to generate a set of visibilities .
each model visibility was compared to measured visibilities that were closest to it ( i.e. , within @xmath69 ) in spatial frequency .
a levenberg - marquardt fitting algorithm was employed to calculate the free parameter set that minimized the global @xmath40 .
the results of the fit for each observing epoch are given in table 3 .
errors for each parameter were calculated using the resulting covariance matrix . in order to confirm the validity of our analyses , we added the visibility and closure phase models of the point sources , that were earlier subtracted from the visibility data from each epoch , to the best - fit model visibilities .
@xmath40 values relating this sum to the original visibility data were calculated for each epoch , and were found to be comparable to the @xmath40 values of the fits to uds and point sources shown in table 2 .
the changes with time of the fitted layer temperatures and outer radii are similar to the changes with time in the effective stellar temperatures and the apparent ud radii respectively .
all our fitted layer radii and temperatures are roughly consistent with the 1.3@xmath21.5@xmath7 ( 0.057  to 0.065diameter ) and 1500@xmath22000k layers modeled by most authors @xcite .
the values we calculate for the 11.15@xmath0 m optical depth of @xmath11 are half of those resulting from the models of @xcite and @xcite .
cccc 2006 & 1.9@xmath180.7 & 0.8@xmath180.1 & 28.6@xmath180.7 + 2007 & 2.5@xmath180.9 & 0.9@xmath180.1 & 25.5@xmath180.9 + 2008 & 2@xmath181 & 1.2@xmath180.2 & 25@xmath181 + 2009 & 2.3@xmath180.8 & 0.99@xmath180.06 & 26.8@xmath180.8 + 2010 & 2.8@xmath180.6 & 1.05@xmath180.03 & 29.8@xmath180.6      the lack of wideband mid - ir spectral information during our observing epochs limits modeling of molecular gas and dust species that contribute to the layer opacity .
we are able , however , to investigate some basic physical characteristics of the layer using the opacity mechanism described by @xcite .
this work showed that , for solar abundance gas at temperatures similar to those in red giant atmospheres , collision processes between neutral hydrogen atoms and free electrons dominate the mid - ir opacity . for temperatures of 2500k to 3000k ,
the primary sources of free electrons were found to be thermally ionized na and al .
the abundances of h@xmath70 , h@xmath12 and h@xmath71 relative to h were found to be too low to contribute to the opacity through collisions with electrons , and photoionization and photodissociation processes were also found to be insignificant .
the absorption coefficient due to electron - hydrogen atom collisions , first derived by @xcite , is given in simplified form as @xmath72 where @xmath73 is the electron number density , @xmath74 is the hydrogen atom number density , and @xmath22 is expressed as a third - order polynomial in @xmath53 @xcite .
this opacity mechanism is used to explain the sizes of the radio photospheres of red giant stars @xcite , at different wavelengths .
for example , @xcite found that the radio photosphere of betelgeuse appeared larger , yet cooler , for longer wavelengths .
we consider the implications of this mechanism dominating the mid - ir continuum opacity .
equation 9 can be used to calculate hydrogen atom number densities from the fitted layer optical depths .
the electron number densities were obtained from table 2 of @xcite , which gives the fractional ionizations , relative to @xmath74 , of various metals at different temperatures .
we further used the hydrogen number densities to calculate the total layer mass , @xmath75 , for each epoch , using a distance of @xmath76pc to betelgeuse @xcite .
our results are given in table 4 .
we find masses for the layer of between @xmath77 and @xmath78 .
our values for the hydrogen number density of between @xmath79@xmath80 and @xmath81@xmath80 , using a hydrogen atomic mass of @xmath82 g , are equivalent to mass densities of @xmath83g@xmath80 .
such values are consistent with the hydrogen density derived for the betelgeuse layer modeled by @xcite , where they derived their hydrogen density by assuming a ratio between the h@xmath12o number density and the hydrogen number density of @xmath84 @xcite and using their calculated values of @xmath85 for the layer .
our values for the hydrogen density are therefore also consistent with the values for the water density derived by @xcite , though they result from very different methods of calculation .
ccc 2006 & 5@xmath181 & 6@xmath181 + 2007 & 4.6@xmath180.9 & 3@xmath181 + 2008 & 9@xmath182 & 4@xmath182 + 2009 & 4.3@xmath180.8 & 3@xmath181 + 2010 & 3.5@xmath180.4 & 5.0@xmath180.9
we now consider a physical interpretation of the point sources used to model the isi interferometric data on betelgeuse .
figure 4 shows images of the intensity distribution of betelgeuse during each observing epoch , calculated from the photosphere and layer models , with the point source locations indicated by filled white disks .
the color brightness of the images is consistent with the same flux density scale , also depicted in figure 4 , and the areas of the disks representing the point sources indicate their flux densities calculated from their fitted fractions of the total intensity .
all point sources are positive in intensity with respect to the star .
non - uniformities on the observed stellar disk of betelgeuse at optical and near - ir wavelengths have been reported by several authors ( e.g. , * ? ? ?
* ; * ? ? ?
these asymmetries are generally modeled either as point sources , or as uds or two - dimensional gaussian profiles of fixed , unresolved size .
the intensities of these ` spots ' are generally found to be between 5% and 20% of the intensity of the stellar disk .
the distribution and strengths of the spots are found to vary significantly across optical and near - ir bands during contemporaneous observations @xcite , with the star sometimes appearing featureless in one band and highly non - uniform in another .
the spots are also found to vary on timescales of less than 8 weeks @xcite .
one such spot was recently resolved in the _ h_-band by @xcite , who reported a diameter of 0.011 for this spot .
these results are consistent with a prediction by @xcite for extremely large - scale convection on red giants .
red supergiants are suggested to have only a few tens of convection cells on their surfaces , as opposed to the thousands present on the surface of the sun .
similar predictions also result from three - dimensional numerical simulations ( e.g. , * ? ? ?
* ; * ? ? ?
* ) , and an initial attempt at matching simulations to interferometric observations was made by @xcite .
@xcite predicted temperature fluctuations of up to @xmath86k with respect to the photospheric temperatures of red giant stars .
this is consistent with the intensities of spots observed in the visible and near - ir .
the point sources fitted to our mid - ir data do not arise directly from large - scale convective features on the surface of betelgeuse .
if we assume a 0.01  diameter for the weakest of our point sources ( see table 2 ) , which contributes 0.7% of the 4400jy flux density of betelgeuse in 2009 , its temperature of @xmath1660k above the effective surface temperature , calculated using equation 2 , is consistent with a photospheric hotspot of @xmath11000k behind the @xmath87 layer .
if we assign a similar diameter to any of the other point sources we observe , even greater temperatures would result for associated photospheric hotspots .
the mid - ir point sources are hence more intense than expected from surface temperature fluctuations caused by giant convection cells @xcite .
furthermore , every one of our point sources is located beyond the edge of the photospheric disk ( see figure 4 ) .
it is clear therefore that the observed asymmetry of betelgeuse in the mid - ir is caused by a non - uniform intensity distribution in the layer .
we do not have sufficient sampling of the @xmath31 plane in any epoch to identify whether the layer is in fact spherically symmetric , but with temperature fluctuations , or whether the layer is extended in particular directions .
the point sources can also be explained in terms of opacity inhomogeneities in the layer . in the rayleigh - jeans approximation to the planck law , the effective radiation temperature , @xmath88 , of a column of material with optical depth @xmath56 and temperature @xmath53
is given by @xmath89 for example , increasing the value of @xmath56 from 1 to 3 increases @xmath88 by a factor of 1.5 .
non - uniform spatial , temperature and opacity distributions of the betelgeuse layer could all explain the observed mid - ir non - uniformities .
molecular layers directly surrounding the photospheres of red supergiant stars , and in particular betelgeuse , have been modeled by a variety of authors under differing assumptions for the layer geometry and the dominant opacity source .
a layer of water with a column density of 10@xmath90cm@xmath49 at a temperature of @xmath11500k was first suggested by @xcite for betelgeuse , based on modeling of near - ir spectra . a combination of spectral and interferometric modeling by @xcite , assuming a layer geometry identical to ours , suggested a layer temperature of 2050k and an outer radius of 1.45@xmath7 , while modeling the mid - ir opacity as being dominated by a continuum of water emission and absorption lines . under the assumption of a layer with a well - defined inner radius
, @xcite modeled a shell beyond 1.3@xmath7 with a temperature of 2250k , again dominated by water opacity .
more advanced spectral modeling by @xcite revealed that water could not play a dominant role in determining the mid - ir continuum opacity , and instead suggested the presence of amorphous alumina in the molecular layer .
this model was consolidated by @xcite to suggest a shell between 1.32@xmath7 and 1.42@xmath7 at a temperature of 1520k .
the properties we derive for a layer surrounding betelgeuse , summarised in table 3 for each observing epoch , are largely consistent with these results .
we find temperatures that vary between 1900k and 2800k and outer radii that vary between 1.14@xmath7 and 1.36@xmath7 .
our modeling technique most closely matched that of @xcite in assuming thermal continuum emission from the layer within our observing band .
this assumption is justified by the featureless spectrum of betelgeuse within , and around , our observing band and by the lack of any significant water lines within and surrounding our band @xcite .
the results are consistent with previous models of the betelgeuse layer .
a consistent picture of the composition of the layer has , however , not been achieved . while spectral modeling has uncovered the presence of numerous molecular species , including h@xmath12o , co , sio and possibly cn , as well as alumina dust @xcite , the opacity of the layer as a function of wavelength is not yet fully determined .
our results for the optical depth of the layer of @xmath11 at 11.15@xmath0 m wavelength are not in agreement with those of @xcite and @xcite , who obtain optical depths of @xmath12 .
the number density of atomic hydrogen , based on the 11.15@xmath0 m opacity arising from electron - hydrogen atom collisions , is @xmath110@xmath92@xmath80 .
this density , given the derived layer size , corresponds to h@xmath12o column densities of 10@xmath93cm@xmath49 ; this result is consistent with h@xmath12o column densities measured by @xcite and those derived from spectral modeling ( e.g. , * ? ? ?
we therefore suggest that electron - hydrogen atom collisions play a significant role in determining the mid - ir continuum opacity of the betelgeuse layer , and that this mechanism must also be accounted for in models for the layer composition .
we have shown that the layer of material surrounding the photosphere of betelgeuse has a non - uniform intensity distribution .
a similar result was recently reported by @xcite based on interferometric measurements of the one - dimensional intensity and velocity profile of co gas in the inner atmosphere of betelgeuse .
@xcite also reported a change in the one - dimensional co intensity and velocity profile between two epochs spaced by a year .
we find that the intensity distribution of the layer varies between each the five observing epochs we present , which are also spaced by approximately one year .
the non - uniformities that we observe in the layer reflect , on much smaller spatial scales , the non - uniformities observed in various components of the circumstellar environment of betelgeuse . a gaseous component at similar temperatures to the betelgeuse layer
was shown to have an asymmetric intensity distribution by @xcite .
our observations and those of @xcite could be linked by the opacities in both wavelength bands being partly caused by electron - hydrogen atom collisions .
the opacity due to this mechanism is proportional to @xmath94 , implying that the observed stellar size decreases with shorter observing wavelengths .
material at chromospheric temperatures is also observed to be asymmetrically distributed @xcite , and observations of a circumstellar envelope by @xcite at near - ir wavelengths also reveal significant asymmetry .
asymmetry in the large - scale ( @xmath11 ) dust distribution surrounding betelgeuse was also observed by @xcite .
molecular layer models have been applied to the few red supergiant stars that have been observed with sufficient spatial resolution , such as @xmath0 cephei @xcite and @xmath5 herculis @xcite .
similar layers are identified around agb stars @xcite .
the properties of individual red supergiant layers vary significantly from star to star , and a greater sample of measurements is required to understand what determines individual layer properties .      the processes by which red supergiant stars form and accelerate dust , while crucial to mass - loss scenarios , are little understood .
existing models of stellar atmospheres ( e.g. , * ? ? ?
* ; * ? ? ?
* ) can not predict locations for dust condensation , and attempts to model dust acceleration through radiation pressure have also so far been unsuccessful ( e.g. , * ? ? ?
our multi - epoch study of a layer directly surrounding betelgeuse provides new insight into the processes involved in red supergiant mass - loss .
m flux density , of betelgeuse over the period 2006@xmath22010 .
the visual magnitude measurements , joined by a solid line , are 30-day averages of data obtained from the aavso online archive .
the 11.15@xmath0 m flux density measurements , also listed in table 1 , are shown here as crosses .
the timespans of each isi observing epoch are bounded by vertical dotted lines . ]
we have shown that the basic properties of the layer , given in table 3 , exhibit striking variability from year to year .
the layer is found to decrease in size and increase in optical depth over the years @xmath95 , and increase in size and decrease in optical depth between @xmath96 .
the layer temperature does not show such systematic variability , although the largest layer size ( in 2010 ) does correspond to the highest temperature .
the variability in the apparent mid - ir size of betelgeuse reported by @xcite can also be interpreted to indicate variations in the layers .
intriguingly , a similar phenomenon was reported by @xcite , who noted a decrease in the size of betelgeuse at visible wavelengths followed by an increase .
figure 5 shows measurements of the visual magnitude of betelgeuse during between 2006 and 2010 , as well as the isi measurements of the total 11.15@xmath0 m flux density ( also listed in table 1 ) .
the timespans of the isi observing epochs are also plotted .
the visual magnitude measurements were obtained from the online database of the american association of variable star observers ( aavso ) .
the data , collected by a multitude of volunteer observers , were averaged over 30-day timespans ; each timespan with less than 10 measurements was discarded .
the visual variability of betelgeuse over our observing epochs was approximately 0.9 magnitudes .
the visual intensity of betelgeuse is not likely to represent that of the photosphere because the apparent stellar diameter at visible wavelengths is greater than that at near - ir wavelengths @xcite .
also , using contemporaneous interferometric observations of betelgeuse at different wavelengths , @xcite found that hotspots evident in the visible were not present in the near - ir .
observations of betelgeuse at visible wavelengths could represent regions within the layer observed at mid - ir wavelengths .
it is interesting that the mid - ir flux density of betelgeuse , as well as the effective temperatures plotted in figure 3 , and the visual magnitude appear somewhat correlated , particularly in their increase during the 2009 and 2010 epochs .
previous authors @xcite have concluded that the visual variability of betelgeuse is probably indicative of quasi - periodic pulsations , and is too great to be explained by variations in the hotspot distribution .
the layer of material surrounding betelgeuse is a candidate location for the onset of dust formation .
a number of authors have identified the presence of alumina dust directly above the photosphere @xcite .
the derived physical parameters of the layer , such as the temperatures and hydrogen mass densities , are suitable for the formation of amorphous alumina @xcite . while some aspects of the modeled variability of the layer might be due to variations in the photosphere , measurements of changes in the apparent mid - ir appearance of betelgeuse must be largely caused by variations in the layer .
if alumina contributes significantly to the mid - ir continuum opacity , these variations would represent variations in the dust content of the layer .
our estimates of changes in the mass of the layer , on the order of @xmath97 ( see table 4 ) , might indicate evolution in the layer .
the point sources we observe in the layer of material surrounding betelgeuse and their changes suggest that the layer dynamics are anisotropic .
asymmetries observed in the circumstellar environment of betelgeuse are generally linked to the presence of giant convection cells on the stellar surface .
such convection cells are thought to shape the circumstellar environment of betelgeuse by elevating cool photospheric material @xcite and driving large - scale chromospheric motions @xcite .
our results indicate that giant convection cells have a significant role in shaping the betelgeuse layer .
this is because the mechanism that elevates the layer material above the photosphere is clearly anisotropic .
the large intensities of the observed point - source components imply the action of large - scale photospheric features , and their variability between observing epochs demonstrates the transient nature of the cause of the point sources .
changes in the distribution of convection cells on the photosphere of betelgeuse could also influence the overall layer properties . while the nature of the non - uniform brightness distribution of the layer is yet to be determined
, our observations represent evidence for the role of giant convection cells in shaping the inner atmosphere of betelgeuse .
observations of a sample of red supergiant stars with high spatial resolution and wide bandwidth , over multiple epochs , are necessary to form a general picture for the role of molecular layers in mass - loss processes .
more detailed interferometric imaging of betelgeuse at mid - ir wavelengths is also required to understand the true role of giant convection cells in shaping the layer .
simultaneous visible , near - ir , and mid - ir interferometry could also provide interesting information on correspondences between different regions of the layer .
interferometric studies with both high spatial resolution and high spectral resolution are planned for the isi @xcite , and have been performed at near - ir wavelengths by @xcite and @xcite .
such investigations are aimed at revealing the molecular compositions , kinematic structures and true extents of red supergiant layers .
more detailed modeling of layer opacities is also required , taking into account contributions from molecular absorption and emission , dust and electron - hydrogen atom collisions .
upcoming optical and infrared interferometer systems , such as the magdalena ridge observatory interferometer and the matisse instrument at the very large telescope interferometer , as well as high spatial resolution observations with the next generation of extremely large telescopes could finally clarify our understanding of red supergiant mass loss .
the isi has been used to conduct observations of betelgeuse over 2006@xmath22010 in the mid - ir continuum with high spatial resolution . by fitting visibility and closure phase measurements with image models ,
the apparent stellar size was measured at each observing epoch , and up to two point sources were sometimes detected at the edge of the stellar disk .
the results were interpreted in terms of an absorbing layer surrounding an assumed 3641k photosphere .
we have shown that :    1 .
observations during every epoch can be modeled in terms of a layer surrounding the photosphere , with optical depths of @xmath11 at 11.15@xmath0 m wavelength , outer sizes between 1.16@xmath7 and 1.36@xmath7 , and temperatures between 1900k and 2800k .
electron - hydrogen atom collisions contribute to the opacity of the layer .
3 .   assuming a non - varying photosphere , the layer exhibits significant variability .
4 .   the layer has a non - uniform intensity distribution that also varies from year to year .
w. fitelson , b. walp , c.s .
ryan , d.d.s .
hale , a.a .
chandler , k. reichl , r.l .
griffith , v. toy , as well as many undergraduate researchers all participated in these observations , and their excellent help is greatly appreciated .
this research made use of the simbad database .
fitting programs used the mpfit idl routines of c. b. markwardt .
we are grateful for support from the gordon and betty moore foundation , the office of naval research , and the national science foundation .
adams , w.  s. , & maccormack , e.  1935 , , 81 , 119 antia , h.  m. , chitre , s.  m. , & narasimha , d.  1984 , , 282 , 574 bester , m. , danchi , w.  c. , hale , d. , townes , c.  h. , degiacomi , c.  g. , mekarnia , d. , & geballe , t.  r.  1996 , , 463 , 336 bloemhof , e.  e. , danchi , w.  c. , & townes , c.  h.  1985 , , 299 , l37 bonneau , d. , & labeyrie , a.  1973 , , 181 , l1 boyajian , t.  s. , et al .
2009 , , 691 , 1243 buscher , d.  f. , baldwin , j.  e. , warner , p.  j. , & haniff , c.  a.  1990 , , 245 , 7p chiavassa , a. , plez , b. , josselin , e. , & freytag , b.  2009 , , 506 , 1351 chiavassa , a. , haubois , x. , young , j.  s. , plez , b. , josselin , e. , perrin , g. , & freytag , b.  2010 , , 515 , a12 dalgarno , a. , & lane , n.  f.  1966 , , 145 , 623 danchi , w.  c. , bester , m. , degiacomi , c.  g. , greenhill , l.  j. , & townes , c.  h.  1994 , , 107 , 1469 de beck , e. , decin , l. , de koter , a. , justtanont , k. , verhoelst , t. , kemper , f. , & menten , k.  m.  2010 , , 523 , a18 dirks , u. , schirrmacher , v. , & sedlmayr , e.  2008 , , 491 , 643 freytag , b. , steffen , m. , & dorch , b.  2002 , astronomische nachrichten , 323 , 213 gehrz , r.  1989 , in interstellar dust , eds .
l.  j.  allamandola & a.  g.  g.  m.  tielens , iau symp .  135 , pg .  445 gilliland , r.  l. , & dupree , a.  k.  1996 , , 463 , l29 gustafsson , b. , edvardsson , b. , eriksson , k. , jrgensen , u.  g. , nordlund ,  . , & plez , b.  2008 , , 486 , 951 hale , d.  d.  s. , et al .  2000 , , 537 , 998 harper , g.  m. , brown , a. , & lim , j.  2001 , , 551 , 1073 harper , g.  m. , brown , a. , & guinan , e.  f.  2008 , , 135 , 1430 harper , g.  m. , richter , m.  j. , ryde , n. , brown , a. , brown , j. , greathouse , t.  k. , & strong , s.  2009 , , 701 , 1464 hartmann , l. , & macgregor , k.  b.  1980 , , 242 , 260 haubois , x. , et al .
2009 , , 508 , 923 hebden , j.  c. , eckart , a. , & hege , e.  k.  1987 , , 314 , 690 heger , a. , & langer , n.  1998 , , 334 , 210 humphreys , r.  m.  1986 , in luminous stars and associations in galaxies , eds .  c.  w.  h.  de loore , a.  j.  willis , & p.  laskarides , iau symp .  116 , pg .  45 jennings , d.  e. , & sada , p.  v.  1998 , science , 279 , 844 jennison , r.  c.  1958 , , 118 , 276 josselin , e. , & plez , b.  2007 , , 469 , 671 kervella , p. , verhoelst , t. , ridgway , s.  t. , perrin , g. , lacour , s. , cami , j. , & haubois , x.  2009 , , 504 , 115 kiss , l.  l. , szab , g.  m. , & bedding , t.  r.  2006 , , 372 , 1721 lim , j. , carilli , c.  l. , white , s.  m. , beasley , a.  j. , & marson , r.  g.  1998 , , 392 , 575 michelson , a.  a. , & pease , f.  g.  1921 , , 53 , 249 monnier , j.  d. , geballe , t.  r. , & danchi , w.  c.  1998 , , 502 , 833 ohnaka , k.  2004 , , 421 , 1149 ohnaka , k. , et al .
2005 , , 429 , 1057 ohnaka , k. , et al .
2009 , , 503 , 183 ohnaka , k. , et al .
2011 , arxiv:1104.0958 pease , f.  g.  1922 , , 34 , 346 perrin , g. , ridgway , s.  t. , coud du foresto , v. , mennesson , b. , traub , w.  a. , & lacasse , m.  g.  2004 , , 418 , 675 perrin , g. , ridgway , s.  t. , verhoelst , t. , schuller , p.  a. , coud du foresto , v. , traub , w.  a. , millan - gabet , r. , & lacasse , m.  g.  2005 , , 436 , 317 perrin , g. , et al .
2007 , , 474 , 599 plez , b.  2010 , in hot and cool : bridging gaps in massive star evolution , eds .  c.  leitherer , p.  bennett , p.  morris , & j.  van loon , asp conf .
ser .  425 , pg .  124 reid , m.  j. , & menten , k.  m.  1997 , , 476 , 327 rosseland , s.  1924 , , 84 , 525 ryde , n. , harper , g.  m. , richter , m.  j. , greathouse , t.  k. , & lacy , j.  h.  2006 , , 637 , 1040 salpeter , e.  e.  1977 , , 15 , 267 schwarzschild , m.  1975 , , 195 , 137 stencel , r.  e. , pesce , j.  e. , & bauer , w.  h.  1989 , , 97 , 1120 tatebe , k. , & townes , c.  h.  2006 , , 644 , 1145 tatebe , k. , chandler , a.  a. , wishnow , e.  h. , hale , d.  d.  s. , & townes , c.  h.  2007 , , 670 , l21 townes , c.  h. , wishnow , e.  h. , hale , d.  d.  s. , & walp , b.  2009 , , 697 , l127 tsuji , t.  2000 , , 538 , 801 tsuji , t.  2000 , , 540 , l99 tsuji , t.  2006 , , 645 , 1448 van belle , g.  t. , et al .
2010 , arxiv:1012.1146 verhoelst , t. , et al .  2006 , , 447 , 311 verhoelst , t. , van der zypen , n. , hony , s. , decin , l. , cami , j. , & eriksson , k.  2009 , , 498 , 127 weiner , j. , danchi , w.  c. , hale , d.  d.  s. , mcmahon , j. , townes , c.  h. , monnier , j.  d. , & tuthill , p.  g.  2000 , , 544 , 1097 weiner , j. , hale , d.  d.  s. , & townes , c.  h.  2003 , , 588 , 1064 weiner , j. , hale , d.  d.  s. , & townes , c.  h.  2003 , , 589 , 976 weiner , j.  2004 , , 611 , l37 wilson , r.  w. , dhillon , v.  s. , & haniff , c.  a.  1997 , , 291 , 819 wishnow , e.  h. , mallard , w. , ravi , v. , lockwood , s. , fitelson , w. , wertheimer , d. , & townes , c.  h.  2010 , , 7734 woitke , p.  2006
, , 452 , 537 woitke , p.  2007 , in why galaxies care about agb stars : their importance as actors and probes , eds .  f.  kerschbaum , c.  charbonnel , & r.  f.  wing , asp conf .
ser .  378 , pg .  156 woodruff , h.  c. , et al .
2009 , , 691 , 1328 yoon , s .- c . , & cantiello , m.  2010 , , 717 , l62 young , j.  s. , et al .
2000 , , 315 , 635