accurate distances to planetary nebulae ( pne ) are critical for calculating the size , mass , luminosity , age , and other properties of pne and their central stars .
distances to pne are also crucial when studying properties of the galaxy such as the galactic rotation , and scale height of their progenitor stars . in addition , distances are required to ascertain the planetary nebula luminosity function in order to use pne as standard candles for cosmological distance studies ( ciardullo _ et al . _
1988 ) .
despite being bright and containing ( dominated by ) a rich line spectrum , pne have distances that are remarkably ill - constrained .
in fact , distances to many extragalactic pne are better known than distances to most galactic pne .
many different distance measurement techniques have been applied to pne in the past .
most of have been statistical in nature , presuming that some property is common to the entire population ( hajian & terzian 1996 ) . the most commonly used pn distance estimation technique is known as the shklovsky method ( shklovsky 1956 ) , and assumes a constant ionized gas mass for all pne .
unfortunately , these statistical techniques yield errors that are large , often a factor of two or more , as reviewed by terzian ( 1997 ) . in the case of the shklovsky method , @xmath1 , where @xmath2 is the nebular distance and @xmath3 the ionized mass . since pne progenitor masses can span a factor of @xmath410 , and assuming that @xmath3 scales with the progenitor mass , distance errors of several hundred percent are typical .    the most reliable distance measurements for pne so far have been provided by measuring their expansion parallaxes ( terzian 1997 ) .
pn shell expansion velocities are typically @xmath5 km s@xmath6 , or @xmath7 mas yr@xmath6 at distances of @xmath8 kpc ( mas = milliarcsecond ) .
expansion parallax distance determinations require two epochs of image data separated by a time baseline @xmath9 yr , which can be used to measure the angular expansion `` parallax '' , @xmath10 , of a fiducial feature .
the physical expansion velocity of the features must be spectroscopically obtained , and converted to a tangential velocity , @xmath11 ( _ i.e. _ , the component of the velocity normal to the line of sight ) using a spatiokinematic model of the nebula .
the distance to a pn can then be computed directly by dividing @xmath11 by @xmath10 ( hajian , terzian , & bignell 1993 ) .
distances to nine pne have been previously measured using radio images from the very large array ( vla ) at different epochs ( masson 1989a ; masson 1989b ; gomez _ et al . _  1993
; hajian , terzian , & bignell 1993 ; hajian , terzian , & bignell 1995 ; hajian & terzian 1996 ; haryadi & seaquist 1998 ) .
one past attempt has been made to measure the expansion of ngc  6543 using the vla .
however , this expansion parallax measurement was unsuccessful due to insufficient image quality in the radio - wavelength data ( hajian & terzian 1996 ) .
the many well - defined , optically bright features of ngc  6543 serve as ideal fiducials in measuring its angular expansion .
the ability to measure the expansions of individual features within this extraordinary nebula provides us with a way to sort through the nebula s recent history , to probe the mass ejection behavior of its nucleus , and to directly determine the distance to this pn .
as discussed by hajian , terzian , & bignell ( 1995 ) , the optical expansion of a pn is particularly detectable and measurable with high accuracy ( @xmath1225% ) using the hubble space telescope ( hst ) . in this paper , we report the detection of the systematic expansion of ngc  6543 between 1994 and 1997 from hst archival observations . in order to quantify the magnitude of the angular expansion rate , we present three separate astrometric analyses of these recorded images . these observations , along with a spatiokinematic model of the nebula derived from independent measurements of the doppler velocity , permit three direct distance determinations , all of which are in excellent mutual agreement .
in addition , the expansion pattern reveals that the kinematic ages of features approximately scale with their angular size .
this provides insight into the highly organized and orchestrated mass loss history of the nucleus , presently a wr - of star ( heap & augensen 1987 ) .
this two - in - one paper discusses both the distance to and the nebular evolution of ngc  6543 .
the bright core of ngc  6543 just fills the higher - resolution planetary camera ( `` pc '' ) of the wide field and planetary camera 2 ( `` wfpc2 '' ) whose projected pixel size is 45.5 mas .
archived unsaturated hst narrowband images from go program 5403 and corresponding calibration images from go program 6943 were used for this study ( table  1 ) .
the data epochs are separated by 35 months .
images through f502n ( hereafter `` [ o  iii ] '' ) , f656n ( `` h@xmath13 '' ) , and f673n ( `` [ s  ii ] '' ) filters were obtained at both epochs .
the use of the same instrumentation for the two epochs , albeit at differing orientations , reduces systematic errors .
measurements of the expansion were based solely on images in the [ o  iii]@xmath14502 nm and [ s  ii]@xmath14673 nm lines . the [ o  iii ] image is optimal for this study since it is characterized by well - defined , bright structures with excellent signal - to - noise ratios in the wfpc2 images .
the h@xmath13 images appear essentially redundant to their [ o  iii ] counterparts .
for this reason the h@xmath13 images were not analyzed .
the [ s  ii ] images were used to confirm the expansion measurements from the [ o  iii ] images , and to investigate the proper motions of specific , bright , low - ionization features .
the pointing and orientation of the hst were different for the 1994 and 1997 observations .
therefore , after cosmic rays were removed and the relative image intensities were calibrated , the images had to be regridded to a common center and orientation .
we used the iraf tasks `` drizzle '' , `` imshift '' , and `` rotate '' to remove geometric distortions from , align , and rotate the images .
the shifts and rotations were dithered until patterns of residuals of the corresponding images from the two epochs showed no trace of alignment or orientation errors .
we estimate the accuracy of the alignment to be @xmath15 0.05 pixel ( 2 mas ) and the orientation to be @xmath16 0.1@xmath17 .
blinking of the 1994 and 1997 images shows clearly that the nebula has changed its scale size ( _ i.e. _ the pattern changes are radial ) .
our strategy is therefore to average shifts in the locations of fiducials on equal and opposite sides of the nucleus ( along radial lines ) using various methods .
the angular expansion rate , @xmath10 , is best measured using fiducials with narrow peaks and/or sharp edges . in our analysis , we concentrated on the edges of several nebular features shown in figure 1 . using the nomenclature of miranda & solf ( 1992 ; hereafter `` ms92 '' ) , we analyzed the `` p.a .
25@xmath17 ellipse '' ( e25 ) and the `` p.a
. 105@xmath17 ellipse '' ( e105 ) . in the outer regions ,
we were able to measure the motion of the bright polar caps and the faint polar condensations , referred to as ( d - d@xmath18 ) and ( f - f@xmath18 ) , respectively , in ms92 .
these structures are identified in figure 1 .
we used three methods to extract @xmath10 from the images .
the most direct approach , called the `` profile method '' , is to compare one - dimensional radial brightness distributions , or `` profiles '' , from images of two epochs taken with the same filter .
the relative positions of sharp nebular features are found by analyzing the flux gradients in the radial profiles in the residual image of the two epochs , as pioneered by masson ( 1986 ) .
this method allows @xmath10 to be estimated along many one - dimensional `` cuts '' through the nebula .
the random errors appear to be dominated by photon statistics .
a summary of the `` profile method '' can be found in hajian , terzian , & bignell ( 1993 , 1995 ) and in hajian & terzian ( 1996 ) , referred therein as the `` expansion parallax algorithm '' .
the analysis for hst images is identical to the analysis for vla data except that hst data are not obtained in the fourier plane .
the method requires high signal - to - noise ratios in the residuals after the images of two epochs are subtracted , and it assumes no change in surface brightness or nebular shape . using this method ,
@xmath10 can be computed for any bright feature intersecting a radial cut : @xmath19 where @xmath20 is the flux gradient of the nebular feature and @xmath21 is the peak value of the difference map near the nebular feature .
we take the uncertainty in expansion measurement along an individual profile to be equal to ( 1 ) the difference between expansion rates inferred from the corresponding features on opposite sides of the nucleus , or ( 2 ) the unweighted mean and standard deviation of multiple cuts intersecting an extended structure , where appropriate .
the `` radial fitting '' method is a similar and somewhat redundant technique .
the positions of bright nebular features along 1-d radial flux profiles are measured by performing gaussian fits , by evaluating the flux - weighted centroid of a fiducial s peak , or by graphically measuring the shift of a sharp edge in each of the two epochs .
the radial fitting method has a somewhat increased potential for yielding expansion estimates for nebular features that are too faint for similar measurements with the profile method .
this is due to the increased noise level in a difference map relative to the epoch maps ( @xmath22 for epoch maps with equal noise levels ) .
the random errors of the measurement are dominated by the uncertainties in the fit , which include the judgement involved in defining a fiducial . for fits to peaks ,
the measurements of a gaussian feature in a complicated environment are sensitive to choices of the edges and base of the fitted gaussian curve .
great care is necessary to use the same procedures for the images of both epochs .
finding the relative positions of sharp edges is straightforward .
it is done by overplotting the 1-d flux profiles of the two epochs , and graphically measuring the radial displacement at multiple positions along the steep flux dropoffs of the edges .
we assume a minimum 20% uncertainty in measured expansion rates since we only employ the radial fitting method once for each fiducial in each filter .
the `` magnification method '' , is the simplest but least accurate of the methods .
a magnification factor is applied to the first - epoch image which minimizes the rms of the difference , or residual , image in the regions of the bright fiducial features .    for each emission line , [ o  iii ] and
[ s  ii ] , we used the iraf task `` geotran '' to magnify the 1994 image by various magnification factors , @xmath23 , in the range of 1.001 to 1.005 . for radial motions with no brightness changes ,
a given feature will disappear in the residual image after the magnified 1994 image is subtracted from the aligned 1997 image .
we can then compute @xmath10 using the following equation : @xmath24 where @xmath25 is the angular distance from the nucleus to the feature .
the advantages of the magnification method are that the entire image is utilized , that the intuitive pattern recognition capability and judgement of the brain can be exploited , and that the residuals may uncover deviations from simple radial expansion in the complex nebula that would not easily be recognized in one - dimensional analyses of the brightness distribution .
the major disadvantage to the magnification method is that errors are difficult to quantify .
we assume an uncertainty of 25% in our magnification method expansion rates .
once @xmath10 has been measured using the above methods , and @xmath11 has been determined based on spectroscopic observations and a spatiokinematic model used to convert radial into tangential velocities , the distance to the nebula in parsecs can be calculated : @xmath26 the tangential component of the expansion velocities of various features were taken from the detailed observations and the nebular models of ms92 .
in addition , the kinematic age , @xmath27 , of a fiducial is given by : @xmath28    we also note that a more accurate distance can be computed by taking all of the measured values of @xmath10 into account ( _ i.e. _ , @xmath10 at multiple p.a.s for multiple features , rather than just the minor axis of e25 , which our distance is based upon ) .
however , since ( 1 ) the only available kinematic information for ngc  6543 was obtained in an emission line , [ n  ii]@xmath29 , that is different from the line(s ) used to determine the astrometric expansion of the nebula , [ o  iii]@xmath30 and [ s  ii]@xmath31 ; and ( 2 ) uncertainty in the adopted expansion velocity is approximately @xmath1610% , we did not feel it was necessary to fit a complex spatiokinematic model to the data
. this will be more appropriate when we obtain positionally resolved velocities in the summer of 1999 .
we next discuss the measurements of @xmath10 using all three methods .
the results are compiled in table  2 .
also shown in the table is the expected value of @xmath10 assuming uniform nebular expansion whose rate is calibrated to the average measurement of @xmath10 along the minor axis of e25 ( see below ) .
for the [ o  iii ] images , we applied the profile analysis to various radial cuts through the nucleus in @xmath32 increments , as shown in figure 2 .
the resulting angular expansion rates are displayed in units of mas yr@xmath6 wherever the profile intersects a bright , thin feature .
for the [ s  ii ] images we measured @xmath10 along radial cuts through the central star at intervals of 2 degrees as shown in figure 3 .
the residuals of the magnification method are shown in figures 4 and 5 , and the corresponding values of @xmath10 are shown in table  2 .
the changes in the nebular structure are seen most clearly in the residual images of the unmagnified 1994 - 1997 image pairs .
see the upper left panels of figures 4 and 5 .
although the residuals are only about 1% of the image brightnesses , the patterns of change are quite clear . none of the residuals from any of these techniques yield evidence of nonradial motions .    the best estimates of the distance rely on sharply defined symmetric structures such as spheres or prolate ellipsoids with sharp edges .
we capitalize on a geometric trick : the measurement of @xmath10 along the minor axis of a prolate ellipsoid is ideal for deriving an expansion distance since no geometric correction for inclination is required .
e25 is the optimum structure for expansion studies because of the prominent and sharp edges along its minor axis and its well - measured expansion velocity .
e105 is somewhat more diffuse than e25 , and ms92 s doppler velocity mapping has shown it to be a tilted circular ring rather than an prolate ellipsoid ; however , it serves as a useful consistency check .
_ e25 : _ for the profile method , the best measured expansion rate of e25 s minor axis in the [ o  iii ] images is @xmath33 mas yr@xmath6 . adopting ms92 s expansion velocity of 16.4 km s@xmath6 from the kinematic model described earlier in this section , and a 10% velocity uncertainty results in a distance d@xmath34 pc . based on a size of @xmath25 = 36 for the minor half - axis of e25 ,
the kinematic age is t@xmath35 yr .    the radial fitting method was used on both the [ o  iii ] and the [ s  ii ] images . using the iraf command `` splot ''
, gaussians were fitted to measure the shift of the fiducial features between the two epochs . for
[ s  ii ] , the radial cuts were offset slightly to avoid a bright cosmic ray @xmath16 15 to the south of the central star .
we measured an angular expansion rate @xmath36 mas yr@xmath6 for the minor axis of e25 using this method .
the result agrees to within 10% with the result of the profile method .
the kinematic model described above implies a distance d@xmath37 pc . the corresponding kinematic age t@xmath38 yr .    the magnification method gives very consistent results . the expansion of most of the inner core ( within the intersection of e25 and e105 ) is well characterized by a single magnification factor midway between 1.0025 and 1.003 .
we adopt a value of 1.00275 and an error of 25% ( _ i.e. _ , @xmath39 mas yr@xmath6 ) .
judging from the darkened residuals along e25 , there may have been a slight decrease in surface brightness as well .
the distance derived from the magnification method results is d@xmath40 pc .
we derive a kinematic age t@xmath41 yr .    combining the results for the minor axis of e25 with equal weights and averaging the errors
, we derive a best - estimate distance @xmath42 pc and a kinematic age @xmath43 yr .
this conservative approach of averaging our uncertainties together is warranted because our measurements all rely on the same images and radial velocities , and are therefore not statistically independent .
next we consider the measurement of @xmath10 along the major axis of e25 .
the structure at the ends of the major axis is much more poorly defined than along the minor axis .
accordingly , the various techniques produce somewhat scattered results .
we measure 4.51 , 1.95 and @xmath165 mas yr@xmath6 using the profile , radial , and magnification methods respectively .
the disagreement in the results is substantial .
the important point is that on average all of the methods fall about a factor of two short of the expected result for uniform ( `` hubble law''-like ) expansion , 7.7 mas yr@xmath6 . in other words , the kinematic age of the major axis of e25
is @xmath44 yr .
similarly large expansion ages are obtained not just along the major axis of e25 but throughout the periphery of the core of ngc  6543 by almost all methods .
whether the peripheral structures are truly older than the material along the minor axis of e25 , or whether the outer material is coeval and decelerated is an issue that will be addressed in  4 .    in principle
, the major axis of e25 can be exploited to give a check on the expansion distance since we are able to measure its expansion rate @xmath10 .
after some consideration we decided this would be a mistake for several reasons : ( 1 ) @xmath10 must be corrected for inclination , and the angle is uncertain ; ( 2 ) the structure of e25 along its major axis is irregular and faint ; ( 3 ) our data show that the ellipse e25 does not expand uniformly , as ms92 assumed , and that the major axis deviates the most from the uniform expansion pattern ; and ( 4 ) there is some doubt that the tips of e25 are actually part of the same physical structure as the minor axis ( balick , wilson , and harrington , in progress ) .
e105 has almost the same projected size and shape as e25 , so it is tempting to use it to derive an independent distance to ngc  6543 . however , e105 is lumpier and more diffuse , and measuring its size and expansion rate is commensurably more uncertain . in addition , its expansion velocity is not as well known .
all three methods again yield consistent results on the expansion rate of e105 along its major and minor axes . moreover , those results concur with the expansion rates derived for e25 .
that is , the values of @xmath10 for e105 follow a uniform hubble - law expansion relation based on the angular size and @xmath10 that we measured for the minor axis of e25 .
this is comforting .
this means that e105 and the minor axis of e25 share the same kinematic ages .
let us now employ the ms92 model wherein e105 is a planar circular ring which surrounds e25 and shares the same symmetry axis .
since e105 is a tilted circular ring , the ring s expansion rate measured along its major axis , @xmath45 mas yr@xmath6 , needs no geometric correction for tilt .
( this value is an average of the results of the various methods in table  2 . )
the measured expansion velocity , is 28 km s@xmath6 ( ms92 ) , and the derived distance @xmath46 is @xmath47 pc , where the errors in the expansion velocity ( which may be substantial for this feature ) have been ignored .
satisfyingly , d@xmath48 is consistent with @xmath49 .
_ caps @xmath50 : _ the caps are lumpy complexes which lie within thin regions in the outermost extremities of the nebular core .
the locations of the caps at the extreme edges of the core suggest that the pair of them lie in or near the plane of the sky .
this means that their doppler shifts should be , in principle , unmeasurable , and no distance can be estimated from them . on the other hand ,
this same geometry guarantees a straightforward computation of the kinematic age once the angular expansion rate of the caps has been measured .
if ngc  6543 expands uniformly at the rate determined from the minor axis of e25 then we would expect the caps to exhibit angular expansion rates @xmath51 10 mas yr@xmath6 .    using the profile method
, we detected the signature of angular expansion from the polar caps in both [ o  iii ] and [ s  ii ] . in the case of [ o  iii ] , the expansion rate of the polar caps averages to 6.70@xmath01.37 mas yr@xmath6 . in the case of
[ s  ii ] , the average expansion of the polar caps is @xmath52 mas yr@xmath6 .
the discrepancy from uniform expansion lies well outside of measurement errors .
these results are verified by the radial method , which gives expansion rates of @xmath53 mas yr@xmath6 and @xmath54 mas yr@xmath6 for [ o  iii ] and [ s  ii ] , respectively .
the residuals of the caps are minimized for values of @xmath55 in both the [ s  ii ] and [ o  iii ] images .
our value of @xmath56 corresponds to @xmath10 of @xmath57 mas yr@xmath6 , which agrees very well with the results from the other two methods .    when all three methods are averaged together , the caps have an expansion rate of @xmath58 mas yr@xmath6 , which is more than a third smaller than the expected uniform expansion rate .
they are characterized by a kinematic age of @xmath59 yr .
the caps , like the major axis of e25 have significantly larger kinematic ages than the inner regions of e25 .    _
condensations @xmath60 , & jets @xmath61 : _ if the condensations follow the expansion rate of the minor axis of e25 then we expect to measure a large angular expansion rate , @xmath62 mas yr@xmath6
. yet their proper motions are barely discernible ( table 2 and figure 5 ) .
this is even more peculiar in light of their relatively large doppler shifts , @xmath6325 km s@xmath6 ( ms92 ) which , after ( presumably ) large correction for inclination suggest true space velocities @xmath64 km s@xmath6 relative to the nucleus . it is therefore quite curious that their proper motions are so small .
unfortunately the roll angle of hst excluded the interesting polar jets j - j@xmath18 from the 1997 images .
future hst observations of ngc  6543 using the 1994 roll angle should allow measurement of their proper motion .    _ decreasing expansion rates . _
all of the methods seem to agree that the apparent expansion rates do not increase linearly with increasing radius , _
i.e. _ the core of ngc  6543 does not expand uniformly .
either the outer nebular gas is decelerating , or it was ejected considerably earlier than was the gas within @xmath65 of the nucleus .
we shall return to these points later .      based primarily on e25
we have derived a distance to ngc  6543 of 1001@xmath0269 pc .
this result is consistent with the previous distance determinations of 1100 , 1170 , and 980 pc by cudworth ( 1974 ) , castor _ et al . _  ( 1981 ) , and cahn _ et al . _  ( 1992 ) , respectively , but is discrepant from the 640 pc distance computed by daub ( 1982 ) by about 1.4@xmath66 .
the various methods that we used to measure angular expansion rates all appear to be feasible .
( the glaring exception is @xmath67 measured for the e25 major axis in [ o  iii ] line . )
indeed , based on the experience that we gained , it seems that proper motions of 0.25 pixels , or about 12 mas with the pc camera , can be measured with errors of order 20% when the signal is strong .
this applies not only to images obtained under ideal conditions ( same guide stars , centering , spacecraft orientation ) but under situations in which substantial regridding of the data and corrections for geometric distortions and cosmic rays are required .    according to the papers by hajian , terzian , & bignell ( 1993 , 1995 ) and hajian & terzian ( 1996 )
, the profile method might even achieve accuracies of 2 mas under better conditions ( sharper features , better control of the observations , ideal exposure durations , etc . )
than those in this study .
however , the ultimate technique would exploit the two dimensions of the image data , perhaps based on the optimum reduction of the residuals along various filamentary features or the nebula as a whole .
of course , at some point expansion parallax studies will be limited by systematic errors in the wfpc2 detectors , but apparently this limit has not been reached for motions as small a few mas .
however , the most severe limitation of any pn expansion distance method remains a combination of assumptions about the three - dimensional geometry of the target and the measurement of expansion velocities , especially with spatial resolutions obtainable from the ground in complex targets such as ngc  6543 .
ms92 built a heuristic ( empirical ) model for the evolution of ngc  6543 from the best ( ground - based ) images and their detailed long - slit spectroscopy available in 1991 .
they assumed that each feature in the core of the nebula was axisymmetric and uniformly expanding . from the data and
this assumption they could derive the dimensions and inclination of each feature .    e25 was modeled as a closed , prolate ellipsoid , and e105 is a circular ring surrounding e25 in its equatorial plane . in their model
a second , or `` outer '' thick ellipsoid surrounds the nebular core .
its inner edges are approximately delineated by the caps @xmath2 and @xmath68 on its major axis and the projected tips of e105 along its minor axis .
ms92 found that both e25 and the outer ellipsoid expand uniformly along nearly the same symmetry axis and inclination .    the newer and more detailed hst images require an entirely new geometric interpretation .
our heuristic model for the geometry of ngc  6543 is shown in figure 6 .
we retain the idea of e25 as a prolate ellipsoid , albeit with numerous small bulges and irregularities near its tips .
we are confident that e25 is a prolate ellipsoid for two reasons , ( 1 ) ms92 s spectroscopic data showed with high certainty that it is a closed ellipsoid , and ( 2 ) the hst images are highly suggestive of a prolate ellipsoidal structure .
the outer thick ellipsoid is replaced with a figure - eight outline of two roughly - spherical bubbles with truncated outer regions where the caps are seen .
e105 is situated at the locus where the bubbles are conjoined at the waist , much like a figure eight with a fat waist , or `` fat eight '' ( 8) .
e25 forms one coherent structure , and e105 and the fat eight together form an ensemble which we denote the `` conjoined bubbles '' .    both e25 and
the conjoined bubbles have sharp edges which are nearly unresolved .
if they are @xmath69 yr old or older , and if these sharp edges ( which are strong nebular pressure gradients ) were free to relax ( expand ) at their internal sound speed of 10 km s@xmath6 then they would be @xmath513 x 10@xmath70 cm , or @xmath512 , in diameter .
however , almost all of the features in the core of ngc  6543 have edges which are @xmath71 thick , suggesting that they are constrained by some pressure , perhaps thermal or ram pressure .
what s more , if the edges of the expanding features have been accreting material with low specific momentum then they will have been decelerated .
the further they have traveled the more they will have slowed , nicely in accord with the expansion patterns discussed in  3.2 . in any event , the growth patterns of the core of ngc  6543 provide a very interesting and useful glimpse into the physics of its formation and evolution , as we discuss next .    _
the standard concept of e25 as a wind heated , thermally expanding bubble moving forward supersonically ( @xmath16 mach 1.5 ) is compatible , on the whole , with its morphology , measured doppler shifts , and pattern of angular expansion .
the gas which the bubble has displaced and swept up is now compressed into the bright rim which outlines it .
the leading shell seems to have developed thin - shell instabilities , or lumps , of small amplitude during its evolution .
this happens because the outward flow of the thermally - driven expanding bubble surface encounters dense regions upstream and bends around them , forming surface lumps .
( however , the scale of the lumps is too small to allow studies of their individual expansion patterns . )
the 1994 and 1997 observations show that although e25 is likely to be a prolate ellipsoid , it is not expanding uniformly . indeed , its expansion rate is smallest at its outermost tips .
the expansion of the conjoined bubbles is more difficult to ascertain ; however , it s angular expansion rate is roughly that of e25 .
therefore , the conjoined bubbles are either older or more decelerated than e25 , or both
.    lumpy ellipsoids like e25 are common in elliptical planetary nebulae such as ngc  6543 .
close counterparts are seen in hst images of ngc  5882 , 6826 , 6884 , and especially ngc  7009 ( balick _ et al._1998 ) . in the last case
the bubble has essentially the same morphology as does e25 .
emanating from the tips of the bubble in ngc  7009 are highly collimated gas jets with bright ansae at their termini .
the symmetry axis of ngc  7009 appears to be oriented close to the plane of the sky .
however , if viewed at an inclination of @xmath72 , the bubble , jets , and ansae would resemble e25 and the jets j - j@xmath18 in ngc  6543 .    the elongated shape of e25 and the bubble of ngc  7009 are generally believed to be the consequence of a high - pressure waist of confining material that forces the thermally - expanding bubble to grow most rapidly along the direction orthogonal to the disk ( _ e.g. _ , icke _ et al . _
thick disks or tori produce prolate ellipsoids because the expansion in a large equatorial region is hindered .
thin disks on the other hand , produce pairs of figure - eight lobes as in the homunculus of @xmath73 carinae ( morse _ et al . _
1998 ) and the conjoined bubbles of ngc  6543 , because the restricting material is confined to a narrow , dense region , while allowing free expansion in the rest of the nebula .
the concept of e25 as a thermally expanding bubble fails two tests .
firstly , the torus that putatively constrains the equatorial growth of e25 should be situated astride the minor axis of e25 and probably interior to e105 .
such a torus is readily discerned in ngc  7009 .
however , no such torus is apparent in ngc  6543 , perhaps owing to the confusion of brighter features nearby .
secondly , the prolateness develops because the wind - driven bubbles expand fastest in the polar directions where the upstream confining pressures are smallest . therefore the angular expansion rate of e25 is predicted to be largest along its major axis .
rather , just the reverse is seen .
the kinematic age of the tips of e25 is roughly twice as large than the kinematic age of its waist .
if the observed trends continue , the e25 ellipsoid will evolve to become increasingly less prolate in time .    to rescue the concept of e25 as a thermally expanding bubble
one might argue that although the tips of the bubble originally expanded much faster than its waist , the tips have now displaced and accreted considerable amounts of gas with lower specific momentum and have been decelerated .
alternately the gas behind the tips may have cooled adiabatically , so the thermal pressure driving the expansion has abated .    _ conjoined bubbles . _ almost certainly the conjoined bubbles represent the projected edges of a pair of bubble - like shells conjoined along their waist . along its waist
a ruffled and slightly irregular series of ragged low - ionization , elongated knots ( e105 ) which seem to point away from the nucleus .
these knots may be the result of an instability formed as the edges of the expanding bubbles pinch and compress the trapped gas in their waist , like the blades of a pair of scissors .
alternately they are ablation flows from dense neutral knots subjected to ionizing radiation ( mellema _ et al._1998 ; redman & dyson 1999 ) .
although e25 has counterparts in many other planetaries , the conjoined bubbles do not .
the closest analogues may be ngc  7027 and some bipolars such as ngc  650 - 1 , 2371 - 2 , and 6309 .
interestingly , the conjoined bubbles also mimics @xmath73 carinae , except that the waist of the latter object looks like a large thin disk with radial streaks , and its figure eight is more pinched .
if @xmath73 carinae follows the pattern of decreasing expansion with radius that presently characterizes the conjoined bubbles , then it is not difficult to imagine that its overall morphology will evolve over time to resemble the conjoined bubbles !
dwarkadas & balick ( 1998 ) showed how this might happen if the hydrodynamic evolution of the system is dominated by thermal expansion and adiabatic cooling .
_ historical evolution .
_ one key fact is that the waist of the conjoined bubbles , e105 , and e25 have the same kinematic age within our error bars . hence e25 and the conjoined bubbles formed at roughly the same time .
based on the increasing kinematic age with radial distance from the nucleus , it appears that e25 is overtaking some of the more distant nebular regions such as the tips of e25 and the caps d - d@xmath18 .
we propose a fairly simple framework for understanding the evolution of the core of ngc  6543 . about 1000 yr ago , the star ejected a pulse of material which now forms the conjoined bubbles .
however , the absence of limb brightening along the edges of the conjoined bubbles ( except at @xmath74 ) suggests that the ejection event filled the interior with low - density gas .
for example , this might have happened if the gas formed a shock which left some of the outflowing gas in a hot state .
starting shortly after the pulse of mass ejection the nucleus has been blowing a high - speed wind , perhaps much like the wind observed today : @xmath75 @xmath76 m@xmath77 yr@xmath6 and @xmath78 km s@xmath6 ( perinotto , cerruti - sola & lamers 1989 ) .
this wind has swept into the interior of the expanding conjoined bubbles .
it has displaced and plowed the upstream gas and has formed the feature e25 : therefore , e25 defines the `` realm '' of the wind - affected portion of the core of the nebula .
e25 is now seen as a prolate balloon pushing into the interior of the conjoined bubbles , fastest near the equator , and decelerating towards the poles . as we noted before , the plowing of slower upstream gas and adiabatic cooling inside the tips may have slowed the expansion of e25 along its major axis .
the upstream portions of the conjoined bubbles have yet to know of the existence of e25 .
thus the morphology of the conjoined bubbles is a result of events that shaped it initially .
this feature `` remembers '' its origins .
its pinched waist suggests that a thin , dense disk shaped it into a bipolar when it was formed .
e105 is the remnant of this disk .
this fanciful scenario fails to account for the symmetrically placed caps , condensations , and jets of ngc  6543 .
we can only conjecture that the caps are dense , neutral material left behind by an even earlier ejection of mass . their low specific momentum and neutral , dense gas form impediments to the polar growth of the conjoined bubbles and , at the same time , cause dense ionization fronts to form .
the origins of the condensations and jets are even less obvious .
we shall return to these issues in a later paper in which the detailed ionization structure of the nebula is discussed .
the excellent agreement of each of the three methods presented confirms that we have , for the first time , measured the optical expansion parallax of a pn , and that our error estimates are adequate .
we have shown that accurate measurements of the expansion of nebular features are readily obtainable with the wfpc2 camera with a reasonable time baseline .
we have further shown that the e25 major axis , the polar caps , the polar condensations , and the conjoined bubbles are all expanding significantly slower than expected of a hubble - like expansion . this result is confirmed in both the [ o  iii ] and [ s  ii ] images , and indicates a region in which ejected gas is either decelerating , or left behind from an earlier epoch , or both .
combining the proper motion of bright fiducials of the nebula with corresponding spectroscopically obtained radial expansion velocities ( ms92 ) has allowed us to determine the distance to ngc  6543 of 1001@xmath0269 pc .
our deduced kinematic age of the brightest inner parts of the core of the nebula , based on three measurements along the minor axis of the tilted prolate ellipsoid e25 , is 1039@xmath0259 yr .
the outer regions , which resemble a pair of conjoined bubbles with a large waist , have longer kinematic ages .
a heuristic and highly conceptual explanation of the evolution of these features is as follows .
a major outburst occurred a thousand years before today s nebula formed .
even earlier outbursts produced the huge halo that surrounds ngc  6543 , the polar caps , and other features beyond about a radius of 20@xmath79 from the nucleus .
these debris form the environment for the latest ejection .
the major outburst formed a dense equatorial disk , the remnants of which are e105 .
in addition , the ejected gas was shaped by this disk into a pair of expanding bubbles which join at e105 .
the `` conjoined bubbles '' were created in such a way that their interiors are uniform and probably in a hot state .
these bubbles are expanding and adiabatically cooling .
stellar winds began to blow at high speeds shortly after the latest massive outburst .
these have created an overpressured hot bubble which expends into the interior of the conjoined bubbles .
the bubble attained a prolate geometry by rapid expansion along its polar direction .
the subsequent evolution of the bubble has formed e25 .    like many other interior bubbles in planetary nebulae ,
e25 is forming lumps on its surface , perhaps the result of the buckling of its thin expanding shell .
the expansion of the tips of the prolate bubble has displaced and accreted material of lower specific momentum , so the tips have decelerated much more than its waist .
e25 is therefore becoming less prolate as it expands .    in the future e25
may overtake the conjoined bubbles .
when this occurs rayleigh - taylor instabilities will set in , shredding the remnants of the nebula , and forming a clumpy halo perhaps much like the larger one that presently surrounds ngc  6543 ( _ e.g. , _ balick _ et  al .
_ 1993 ) .    the ongoing go
7501 program has selected 30 pne for wfpc2 observations for its ongoing multi - epoch expansion parallax program .
first epoch observations are nearly complete .
subsequent observations will allow accurate expansion parallax distance measurements of these pne .
third epoch observations of ngc  6543 are planned for late in the year 2000 during which the observing configuration of the 1994 observations by harrington & borkowski ( go 5403 ) will be replicated in the bright , complementary nebular lines of [ o  iii ] and [ n  ii ] as closely as possible .
this observation should provide excellent data spanning six years under highly controlled conditions , allowing even more accurate measurements of the expansion patterns and the distance to ngc  6543 .
we are extremely grateful to virginia player for her diligent efforts in making astrometric measurements for our analyses .
we are most grateful to the space telescope science institute for providing archived data used in this study .
support for this work was provided by nasa through grant number go 7501 from the space telescope science institute , which is operated by aura , inc .
, under nasa contract nas 5 - 26555 .
[ o  iii ] ( f502n ) & 2x200s , 2x600s & 17@xmath80 58@xmath81 28@xmath82 & 66@xmath8337@xmath84 & 274.105@xmath83 & 9 - 18 - 94 & 5403 & harrington ( f502n ) & 200s & @xmath85 58@xmath81 28@xmath86 & 66@xmath833814@xmath87 & 305.434@xmath83 & 8 - 17 - 97 & 6943 & casertano              e25 minor & 72 ( in [ o  iii ] ) & profile & 3.63@xmath00.82 & & & rad .
fit & 3.36@xmath00.86 & & & mag . &
3.4@xmath00.9 & & & & & e25 major & 161 ( in [ o  iii ] ) & profile & 4.51@xmath00.85 & & & rad .
fit & 1.95@xmath01.01 & & & mag . &
4.8@xmath01.2 & & & linear ( hubble ) expectation & 7.7 & & & & & e105 minor & 90 ( in [ s  ii ] ) & profile & 5.8@xmath01.5 & & & rad .
fit & 3.98@xmath01.17 & 3.28@xmath00.7 & & mag .
& 4.6@xmath01.2 & 4.6@xmath01.2 & & linear ( hubble ) expectation & 4.3 & & & & & e105 major & 146 ( in [ s  ii ] ) & profile & 7.2@xmath01.5 & & & rad .
fit & 4.68@xmath02.34 & 8.59@xmath01.7 & & mag . &
7.5@xmath01.9 & 7.5@xmath01.9 & & linear ( hubble ) expectation & 7.0 & & & & & polar caps ( dd@xmath89 ) & 212 ( in [ s  ii ] ) & profile & 6.70@xmath01.37 & 6.94@xmath01.42 & & rad .
fit & 6.56@xmath01.3 & 5.86 @xmath01.2 & & mag . &
6.4@xmath01.6 & 6.4@xmath01.6 & & linear ( hubble ) expectation & 10.2 & & & & & condensations ( ff@xmath89 ) & 247 ( in [ s  ii ] ) & rad .
fit & & 2.73@xmath01.64 & & mag .
& & 4.2@xmath01.1 & & linear ( hubble ) expectation & 11.9 &                                                            * figure 1 : * images of ngc  6543
. left : a linear representation of the [ o  iii ] surface brightness .
various features are identified ( ms92 terminology shown in italics ) .
white dots show the locations of the radial fits of angular expansion rates of table 2 .
center : a negative linear representation of the [ o  iii ] surface brightness superimposed on a logarithmic display of the faint nebular background in the same hst image .
right : like the center panel , but for the [ n  ii ] line .
the [ s  ii ] image is noisier but otherwise similar in appearance to this panel .    * figure 2 : * a linear representation of the [ o  iii ] surface brightness of ngc  6543 showing the location of the profiles analyzed in  3.2 .
values of the angular expansion rate in mas yr@xmath6 are shown corresponding to the locations ( black dots ) where the profile cuts intersect bright regions of the nebula .
* figure 3 : * same as for figure 2 , except for the [ s  ii ] surface brightness of ngc  6543 .
the values listed correspond to the measured angular expansion rates in mas yr@xmath6 where the profiles intersect the bright polar caps .
* figure 4 : * a montage of f502n ( [ o  iii ] ) residual images after the 1997 images were subtracted from magnified and aligned 1994 images .
the magnification factors that were applied to the 1994 images are shown in the upper right corners of each panel .
* figure 6 : * schematic of a conceptual geometric model of ngc  6543 . a plan view of the model is shown on the right side .
how it might appear if tilted at about 30@xmath17 is shown in projection on a combined [ n  ii ] and [ o  iii ] image .


chromospheric evaporation was first suggested by @xcite to explain the origin of the hot , dense , soft x - ray ( sxr ) emitting plasma confined in the coronal loops during solar flares .
the basic scenario is as follows .
magnetic reconnection , believed to be the primary energy release mechanism , heats the plasma and accelerates particles high in the corona .
the released energy is transported downward along the newly reconnected closed flaring loop by non - thermal particles and/or thermal conduction , heating the chromospheric material rapidly ( at a rate faster than the radiative and conductive cooling rates ) up to a temperature of @xmath0@xmath1 k. the resulting overpressure drives a mass flow upward along the loop at a speed of a few hundreds of km s@xmath2 , which fills the flaring loop with a hot plasma giving rise to the gradual evolution of sxr emission .
this process should also result in a derivative of the sxr light curve in its rising portion that closely matches the hard x - ray ( hxr ) light curve , which is called neupert effect and is observed in some ( but not all ) flares @xcite .
hydrodynamic ( hd ) simulations of chromospheric evaporation have been carried out with an assumed energy transport mechanism ( e.g. , electron `` beam '' or conductive heating ) @xcite leading to various predictions on the uv - sxr spectral lines produced by the evaporated plasma , as well as the density and temperature profiles along the flaring loop .
most of the observational tests of these predictions rely on the blue - shifted components of sxr emission lines produced by the up - flowing plasma , first reported by @xcite and @xcite using spectra obtained from the _ p78 - 1 _ spacecraft .
similar observations were subsequently obtained from x - ray spectrometers on the solar maximum mission ( _ smm _ ; antonucci et al .
1982 , 1984 ) , the _ hinotori _ spacecraft @xcite , the _ yohkoh _ spacecraft @xcite , and the solar and heliospheric observatory ( _ soho _ ; brosius 2003 ; brosius & philips 2004 ) .
@xcite , on the other hand , observed co - spatial sxr blue - shifts ( up - flows ) and h@xmath3 red - shifts ( down - flows ) as expected from hd simulations @xcite .
a summary of relevant observations from _ smm _ can be found in @xcite .
all the aforementioned observations , however , were indirect evidence in the sense that the evaporation process was not imaged directly . based on hd simulations , @xcite derived the expected x - ray brightness profile across the evaporation front and suggested that _
yohkoh_/sxt or x - ray imagers with equivalent or better spatial and temporal resolution should be able to detect the front .
indeed , @xcite found that the hxr and sxr sources of the 1994 june 30 flare moved toward the loop top ( lt ) during the impulsive phase . since the flare was located near the center of the solar disk , they identified such motions as the horizontal counterpart of the line - of - sight motion revealed by the blue - shifted emission lines observed simultaneously by _
yohkoh_/bcs .
the reuven ramaty high energy solar spectroscopic imager ( _ rhessi _ ) with its superior spatial , temporal , and spectral resolution @xcite provides us with opportunities to study the chromospheric evaporation process in unprecedented detail .
we report in this paper our analyses of the spatial and spectral evolution of a simple flare on 2003 november 13 with excellent _ rhessi _ coverage . because the flare occurred near the solar limb it presented minimum projection effects and a well - defined loop geometry that allows direct imaging of the hxr brightness profile along the loop .
the observations and data analyses are presented in
[ obs ] , followed by a derivation of the evolution of the density profile along the flaring loop in  [ dens ] .
we summarize the major findings of this paper and draw conclusions in  [ conclusion ] .
the flare under study is a _ goes _
m1.7-class flare that occurred on 2003 november 13 in the active region 0501 after it appeared on the east limb .
this event followed a period of extremely high solar activities in late october and early november when a series of x - class flares including the record setting x28 flare of 2003 november 4 took place @xcite .
_ rhessi _ had an excellent coverage of this flare .
figure [ lightcv ] shows the _ rhessi _ and _ goes_-10 light curves .
the _ goes _ 8 - 1  ( 1.6 - 12.4 kev ) and 4.0 - 0.5  ( 3.1 - 24.8 kev ) fluxes rise gradually and peak at 05:00:51 and 05:00:15 ut , respectively .
the _ rhessi _
high energy ( @xmath4 kev ) count rates , on the other hand , exhibit two pulses peaking at 04:58:46 and 05:00:34 ut with the first one stronger .
the steps in the _ rhessi _ light curves are due to the attenuator ( shutter ) movements @xcite . before 04:57:57 ut and after 05:08:59 ut
, there were no attenuators in and between the two times the thin attenuator was in except for a short period near 05:05 ut when the attenuator briefly moved out .
figure [ mosaic ] shows the evolution of the flare at different energies , which may be divided into three phases : ( 1 ) before 04:57:57 ut , corresponding to the rising phase when the emission mainly comes from a flaring loop to the south .
( 2 ) between 04:57:57 and 05:08:59 ut , the impulsive phase , during which another loop to the north dominates the emission .
this loop appears to share its southern footpoint ( fp ) with the loop to the south , which is barely visible because of its faintness as compared with the northern loop and _ rhessi _ s limited dynamic range of @xmath010 .
( 3 ) after 05:08:59 ut , the decay phase , when the shutters are out and two off - limb sources ( identified as the lts of the two loops ) dominate .
the relatively higher altitudes compared with earlier lt positions are consequences of the preceding magnetic reconnection , as seen in several other _ rhessi _ flares @xcite .
clearly the southern loop , which extends to a relatively higher altitude , evolves slower and is less energetic than the northern one .
we shall focus on the evolution of the northern loop during the first hxr pulse ( 04:58 - 05:00 ut ) in this paper .
it is necessary to check if pulse pileup is important in this flare , before we can make more quantitative interpretation of the data .
the reason is that although we have applied the first order pileup correction @xcite in our spectral analysis , such a correction is challenging for images and not available at present .
there are several ways to do the check , the detector livetime being the first and simplest indicator .
we first accumulated spatially integrated spectra for every 1 s time bins during the interval of 04:58:01 - 04:59:49  ut , ( see text about figure [ prof_3e ] ) , which covers the bulk of the first hxr pulse . ] using the front segments of all the nine detectors except detectors 2 and 7 which have degraded energy resolution @xcite .
we then obtained the livetime ( between data gaps ) from the spectrum object data and averaged it over the seven detectors being used .
the resulting livetime generally decreases with time , ranging from 96% to 89% , with a small modulation produced by the spacecraft spin . in this m1.7 flare ,
such a livetime is comparably high ( cf .
livetime of @xmath5 during the 2002 july 23 x4.8 flare and @xmath6 during the 2002 february 20 c7.5 flare ) and indicates minor pileup severity .
another approach involves inspecting the change of the spectrum due to pileup .
we accumulated spectra over every 1 spacecraft spin period ( @xmath7 s , with the same set of detectors mentioned above ) and used the pileup correction to obtain the relative fraction of the pileup counts among the total counts as a function of energy @xcite .
we find that the pileup counts amount to less than @xmath8 of the total counts at all the energies until 04:59:01  ut when the livetime drops to 91% .
after that , the relative importance of the pileup counts continues increasing , but remains below @xmath9 of the total counts before 04:59:17  ut . toward the end of the first hxr pulse ( 04:59:45 - 04:59:49  ut , livetime @xmath10 ) , the pileup counts to total counts ratio exceeds @xmath11 in the entire 20 - 40  kev range and humps up to @xmath12 near 28 kev .
we respectively integrate the pileup counts and total counts over the 20 - 40  kev band , and plot their ratio versus time , as a general indicator of pileup severity ( see figure [ pileup ] ) .
clearly this ratio is @xmath13 during the first 2/3 of the interval shown and does not reach the moderate @xmath14 level until the very end .
we therefore conclude that pileup effects are generally not very significant for this flare , especially during the first minute of the impulsive phase , because the count rate is not too high and the thin shutter is in at times of interest which further attenuates the count rate .
it should be noted that the two piled up photons ( that result in a single photon seen in the image ) most probably originate from the same location on the sun and pileup of photons across difference sources is relatively unimportant ( g. hurford , private communication ) .
therefore the source geometry would not be significantly affected by pileup , except that there could be a  ghost " of a low - energy source appearing in a high - energy image for very large ( e.g. x - class ) flares . however , the spectra of individual sources derived from images are distorted , which is relatively more significant at the lt than at the fps .
this is because ample low - energy photons dominate in population over high - energy photons , having the highest probability to produce pileup , and generally most of the low - energy photons are emitted by the lt source .
we now examine the images in greater detail .
top left panel _ of figure [ mdi ] shows _ rhessi _ clean @xcite images of the northern loop at 9 - 12 , 12 - 18 , and 28 - 43 kev for 04:58:22 - 04:58:26 ut .
( although the 4-second integration time is rather short , the image quality is reliable with a well - defined source structure . ) at 9 - 12 kev the lt dominates and the emission extends towards the two fps , which dominate the emission at 28 - 43 kev and above with the northern fp ( n - fp ) much brighter than the southern one ( s - fp ) .
one of the most interesting features of the source structure is that emission from the legs of the loop dominates at the intermediate energy ( 12 - 18 kev ) .
similar structures are also observed for several other time intervals during the first hxr pulse ( see discussions below ) .
we find emission from the legs is a transient phenomenon at intermediate energies , because when integrating over a long period and/or a broad energy band , the lt and/or fp sources become dominant . to our knowledge
no images like this have been reported before .
we attribute this in part to the relatively short integration time and _ rhessi _ s high energy resolution .    to be compared with observations at other wavelengths , the same images at 9 - 12 and 28 - 43 kev ( solid contours )
are shown with the _
soho_/eit , mdi magnetogram , and mdi white light maps in the other three panels of figure [ mdi ] , where the dashed contours depict the southern loop at 6 - 9 kev for 04:57:40 - 04:57:52 ut . the eit image at 04:59:01 ut ( _ upper right panel _ ) shows emission at 195  co - spatial with the sxr emission from the northern loop .
the brightest 195  emission , an indicator of the highest differential emission measure ( thus the highest density ) at @xmath15 k , appears to be close to the n - fp which is also the strongest fp in hxrs .
k ( see * ? ? ?
12 ) , and emission intensity would be lower for both higher and lower temperatures . ]
the _ lower left panel _ displays the x - ray emission along with the post - flare ( 05:57 ut ) mdi magnetogram .
this clearly shows that the northern loop straddles across a polarity reversal with the brighter n - fp associated with a stronger magnetic field . the southern loop ( dashed contours )
is associated with a even weaker magnetic field . here
we show the mdi magnetogram recorded one hour after the flare s impulsive phase because during a flare there are many uncertainties in the magnetic field measurement .
the _ lower right panel _ shows the mdi continuum map at 12:47 ut ( about 8 hours after the flare ) , suggesting that the flare occurred above the lower sunspot region ( dark area ) .
note that during this interval the sunspot has moved westward for about @xmath16 in heliographic longitude .
we do not plot the mdi white light map at the time of the flare because then the sunspot was nearly on the limb and barely visible .
next we consider the evolution of the northern loop .
we notice that , as shown in the four columns for 04:58:00 - 04:59:20 ut ( boxed by the dotted line ) in figure [ mosaic ] , the fps initially appear at all energies but later on dominate only in the high energy bands , while the lt is first evident at low energies and becomes more and more prominent at relatively higher energies as indicated by the dashed diagonal line .
the emission from the lt also extends towards the legs at intermediate energies and in a given energy band the emission concentrates more and more at the lt with time .
these are expected to be common features of flares with a single loop because of chromospheric evaporation that can increase the plasma density in the loop making the lt dominate at progressively higher energies .
however , because the 20 second integration time is relatively long , these images do not uncover the details of the evaporation process . to remedy this
we have carried out three different but complementary analyses of the images with higher temporal or energy resolution .
\1 . to study the source morphology change over short time intervals , we model the loop geometry and study the evolution of the hxr brightness profile along the loop .
we first made clean images in two energy bands of 6 - 9 kev ) photons by the attenuator and escape of the germanium k - shell fluorescence photons that are produced by photoelectrical absorption of higher energy ( 10 - 20 kev ) photons in the germanium detector ( see * ? ? ?
. however , for the flare under study , the 6 - 9 kev image most likely reveals the real lt morphology , because there are ample thermal photons at lower energies originating from the lt source and photons at slightly higher energies seem to come from the same location . ] and 50 - 100 kev over the time interval of 04:58:12 - 04:58:53 ut that covers the plateau portion of the first hxr pulse . from these two images we obtained the centroids ( indicated by the white crosses in figure [ loop - model]_a _ ) of the sources identified as the lt ( 6 - 9 kev ) and two fps ( 50 - 100 kev ) , respectively . assuming a semi - circular loop that connects the three centroids , we located the center of the circle marked by the plus sign in figure [ loop - model]_a_. the grey scale in figure [ loop - model]_a _ was obtained by superposition of 30 images ( six 8-second intervals from 04:58:08 to 04:58:56 ut in five energy bands : 9 - 12
, 12 - 15 , 15 - 20 , 20 - 30 , and 30 - 50 kev ) reconstructed with the pixon algorithm @xcite . figures [ loop - model]_b _ and [ loop - model]_c _ respectively show the intensity profiles perpendicular to and along the loop ( averaged over the respective orthogonal directions ) .
the inner and outer circles ( at @xmath17 and @xmath18 ) in figure [ loop - model]_a _ show the positions of the 50% of the maximum intensity in figure [ loop - model]_b_. however , for inferring the intensity profile along the loop we use radially integrated flux down to the 5% level .
this enables us to include as much source flux as possible ( with little contamination from the southern loop ) .
we define the mean of the radii at the 5% level as the radius of the central arc of the loop ( the white dot - dashed line in figure [ loop - model]_a _ ) .    with the above procedure , one can study the evolution of the brightness profile along the loop at different energies .
figure [ prof_3e ] shows the results obtained from pixon images with an integration time of 1 spacecraft spin period ( @xmath04 seconds ) from 04:58:01 to 04:59:49 ut for three energy bands ( 20 - 30 , 15 - 20 , and 12 - 15 kev ) . using a simple algorithm we determine the local maxima whose slopes on both sides exceed some threshold value and mark them by filled circles .
we compare each profile with its counterpart obtained from the clean image ( with the same imaging parameters ) , and use the rms of their difference to estimate the uncertainty as indicated by the error bar near the right - hand end of the corresponding profile . for each panel ,
the rms difference of all the profiles , as a measure of the overall uncertainty , is shown as the error bar in the upper - right corner .
this uncertainty is about 10% for the three energy bands ; as expected , it increases slightly at higher energies which have lower counts .
figure [ prof_3e]_a _ displays the profile at 20 - 30  kev , which as expected ( see figure [ mosaic ] ) shows emission from the two fps with fairly constant positions until the very last stage when the lt emission becomes dominant . at this stage ,
the s - fp becomes undetectable and the n - fp has moved very close to the lt . at the lower energy ( 15 - 20 kev , figure [ prof_3e]_b _ ) the maxima tend to drift toward the lt
gradually and eventually merge into a single lt source . at the even lower energy ( figure [ prof_3e]_c _ )
this trend becomes even more pronounced and the drift starts earlier except here the shift is not monotonic and there seems to be a lot of fluctuations .
we also repeated the same analysis at a higher cadence ( every 1 second , @xmath7 second integration interval ) with both pixon and clean algorithms .
the evolution of the resulting profiles ( although oversampled and thus not independent for neighboring profiles ) appears to be in line with that shown here at a 4-second cadence obtained with pixon .
the general trends of these results indicate that high energy hxr producing electrons lose their energy and emit bremsstrahlung photons higher and higher up in the loop as the flare progresses .
this can come about simply by a gradual increase of the density in the loop , presumably due to evaporation of chromospheric plasma . from the general drift of the maxima
we obtain a time scale ( @xmath010 s of seconds ) and a velocity of a few hundred km s@xmath2 consistent with the sound speed or speed of slow magnetosonic waves .
as stated above at low energies we see some deviations from the general trend , some of which do not appear to be random fluctuations . if so and if we take one of the evident shorter time scale trends shown by the dashed line in figure [ prof_3e]_c _ we obtain a large velocity km s@xmath2 @xcite and 800 km s@xmath2 @xcite , obtained from blue - shifted fe xxv spectra . ]
( @xmath0@xmath19 km s@xmath2 ) which is comparable to the alfvn or fast magnetosonic wave speed .
this may indicate that another outcome of energy deposition by non - thermal particles is the excitation of such modes which then propagate from the fps to the lt and might be responsible for the circularly polarized zebra pattern observed in the radio band @xcite .
this , however , is highly speculative because the spatial resolution ( @xmath0@xmath20 ) is not sufficiently high for us to trust the shorter time scale variation .
the longer time scale general trend , however , is a fairly robust result .
\2 . instead of examining the source structure with high time resolution
, we can investigate it with higher energy resolution at longer integration intervals as tradeoff for good count statistics and image quality .
to this end , we have made pixon images during three consecutive 24 s intervals starting from 04:58:00 in 20 energy bins within the 6 - 100 kev range .
figure [ box ] shows a sample of these images at 04:58:24 - 04:58:48 ut .
figures [ prof_24s]_a_-[prof_24s]_c _ show the x - ray emission profile along the loop at different energies for the three intervals . ) . ] as in figures [ prof_3e]_a_-[prof_3e]_c _ the high energy emission is dominated by the fps but there is a decrease of the separation of the fps with decreasing energies and with time . again at later stages the lt dominates and the profile becomes a single hump .
the general trend again suggests an increase of the gas density in the loop . at lower energies ( @xmath21 kev ) , the profile is more complicated presumably due to many physical processes ( in addition to chromospheric evaporation ) , such as thermal conduction and transport of high energy particles , thermal and non - thermal bremsstrahlung , wave excitation and propagation , wave - particle coupling , and even particle acceleration , which may be involved .
we believe that a unified treatment of acceleration and hd processes with physical conditions close to the flare is required for interpreting these results to uncover the details . to quantify this aspect of the source structure evolution , we divided the loop into two halves as shown by the boxes in figure [ box ] and calculated their emission centroids .
the resulting centroids at the three times together with the central arc of the model loop are plotted in figure [ centroids ] . as can be seen , for each time interval
the centroids are distributed along the loop with those at higher energies being further away from the lt , and the entire pattern shifts toward the lt with time .
figure [ centr - dist ] shows the centroid positions of the northern half of the loop ( where the source motions are more evident ) along and perpendicular to the loop during the three intervals .
this again shows that higher energy emission is farther away from the lt and the centroids shifted towards the lt with time , but similarly there are some complicated patterns at low and intermediated energies .
all these are consistent with the general picture proposed above for the chromospheric evaporation process .
\3 . to further quantify the source motions
, we obtained the brightness - weighted standard deviation or the 2nd moment of the profiles .
in general the moment measures the compactness of the overall emission , but does not yield the sizes of individual sources whose measurement is still challenging for _ rhessi _ @xcite . hence our attention
should be paid to the general trend of the moment rather than its absolute values , which may be subject to large uncertainties and thus less meaningful .
the moments of the profiles resulting from clean images ( in three energy bands over 8 s intervals ) are plotted in figure [ time]_b_. there is a general decrease of the moment with the decline starting earlier at lower energies .
such a decrease is expected if the two fps move closer to each other . however , caution is required here because a decrease of this quantity could also come about by other causes , say by an increasing dominance of the brightest source .
we therefore checked the original images and the corresponding profiles when interpreting our results . to estimate the uncertainty of the moment , for each energy band we repeated the calculation with different integration time ( _ panel c _ ) .
the resulting moments remain essentially unchanged and as excepted the fluctuations of the moment decrease with increasing integration time .
we also plot in _ panel c _ the moment ( solid line ) obtained from pixon images with an integration time interval of two spin - period ( @xmath22 s ) , which basically agrees with its clean counterpart in the general trend .
the gradual decrease of the moment is consistent with the motion of the centroids of sources up along the legs of the loop , which can take place by a continuous increase of the gas density in the loop due to evaporation .
spectral analysis can be used to study the evaporation process as well . with an isothermal plus a power law model
, we fitted the spatially integrated _ rhessi _ spectra down to 6 kev @xcite for every 8 second interval during the impulsive phase . the emission measure ( @xmath23 ) and temperature of the isothermal component ( asterisk symbols )
are plotted in _
panels d _ and _ e _ of figure [ time ] , respectively .
the @xmath23 rises almost monotonically with time from 0.6 to 14.2 @xmath24 @xmath25 .
this translates into an increase of the plasma density ( @xmath26 ) by a factor of @xmath05 assuming a constant volume @xmath27 .
the temperature remains almost constant with a trend of slight decrease with time .
the @xmath23 and temperature derived from the _ goes _ data ( plus symbols ) are also shown for comparison . in general the _ goes _ results are smoother and the temperature increases monotonically but remains below that of the _ rhessi _ , consistent with previous results @xcite .
this is expected because _ rhessi _
is more sensitive to higher temperatures than _
goes_. however , surprisingly , the _ goes _ emission measure is also lower than that of _ rhessi _ as opposed to what is the case more generally ( see holman et al .
it is not clear whether or not this is due to a problem related to the _ rhessi _ calibration at low energies .
nevertheless , the continuous increase of the @xmath23 s at comparable rates does suggest a gradual increase of the plasma density .
the best fit parameters of the power law component with a low energy cutoff are plotted in figure [ time]_f_. the power - law index @xmath28 ( plus symbols ) is anti - correlated with the high energy light curves ( see figure [ time]_a _ ) and shows a soft - hard - soft behavior . it starts with 4.43 at 04:58:02 ut , drops to 3.82 at the impulsive peak ( 04:58:26 ut ) , and rises up to 7.12 at 04:59:46 ut . the high indexes ( @xmath29 ) may be an indicator of high - temperature thermal rather than non - thermal emission .
thus in what follows we limit our analysis to time up to 04:59:20 ut . the low - energy cutoff ( asterisks ) of the power law is about @xmath30 kev and is near the intersection of the isothermal ( exponential ) and power - law components .
the neupert effect is commonly quoted as a manifestation of chromospheric evaporation @xcite and a simple energy argument ( e.g. li et al . 1993 ) is often used to account for the relationship between sxr and hxr fluxes ( @xmath31 and @xmath32 ) . in the thick - target flare model , the non - thermal @xmath32 represents the _ instantaneous _ energy deposition rate ( @xmath33 ) by the electron beam precipitating to the chromosphere , but the thermal @xmath31 is proportional to the _ cumulative _ energy deposited , i.e. , the time - integral of @xmath33 .
it naturally follows that the temporal derivative of the sxr flux , @xmath34 , should be related to @xmath32 .
the simplest test of the neupert effect is usually carried out by plotting @xmath34 and @xmath32 in some energy band .
there are many reasons why a simple linear relationship would not be the case here .
the first and most important is that @xmath33 is related to @xmath32 through the bremsstrahlung yield function @xmath35 ( @xmath36 ) which is not a constant and depends on the spectrum of the electrons or hxrs ( see e.g. petrosian 1973 ) . here
the most crucial factor is the low energy cutoff ( @xmath37 ) of the non - thermal electrons , but the spectral index also plays some role .
the total yield of all the bremsstrahlung photons produced by a power law spectrum of electrons with energies above @xmath37 ( in units of 511 kev ) is @xmath38 and the yield of the photons whose energies are greater than @xmath37 is @xmath39 where @xmath40 , @xmath41 is the coulomb logarithm , and @xmath42 is the spectral index of the power - law electron flux .
as shown in figure [ time]_f _ both the low energy cutoff and spectral index of the non - thermal emission vary during the pulse , indicating variations in the electron spectrum and thus breaking the linearity of the sxr - hxr relationship .
other factors which can also produce further deviations are energy deposition by protons ( and other ions ) , by conduction and other possible ways of dissipation of energy than simply heating and evaporating the chromospheric plasma by non - thermal electrons .
a detailed treatment of the problem requires solutions of the combined transport and hd equations , which is beyond the scope of this paper .
@xcite by inclusion of some of these effects in an approximate way found that the expected relationship was mostly not present in several _ rhessi _ flares .
finally one must include the fact that the chromospheric response of sxr emission will be delayed by tens of seconds depending on the sound travel time ( and its variation ) and other factors .
the flare under study has shown no indication of gamma - ray line emission which means that the contribution of protons most probably is small . in the currently most favorable model
where the electrons are accelerated stochastically by turbulence ( see e.g. petrosian & liu 2004 ) the turbulence can suppress heat conduction during the impulsive phase and possibly also during the decay phase @xcite . because there does not appear to be large changes in the shape of the loop during the impulsive phase other energy dissipation processes such as cooling by expansion may also be negligible . assuming these to be the case we have performed the neupert effect test in two ways , the first being the common practice of examining the relation between @xmath34 and @xmath32 .
we then examine the relation between @xmath33 and @xmath34 by taking into account the variation of the bremsstrahlung yield .
the temporal derivatives of the fluxes of the two _ goes _ channels are shown in the _ bottom panel _ of figure [ lightcv ] . as evident , during the rising portion of the _ goes _ fluxes both channels derivatives indeed match the first pulse of the _ rhessi _ hxr light curves ( @xmath43 kev ) , but not during the second weaker pulse ( where the 1 - 8  derivative shows some instrumental artifacts ) .
this may be due to the fact that the neupert effect of the second pulse is overwhelmed by the cooling of the hot plasma produced during the first stronger pulse .
nevertheless , the sxr light curves ( of both _ goes _ and _ rhessi _ ) exhibit slightly slower decay rate than expected from the first pulse alone .
this most likely is the signature of the energy input by the second pulse , which slows down the decay of the first pulse alone .
we note in passing that the sxr light curves start rising several minutes prior to the onset of the hxr impulsive phase .
this is an indication of preheating of the plasma before production of a significant number of suprathermal electrons .
the 6 - 12 kev curve rises faster than the _ goes _ curves at lower photon energies , which is consistent with the picture that the primary energy release by reconnection occurs high in the corona where the relatively hotter plasma is heated before significant acceleration of electrons ( as suggested in petrosian & liu 2004 ) , and before transport of energy ( by accelerated electrons or conduction ) down the flare loop to lower atmosphere where cooler plasmas are heated subsequently and produce the _ goes _ flux . on the other hand , the increase of the sxr flux at the beginning is dominated by the southern loop , which shows little evidence of chromospheric evaporation . the phenomenon therefore may be a unique feature of this flare .    to quantify the sxr - hxr relationship ,
we cross - correlated the _ rhessi _
30 - 50 kev photon energy flux ( @xmath44 , figure [ neupert]_a _ ) and the derivative of the _ goes _ low energy channel flux ( @xmath34 , figure [ neupert]_c _ ) in the sxr rising phase ( 04:58:00 - 04:59:51 ut ) .
the resulting spearman rank - order correlation coefficient ( see figure [ neupert]_f _ ) , an indicator of an either linear or nonlinear correlation , shows a single hump with a maximum value of 0.91 ( corresponding to a significance of @xmath0@xmath45 ) at a time lag of 12 s. this suggests a delay of @xmath34 relative to @xmath44 , which is expected given the finite hydrodynamic response time ( on the order of the sound travel time of @xmath020 s for a loop size of @xmath0@xmath46 cm and @xmath47 k ) required for redistribution of the deposited energy .
such a delay is evident in the numerical simulations of @xcite who in addition found that the density enhancement contributes more to the total sxr emissivity than the temperature increase for longer duration ( @xmath48 s ) hxr bursts during the decay phase . in figure
[ neupert]_d _ , we plot the two quantities with the _ goes _ derivative shifted backward by 12 s to compensate the lag of their correlation .
a linear regression ( dotted line ) gives @xmath49 with an adjusted coefficient of determination ( or so - called r squared ) @xmath50 close to 1 suggesting a good linear correlation .
we also carried out the same analysis for the electron energy power @xmath33 , assuming a thick - target model of power - law electrons with a low - energy cutoff of @xmath51 kev .
we first obtained the energy flux of all the photons with energies greater than @xmath37 , @xmath52 , from the 30 - 50 kev photon energy flux @xmath44 : @xmath53 where @xmath54 is the photon flux distribution at the sun ( photons kev@xmath2 s@xmath2 ) which is obtained from spectrum fitting ( see  [ spec ] ) and assumed to extend to infinity in energy space .
we then calculated the power of the electrons by @xmath55 where the bremsstrahlung yield @xmath56 is given by equation ( [ y ] ) . ) here .
however , one can still use equation ( [ y ] ) with a simple correction factor of @xmath57 in the range @xmath58 to achieve an accuracy of @xmath59 . ]
the resulting @xmath33 is plotted versus time and _ goes _ derivative in _ panels b _ and _ e _ of figure [ neupert ] respectively .
the dotted line in _ panel e _ shows a linear fit ( @xmath60 ) to the data : @xmath61 the corresponding spearman correlation coefficient has a peak value of 0.78 ( significance @xmath0@xmath62 ) at a time lag of 3 s ( figure [ neupert]_f _ ) . as evident ,
@xmath33 yields no better correlation with @xmath34 than @xmath44 , which is similar to the conclusion reached by @xcite . during the hxr decay phase ( after 04:59:20 ut )
the spectrum becomes softer ( @xmath63 ) and @xmath64 decreases much slower than @xmath44 since the bremsstrahlung yield ( equation [ [ y ] ] ) decreases with the spectral index . as noted above for these high spectral indexes , the emission might be thermal rather than non - thermal .
the inferred electron power is thus highly uncertain for these times .
as stated earlier , the total energy of the non - thermal electrons is very sensitive to the low - energy cutoff @xmath37 which is generally not well determined ( cf .
sui et al .
we thus set @xmath37 as a free parameter and repeated the above calculation for different values of @xmath37 ( ranging from 15 to 28 kev ) .
we find that , as expected , the temporal @xmath33-@xmath34 relationship highly depends on @xmath37 . for a small @xmath37 ( @xmath65 kev ) , @xmath33 keeps rising till @xmath004:59:50 ut ( near the bottom of the @xmath44 light curve ) , which makes the @xmath33-@xmath34 correlation completely disappear . on the other hand , for a large @xmath37 ( @xmath66 kev )
the correlation is generally good during the impulsive pulse ( through 04:59:10 ut ) and the larger @xmath37 the better the correlation .
this is because the conversion factor @xmath67 in equation ( [ factor ] ) is an increasing ( decreasing ) function of the photon spectral index @xmath28 if @xmath37 is sufficiently small ( large ) . for a small @xmath37 , for example
, the photon energy flux @xmath52 may have a somewhat large value in the valley of the @xmath44 light curve when @xmath28 is high .
in addition , during this time interval the bremsstrahlung yield @xmath56 becomes small since @xmath42 is large ( see equation [ [ y ] ] ) and consequently this may result in a very large @xmath33 by equation ( [ edot ] ) .    as to the magnitude of the energy flux of non - thermal electrons , @xcite in their hd simulations found that the dynamics of the flare loop plasma is very sensitive to its value .
for a low energyflux ( @xmath68 ergs @xmath69 s@xmath2 ) , the up - flow velocity of the evaporating plasma is @xmath010 s of km s@xmath2 ; for a high energy flux ( @xmath70 ergs @xmath69 s@xmath2 ) , a maximum up - flow velocity of @xmath0100 s of km s@xmath2 can be produced .
for the flare under study , we estimate the area of the cross - section of the loop to be @xmath71 @xmath72 , where the upper limit corresponds to the loop width determined by the 5% level in figure [ loop - model]_b_. we read the maximum electron power of @xmath73 ergs s@xmath2 from figure [ neupert]_b _ , which is then divided by @xmath74 ( assuming a filling factor of unity ) to yield the corresponding electron energy flux : @xmath75 ergs @xmath69 s@xmath2 . the source velocity on the order of a few hundreds of km s@xmath2 estimated in  [ morp ] is consistent with that predicted by @xcite . for comparison
, we note that @xcite from _ rhessi _ data also obtained an energy flux of @xmath76 ergs @xmath69 s@xmath2 for an m2.2 flare during which an up - flow velocity of @xmath0230 km s@xmath2 was inferred from simultaneous co - spatial _
soho_/cds doppler observations . in summary ,
the _ goes _ sxr flux derivative @xmath34 exhibits a neupert - type linear correlation with the _ rhessi _ hxr flux @xmath44 during the first hxr pulse .
however , unexpectedly , the correlation between the electron power @xmath33 and @xmath34 is not well established based on the simple analysis presented here , which suggests that a full hd treatment is needed to investigate the chromospheric evaporation phenomenon ( see discussions in  [ conclusion ] ) .
for the 1994 june 20 disk flare , @xcite interpreted the moving sxr sources as thermal emission from the hot ( @xmath77 mk ) plasma evaporated from the chromosphere based on the good agreement of the emission measure of the blue - shifted component and that of the sxr from the fps . for the limb flare under study here ,
doppler shift measurements are not available . meanwhile
, a purely thermal scenario would have difficulties in explaining the systematic shift of the centroids towards the fps with increasing energies up to @xmath070 kev as shown in figure [ centr - dist ] .
a non - thermal scenario appears more appropriate .
that is , the apparent hxr fp structure and motions can result from a decrease in the stopping distance of the non - thermal electrons with decreasing energy and/or increasing ambient plasma density caused by the chromospheric evaporation ( as noted earlier in  [ morp ] ) .
one can therefore derive the density distribution along the loop from the corresponding x - ray emission distributions ( e.g. , figure [ prof_24s ] ) , without any pre - assumed density model ( c.f . * ? ? ?
this approach is described as follows .    for a power - law x - ray spectrum produced by an injected power - law electron spectrum
, @xcite obtained a simple empirical relation ( also see * ? ? ? * ) for the x - ray intensity @xmath78 per unit photon energy @xmath79 ( in units of 511 kev ) and column depth @xmath80 ( in units of @xmath81 = 5 \times 10^{22}$ ] @xmath69 for @xmath82 cm and @xmath83 ) : @xmath84 where @xmath28 and @xmath42 @xmath85 are the photon and electron spectral indexes , respectively , @xmath86 is a constant normalization factor , and @xmath87 , where @xmath88 is the distance measured from the injection site .
this equation quantifies the dependence of the emission profile ( or source morphology ) on the electron spectral index and column depth . in general , when @xmath42 decreases ( spectrum hardening ) , the intensity at a given photon energy rises ( drops ) at large ( small ) @xmath80 s and thus the emission centroid shifts to larger @xmath80 s .
this is expected because for a harder spectrum , there are relatively more high energy electrons that can penetrate to larger column depths and produce relatively more bremsstrahlung photons there .
the opposite will happen when the spectrum becomes softer . during the impulsive peak showing a soft - hard - soft behavior ( see  [ spec ] ) , one would expect that the emission centroids shift first away from and then back toward the lt ( if the density in the loop stays constant ) .
knowing the spectral index , the emission profile can therefore yield critical information about the density variation in both space and time .    to compare the above empirical relation with observations , we first integrate @xmath78 over an energy range @xmath89 $ ] , @xmath90 and then integrate @xmath91 over @xmath80 to obtain the cumulative emission , @xmath92              k^{-\gamma } dk ,    \label{ftau}\ ] ] where we have chosen @xmath93
so that @xmath94 . comparison of @xmath95 with the observed emission profiles gives the column depth @xmath96 whose derivative with respect to @xmath88 then gives the density profile along the loop .    specifically for this flare , we assume that the non - thermal electrons are injected at the lt indicated by the middle vertical dotted line in figure [ prof_24s ] and denote the profile to the right - hand side of this line ( i.e. , along the northern half of the loop ) as @xmath97 , where @xmath89 $ ] is the energy band of the profile . the observed cumulative emission is then given by , @xmath98 where @xmath99 ( corresponding to @xmath100 ) is the maximum distance considered and @xmath101 has been properly normalized .
then @xmath102 can be obtained by inverting @xmath103 where the integration over @xmath79 in equation ( [ ftau ] ) can be calculated numerically .
it should noted that , however , not all the profiles in figure [ prof_24s ] are suitable for this calculation , because low energy emission is dominated by a thermal component especially in the lt region and at later times .
we thus restrict ourselves to the energy ranges of 12 - 72 , 13 - 72 , and 17 - 72 kev , respectively for the three 24 s intervals .
the lower bound is the energy above which the power - law component dominates over the thermal component , determined from fits to the spatially integrated spectrum for each interval as shown in figure [ spec_24s ] . within these energy ranges , separate leg or fp sources rather than
a single lt source can be identified in the corresponding image , which is morphologically consistent with the non - thermal nature of emission assumed here . to further minimize the contamination of the thermal emission in our analysis ,
we have excluded the lt portion of the emission profile in excess of the lowest local minimum ( if it exists ) between the lt and leg ( or fp ) sources .
an example of this exclusion is illustrated with the hatched region in figure [ prof_24s]_c _ for the 19 - 21  kev profile .
this was done by simply replacing the profile values between the lt and the local minimum positions with the value at the minimum .
we calculated @xmath104 for every emission profile within the energy ranges mentioned above for the three intervals in figure [ prof_24s ] , with the photon indexes , @xmath105 4.46 , 3.97 , and 4.23 respectively .
from the geometric mean of the column depths obtained at different energies , @xmath106 , we derived the density profile @xmath107 for each time interval .
the results are shown in figure [ dens - fig ] , where we shall bear in mind that attention should be paid to the overall trend rather than the details of the density profile and its variation , because the profile here only spans about 3 times the resolution ( @xmath108 ) and thus is smoothed , making neighboring points not independent .
as can be seen , between the 1st and 2nd intervals , the density increases dramatically in the lower part of the loop , while the density near the lt remains essentially unchanged .
the density enhancement then shifts to the lt from the 2nd to the 3rd interval .
this indicates a mass flow from the chromosphere to the lt .
the density in the whole loop is about doubled over the three intervals , which is roughly consistent with the density change inferred from the emission measure ( see figure [ time]_d _ ) .
these results are again compatible with the chromospheric evaporation picture discussed in  [ morp ] .
we have presented in this paper a study of _ rhessi _ images and spectra of the 2003 november 13 m1.7 flare .
_ rhessi _ s superior capabilities reveal great details of the hxr source morphology at different energies and its evolution during the impulsive phase .
the main findings of this paper are as follows .
( 1 ) the energy dependent source morphology in general shows a gradual shift of emission from the lt to the fps with increasing energies . over some short integration intervals emission from the loop legs
may dominate at intermediate energies .
( 2 ) the emission centroids move toward the lt along the loop during the rising and plateau portions of the impulsive phase .
this motion starts at low energies and proceeds to high energies .
we estimate the mean velocity of the motion to be hundreds of km s@xmath2 , which agrees with the prediction of the hydrodynamic simulations by @xcite .
there are also shorter time scale variations that imply much higher velocities ( @xmath0@xmath19 km s@xmath2 ) but we are not certain if they are real because of instrumental limitations .
( 3 ) fits to the spatially integrated _ rhessi _ spectra with a thermal plus a power law model reveal a continuous increase of the emission measure while the temperature does not change significantly .
the _ goes _ data show a similar trend of the @xmath23 but a gradual increase of the temperature .
( 4 ) the time derivative of the _ goes _ sxr flux is correlated with the _ rhessi _ hxr flux with a peak correlation coefficient of 0.91 at a delay of 12 s in agreement with the general trend expected from the neupert effect .
however , the correlation between the electron power and the _ goes _ derivative is no better than the sxr - hxr correlation .
( 5 ) from the observed brightness profiles we derive the spatial and temporal variation of the plasma density in the loop , assuming a non - thermal thick - target bremsstrahlung model .
we find a continuous increase of the density , starting first at the fps and legs and then reaching to the lt .
all these results fit into a picture of continuous chromospheric evaporation caused by the deposition of energy of electrons accelerated during the impulsive phase .
several of the new features of this event ( such as the leg emission at intermediate energies ) may be common to many solar flares .
expanding the sample of flares of this kind will be very helpful to understand the underlying physical processes .
the new findings reside near the limit of _ rhessi _ s current temporal , spatial , and spectral resolution . as
advanced imaging spectroscopy capabilities are being developed and spatial resolution is being improved in the _ rhessi _ software @xcite , it will be critical to obtain the spatially resolved photon spectrum along the loop .
this will yield incisive clues to the nature of the moving x - ray sources and relevant energy transport mechanisms and will be useful to check the reality of the short time scale variations .
there are several important questions that need to be further addressed in future observational and theoretical investigations : ( 1 ) what is the nature of the moving x - ray sources ?
could they be characterized as thermal emission from the evaporated hot plasma or non - thermal emission from the precipitating electrons , or a mixture of both ? could they be related to mhd waves or evaporation fronts ?
( 2 ) what are the roles of different heating agents of the chromosphere , i.e. , electron beams , thermal conduction , and/or direct heating by turbulence or plasma waves during the impulsive phase ?
we have pointed out to some of the many physical processes that come into play in answering such questions . here
we describe possible directions for future theoretical studies .
we have shown that a more physical based test of the neupert effect between the electron power and sxr flux derivative does not reveal a better correlation than the usual hxr vs sxr derivative correlation .
although the observed source velocity agrees with those of hd simulations , there are some features that current simulations have not addressed . to answer these questions requires an updated numerical calculation where one combines the model of particle acceleration and transport with the hd simulation of the atmospheric response to energy deposition to form a unified picture of solar flares . for example
, one can use the output electron spectrum from the stochastic particle acceleration model @xcite as the input to the transport and hd codes , rather than simply assume a power - law electron spectrum as previous hd simulations .
such a study can shed light on the relative importance of particle beams and thermal conduction in evaporating chromospheric plasma and the roles that mhd waves may play in heating the flaring plasma , in particular , addressing our tentative observation of the fast source motion which suggests possible presence of mhd waves in the flare loop
. a better understanding of their propagation , damping , and excitation mechanisms is necessary for uncovering the energy release process during flares .
this is particularly true in the context of the stochastic particle acceleration model .