Instruction: Write an abstract for the following acadamic paper.

Input:
Paper: we have recorded the flight of a fly during take off and landing using digital high speed photography . it is shown that the dynamics of flexible wings are different for these procedures . during this observation fly flew freely in a big box and it was not tethered .
Abstract: in this fluid dynamics video , we demonstrated take off and landing of a fly . the deformation of wings is in focus in this video .

Paper: the powerful infrared source and oh maser vy canis majoris is one of a few very luminous cool hypergiant stars that define the empirical upper luminosity boundary in the hr diagram @xcite . it is thought to represent a very rapid stage of evolution characterized by a high mass loss rate . the red supergiant stage is a few hundred thousand years long , or less for the most massive stars , but the transition time from main sequence to rsg and later from rsg back to a yellow hypergiant may be only a few thousand years @xcite . these evolved late - type stars are major sources of gas , dust and molecules that contribute to the replenishment of the interstellar medium . vy cma , in particular , is ejecting large amounts of gas and dust at a prodigious rate . with its extended and complex circumstellar nebula , it offers a unique opportunity to study its high mass loss episodes and the detailed kinematics of this process . vy cma is a special case even among the rare hypergiants . it is associated with a region of recent star formation , the young cluster ngc 2362 and h ii region sharpless 310 , and the dark cloud l1667 , and has been considered a pre  main sequence star ( herbig 1970a , b ; lada & reid 1978 ) as well as an evolved massive star @xcite . the lack of li i @xmath26707 @xcite plus its powerful double  peaked maser emission , like other evolved stars , is strong evidence that it is a massive post  main sequence object . the measured @xmath3c/ @xmath4c ratio @xcite however is intermediate between what is measured for normal m supergiants and what is expected for the ism . at a distance of 1.5 kpc @xcite it is one of the most luminous red supergiants known ( @xmath5 ) . with its very visible asymmetric nebula 10@xmath6 across and high mass loss rate ( @xmath7 ) @xcite , vy cma is also one of the most important for understanding the high mass loss episodes in massive star evolution . interferometric maps of the oh and h@xmath8o emission reveal a complex structure interpreted as either an expanding disk or bipolar outflow ( bowers , johnston , & spencer 1983 ; richards , yates & cohen 1998 ) with expansion velocities of 30  40 km s@xmath1 . speckle imaging @xcite revealed an asymmetric resolved component 0@xmath91 across , elongated approximately north - south . a geometry consistent with the maser maps and the infrared images @xcite is an axially symmetric model with a slightly flattened optically thick envelope or an equatorial disk and a bipolar flow / rotational axis aligned roughly northeast - southwest and inclined 15 to 30@xmath10 to our line of sight @xcite . however , our multi - wavelength hst / wfpc2 images of vy cma revealed a complex asymmetric distribution of material with no clear axisymmetric structures ( see figures 1,2 , & 3 in smith et al . ( 2001 ) , and with no evident disk - like structure in the reflected optical light extending out to 2@xmath6 or more from the star as proposed by herbig ( 1970a , b ) , efstathious & rowan - robinson ( 1990 ) , and monnier et al ( 1999 ) . instead , the circumstellar nebula is dominated by a prominent , nebulous arc to the northwest and two bright filamentary loops or arcs to the southwest of the star , plus some relatively bright condensations or clumps of dusty knots near the star , all of which are evidence for multiple and asymmetric mass loss events . the apparent random orientations of the arcs suggest that they were produced by localized events not strongly aligned to the star s axis or equator . the models for disks are based on modelling the spectral energy distribution and therefore may be due reflected light from the prominent tail and the bright condensations near the star rather than a true disk . the arcs and loops also do not have the appearance of continuous outflow streams as suggested by monnier et al . ( 1999 ) , and are too massive to have been ejected by radiation pressure on dust . the initial ejection must have been induced by other processes . we therefore speculated that the arcs may be expanding loops caused by localized activity on the star s ill - defined surface , involving magnetic fields and convection . although this suggestion was conjectural , large starspots do occur on the surface of @xmath11 ori ( gilliland & dupree 1996 ; kluckers et al . 1997 ; uitenbroek , dupree & gilliland 1998 ; dupree , lobel & gilliland 1999 ) with accompanying outflows and chromospheric activity @xcite . non - radial pulsational instability is potentially an alternative ejection mechanism , but the distinction between pulsation and convection may be blurred in vy cma .    to learn more about the morphology , kinematics , and possible origin of vy cma s complex ejecta , we have obtained long - slit spectra to map the emission and absorption lines across several structures , including the prominent nw arc and the two bright filamentary loops to the southwest . our results reveal the presence of a complex pattern of velocities in the the strong k i emission line and a significant velocity gradient across the dusty northwest arc . we find that ; + @xmath12 most of the gas and dust in the nebula , however , is nearly stationary or moving very slowly with respect to the star . this is not only surprising but may be difficult to explain for an evolved , massive star . + @xmath12 embedded in this nebula are several discrete structures that are kinematically distinct from the surrounding stationary nebulosity . the nw arc and the two prominence  like arcs are expanding at 40  50 km s@xmath1 relative to the star . + @xmath12 they were ejected in separate events over several hundred years from different locations on the star . + @xmath12 a strong absorption line in the reflected spectrum in the ejecta is identified with ca i ] . although absorption from a semi - forbidden line seems unlikely , we demonstrate that the identification is probably correct . + the strong k i emission in the star and the ejecta has never been adequately explained . in the appendix , we discuss a possible mechanism for forming the k i emission and other rare emission lines , including ca i ] .    in the next section we briefly describe the observations and the data reduction procedures . the spectrum of the central star and the origin of its peculiar emission lines are discussed in  3 and in the accompanying appendix . the kinematics of the gas and dust , the expansion of the arcs and loops , and the evidence for stationary ejecta are described in  4 . in the final two sections we summarize the evidence for localized high mass loss events and comment on the perplexing nature of vy cma . spectra were obtained at four slit positions across the nebula ( figures 1a and 1b ) with the hires echelle spectrograph on the keck1 telescope on mauna kea in december , 2002 . the resolution was 40000 with a dispersion of .06 per pixel and a spatial scale of 0@xmath919 per pixel . our spectral coverage from @xmath136410 to 8727 in 15 orders with inter  order gaps was chosen to measure the velocities of known emission lines plus several absorption lines in the 8400  8700 region . we used the b3 decker which gives a 14@xmath6 long slit with a 0@xmath957 slit width . three of the slit positions are parallel across the nebula and cross clumps of knots , filaments , and the nw arc . one of the slit positions(i ) includes the star . the fourth slit ( v ) crosses the two prominence  like arcs , 1 and 2 , in the outer parts of the ejecta to the southwest and the clump of knots south of the star . the seeing was estimated as 0@xmath98 or better throughout the night . the observations are summarized in table 1 . all of the data reductions and extractions were performed using iraf . the frames were trimmed , bias subtracted and flatfielded . calibration and object exposures were combined for each slit position . using the quartz exposures for reference , we extracted each echelle order from the combined frames using the ` strip ' option in the iraf task _ apsum_. the dispersion corrections were determined for each order separately . no sky subtraction was done because emission and absorption lines are still apparent in the most distant apertures , beyond the visible nebulosity . the spectra were then extracted in 1 long apertures deliberately placed across prominent features including the nw arc , the filamentary arcs , the knots and the embedded star . additional extractions were then made along the slit keeping 0@xmath94 separation between apertures . the positions of the apertures are summarized in table 2 ; the positions are given along each slit with respect to the point nearest the star , the radial distance of each aperture from the star , and the position angle from the star . figures 1a and 1b show the positions of the extraction apertures for slits i and ii , and for slits iii and v , respectively . four extractions were made along slit iii , across specific features listed in table 2 plus four additional apertures to the end of the slit . the wavelength calibration was then applied to each extracted spectrum . the separate exposures were averaged to remove cosmic rays , except for those along slit i which were summed because only two exposures were available . vy cma probably represents the most extreme known case of stellar activity , in the same sense as , e.g. , @xmath11 ori . its inner wind ( chromosphere ? ) has a density comparable to that of @xmath11 ori , but the volume , mass , and energy are larger by almost two orders of magnitude . the photosphere is appreciably cooler than betelgeuse and is about one order of magnitude more luminous . the observed emission lines  particularly the strong k  i and the ca  i ] emission features  differ from @xmath11 ori and have not been explained . in the next few paragraphs we review what one _ expects _ to see in this situation . the photospheric temperature , luminosity , and radius are thought to be , respectively , around 3000 k characteristic of an m5 type supergiant , @xmath14 , and @xmath15 @xmath0 12 au . ( if the wind is opaque , then @xmath16 is larger than the true stellar radius @xmath17 , the sonic point in the flow . ) the usually quoted mass - loss rate and wind speed , @xmath18 @xmath19 yr@xmath1 and 35 km s@xmath1 , imply that the wind density at @xmath20 au is @xmath21 @xmath22 or @xmath23 g @xmath22 .    given those parameters , we expect dust grains to condense near radius @xmath24 au , at temperatures around 1000 k. if the resulting dust - to - gas ratio is `` normal '' , then , at red wavelengths @xmath25  , the opacity is @xmath26 @xmath27 g@xmath1 and the average optical thickness due to dust outside @xmath28 is @xmath29 . observations of scattered light at radii @xmath30 2000 au seem consistent with this amount of dust @xcite . with such a large average @xmath31 , one might naively expect the object to be seen only as a mid - ir source . in reality , however , the dust should be very inhomogeneous . even if the inner wind is smoothly isotropic ( which we doubt ) , radiative acceleration acting on the dust near @xmath28 is unstable in a rayleigh - taylor - like manner ; the ensuing flow must be clumpy with a fairly small size scale . = wind speed , @xmath32 = sound speed in the wind , and @xmath33 = acceleration due to radiation acting on the dust ; then the three most obvious characteristic size scales for instabilities near @xmath28 are @xmath34 , @xmath35 , and @xmath36 . for vy cma , the last of these is of the order of 510 au while the others are much less ; so the expected condensations in the region just outside @xmath28 should be , typically , no larger than 10 au , and very likely smaller . ] therefore we should not be surprised to see visual- or red - wavelength light from the star , escaping along low - density paths through a very inhomogeneous dusty shell . since the material is expanding , and a typical photon may be scattered one or more times before it escapes , the emergent spectrum is likely to be red - shifted by an amount comparable to the wind speed @xmath37 , and somewhat blurred in wavelength . van blerkom and van blerkom ( 1978 ) described this process for a constant wind speed , and the effect may be supplemented by radiative acceleration of material in the regions most favorable for photon escape , i.e. , in relatively low - density gas between condensations . after escaping from the inner dusty region , some photons can be scattered sideways by dust in the outer ejecta , causing us to see an extended , inhomogeneous halo around the bright central region . any photon that is scattered toward an appreciably non - radial direction is doppler shifted by a net amount @xmath38 , where @xmath37 is the outward wind speed and @xmath39 is ( @xmath40scattering angle ) . thus we expect the spectrum reflected in the halo to be more red - shifted than that seen when we observe the central source . since a line of sight through the halo typically samples a range of @xmath39-values , the reflected spectrum should be blurred more than the central spectrum .    generally speaking , observations of vy cma confirm the above expectations . at wavelengths around 7000  , very roughly 2 percent of the star s light escapes through the dust ( see measurements reported in smith et al . 2001 ) , suggesting typical optical depths of the order of 4 or 5 , not 50 , along the photon escape routes . absorption lines in the central emergent spectrum are doppler - shifted by about + 25 km s@xmath1 relative to the supposed systemic velocity , and are somewhat broader than one expects for such a low - gravity star . when reflected in the halo , the same features have larger doppler shifts and are more blurred ( fig . 6 ) . throughout this paper we shall note additional details that do _ not _ fit the above picture in straightforward ways . the k i resonance lines ( @xmath13 7665 , 7699 ) are not only the strongest emission lines in the spectrum of vy cma @xcite , but these lines , rarely observed in emission , are stronger in vy cma than in any other known late - type star ( see bernat & lambert 1976 , guilain & mauron 1996 ) . in addition to k i , strong emission is also present in the na i `` d '' lines , ca i ] @xmath26572 , and rb i ( @xmath137800 , 7948)@xcite . narrow band heads of sco ( wallerstein in hyland et al 1969 and herbig 1974 ) , and band heads of tio and vo , plus lines of fe i , ti i and cr i @xcite are also in emission . our wavelength coverage includes k i @xmath27699 , ca i ] @xmath26572 and rb i @xmath27948 , but not the na i `` d '' lines . in addition to the emission features , there are strong absorption lines of ca ii , fe i , and ti i in the wavelength region from approximately 8400 to 8700 that are relatively free from molecular absorption bands so they can be confidently used for radial velocity measurements . in our data , the central star is in aperture 2 on slit i ( fig . 1a ) , and examples of the emission and absorption lines in the ejecta are shown in figures throughout the paper . the spectrum at the star does not show any significant differences from those published earlier and the velocity measurements for the star summarized in table 3 are consistent with previously published results . all velocities in this paper are heliocentric .    in figure 2 we show the k i and na i line profiles on the star from a hires spectrum obtained by george herbig in 2000 that includes the d lines . both the k i and na i emission have peak fluxes several times as high as the continuum , and not surprisingly , their doublet ratios show that they are optically thick . the minima in the strong p cygni profiles have the same velocity at @xmath0 15 km s@xmath1 . for the k i lines , this is a shift of @xmath0 30 km s@xmath1 from the emission peaks which have a mean velocity of 45 km s@xmath1 , consistent with our measured k i and ca i ] emission velocities . the na i emission peaks however occur at @xmath41 60 km s@xmath1 which agrees with our mean absorption line velocity at the star ( table 3 ) . it has been recognized for some time @xcite that the absorption line velocities are not only different from the k i and ca i ] emission lines , but also deviate significantly from the velocity of the system inferred from the oh and h@xmath8o maser observations , v@xmath42 20  22 km s@xmath1 ( v@xmath43 @xmath41 35  37 km s@xmath1)(bowers , johnston & spencer 1983 ) . herbig ( 1970a ) first suggested that the stellar absorption lines are red - shifted and broadened due to scattering by an expanding dust shell or equatorial disk , and as mentioned above , this was modelled by kwok ( 1976 ) and van blerkom and van blerkom ( 1978 ) . note that our velocities for the k i and ca i ] emission lines at @xmath41 41 km s@xmath1 are significantly less than the mean absorption line velocity . this difference is usually attributed to the location where the emission lines are formed , outside the inner radius of the dusty shell or disk .    the extremely strong k i emission and other rare emission lines such as semi  forbidden ca i ] are difficult to explain . in the appendix we discuss problems associated with the formation of the k i lines in vy cma . some obvious excitation processes such as collisional and radiative recombination are eliminated , and while resonant scattering is the most promising explanation , it requires a special geometry to cause strong net emission . we propose a model for the production of the very strong k i emission observed at the star , but problems associated with resonant scattering remain , especially in the ejecta . the strong k i emission line is our best tracer of the gas in vy cma . it shows considerable variability with multiple emission peaks throughout the ejecta ( see also smith 2004 ) . examples of the line profiles along slits i , ii and iii are shown in figure 3 , with their velocities summarized in table 4 . the emission profile on the star is single ( slit i , ap 2 ) with a strong p cygni absorption feature formed in its expanding wind . the p cygni profile is present in the reflected spectrum near the star and is gradually red - shifted due to scattering by the dust expanding outwards . off the star , the k i profiles along the slits are much more complex , often with two or more recognizable emission peaks . in the outermost apertures the profiles become single and narrow . an emission peak with the expected systemic velocity is observed throughout much of the ejecta to the west and northwest of the star . it has a nearly constant velocity at increasing distances and is essentially motionless relative to the star . this and other evidence for stationary gas and dust in vy cma s ejecta are discussed in  4.3 . interestingly , a broad and shallow k i absorption feature to the red of the emission profile is present at positions where the slits cross the dusty condensations and nebulous arcs primarily to the west and nw of the star . this may be an inverse p cygni profile , but the broad absorption minimum has a red - shifted velocity , typically @xmath41 86 km s@xmath1 ( table 5 ) like the other reflected absorption lines at the same positions ( see  4.1 ) . note that on average , the k i emission greatly exceeds the strength of the absorption in our data . one would normally expect emission to balance absorption for a resonance  scattered line ( see the appendix ) . the semi - forbidden ca i ] @xmath26572 emission line ( figs . 4 and 5 ) is visible only in the extractions on and near the star . the line is very narrow , does not show the multiple peaks observed in the k i line , and has essentially a constant doppler velocity at @xmath0 41 km s@xmath1 like that for the k i emission on the star ( table 4 ) . a strong , reflected absorption line which we attribute to ca i ] absorption is present in the ejecta . although absorption in a semi - forbidden line seems unlikely , we demonstrate later ,  4.1 , that this identification is almost certainly correct . another rare emission line , rb i at @xmath2 7948 , is present in the reflected spectra , but is relatively weak and broad . h@xmath11 emission in vy cma is usually reported as very weak and variable . we find that the relative strength of the h@xmath11 emission depends on position in the nebula ( see figs . 4 and 5 ) . in the inner region of the ejecta it is weak and broad with an associated broad shallow absorption . the narrow h@xmath11 emission line visible in the outermost extractions along all of the slit positions is nebular from the nearby h ii region sharpless 310 . vy cma lies just to the west of an arc of nebular emission . weak nebular emission from [ n ii ] and [ s ii ] is also present in the spectra at the same positions . the velocity of this h@xmath11 line at @xmath0 35  37 km s@xmath1 agrees with the systemic velocity and the velocities of the k i lines measured at the same positions , thus confirming the close association of vy cma with the cluster and h ii region . the reflected and scattered spectrum in the circumstellar ejecta reveals a remarkably complex behavior of both the emission and absorption lines . the k i emission line profiles and their doppler velocities vary with position in the nebula as the slits cross recognizable knots , arcs , and filamentary features in the image . they also indicate apparently stationary gas in the same projected locations , which is difficult to understand for the reasons noted below . the reflected absorption lines show a strong velocity gradient across the prominent nebulous arc to the nw of the star , while the same absorption lines along slit v show virtually no variation in doppler velocity . these phenomena are discussed in the following subsections . the absorption lines from the obscured star and the dusty ejecta are reflected or scattered by the embedded dust grains and reveal both the relative motions of the dusty clouds and the velocity dispersions within them . the mean velocities of the strong absorption lines of fe i , ti i , and ca i in the far  red spectrum are summarized in table 5 and examples of their profiles are shown in figure 6 . a few of the lines in apertures on and near the nw arc are double or slightly asymmetric . for example , four of the absorption lines in i ap 3 are double so we used the blue components , which agree with the velocities of the single lines , for the mean velocity in table 5 . in the extractions immediately on either side of the nw arc , a few of the absorption lines also have asymmetric profiles . only the single lines and well defined absorption minima of the some of the asymmetric lines were used for the mean velocity . the asymmetries may be due to additional dusty clouds along the line of sight or to the expected gradient along the nw arc . the absorption line profiles from the dusty arcs are broad and shallow compared to the narrower profiles on and near the star . along slit i , for example , the equivalent widths of the absorption lines show only a small decrease from the star to ap 3 and ap 4 ( nw arc ) , but the fwhm of the same profiles increase by a factor of 2 to 3 due to the scattering and motions of the grains . to determine the velocity dispersion of the particles responsible for the broadening , we smeared the stellar profiles of the fe i and ti i absorption lines with a gaussian . we used a range of values of @xmath44 from 0.5 to 1.6 and compared the smeared stellar lines with the observed profiles along the slits . the best fit varied slightly from line to line and of course for the different apertures , but for those profiles across the nw arc , the required smearing corresponds to a fwhm of 85 to 120 km s@xmath1 . the total velocity range will be even larger . the dusty scattering regions in this visible feature thus have a relatively high velocity dispersion in addition to their net motion with respect to the star . the mean velocities for the absorption lines are plotted in figure 7 with respect to the distance along each slit from the point nearest the star . they show a large redshift as the three parallel slits ( i , ii and iii ) cross the dusty arcs to the west and nw of the star , and then decrease in the outer nebulosity beyond the prominent nw arc . _ the profiles across the nw arc are significantly red - shifted with a velocity difference of @xmath41 + 50 km s@xmath1 relative to the star , suggesting a net motion of the expanding arc of this order . _ this difference is presumably the `` moving mirror '' redshift , @xmath38 , where u is the outward expansion of the reflecting material relative to the star and @xmath39 represents the direction of the flow relative to the plane of the sky ( @xmath45 for motion toward us ) .    if the nw arc represents the material from a single ejection event we expect the expansion velocity to be proportional to distance along the arc . this expected trend is not immediately obvious in the observed velocities . it is probably blurred in our data because of the size of the extraction apertures , although some subtle evidence for this effect may be present in the velocities between slits ii and iii in apertures directly on the arc ; iii ap 4 is slightly further from the star and has a higher velocity than ii ap 5 . moreover , in figure 6 , one can see a velocity component on the short wavelength side of the lines at position ii ap 4 and on the long wavelength side at ii ap 6 . the velocity difference is consistent with the uniform expansion expected for the nw arc .    unlike most other locations in the ejecta , appreciable k i _ absorption _ is observed in the extractions on and near the nw arc . it also shows an appreciable redshift at @xmath0 85 km s@xmath1 ( table 5 ) , and the line is also quite broad and asymmetric . unlike the other absorption lines , the k i velocity is approximately constant across the arc . this is probably due to the strong emission line covering or overlapping part of the absorption feature . the red wing of the k i absorption , however , does show a redshift as the slits cross the arc , from 135 to 151 km s@xmath1 along slit ii .    using the velocities from the absorption lines at position ii ap 5 , directly on the nw arc , we apply its  moving mirror " redshift of 54 km s@xmath1 to the k i velocity at this position , and get 31 km s@xmath1 , close to the 35 km s@xmath1 systemic velocity of the star . we therefore assume that the k i absorption in some sense represents the spectrum of the star as seen from the nw arc . in figure 8 we show the k i and ca i ] line profiles shifted by ( 35 + 54 ) km s@xmath1 , one of the ordinary absorption lines , and the k i emission at the star s position , i ap 2 . the k i and ca i ] absorption features are centered close to zero and are thus nearly in the rest frame of the star . the ti i absorption line is shifted an additional 25 km s@xmath1 in figure 8 , because as explained in  1 , it arises inside the circumstellar dust formation region . since the k i and ca i ] absorption features do not need the extra 25 km s@xmath1 shift , we conclude that they occur _ more than 100 au from the star , outside the ring or shell of dust . _ @xmath26573 is a semi - forbidden line , one would not have expected it to appear in absorption as observed . but the wavelength coincidence is essentially perfect , and we have found no other likely identification for the observed feature . thus we must conjecture that some atomic - physics trick makes the oscillator strength much larger than one would expect . ] the red wing of the k i absorption , at @xmath0 60 km s@xmath1 with respect to the absorption minimum , indicates the presence of a substantial dispersion within the expanding gas of 100 km s@xmath1 or more , about the same as observed for the other absorption lines at the same locations . this means that either the nw arc has a large disperison in its expansion velocity or a large extent along the line of sight , or both . note that the `` terminal velocity '' of the k i p cyg profile as seen at the star is also about 45  50 km s@xmath1 . thus , the moving - mirror redshift at the nw arc , the p cyg `` terminal velocity '' , the red wing of the k i absorption , and the reflected line widths _ all suggest flow speeds of about 50 km s@xmath1_. expansion velocities greater than 35 km s@xmath1 are required by the absorption line widths and there is no indication of velocities much greater than 70 km s@xmath1 relative to the star . an expansion speed of 50 km s@xmath1 implies that @xmath39 must be near zero so the nw arc is basically moving across our line of sight . this is not a surprising result given its orientation and appearance within the nebula . furthermore , the limits on the expected expansion velocity discussed above , 35 to 70 km s@xmath1 , would restrict @xmath39 to between -25@xmath10 and + 17@xmath10 . assuming mostly transverse motion for the arc , at its distance from the star of about 3@xmath6 or 4500 au at 1.5 kpc from the sun , the material in the _ nw arc was ejected about 400 years ago _ ( -100 , + 200 yrs . ) . we get a somewhat different result for slit i which crosses the nw arc at a section where it appears to be bending back toward the star . at i ap 4 the doppler velocities of the absorption lines are slightly lower and the  moving mirror " redshift relative to the star is 40 km s@xmath1 . assuming 50 km s@xmath1 for @xmath37 , @xmath39 is roughly @xmath46 for this position . thus , this section of the nw arc may actually be the nearer side with a small component of motion toward us . the corresponding k i emission , positions ii ap 5 , ii ap 6 , iii ap 4 and iii ap 5 , has two definite velocity peaks ( figure 3 ) which can also be seen as two separate emission features in the image of slit iii in figure 9 . the shorter wavelength one presents some problems that are discussed in the next section . the red - shifted component with velocities of 60 to 70 km s@xmath1 exhibits a small velocity gradient across the arc as one would expect for an ejection . recognizing this trend , we can identify a corresponding bump in the line profiles at ii ap 4 and iii ap 3 at @xmath41 52 km s@xmath1 that probably corresponds to the same gas . since this feature spatially coincides with the nw arc and is expanding , we presume that this is the associated gas . however , the doppler velocity , at ii ap 5 , for example , corresponds to a shift of 25 km s@xmath1 relative to the star , half that for the absorption lines . the slope of the velocity gradient , 5 km s@xmath1 per 1000 au yields an expansion age of 950 years . if this emission is due to resonant scattering , then it is moving away from us at an angle of @xmath41 30@xmath10 which sems unlikely given the results from the absorption lines in the nw arc . if , on the other hand , it is reflected emission , it will be directed toward us at @xmath41 30@xmath10 since the redshift is smaller than for the absorption lines . the expansion age for the emitting gas is then about 550 years , comparable to the above results for the nw arc . slit positions iii ap 3 and ii ap 4 cross a filamentary , somewhat twisted arc just west of the star , between it and the nw arc . if we apply the same procedure described above for the nw arc with the k i and ca i ] absorption lines , we find that the moving mirror velocity for the dust is 45 km s@xmath1 and the flow speed is 55 km s@xmath1 with @xmath39 of @xmath47 . at its distance from the star , this material was ejected only 200 years ago . in contrast with the results for the ejecta to the west and nw of the star , the doppler velocities of the reflected absorption lines along slit v are nearly constant ( table 5 ) , show no gradient with position ( figure 7 ) , and no significant velocity shift with respect to the velocity at the star . this suggests that the reflecting material is basically stationary relative to the star . this and other evidence for surprisingly stationary gas and dust in the circumstellar nebula is discussed in  4.3 .    in the next section we discuss the significant velocity shifts in the scattered light from some of the k i emission lines in the ejecta and additional evidence for ejection of discrete features in the nebula . the k i emission line often shows a complex profile with multiple peaks that vary with position in the ejecta . if the k i emission lines in the ejecta are due to resonant scattering , the doppler shifted and multiple peaks observed at different places must be due to gas moving at these different velocities thus permitting us to determine the expansion of the nebula and motions within it . however reflection by the dust can also contribute to the observed profiles . in this case we would expect to see either a single profile like that on the star , or an emission peak , with a velocity shift between the emission velocity at the star and the absorption lines at the same position . a few of the k i profiles closest to the star may be due to reflection or have a contribution from reflection . these appear to be either single or lack obvious multiple peaks even though the lines are quite broad . some possible examples are the k i emission at i ap 1 , v aps . 1 and 2 and iii ap 2 in figures 3 and 11 . for the multi - peaked profiles further from the star , it is difficult to determine if the emission is due to resonant scattering or reflection . both processes are probably occurring .    in figure 10 we show the variation of the emission line velocities along the different slits as we did in figure 7 for the absorption lines . slit v is shown separately . along slits i , ii and iii we see a pattern to the west and northwest of the star in which most of the k i profiles have an emission peak at a velocity near 35 km s@xmath1 ( see table 5 ) . this is very apparent in figure 9 as the strong , straight emission feature in the two - dimensional image of slit iii . whether these emission peaks and single profiles are due to reflection or to resonant scattering , their velocities are consistent with stationary gas , near zero velocity with respect to the star . this possibility is discussed more fully in  4.3 . a second , somewhat weaker emission peak present in extractions on the nw arc was discussed in the preceding section . the changing profile and multiple peaks of the k i emission line along slit v ( figure 11 ) show the most dramatic shifts in doppler velocity with position in the nebula , and if these are due to resonant scattering in the gas , imply the existence of clouds or filaments of gas with a wide range of velocities in the ejecta . the k i line varies from a single profile in ap 1 with a strong p cygni absorption , similar to that viewed on the star , to multi - peaked profiles with large velocity shifts where the slit crosses knots and filaments in the outer parts of the nebula . the k i velocity to the north and east of the star , v ap 1 and ap 2 , is consistent with that observed on the star and may simply be the reflected profile from the star as we mentioned earlier . in apertures south of the star across the dusty knots and arcs , one of the strong emission peaks in apertures 3 to 6 has a velocity like that of the reflected absorption lines ( 60  65 km s@xmath1 ) . this emission component may be reflected by the dust or due to resonant scattering by gas moving @xmath41 30 km s@xmath1 relative to the star . likewise , the blueward peak near 25 km s@xmath1 in the same apertures , may be either reflected or scattered emission and represents a separate flow of emitting gas along the same line of sight . the most striking change occurs across the two outer filamentary arcs ( apertures 5 to 7 ) with the appearance of a strong blue - shifted emission feature . the sharp transition from an emission peak near 20 km s@xmath1 in aperture 5 on arc 2 to the blue - shifted emission across arc 1 at @xmath04 km s@xmath1 , suggests that this is resonant scattering in kinematically and spatially separate gaseous filaments or gas flows moving very differently from the reflecting dust . this shift can be clearly seen in the two - dimensional image in figure 12 where narrow emission appears in the p cyg absorption core and broadens out across arc 1 . the p cyg absorption is either no longer present or has been blocked by the strong blue - shifted emission . this is not a continuous outflow from a single ejection as suggested by smith ( 2004 ) , because these two emission features are resolved in the image of the slit . the velocities of the emission peaks associated with these three features are shown connected by a lines in figure 10 . we can trace this prominent blue - shifted emission feature in the two dimensional image over @xmath0 3@xmath6 in the ejecta and directly across arc 1 in aperture 7 . it corresponds to the strong emission peaks in apertures 6 and 7 and may also be weakly present in ap 8 and possibly in ap 5 as a small bump on the blue side of the profile . in ap 6 broad , weaker emission scattered or reflected in gas at a range of velocities up to 60 km s@xmath1 is also present , but is not observed in ap 7 across the arc . we therefore identify this strong blue - shifted emission feature with resonant scattering by gas in a stream or flow associated with the material in arc 1 . the doppler velocities at 4  8 km s@xmath1 indicate that this material is moving relative to the star at @xmath41 30 km s@xmath1 toward the observer . this of course contrasts sharply with the nearly constant 60 km s@xmath1 velocity of the absorption lines , and suggests a significant kinematic difference between the gas in the arc and the surrounding circumstellar material similar to what we found for the nw arc . although arc 1 has a significant velocity toward us relative to the star , its motion across our line of sight is not known . if we adopt 30  35 km s@xmath1 from the width of the broadened absorption lines ( see  4.3 ) at the same positions , as representative of the motions in the surrounding material for the transverse velocity , then the total motion relative to the star may be as much as 45 km s@xmath1 . the material associated with arc 1 would then have been ejected @xmath41 1000 years ago assuming it is moving in the plane of the sky .    the two prominent emission peaks at 22 and 65 km s@xmath1 across arc 2 , in aperture 5 , are well separated kinematically in the two dimensional image of the slit , but a spatial separation or identification with distinct features is not obvious in the images at these positions . the 65 km s@xmath1 feature agrees with the velocity of the absorption lines in the same aperture , so we assume that this gas is associated with the reflecting dust , and since it is more red - shifted that it comes from background material . the second emission peak may then be due to reflection or resonant scattering by gas in arc 2 . this implies a motion of @xmath0 13 km s@xmath1 relative to the star for the emitting gas along the line of sight . the corresponding motion of the arc across our line of sight is not known and may be different from the assumed 30  35 km s@xmath1 , but with this motion , arc 2 is moving at @xmath41 37 km s@xmath1 and was ejected about 800 years ago . given the uncertainty in these estimates , it is possible that arcs 1 and 2 were ejected at about the same time . however the slit image shows that they are kinematically separate emission features which are also spatially separate in the wfpc2 image . therefore assuming comparable outflow velocities relative to the star , arc 1 would have been ejected earlier , but the ejection times also depend on the projection angle for these structures which is not known . other distinct features in the ejecta such as the two clumps of dusty knots , s and sw , are much closer to the star and presumably represent material ejected much more recently . the evidence for localized mass loss events and stellar activity are discussed in  5 . the presence of gas and dust throughout the nebula , apparently at the same velocity of the star , presents a serious problem because it implies that much of the circumstellar material is nearly stationary or moving very slowly relative to the star . this certainly conflicts with what one expects for a red supergiant with a high mass loss rate .    in the ejecta to the west and nw of the star , the k i emission line has a velocity component consistent with relatively motionless gas , at the systemic velocity of the star , extending over @xmath0 6@xmath6 or @xmath41 9000 au in the nebula ( figure 10 ) . if this emission is due to reflection by dust , we should see a shift in velocity with respect to the rest frame of the star even if the motion is transverse from a uniform flow of gas outwards . this is not observed . if it is resonant scattered emission , and if the emission comes from widespread transparent gas , we would expect to see a broad profile centered at the average velocity of the system . however , much of the nebula where this k i emission is observed is not transparent and , more importantly , the emission components are quite narrow . the fwhm s range from 0.25 to 0.4  or only 10  15 km s@xmath1 wide , compared to @xmath41 1 for the single k i profiles on and near the star . this is much less than the expected expansion speed of the ejecta and suggests that there is very little doppler broadening in the lines . _ the emitting gas is therefore either not moving relative to the star , or all of the material must be moving in the plane of the sky _ i.e. in a narrow plane viewed perpendicularly . the latter seems highly unlikely over such a large area . furthermore , some of the profiles show two emission peaks , one near 35 km s@xmath1 and the other at 41 km s@xmath1 , the velocity at the star , suggesting that the emission is due to both reflection and resonant scattering . the narrow k i emission lines in the outer ejecta also have very broad emission wings that extend to 100  200 km s @xmath1 from the center of the emission line ( figure 13 ) . they are formed too far from the star to be electron scattering wings . it is possible that these emission wings are from more distant gas and dust , expanding outwards from vy cma perhaps in a more spherical distribution from a much earlier unstable period in the star s evolution . there is also some evidence for constant velocities , or zero motion relative to the star , in the reflected light to the northwest of the star . after crossing the near side of the nw arc , the absorption line velocities measured in the outermost apertures along slit i have a velocity like that observed at the star ( see figure 7 ) . along slits ii and iii , however , the absorption line velocities decrease , but stay red - shifted relative to the star by @xmath0 20 km s@xmath1 perhaps because the material , in the same plane , is being pushed by the expanding arc . the ca i ] emission has a nearly constant velocity at 40  41 km s@xmath1 throughout the ejecta . since this is the same as its velocity on the star , we assume that this is reflected emission . the cai ] emission profile does not show the variability with position and multiple peaks readily apparent in the k i emission line . as mentioned in the appendix , the presence of this semi - forbidden line in both emission and absorption is not understood . along slit v we observe absorption lines , including ca i ] , and some k i emission with nearly constant velocities extending over much of the length of the slit and across the nebula . all of these indicators have apparent doppler velocities near 60  65 km s@xmath1 which is also the red - shifted velocity of the absorption lines scattered by the dust envelope around the star . for the absorption lines , this implies that the reflecting dust is not only barely moving but also has little motion with respect to the star . unlike resonant scattered emission , reflected absorption lines are expected to show a redshift relative to the star even if the reflecting dust is moving perpendicular to our line of sight . admittedly , the problem is reduced if the material is moving toward us ; in that case , the `` moving - mirror '' effect is less , but there should still be a measurable redshift . two selection effects , forward scattering by the dust plus extinction , can also serve to reduce the measured redshft because we preferentially see the side of the ejecta toward us . but these considerations do not explain the near constancy of the velocities with distance from the star and the lack of any additional shift in the apertures across arcs 1 and 2 , which are quite red and dusty in the images . the absorption line profiles along slit v are also broadened compared to those on the star , although not as much as across the nw arc ( see figure 6 ) . we compared the gaussian smeared stellar profile with the observed profiles and obtained a range of results which varied with position along the slit . the reflected profiles from apertures 2 and 3 closest to the star , showed only marginal broadening with respect to the star , even though aperture 3 crosses the relatively bright dusty knot or condensation just south of the star . the smeared profiles that best matched those observed from aperture 4 , which crosses several small arcs , have fwhm s corresponding to 20 - 40 km s@xmath1 , and for apertures 5 through 7 across the outer filamentary arcs 1 and 2 , the fwhm s are 65 km s@xmath1 . although this is smaller than measured for the nw arc , it represents significant motions along the line of sight of @xmath48 30  35 km s@xmath1 for the material within these arcs , in contrast with the apparent lack of motion with respect to the star . consequently , we would expect to see a redshift relative to the star of at least 10 km s@xmath1 or so , especially across the arcs even with the above considerations for a reduced redshift . instead , the motion relative to the star is at most @xmath41 5 km s@xmath1 assuming that the material along slit v sees the same reflected and red - shifted spectrum of the star that we see , i ap 2 . it is possible that at different positions in the nebula , we may be observing radiation that escaped from different regions on the star , perhaps through holes or gaps in the dusty envelope . consequently , the gas and dust at these locations may be getting a different view or spectrum of the star . for example , if there is a ne - sw axis of symmetry with a bipolar flow , as has been suggested , then along slit v we may be observing gas and dust that views the star from near the polar direction which would be less obscured if there is a disk or flattened dusty envelope around the star . if gas is escaping through a low density region in the dust envelope , then the reference velocity for the reflected absorption lines would more likely be closer to the systemic velocity than the red - shifted velocity observed at the star . if this is the case , the reflected material may be moving at 30 km s@xmath1 relative to the star . this is just a hypothesis , but it might solve the problem of the zero - velocity gas and dust along slit v. however , it would also require that along most of slit v , we see the star from exactly the same direction or angle which then poses problems for the assumed geometry .    a person unfamiliar with the nature of vy cma would conclude from the above description that a substantial amount of nearly stationary or very slow moving gas and dust exists in its ejecta . this is obviously a paradox given the star s strong wind and probable mass loss history . at this time we have no satisfactory explanation . it is possible that the general flow of the diffuse gas in the extended circumstellar material is extremely slow , but the escape velocity for vy cma is @xmath41 60 km s@xmath1 at its photosphere and 20 km s@xmath1 at the dust formation radius at 100 a.u . conceivably this material has been slowed by some process or it did not have quite the escape velocity and is now falling back . but there is no evident reason for this to happen , other than it would help solve our paradox . the results of our high resolution spectroscopy demonstrate that the three primary arcs in the visual images of vy cma are kinematically distinct from the surrounding stationary circumstellar material . the nw arc and arcs 1 and 2 have expansion speeds of 50 km s@xmath1 to @xmath41 40 km s@xmath1 , respectively . the gas and dust associated with these structures was thus ejected in specific mass loss episodes ; @xmath41 400 years ago for the nw arc and possibly 1000 to 800 years ago for arcs 1 and 2 . for the `` western arc '' we found an expansion age of only slightly more than 200 years . thus we have evidence for a period of considerable activity during the past 1000 years , including some recent ejections . in addition to the arcs , numerous clumps of small knots or condensations are observed in the images relatively close to the star . the dusty clumps to the south(s ) and southwest(sw ) of the star are the brightest , but several very small knots are observed just to the west and northwest of the star . these are clustered quite close , within 0@xmath95 of the star , and must have been ejected quite recently , only @xmath41 70 years years ago assuming an outflow speed of 50 km s@xmath1 .    with the possibility of ejection episodes this recently , there may be some hint of activity in the star s observational record . numerous ground - based observers in the past have often described the appearance of knots of nebulosity that come and go and may vary in position , but the timescales for the motions and variablility seem too short to be related to the expanding ejecta . ( see smith et al 2001 for a discussion and the references ) these variations could be due to changing illumination of the dusty ejecta by radiation escaping from gaps or low density regions in the envelope @xcite . vy cma s historical light curve @xcite shows some very intriguing variability . the most significant change is its decline from m@xmath49 @xmath41 6.5 mag . beginning in 1872 to @xmath41 8 mag . by 1880 with some brief brightening episodes superposed . whether this 1.5 mag decrease was due to a change in bolometric correction or to the formation of dust , is not known , but it may correspond to the appearance of several knots and small arcs about 100 to 200 years ago . the star apparently then slowly faded to about 8.5  9th magnitude , and although the record is spotty , it may have had periods of greater variability in the 1910 s and again in the 1930 s . the more modern record based on the aavso light curve and published photometry , shows that vy cma displayed considerable variability of 1 to 2 magnitudes from about 1984 to 1995 with oscillations lasting several 100 to a thousand days that resemble the pulsations in irregular variables . since then it has been relatively steady at about 8.5 mag . the nw arc is the most clearly defined structure or mass concentration within the ejecta and is therefore the best candidate for a mass estimate . in our imaging paper , we estimated its total mass at @xmath50 m@xmath51 from the surface brightness in a well resolved 2 arcsec@xmath52 section . in this estimate we assumed an optical depth of unity . this assumption plus the evidence from the velocity gradient across the arc that more material may be participating in the outflow than is apparent in the images , suggests that this may be an underestimate for the arc s total mass . the dynamical timescale for vy cma is on the order of three years as are the episodes of large variability in its light curve . so whether the observed variability and the ejection events are due to non  radial pulsations or large scale convection and activity on the star , the timescales will be comparable . therefore the short term mass loss rate associated with the ejection of the nw arc was at least 10@xmath53 m@xmath51 per year . the corresponding kinetic energy of the nw arc is @xmath41 10@xmath54 ergs which is not unrealistic for a star like vy cma with a total energy output of 6@xmath55 ergs per year .    given the evidence we have for a relatively stationary nebula , most of vy cma s mass loss may be occurring in these separate high mass loss events perhaps due to activity on its surface . if we associate the ejections and surface activity with vy cma s photometric variability then they may occur approximately every 50 years or so . the nw arc may be an especially massive outflow , but if we assume @xmath0 10@xmath53 m@xmath51 per ejection , then over a thousand years the mass loss rate would average @xmath56 m@xmath51 per year which is close to what we estimated in smith et al ( 2001 ) for the entire circumstellar nebula and to previous independent mass loss measurements @xcite . the random orientations of the three primary arcs relative to the star together with their different ages , suggest that they were ejected not only at different times , but by localized events or disturbances at different positions on the star . our results show that the nw arc is moving primarily across our line of sight with little projection with respect to the plane of the sky . vy cma s geometry is unknown , but if we assume the ne - sw rotation / symmetry axis from the maser model , then the nw arc is very likely within 10@xmath10 to 15@xmath10 of the expected equatorial plane . given the relative positions of arcs 1 and 2 within the nebula and with respect to the presumed axis , each has likely been ejected from different locations on the star , at higher latitudes relative to the nw arc , and possibly from the polar region , but that depends on the unknown projection angle . one of the uncertainties in the above discussion is our lack of information on the transverse motions of the ejecta and the geometry of the nebula . we will soon be obtaining second epoch images and imaging polarimetry of vy cma with hst to measure its expansion and map the morphology and relative distribution of the embedded structures . our observations confirm that vy cma has experienced several mass loss events , presumably from localized active regions . these clouds or loops of gas and dust are now expanding outward into what appears to be a surrounding nebula of surprisingly stationary gas and dust . this stationary or slow moving material material is difficult to explain around a star that has had a long history of high mass loss . we can only speculate that the more diffuse , uniformly distributed gas has been created by a continuous , but slow general outflow of material from the overall surface . it therefore seems likely that the high mass loss for rate for vy cma is really a measure of the mass carried out by these specific ejections perhaps accompanied by streams or flows of emitting gas through holes or gaps in the dust envelope . vy cma is the only m  type supergiant known to have a large , visible reflection nebula . schuster & humphreys ( 2005 ) have obtained hst / wfpc2 images of several cool hypergiant stars including the oh / ir m supergiants nml cyg , s per and vx sgr to search for associated nebulosity . the results are basically negative or yield only marginal detections , except for nml cyg which has a small bean - shaped nebula that has probably been shaped by the strong ionizing winds from the nearby hot stars in the cygnus x complex @xcite . in many ways vy cma is most like irc + 10420 , another high mass  loss , strong oh / ir source believed to be in a post  red supergiant stage and with a complex circumstellar environment @xcite . its extended reflection nebula has a variety of structures including knots and condensations , ray - like features and numerous arcs and scalloped loops which in the outer parts closely resemble the arcs observed in vy cma s ejecta . spectroscopy of its ejecta though showed a uniformly expanding outflow of gas @xcite ; there was no evidence for stationary gas or dust in its ejecta . the referee has suggested that the slow moving gas is due to gravitationally bound material in circular orbits ( jura and kahane 1999 ) . this may be the case but the examples of orbiting material about evolved stars come from a very different class of objects ; carbon - rich red giants and post agb stars which are binaries or presumed binaries ( jura , bahm , and kahane 1995 , bujarrabel et al . 2003 ) with `` reservoirs '' of long  lasting molecular material . at the distances where the motionless k i emission is observed ( 4000 - 10000 au from the star ) , the orbital period will be several hundred thousand to a million years , exceeding the expected duration of the red supergiant stage for a star this luminous . furthermore , the oh and h@xmath8o maser observations show that the clouds are expanding at 30  40 km s@xmath1 , greater than the expected escape velocity from the dust shell . therefore we think this explanation is less likely for vy cma if it is an evolved red supergiant . vy cma is a very perplexing star experiencing discrete mass loss episodes that may be the primary mechanism by which it is losing so much mass . we have also uncovered some puzzling questions about the origin of its extended diffuse nebula with the slow moving gas and dust . if it is an evolved massive star , it is apparently in a unique stage , with no obvious counterparts . alternatively , the apparently `` stationary '' ejecta with embedded streams of gas and massive outflows is also consistent with a pre - main sequence state for vy cma . we are especially grateful to george herbig for his continuing interest in vy cma , his comments and advice on this study , and for obtaining the telescope time for these observations . it is also a pleasure to thank michael koppelman for assistance with the data reduction . this work was supported in part by the university of minnesota . the k  i emission line is difficult to assess , because this atom has an extremely small ionization potential ( 4.34 ev ) and because some critical parameters are inaccessible for vy cma . since k  i @xmath577665,7699 is a low - excitation resonance doublet , one is tempted to think `` it s merely resonance scattering ''  which , however , does not answer the most important questions . note , for instance , the strength of @xmath27699 emission relative to the continuum when we look directly toward the star , position i ap 2 in fig .  3 . in order to achieve the observed equivalent width of @xmath41 5  in emission simply by resonance scattering , a substantial amount of radiation must be removed from the continuum along some ( indeed most ) lines of sight . on average , _ resonance scattering in itself neither augments nor diminishes the net emergent flux of @xmath27699 photons , _ excepting the small fraction that are scattered back into the photosphere . however , there is no hint of sufficient `` absorption '' to balance the observed `` emission '' anywhere along our slit positions , including reflected spectra that represent the star as seen from various directions in space . other puzzles involve the apparent velocity and the line profiles seen in the extended halo ,  4.2 above .    as a preliminary detail , note the optical thickness for k  i @xmath27699 . ( @xmath27665 has a larger oscillator strength but is not included in our data . ) assuming a normal potassium abundance @xmath58 and a local line width of 4 km s@xmath1 , the line - center opacity is @xmath59 @xmath27 g@xmath1 . at least in the dusty zones , compare this to the continuum opacity for scattering by dust grains , @xmath60 @xmath27 g@xmath1 at the same wavelength . if k  i resonance scattering is to be at least comparable to reflection by dust , we require @xmath61 , i.e. , @xmath62 . given potassium s low ionization potential , this is a stronger constraint than one might guess at first sight .      can k  i emission arise within @xmath63 au , before dust has formed ? first consider the @xmath64 ionization ratio , which , in a simple model , depends on the star s photon supply at energies above 4.3 ev and on the electron density in the wind . suppose that free electrons are provided by the elements with low ionization potentials , giving @xmath65 . at @xmath20 au , for example , in a simple model we expect @xmath66 @xmath22 . if the continuum spectrum of vy cma were a 3000 k planckian , then it would produce about 10@xmath67 potassium - ionizing photons per second ; but this is embarrassingly sensitive to temperature and other parameters , and the star s uv brightness may considerably exceed a simple black - body value ( compare @xmath11 ori , glassgold & huggins 1986 ) . anyway , using reasonable estimates of the photoionization cross - section and recombination coefficient , the values quoted above imply @xmath68 which applies throughout the inner wind , 12 au @xmath69 @xmath70 @xmath69 100 au . this ratio would be decreased by any additional supply of free electrons , e.g. , by some process that can ionize a tiny fraction of the hydrogen and/or an appreciable fraction of carbon ; but it is difficult to reduce the k@xmath71/k@xmath72 ratio by a factor much smaller than 0.2 in this way . most other likely complications tend to _ increase _ the ionization ratio , see glassgold & huggins ( 1986 ) . altogether , at this level of reasoning , one expects the line - center opacity in the inner wind to be rather small for a resonance line , @xmath73 @xmath27 g@xmath1 . the total line - center optical depth along a radial path may exceed 10 , but if so the effect on the total emergent spectrum is merely a fractional -  wavelength redistribution of photons near 7699   i.e. , quite different from the observed strong emission feature which is not accompanied by a comparable amount of absorption at nearby wavelengths .    if there are large - scale inhomogeneities in the local scattering coefficient @xmath74 , the emergent spectrum can depend on viewing direction ; a net surplus of 7699  photons may escape in some directions at the cost of a net deficit elsewhere . in order to obtain the observed strong feature in this way , however , one must view the star from a statistically unusual direction . a fundamentally non - spherical model in this vein  for instance , imagine a circumstellar torus seen from a near - axis viewpoint  is difficult to construct for the following reason : in order to obtain the large observed ratio ( @xmath27699 peak height)/(continuum level ) @xmath41 5 , the solid angle of `` resonance photon escape '' directions must be far smaller than that of the `` opaque '' directions . ( admittedly this may be the case south of the star , along slit v. ) in summary , this general type of scenario can not be rejected outright but we hope to find alternatives that are quantitatively more plausible .    can the observed k  i emission be excited by collisional processes in the inner wind ? in such a model the required luminosity in the @xmath27699 line is of the order of 150 @xmath75 or @xmath76 erg s@xmath1 . electron collisions are rather ineffective because @xmath77 is insufficient . if the inner - wind temperature is , say , 4000 k , and we adopt the relative electron density mentioned above , then collisional excitation of k  i @xmath137665,7699 proves too weak by a factor of @xmath78 , i.e. , by 5 to 9 orders of magnitude . radiative recombination is too weak by about 8 orders of magnitude . unusual effects such as collisional excitation by fast non - thermal electrons , collisional excitation by hydrogen atoms , etc . , are likewise unpromising , though we omit details here . should we consider fluorescence , radiative excitation to higher levels of k  i ? the upper level for @xmath137665,7699 is 4p @xmath79p@xmath80 , and the next resonance level , 5p @xmath79p@xmath80 , can be excited from the ground level by photons near 4044  . an ensuing cascade , 5p @xmath81 ( 3d or 5s ) @xmath81 4p , would lead to @xmath27699 emission . since the photosphere of vy cma is relatively faint at violet wavelengths , this process seems worthwhile only if it is pumped by some other strong emission line near 4044  . so far as we know , there is no likely candidate for this role , nor for pumping k  i levels above 5p .    in summary , we can almost eliminate all the obvious excitation processes other than resonance scattering , and even the latter does not seem promising in the inner wind . therefore let us focus on resonance scattering farther out , @xmath82 au , where inhomogeneous dust may help to solve the puzzle . once dust has begun to form near @xmath83 au in the outward mass flow , the k@xmath71/k@xmath72 ionization ratio becomes even more difficult to predict . the uv flux is greatly reduced by even a small column density of dust , while , on the other hand , some ion species that help provide @xmath77 become depleted . meanwhile the recombination time is several years , long enough for the gas to move a distance comparable to @xmath28 . since the data indicate that the k  i column density is indeed substantial in dust farther out , let us _ assume _ that @xmath84 @xmath27 g@xmath1 , i.e. , at least comparable to scattering by dust , in the inner parts of the dusty region . this assumption is questionable but it is not implausible . figure 14 shows , in an idealized way , how resonance scattering can result in a net emission line if large - scale dusty condensations are present . suppose , for the sake of illustration , that  + @xmath12 each localized dense blob in the figure contains enough dust to be practically opaque at all relevant wavelengths ; + @xmath12 most photons from the star encounter such an obstacle ; + @xmath12 along some escape paths  not necessarily straight  the optical thickness of the lower - density medium is less than 3 ; + @xmath12 the parameters allow about one or two percent of the red - wavelength continuum radiation to escape , as observed for vy cma . often the escape process involves an encounter with one of the dense blobs , scattering by a dust grain , and perhaps further scattering in the less - dense medium . the average dust - grain albedo is obviously relevant .    in fig . 14 , a 7699  photon from the star encounters a dusty condensation , is resonance - scattered by a potassium atom , and then escapes after one more scattering event which may be due to either a k  i atom or a dust grain . this type of behavior , averaged over all photons , leads to _ an apparent emission line in the emergent spectrum . _ there are at least two equivalent ways to see why : + ( 1 ) temporarily imagine , for simplicity , that the dust absorbs but does not reflect radiation . in that case , if the photon in fig . 14 had not been scattered by k  i , then it would have been absorbed by dust in the dense condensation .. this answers a question we asked above , namely , why does the observed spectrum not show k  i @xmath27699 absorption lines to balance the apparent emission ? the answer is that most @xmath27699 scattering events occur in locations where the continuum radiation would not otherwise escape anyway , being absorbed by dust if it is nt scattered by k  i. + ( 2 ) of course , real dust grains reflect radiation as well as absorbing it . in effect , resonance scattering by k  i enhances the overall albedo at a wavelength of 7699  . in other words , resonance scattering provides an extra way for 7699  photons to escape from a dusty condensation . once they escape into a less - dense channel , they have an appreciable probability of escaping from the configuration . ( moreover , k  i scattering may be weak in the low - density regions , so it does not necessarily impede the final escape process . ) the main point is that the k  i emission line may be `` formed '' by an indirect process at radii beyond 100 au or so . our observations indicate that expansion of the dusty medium does not red - shift the emergent k  i emission line as much as it does the stellar absorption spectrum ( see  3 ) ; but a quantitative model is necessary to judge whether this is to be expected . unfortunately , a model of this type has many adjustable parameters and probably requires monte - carlo - style radiative transfer calculations , far beyond the scope of this paper . similar reasoning can probably be applied to the ca  i ] emission . however , as we noted in  4 , the existence of ca  i ] in both emission and absorption seems to require an oscillator strength far larger than one would expect for a semi - forbidden transition . if this statement is wrong due to some geometrical or radiative - transfer fallacy , we have not been able to identify the error . bernat , a. p. & lambert , d. l. 1976 , , 210 , 395 bowers , p. f. , johnston , k. j. , & spencer , j. h. 1983 , , 274 , 733 bujarrabel , v. , neri , r. , alcolea , j. & kahane , c. 2003 , , 409 , 573 danchi , w.c . , bester , m. , degiacomi , c.g . , greenhill , l.j . , & townes , c.h . 1994 , , 107 , 1469 de jager , c. 1998 , , 8 , 145 dupree , a.k . , lobel , a. & gilliland , r.l . 1999 , , 194 , 6605 efstathiou , a. & rowan - robinson , m. 1990 , , 245 , 275 . gilliland , r.l . & dupree , a.k . 1996 , , 463 , l29 glassgold , a. e. & huggins , p. j. 1986 , , 306 , 605 guilain , g. & mauron , n. 1996 , , 314 , 585 herbig , g.h . 1970a , mem . soc . liege , 19 , 13 herbig , g.h . 1970b , , 162 , 557 herbig , g.h . 1972 , , 172 , 375 herbig , g.h . 1974 , , 188 , 533 humphreys , r. m. 1970 , , 82 , 1158 humphreys , r. m. 1975 , , 87 , 433 humphreys , r.m . & davidson , k. 1994 , , 106 , 1025 humphreys , r.m . , et al . 1997 , , 114 , 2778 humphreys , r.m . , davidson , k. , & smith , n. 2002 , , 124 , 1026 hyland , a.r . , becklin , e.e . , neugebauer , g. & wallerstein , g. 1969 , , 15 8 , 619 jones , t. j. et al . 1993 , , 411 , 323 jura , m. & kahane , c. 1999 , , 521 , 302 jura , m. , balm , s. p. & kahane , c. 1995 , , 453 , 721 kluckers , v.a . , edmunds , m.g . , morris , r.h . & wooder , n. 1997 , , 284 , 711 lada , c.j . & reid , m.j . 1978 , , 219 , 95 lobel , a. & dupree , a. k. 2000 , , 545 , 454 lobel , a. & dupree , a. k. 2001 , , 558 , 815 marvel , k.b . 1997 , , 109 , 1286 monnier , j.d . , tuthill , p.g . , lopez , b. , cruzalebes , p. , danchi , w.c . , & haniff , c.a . 1999 , , 512 , 351 morris , m. & bowers , p. f. , 1980 , , 85 , 724 morris , m. & jura , m. 1983 , , 267 , 179 nercessian , e. , guilloteau , s. , omont , a. & benayoun , j. j. 1989 , , 210 , 225 richards , a. m. s. , yates , j. a. , & cohen , r. j. 1998 , , 299 , 319 robinson , l. j. 1971 , ibvs , 599 schaller , g. , schaerer , d. , meynet , g. & maeder , a. 1992 , , 96 , 269 schaerer , d. , meynet , g. , maeder , a. & schaller , g. 1993 , , 98 , 523 schuster , m. t. & humphreys , r. m. 2005 , in preparation smith , n. 2004 , , smith , n. , humphreys , r. m. , davidson , k. , gehrz , r. d. , and schuster , m. t. & krautter , j. 2001 , , 121 , 1111 uitenbroek , h. , dupree , a.k . , & gilliland , r.l . , 1998 , , 116 , 2501 wallerstein , g. 1958 , , 70 , 479 wallerstein , g. 1971 , , 169 , 195 wallerstein , g. 1977 , , 211 , 170 wallerstein , g. 1978 , observatory , 98 , 224 wallerstein , g. 1986 , , 164 , 101 wallerstein , g. & gonzalez , g. 2001 , , 113 , 954 wittkowski , m. , langer , n. & weigelt , g. 1998 , . , 340 , l39      lcccc i & 1 & 1.40 & 1.40 & 122 + & 2 & n.a . & 0.00 & 0 + & 3 & 1.75 & 1.75 & -58 + & 4 & 3.50 & 3.50 & -58 + & 5 & 5.00 & 5.00 & -58 + & 6 & 6.50 & 6.50 & -58 + & 7 & 8.00 & 8.00 & -58 + & 8 & 9.50 & 9.50 & -58 + ii & 1 & -2.80 & 3.23 & 145 + & 2 & -1.40 & 2.02 & 161 + & 3 & 0.00 & 1.28 & -156 + & 4 & 1.40 & 1.77 & -104 + & 5 & 2.80 & 2.93 & -84 + & 6 & 4.20 & 4.23 & -76 + & 7 & 5.60 & 5.58 & -71 + & 8 & 7.00 & 6.95 & -69 + & 9 & 8.40 & 8.34 & -67 + iii & 1 & -0.60 & 1.21 & 178 + & 2 & 0.60 & 1.12 & -121 + & 3 & 1.60 & 1.81 & -91 + & 4 & 3.00 & 3.13 & -76 + & 5 & 4.50 & 4.48 & -71 + & 6 & 5.90 & 5.58 & -68 + & 7 & 7.30 & 7.23 & -66 + & 8 & 8.70 & 8.62 & -65 + v & 1 & -2.10 & 2.28 & 54 + & 2 & -0.70 & 0.28 & 18 + & 3 & 0.70 & 0.81 & 177 + & 4 & 2.10 & 2.07 & -155 + & 5 & 3.50 & 3.44 & -149 + & 6 & 4.90 & 4.82 & -147 + & 7 & 6.30 & 6.22 & -145 + & 8 & 7.70 & 7.61 & -144 + & 9 & 8.10 & 9.01 & -143 +      lcccccc 1 & 31.6 , 40.2 & 6.7 & 16.8 , 29.3 & & 25.4 , 60.4 & 6.3 + 2 & 41.0 & 7.8 & 25.0 & 6.3 & 27.3 , 44.5 & 8.2 + 3 & 23.4 , 45.6 , 52.6 & 10.2 & 25.8 , 42.5 & 6.3 & 26.9 , 41.7 & 8.2 + 4 & 30.8 , 56.5 & 10.2 & 42.1 & 8.2 & 34.3 , 60.4 & 13.3 + 5 & 33.9 & & 33.9 , 59.7 & 12.9 & 35.0 , 70.0 & 12.0 + 6 & 35.1 & & 34.3 , 67.8 & 12.5 & 34.0 , 41.0 & + 7 & 38.6 & & 33.5 , 41.0 & & 36.0 & + 8 & 35.1 & & 34.7 , 41.0 & & 32 , 39 & + 9 & & & 32.0 & & & + 1 & 41.5 & & 42.9 & 24.6 & 13.2 , 40.1 & + 2 & 41.9 & & 41.9 & 31.0 & 39.7 & + 3 & 41.0 & & 41.0 & & 41.0 & + 4 & & & 40.6 & & 38.8 & + 5 & & & 41.5 & & 39 & + 6 & & & & & 41 & + 7 & & & & & 42 & + 1 & -14 & & 34.0 & & -8.1 & + 2 & & & -7.6 & & -5.3 & + 3 & 10 & & -4.4 & & & + 4 & & & 8.4 & & 34.0 & + 5 & & & 35.3 & & 36 & + 6 & 37 & & 34.0 & & 36 & + 7 & 37 & & 37.2 & & 36 & + 8 & 37 & & 38.5 & & 36 & + 9 & & & 37.6 & & & +    lcccccccc 1 & 63.2 & 9 & 63.2 & 9 & 69.8 & 8 & 62.6 & 9 + 2 & 59.9 & 9 & 69.0 & 8 & 75.9 & 8 & 59.3 & 9 + 3 & 76.7 & 8 & 73.6 & 8 & 107.9 & 5 & 64.7 & 9 + 4 & 100.4 & 8 & 101.8 & 4 & 116.0 & 8 & 64.4 & 8 + 5 & 76.7 & 5 & 113.6 & 8 & 102.3 & 7 & 64.6 & 8 + 6 & 68.1 & 3 & 107.4 & 6 & 82.4 & 8 & 64.8 & 8 + 7 & 58.5 & 2 & 85.9 & 7 & 78.0 & 6 & 62.7 & 8 + 8 & 56.5 & 1 & 78.9 & 6 & 81.3 & 3 & 64.2 & 4 + 9 & & & 83.9 & 3 & & & 59.5 & 3 + 1 & & & 63.4 & & 57.5 & & & + 2 & & & 58.8 & & 85.3 & & 64.8 & + 3 & 89.9 & & 84.4 & & 93.1 & & 57.9 & + 4 & 103.1 & & 92.6 & & 92.6 & & 59.7 & + 5 & 80.8 & & 91.2 & & 99.0 & & 62.9 & + 6 & & & 89.0 & & 81.0 & & 68.9 & + 7 & & & 79.8 & & & & 64.8 & + 8 & & & 76.2 & & & & & + 1 & & & & & & & 92.0 + 2 & & & & & 88.1 & & + 3 & & & 89.7 & & 86.5 & & + 4 & 86.9 & & 89.7 & & 84.6 & & + 5 & & & 84.6 & & 87.0 & & + 6 & & & 84.6 & & & & + 1 & 48 , 74 & & 76 & & 73.7 & & 92.0 & + 2 & 48 & & 73 & & 86.5 & & 71.9 & + 3 & 114 & & 76 & & 121.7 & & 73.7 & + 4 & 132 & & 118 & & 118.1 & & 68.2 & + 5 & & & 117 & & 126 & & 79.7 & + 6 & & & 119 & & 109 & & 77.8 & + 7 & & & 111 & & & & 68.7 & +    lccccc 1 & 43.3 & 4.3 & 41.9 & 20.0 & -42.8 , -2.1 , 31.7 + 2 & 37.1 & 6.3 & 41.9 & 20.5 & -3.5 , 30.3 + 3 & 24.1 , 60.4 & 6.7 & 40.6 & & -3.5 + 4 & 28.1 , 52.6 , 67.1 & 7.8 & & & -14.5 + 5 & ( 8.6),21.9 , 65.5 & & 41.5 & & + 6 & 8.6 , 42.5 , 56.2 & & 40.6 & & 35.3 + 7 & 3.9 & & 38.3 & & 34.9 + 8 & ( 3.9 ) , 39.0 & & & & 35.8 + 9 & 40.2 & & & & 37.6 +
Abstract:

Answer: high spatial and spectral resolution spectroscopy of the oh / ir supergiant vy cma and its circumstellar ejecta reveals evidence for high mass loss events from localized regions on the star occurring over the past 1000 years . the reflected absorption lines and the extremely strong k i emission lines show a complex pattern of velocities in the ejecta . we show that the large , dusty nw arc , expanding at @xmath0 50 km s@xmath1 with respect to the embedded star , is kinematically distinct from the surrounding nebulosity and was ejected about 400 years ago . other large , more filamentary loops were probably expelled as much as 800 to 1000 years ago while knots and small arcs close to the star resulted from more recent events 100 to 200 years ago . _ the more diffuse , uniformly distributed gas and dust is surprisingly stationary with little or no velocity relative to the star . _ this is not what we would expect for the circumstellar material from an evolved red supergiant with a long history of mass loss . we therefore suggest that the high mass loss rate for vy cma is a measure of the mass carried out by these specific ejections accompanied by streams or flows of gas through low density regions in the dust envelope . vy cma may thus be our most extreme example of stellar activity , but our results also bring into question the evolutionary state of this famous star . in a separate appendix , we discuss the origin of the very strong k i and other rare emission lines in its spectrum .