the present paper is a sequel of a previous work@xcite concerning the evaluation of saddle - point integrals of the form @xmath0\,\mathrm{d}s , \label{sd.1}\ ] ] where @xmath1 is a suitable integration path in the complex @xmath2-plane , @xmath3 and @xmath4 are functions which , for simplicity , will be assumed to be nonsingular , and @xmath5 will be intended as a `` large '' ( in modulus ) complex parameter . as is well known , the numerical evaluation of integrals of the kind in eq .
( [ sd.1 ] ) is customarily required for solving several classes of physical problems , occurring in optics , quantum mechanics , statistical physics , fluid mechanics , and so on . in optics ,
the evaluation of several diffraction integrals is customarily carried out asymptotically by identifying the parameter @xmath5 as the wavenumber of the radiation . in quantum mechanics , the same role
is played by the inverse of the planck s constant , while in fluid mechanics by the reynold s number@xcite .    in the stationary phase treatment of diffraction
integrals the values of the associated complex wavefield are asymptotically evaluated by taking the contributions coming from the stationary points of @xmath4 , each of them associated to a  ray " in the corresponding geometrical picture .
of particular importance is the birth and the death , as the spatial parameters ruling the  phase " function @xmath6 vary , of  evanescent " rays across sets of codimension 1 , named  stokes sets " @xcite .
the @xmath7- , or weniger , transformation@xcite ( wt for short ) , is particularly efficient for resumming the factorial divergent asymptotic series well away from stokes sets , as well as sets where two or more saddles are symmetrically placed in the complex singulant space@xcite .
unfortunately , as with other resummation techniques@xcite , the wt fails to perform across stokes sets .
the reason of such a failure stems to the extreme `` specialization '' of the transformation itself , which requires , for a successful resummation , an alternating sign pattern of the sequence of the single terms of the series@xcite .
several methods have been conceived for resumming nonalternating , slowly convergent or divergent , sequences @xcite , some of them being based on the serial combination of various resummation techniques@xcite . for the class of saddle - point integrals in eq .
( [ sd.1 ] ) , the marriage between hyperasymptotics@xcite ( h for short ) and the wt @xcite , generating the so - called h - wt ( which stands for hyperasymptotic - weniger transformation ) , allows the wt to successfully operate also across stokes sets .    basically , the h - wt consists in the sequential application , to the integral in eq .
( [ sd.1 ] ) , of a classical hyperasymptotic treatment , as described in ref .
@xcite , followed by the action of the wt on all asymptotic divergent series generated by h. in particular , the results obtained have shown how the 1st - order h - wt , for which only the first - stage of h is anticipated to the wt , is able to provide relative errors several orders of magnitude smaller than those achievable via the use of full hyperasymptotic treatments and with considerably lighter computational complexity and effort .
a key aspect is that , differently from h , the first truncation operated on the starting asymptotic series has not to be an optimal , in the sense of superasymptotics ( i.e. , at the least term ) one , but rather the corresponding truncation order , say @xmath8 , must be used as a free parameter for the subsequent application of the wt .    a question which was not addressed in ref .
@xcite , but mentioned only in the last sentence , is whether wt and h can be combined to higher orders in h , and if so , how the accuracy improves with order .
the present paper is aimed at giving a first answer to such a question .
we shall limit our analysis only to the second stage of h. further increasing of the h - wt order would be achievable along the same guidelines outlined here . on the other hand , it should be noted how , on increasing the order of h , the number of asymptotic series associated to the corresponding remainder that have to be resummed exponentially grows for topologies involving more than two saddles and , at the same time , the number of free parameters ( i.e. , the truncation orders at each h step ) linearly increases .
accordingly , from a mere computational viewpoint , it is mandatory to find a compromise between the h - wt order and the computational effort .
some limits of the 1st - order h - wt have already been emphasized in ref .
@xcite , where asymptotic evaluations of saddle - point integrals for  small " values of the asymptotic parameter were considered . in such cases , in fact , to increase the parameter @xmath8 of the 1st - order h - wt does not necessarily lead to an improvement of the reached accuracy , but often results in the opposite , i.e. , a worsening of it .
it is just the above scenario that we are interested in when the 2nd - level h - wt will be developed .
numerical experiments will be carried out on the class of saddle - point integrals already considered in the numerical sections of the first paper @xcite . moreover
, asymptotic evaluations of the so - called swallowtail diffraction catastrophe @xcite will be proposed as a new numerical experiment .
the swallowtail function is defined via eq .
( [ sd.1 ] ) with @xmath4 being a 5th - order polynomial with respect to @xmath2 , thus involving a four saddle network .
we will present a study of the accuracy achievable via h - wt asymptotic evaluations of the swallowtail diffraction catastrophe for points placed at the stokes set , following the prescriptions by berry and howls@xcite . in doing so
, we will find that the corresponding asymptotic expanding coefficients can be expressed in closed - form terms .
for any practical implementation of the h - wt , a key role is played by the numerical evaluation of the corresponding hyperterminants@xcite which are defined through suitable multiple integrals .
for the lowest - order hyperterminant the exact analytical expression is available from literature@xcite , but unfortunately this is not true for higher - order hyperterminants , included those involved in the 2nd - level h - wt . in the present paper
we solve the problem of the 2nd - level hyperterminant exact evaluation for a particular , but very important , choice of the hyperterminant parameters , which often occurs in the implementation of h for evaluating a wide class of saddle - point integrals . up to our knowledge
, this is a new result which also provides an interesting connection of such hypeterminants to the meijer - g functions @xcite .
moreover , although the closed - form evaluation of 2nd - level hypeterminants for arbitrary choices of their parameters seems to remain an open problem , in the present paper we find a semi - analytical representation which turns out to be suitable for numerical calculations via standard integration packages . similarly as done in ref .
@xcite , one of our aims is to keep the paper reasonably self - consistent .
accordingly , in the next section a brief review of h , up to the 2nd - level , is given .
as far as the wt is concerned , we believe that what is contained in ref . @xcite , together with the extensive bibliography , should be enough also for a nonexpert reader . for this reason , we do not repeat it in the present paper .
for simplicity , we shall refer to the asymptotic evaluation of saddle - point integrals of the type in eq .
( [ sd.1 ] ) where the set of saddle points of @xmath4 will be denoted @xmath9 and the integration path @xmath1 will be thought of as the union of a finite number of steepest descent arcs each of them , say @xmath10 , passing through the contributive saddle point @xmath11 , which will be supposed to be a simple one .
accordingly , the quantity @xmath12 can generally be written as @xmath13\,\mathrm{d}s= \displaystyle\sum_{n \in \mathcal{s}'}\,\mathcal{i}^{(n)}(k ) , \label{sdreview.5}\ ] ] where @xmath14 denotes the subset of @xmath9 containing all the contributive saddles , and @xmath15\,\mathrm{d}s .
\label{sdreview.5.1}\ ] ] the last integral can be written as@xcite @xmath16 where @xmath17 , and where @xmath18 can _ formally _ be written through the following asymptotic series expansion : @xmath19 the expanding coefficients @xmath20 being expressed via the integral representation @xcite @xmath21^{r+1/2}}\,\mathrm{d}s , \label{sdreview.5.0}\ ] ] where the subscript @xmath22 denotes a small positive loop around the saddle @xmath11 .
h starts by writing eq .
( [ sdreview.4 ] ) in the form @xmath23 where @xmath8 represents a positive integer and @xmath24 , denotes the corresponding remainder which , due to the diverging character of the asymptotic series , turns out to be a _ diverging _
quantity too .
h is based on a couple of fundamental results , found via a nontrivial analysis in ref .
the first is that the value of the expanding coefficients @xmath20 at the saddle @xmath11 is strictly related to the values of the expanding coefficients @xmath25 at all those saddles , say @xmath26 , which are _ adjacent _ to @xmath11 , via the following formal resurgence relation@xcite : @xmath27 where @xmath28 denotes the set containing the indexes pertinent to all saddles adjacent to @xmath11 , the quantities @xmath29 , called _
singulants _ , are defined by @xmath30 and the binary quantities @xmath31 are obtained through a topological rule@xcite .
the other fundamental tool of h is the following integral representation of the remainder @xmath32@xcite : @xmath33 equations ( [ sdreview.4])-([remainder.2 ] ) allow hyperasymptotic expansions for the saddle integral in eq .
( [ sdreview.5.1 ] ) to be built up in principle to any order@xcite .
for instance , the direct substitution of eq .
( [ sdreview.4 ] ) into eq .
( [ remainder.2 ] ) leads to the 1st - stage hyperasymptotic expansion @xcite , @xmath34 where the function @xmath35 , called _ hyperterminant _ of order 1@xcite , is defined through the integral @xmath36 where , in order for it to converge , @xmath37 .
moreover , it can be shown that@xcite @xmath38 where @xmath39 denotes the exponential integral function @xcite , while @xmath40 equals 1 if @xmath41 and zero otherwise .
the presence of the term containing @xmath40 has to be ascribed to the evaluation of the integral in eq .
( [ remainder.3.1 ] ) , when @xmath41 , in the cauchy principal - value sense
. equation ( [ istage.complete ] ) represents the first hyperasymptotic stage , at which the divergence of the asymptotic series in eq .
( [ sdreview.4 ] ) is led back to the presence of adjacent saddles @xcite .
furthermore , the asymptotic series in eq .
( [ istage.complete ] ) are only formal , since for @xmath42 the terminant @xmath43 diverges . in ref .
@xcite eq .
( [ istage.complete ] ) was taken as the starting point for introducing the h - wt .
in particular , instead of using the wt directly on the single terms of the series in eq .
( [ sdreview.4 ] ) , it is employed for resumming the asymptotic series associated to all saddles @xmath44 , with @xmath45 , which appear in eq .
( [ istage.complete ] ) . of course , due to the fact that @xmath46 in eq .
( [ istage.complete ] ) , it is mandatory that @xmath8 be left as a free parameter , in order for the wt to be able in decoding the above asymptotic series .
the 2nd - level h can be derived by truncating each of the asymptotic series in eq .
( [ istage.complete ] ) at an order , say @xmath47 , and by generating , for _ each _ adjacent saddle @xmath44 , with @xmath48 , a list of asymptotic series associated to all saddles , say @xmath49 , such that @xmath50 . in appendix
[ iistage ] , only for the reader convenience , the derivation of the 2nd - level hyperasymptotic expansion of the integral in eq .
( [ sd.1 ] ) is briefly recalled according to the formalism of ref .
in particular , it is found that @xmath51 where @xmath52 , the hyperterminant of order 2 , is now defined through the double integral @xmath53 similarly as we done in ref .
@xcite , eq .
( [ iistage.complete ] ) can be used to give estimates of @xmath18 , as functions of the two ( free ) parameters @xmath8 and @xmath47 , by resumming , via the wt , all asymptotic series generated at the second stage of h which are inside the double sum with respect to @xmath54 and @xmath55 .
the numerical evaluation of the hyperterminants represents a step of fundamental importance for any practical implementation of the h - wt algorithm .
unfortunately , differently from the lowest - order h - wt , for which the corresponding hyperterminants are achievable via the closed - form expression in eq .
( [ remainder.3.1.1 ] ) , there are no analytical expressions available for higher - order hyperterminants . in a series of important papers ,
olde daalhuis  @xcite addressed the general problem of the hyperterminants evaluation , up to arbitrary precisions , through the use of convergent series representations based on hypergeometric functions .
however , for the particular case of the 2nd - level hyperterminant , it seems that some new , at least up to our knowledge , results could be established .    from eq .
( [ terminant2 ] ) where , in order for it to converge , @xmath37 and @xmath56 , the hyperterminant can be written as @xmath57 and , by formally expanding the factor @xmath58 as a geometric series , after some algebra takes the form @xmath59 where @xmath60 denotes the hypergeometric function@xcite . although the series in eq .
( [ k2.0.2 ] ) is divergent , it can be decoded via borel summation , i.e. , by replacing the term @xmath61 by its integral representation , i.e. , @xmath62 which , once substituted into eq .
( [ k2.0.2 ] ) , leads to @xmath63 although , as we shall in a moment , it is possible to express the series inside the last equation through a closed form , it is better to carry out the evaluations for the case @xmath64 and @xmath65 separately .    on letting into eq .
( [ k2.0.3.1 ] ) @xmath64 we have @xmath66 the integral in eq .
( [ k2.0.3.2 ] ) can be evaluated by using the representation of the hypergeometric function given by formula 9.34.7 of ref .
in particular , it turns out that @xmath67 where @xmath68 denotes the meijer function@xcite . finally , by using formulas 9.31.5 , 7.813.1 , and 9.31.2 of @xcite , after some algebra it is found that @xmath69 equation  ( [ terminant2.1 ] ) represents one of the main results of the present paper .
as we shall see in the numerical section , in applying the 2nd - level h - wt the evaluation of the hyperterminants @xmath52 is often required for @xmath64 .
this happens whenever the contributive saddle @xmath70 turns out to be adjacent to itself after two hyperasymptotic stages , i.e. , when @xmath71 into eq .
( [ iistage.complete])@xcite .    for @xmath65 , the series inside the integral in eq .
( [ k2.0.2 ] ) can still be expressed through a closed form , although the subsequent integral unfortunately not . however , a semi - analytical expression , which turns out to be suitable for being evaluated via standard numerical integration packages can be derived . in doing this , eq .
( [ k2.0.2 ] ) is first rewritten as @xmath72 , \end{array } \label{gn1.1}\ ] ] where @xmath73 so that the task is to evaluate the series in eq .
( [ gn1.2 ] ) for @xmath65 . on substituting from eq .
( [ k2.0.3 ] ) into eq .
( [ gn1.2 ] ) we have @xmath74 where use has been made of the relation[see ref .  @xcite p.  347
] @xmath75 finally , on writing eq .  ( [ gn1.3 ] ) as @xmath76 where @xmath77 denotes the pochhammer symbol , formula 6.7.1.8 of ref .
@xcite gives at once @xmath78 notice that , when @xmath79>0 $ ] , the function @xmath80 can also be evaluated through the alternative form @xmath81}{p}\ , f(1,1;n+1;1-\gamma p ) .
\end{array } \label{gn1.7}\ ] ] although it seems that the above expressions can not be further simplified , the numerical evaluation of the function @xmath80 can be done with high accuracies by using standard integration packages .
finally , it should be stressed that , for @xmath41 , the evaluation of the double integral defining @xmath52 has to be done , with respect to the @xmath82-variable , in the cauchy principal value sense , in order to overcome the singularity placed at @xmath83 .
this , in turn , implies that an extra term must be added to the result . in appendix
[ appc ] such term is analytically evaluated starting from the definition in eq .
( [ terminant2 ] ) , and turns out to be @xmath84 } { ( -\beta)^{n-1 } }          \,e_{n}(-\beta\gamma ) .
\end{array } \label{extraterm}\ ] ] + all subsequent numerical experiments will be done within the _ mathematica _ language .
consider the evaluation of the airy function , defined as @xmath85\ , \mathrm{d}s , \label{airy.1}\ ] ] which is of the form given in eq .
( [ sd.1 ] ) with @xmath86 , @xmath87 , and @xmath88 .
the detailed analysis of the saddle topology , as well as the expanding coefficients @xmath20 has been summarized in ref .
@xcite , so that here it will not be given again .
we only recall that the two saddles are @xmath89 and @xmath90 , and that @xmath91 , @xmath92 , @xmath93 .
we focus our attention on the evaluation of the airy function across the stokes line@xcite , i.e. , for @xmath94 , in order to compare the performances of the 2nd - level h - wt with respect those displayed by the 1st - order h - wt in the same situation .
more precisely , we write the argument of the airy function as @xmath95 , where @xmath96 is a real positive parameter , whose value coincides with the singulant @xmath97 .
the study of the asymptotic evaluation of the airy function across its stokes line has played a pivotal role in the development of several asymptotic technique , mainly in light of the relative simplicity of the involved saddle topology .
such a simplicity could help in grasping , whereas possible , some conceptual aspects related to the use of the h - wt . differently from what done in ref .
@xcite , where the relative error values were displayed via the use of tables , in the present paper we are going to resort to graphical visualizations , due to the presence of the two `` free '' parameters @xmath8 and @xmath47 . in the first experiment ,
whose results are shown in figure  [ figerrorairyf16 ] , the airy function is evaluated for @xmath98 .
note that the same experiment was carried out in ref .
@xcite via the use of h. the values of the relative error , obtained through the 1st - level h - wt , are shown , as black dots , versus the values of @xmath8 , reported on the abscissa axis .
for each value of @xmath8 , the values of the relative error obtained via the 2nd - level h - wt , with @xmath99 $ ] , are also plotted and , for the sake of clarity , are joined with lines of different color each of them , which corresponds to a different value of @xmath8 , departing from @xmath8 itself .
this can be noted from the figure , where it is immediately seen how the higher the @xmath8 , the longer the corresponding coloured `` leg '' is . from a first look to the figure , it appears that the relative error , obtained with both the 1st- and the 2nd - level h - wt , is lower bounded .
we shall find that all subsequent numerical experiments present the same characteristic . as a first remark
, it should be noting that the improvement of the estimate accuracy induced by the 2nd - level h - wt with respect to that obtained via the 1st - level one of the same order appears , at least for values of @xmath8 not too large , not to adequately refund the unavoidable increase of the computational complexity required by the application of the 2nd - level transformation .
for example , it is seen from fig .
[ figerrorairyf16 ] that a relative error of the order of @xmath100 , achieved through the 2nd - level h - wt with @xmath101 and @xmath102 , would be reached via a 1st - order h - wt by letting @xmath103 , but with a considerable saving of computational effort .
the above example clearly suggests the use of the 2nd - level h - wt only in those cases where the best accuracy attainable via the 1st - level h - wt turns out to be not adequate .
this happens , for instance , when the integral is attempted to be evaluated beyond the asymptotic realm . to put into evidence this aspect , figure  [ figerrorairyf14to2 ]
shows the same as in fig .
[ figerrorairyf16 ] , but for a decreasing sequence of values of @xmath96 , namely 14 ( a ) , 10 ( b ) , 6 ( c ) , and 2 ( d ) . in particular , in fig .
[ figerrorairyf14to2](d ) , where the airy function argument is located at a distance @xmath104 from the origin of the complex plane , the 1st - level h - wt provides a best error of the order of @xmath105 , achieved for @xmath106 .
higher accuracies are not allowed because the information gained at the first h stage turns out to be no longer sufficient to generate wt - resummable sequences .
the 2nd - level h - wt , on the other hand , provides a best error of the order of @xmath107 , which is attained for @xmath108 .
some intuitive insights about the resummation process associated to the h - wt could be grasped by noting that a lower bounded error is an intrinsic imprint of superasymptotic and hyperasymptotic resummations , whereas it is not generally featured by the application of the wt to alternating factorial divergent series@xcite .
accordingly , one should be inclined to think that such an error behavior could be ascribed to the presence of the  regularization " step operated by h on the raw input data . speaking within a more general context
, this should be somewhat related to the possible presence on nonanalytic , nonperturbative correction terms which can not be grasped simply by resummation processes , but rather require the use of  generalized nonanalytic expansions " @xcite .
+ in a second experiment concerning the airy function , the asymptotic parameter @xmath96 is let running within the interval @xmath109 $ ] and , for each value of @xmath96 , an exhaustive search of the optimal values of the truncations @xmath8 and @xmath110 , which minimize the 1st- and 2nd - level relative errors , respectively , is operated .
the results are shown in fig .
[ figoptimalerrorairy ] , where the optimal relative errors obtained via the 1st- ( open circles ) and the 2nd - level ( dots ) h - wt are shown as functions of @xmath96 .
the values of the optimal truncation @xmath8 for the 1st - level h - wt are also reported , versus @xmath96 , in fig .
[ figoptimalnairy ] , while those of @xmath8 and @xmath47 , for the 2nd - level h - wt , in fig .
[ figoptimalnmairy](a ) and ( b ) , respectively .
we will come back later on the above results .
+ in concluding the present section , however , we want to provide a table of explicit values obtained through the use of the 2nd - level h - wt .
we choose to evaluate the airy function for @xmath111 , for which the optimal setting of the truncation parameters turns out to be @xmath108 .
the preliminary step is the evaluation , through a simple wt , of the contribution , to the airy integral , coming from the saddle @xmath112 .
the result is shown in table [ table.1 ] .
furthermore , the subsequent action of the 2nd - level h - wt on the saddle @xmath113 is shown in table [ table.2 ] , where the complete estimates of the airy function provided are reported together with the corresponding values of the truncation @xmath114 $ ] , with @xmath115 .
the second numerical experiment concerns the evaluation of the instanton integral @xmath116\,\mathrm{d}s , \label{inst.1}\ ] ] with @xmath117 , already considered in ref .
@xcite as a simplified prototype for the modeling of instanton tunneling between symmetric double wells .
it was shown in ref .
@xcite that the integral in eq .
( [ inst.1 ] ) can be written as@xcite @xmath118 where , by referring to eq .
( [ sd.1 ] ) , @xmath119 , @xmath120 , and where @xmath1 is the steepest descent path connecting the points @xmath121 and @xmath122 via the lines @xmath123 and @xmath124 . the complete saddle topology , as well as the expressions of the expanding coefficients associated to all saddles have been described in ref .
in particular , there are three saddles , @xmath125 , @xmath126 , and @xmath127 , with @xmath128 , @xmath129 , and @xmath130 .
the saddles involved in the evaluation of @xmath12 are @xmath112 and @xmath131 , but only the latter requires a h - wt treatment , since the associated singulant is @xmath132 , and the corresponding asymptotic series turns out to be nonalternating .
furthermore , @xmath133 , @xmath134 , and @xmath135 , while we recall that the integral in eq .
( [ inst.1 ] ) can be expressed in closed form via @xmath136 , \label{inst.4}\ ] ] where @xmath137 denotes the @xmath22th - order modified bessel function of the first kind .
the first experiment concerns the evaluation of @xmath138 via the 1st- and the 2nd - level h - wt . in fig .
[ figinstantonek1ov2 ] it is seen how the 2nd - level relative error is bounded , with a minimum value of the order of @xmath139 , achieved for @xmath140 . on the opposite
, the 1st - level h - wt turns out to be completely inadequate to provide a reasonably accurate estimate of the function , due to the very low value of @xmath5 .
the searching for optimal values has also been carried out in the present case , but using , as the varying asymptotic parameter , @xmath141 $ ] , i.e. , where the 1st - order h - wt displays the worst results in terms of accuracy , as shown in fig .
2a of ref .
the error values are shown in fig .
[ figoptimalerrorinstanton ] , versus @xmath5 , while the optimal settings of @xmath8 and of @xmath110 are plotted , against @xmath5 , in fig .
[ figoptimalninstanton ] and  [ figoptimalnminstanton ] , respectively .
it is worth now comparing the results pertinent to the airy and the @xmath142 functions .
what we are going to show can seem at first sight somewhat surprising , but gives a possible first hint toward the understanding of the h - wt mechanisms .
for simplicity we shall refer to the 1st - level transformation , but the results will apply also to higher - order levels . in fig .
[ figcomparisonairyinst ] , the values of the relative error obtained for the airy function ( dots ) and for the instanton function ( solid curve ) are plotted , versus @xmath8 , when the values of the parameter @xmath96 and the parameter @xmath5 are numerically equal . in particular , figure ( a ) corresponds to @xmath143 , ( b ) to 7 , ( c ) to 12 , and ( d ) to 20 .
it is clearly seen that the behavior of the relative error follows basically the same law . to give a possible explanation of this , in fig .
[ figcompairyinst ] a pictorial representation of the complete saddle network and the complex integration path involved in the evaluation of the airy ( a ) and of the instanton ( b ) functions is plotted . in both pictures ,
the black dot denotes the saddle for which the h - wt is required . although the two saddle distributions are clearly different , they present some common features that , together with eq .
( [ istage.complete ] ) , are enough to justify what happens in fig .  [ figcomparisonairyinst ] .
each of the `` black '' saddles is adjacent to a single saddle . for the airy function @xmath144
is adjacent to @xmath145 , while for the instanton function @xmath146 is adjacent to @xmath145 .
the values of the corresponding singulants are @xmath96 and 1 , respectively . the use of the resurgence relation in eq .
( [ resurgence.1 ] ) now gives , for the two `` black '' saddles , @xmath147 for the airy function and @xmath148 for the instanton function . from the above equation it is seen that the behavior of the expanding coefficients follows the same asymptotic law as soon as @xmath149 . at the same time , however , the above equality guarantees that also the asymptotic laws for the expanding coefficients corresponding to the adjacent saddles is identical .
in fact , for the airy function the saddle adjacent to @xmath145 is @xmath144 itself , with a singulant value of @xmath150 .
as far as the instanton function is concerned , the saddles adjacent to @xmath145 are @xmath144 and @xmath146 , but for both of them the singulant values are -1 .
accordingly , the use of eq .
( [ resurgence.1 ] ) , together with the condition @xmath149 , provides again an equivalence between the asymptotic laws of @xmath151 for the airy function and @xmath151 for the instanton function . finally , on using eq .
( [ istage.complete ] ) it is not difficult to convince that the retrieving process is the same for the two functions at the 1st - level@xcite . leaving a deeper understanding of this phenomenon to future investigations , it is here worthwhile to point out how an immediate consequence of the above described  topological equivalence " could be the restriction of the study of the h - wt retrieving performances to a few classes of prototype test cases .      as a last numerical experiment we consider asymptotic evaluations of the so - called swallowtail diffraction catastrophe @xcite , which is defined through the following integral : @xmath152\ , \mathrm{d}s , \label{sw.1}\ ] ] which is of the form given in eq .
( [ sd.1 ] ) with @xmath119 , @xmath153 , and @xmath88 .
the integration path @xmath1 can be thought of as the union of steepest descent paths approaching , for @xmath154 , the directions @xmath155 , with @xmath156 .
although a systematic treatment of the swallowtail asymptotics , along the general classical rules recalled in sec .
[ notations ] , can be derived by paralleling the analysis carried out , for the pearcey function , in ref .
@xcite , up to our knowledge it is not present in the current literature . as shown in appendix [ swt ] , all expanding coefficients @xmath20
are given by @xmath157 where @xmath158 and where the polynomials @xmath159 are defined via the generating function formula @xmath160 it is not difficult to prove that eq .
( [ generatingc.3 ] ) allows the numerical evaluation of the polynomials @xmath161 to be efficiently performed via the use of the following recurrence rule , whose derivation is outlined in appendix [ recurrencerule ] : @xmath162 with the triggering values @xmath163 , @xmath164 , and @xmath165 .
the numerical experiments we are going to illustrate concern asymptotic evaluations of @xmath166 at points belonging to the corresponding stokes set , which has been estensively studied in ref . @xcite ( see , in particular , fig . 3 of this reference ) .
accordingly , the triplets @xmath167 have been chosen following the prescriptions given in ref .
@xcite , in order to investigate points at the intersection between the stokes surface and the plane @xmath168 , along the branch corresponding to @xmath169 and @xmath170 .
this leads to triplets of the form @xmath171 , with @xmath172 being a positive parameter@xcite .
we start by considering the case @xmath173 .
the saddle topology is constituted by four saddles , which are listed below together with the corresponding list of adjacent ones : @xmath174 with the orientation matrix being @xmath175 , \label{sw.1.2}\ ] ] where the values of @xmath176 which are not relevant for the present experiment have been replaced by dots .
three of the four saddles , namely @xmath145 , @xmath146 and @xmath177 , do contribute to the integral .
in particular , the integration path consists in the union of three steepest - descent arcs , the first connecting the points @xmath178 and @xmath179 passing through @xmath131 , the second connecting the point @xmath179 to the saddle @xmath112 , passing through @xmath180 , and the third connecting the saddle @xmath112 to the point @xmath181 .
accordingly , the stokes phenomenon occurs via the so - called _ saddle connection _ between saddles @xmath180 and @xmath112@xcite .
a pictorial representation of the topology above described is given , for the reader convenience , in fig .
[ figswpaths ] .    before showing the numerical results about the performances of the 1st- and the 2nd - level h - wt ,
it is worth giving some details about the way the saddle topology of the swallowtail integral influences the retrieving capabilities of the wt .
we refer , in particular , to the contribution , to the swallowtail integral , of @xmath146 and @xmath177 . as far as the former is concerned , since there is only one adjacent saddle , @xmath177 , due to the fact that the corresponding singulant @xmath182 is purely imaginary , it turns out that @xmath183 , thus allowing the wt to operate the resummation .
similar considerations can be done for the contributive saddle @xmath145 .
the situation is somewhat different for the saddle @xmath177 , which is connected to @xmath145 .
in fact , from the above described topology , it turns out that the adjacent saddle @xmath146 is dominant , i.e. , presents the minimum value of @xmath184 , with @xmath185 .
accordingly , one should conclude also for @xmath186 an asymptotic `` factorial divided by power '' law , similar to that corresponding to @xmath187 and , due to the fact that @xmath188 , one should expect the wt to be able in resumming the corresponding asymptotic series .
but this , on the contrary , does not happen , because of the presence of the other two , nondominant , saddles @xmath144 and @xmath145 , symmetrically placed in the complex singulant space . this can be explained by expliciting eq .
( [ resurgence.1 ] ) as @xmath189 , \label{resurgence.1bis}\ ] ] where @xmath190^{*}$ ] .
this example shows how the divergent asymptotic series generated on stokes sets do not necessarily display a strictly nonalternating sign pattern as , for example , happened for the airy function , but rather how the asymptotic behavior of their single terms can display more complex patterns , depending on the whole saddle topology .
figure  [ figswk2 ] shows the relative errors obtained through the 1st- and the 2nd - level h - wts in the case of swallowtail evaluation across the stokes set defined above , when @xmath173 . as for figs .
[ figerrorairyf16 ] ,  [ figerrorairyf14to2 ] , and   [ figinstantonek1ov2 ] , the errors are plotted versus the parameters @xmath8 and @xmath47 . in this case
their optimal values turn out to be @xmath106 ( for the 1st - level ) and @xmath191 ( for the 2nd - level ) , with corresponding error values of @xmath192 and @xmath193 , respectively , evaluated with respect to the  exact " value of @xmath194 , obtained via the method recently proposed in ref .
finally , an experiment about optimal resummation of swallowtail functions has been carried out , by using @xmath195 $ ] as the parameter representative of the  asymptoticity " features .
the error values are shown , versus @xmath172 , in fig .
[ figoptimalerrorsw ] , while the optimal settings of @xmath8 and of @xmath110 are plotted , against @xmath172 , in fig .
[ figoptimalnsw ] and  [ figoptimalnmsw ] , respectively .
the h - wt was introduced in ref .
@xcite as a powerful and easily implementable refinement of h aimed at allowing the wt to successfully decode those divergent asymptotic series generated through the application of the steepest descent method to saddle point integrals evaluated across stokes sets , for which their single terms do not display a strictly alternating sign pattern .
the scheme proposed in @xcite employed the wt only on the asymptotic series generated by the first - stage hyperasymptotic treatment of the corresponding diverging remainder . in the present sequel we reported on the possibility of combining h and the wt to higher orders in h. in particular
, the full development of the 2nd - level h - wt has been detailed within the classical framework of h. the results obtained from the application of the 2nd - level h - wt , also in comparison to those obtained via the 1st - level one , on the different types of saddle - point integrals considered , showed how the increase of complexity and computational effort required by the new transformation be adequately refunded , in terms of accuracy of the estimate , particularly when the integrals are evaluated for values of their parameters which are beyond the asymptotic regime , whereas h turns out to be inapplicable and the 1st - order h - wt unavoidably lacks precision . at the same time , however , it should be noted how , for  ordinary " asymptotic evaluations , at least in the cases considered in the present work , the performances of the 1st- and the 2nd - level h - wts seem to be comparable in terms of the estimate accuracy , against a considerably difference in the computational efforts required by the two transformations .
although the h - wt has been developed , here and in ref .
@xcite , with reference to the evaluation of saddle - point integrals of the type in eq .
( [ sd.1 ] ) , we believe it could be useful also in dealing with problems of different nature like , for instance , the hyperasymptotic treatment of a wide class of linear and nonlinear ordinary differential equations , which has been recently considered  @xcite . the semi - analytical algorithms proposed for the numerical evaluation of the 2nd - level hyperterminants would reveal useful in this perspective .    before concluding , it is worth pointing out , as an important open problem , the need for an _ a priori _ algorithm for estimating the values of @xmath8 and @xmath47 that lead to optimal results . especially in cases where it is not convenient ( or possible ) to evaluate the original function , such an algorithm would certainly be of great help to a typical user .
unfortunately , differently from h , for which the optimal settings of the hyperseries truncations are directly extracted from the singulant values@xcite , at present it does not seem possible to provide similar information for the h - wt . the difficulties in giving practical guidelines for the choice of @xmath8 and @xmath110 can also be appreciated from the results presented in sec . [ numericalexperiments ] and especially from figs .
[ figoptimalnairy ] , [ figoptimalnmairy ] , [ figoptimalninstanton ] , [ figoptimalnminstanton ] , [ figoptimalnsw ] , and [ figoptimalnmsw ] , where it seems quite difficult to obtain general rules for the optimal settings of them .
a possible hint , grasped from a quantitative comparison between the results obtained for the airy and the instanton functions , seems to be given by the strong connection between the h - wt retrieving performances and the saddle topology associated to the integral under consideration .
what we found is that different saddle networks can share a sort of  topological equivalence " property , which is related to the set of the saddles adjacent to that under consideration and to the values of the relevant singulants . if two networks turn out to be equivalent at a certain hyperasymptotic level , this would result in the same computational effort , in terms of relative error , as far as the corresponding h - wt retrieved estimates are concerned .
this , in particular , would imply that the study of the h - wt retrieving performances could be , in principle , carried out only on a restricted class of prototype functions .
i would like to thank all anonymous reviewers for their constructive criticisms and suggestions .
i am also grateful to turi maria spinozzi for his invaluable help during the preparation of the present work .
the starting point is again the integral representation of the remainder @xmath32 given in eq .
( [ remainder.2 ] ) , where we let @xmath196 with @xmath47 being a _ new _ truncation order .
substitution of eq .
( [ iistage.1 ] ) into eq .
( [ remainder.2 ] ) gives @xmath197 on using resurgence , i.e. , eq .
( [ resurgence.1 ] ) , on @xmath25 , namely @xmath198 and @xmath199 after straightforward algebra the second term in eq .
( [ iistage.2 ] ) becomes @xmath200 or , by interchanging the integral with the @xmath201-sum , @xmath202 eventually , on making the substitution @xmath203 , after some algebra the quantity in eq .
( [ iistage.3.2 ] ) can be written as @xmath204 where the 2nd - level hyperterminant @xmath52 , defined by @xmath205 has been introduced .
accordingly , the complete expression of the 2nd - level hyperasymptotic expansion is @xmath206 which coincides with eq .
( [ iistage.complete ] ) .
we first recast the double integral in eq .
( [ terminant2 ] ) as follows : @xmath208 which , if @xmath209 were not negative , would take the form @xmath210 ,    \end{array } \label{betanegativo.2}\ ] ] where the integrand turns out to be continuous even for real , positive values of @xmath209 , since @xmath211 if @xmath207 , eq .
( [ betanegativo.1 ] ) must be written as @xmath212 where @xmath213 denotes the principal part operator .
on evaluating the integral in eq .
( [ betanegativo.4 ] ) we have @xmath214\\      \\- \mathrm{i}\pi\exp(\beta)\,\displaystyle\frac{(-\beta)^{m}}{\beta - v}.      \end{array } \label{betanegativo.5}\ ] ] finally , the direct substitution of eq .
( [ betanegativo.5 ] ) into eq .
( [ betanegativo.4 ] ) leads to the closed - form of the extra term as @xmath215 } { ( -\beta)^{n-1 } }          \,e_{n}(-\beta\gamma ) .
\end{array } \label{betanegativo.6.1}\ ] ]
the swallowtail function in eq .
( [ sd.1 ] ) corresponds to let @xmath216 , @xmath88 , and @xmath217 in eq .
( [ sw.1 ] ) .
the ( four ) saddle points @xmath218 are solutions of the algebraic equation @xmath219 which can be exactly solved by using , for instance , cardano s formula .
the evaluation of the expanding coefficients @xmath20 requires to evaluate the integral in eq .
( [ sdreview.5.0 ] ) .
the first step is to expand @xmath220 around @xmath11 , i.e. , @xmath221=\\ \\
= \displaystyle\frac{(s - s_n)^2}{5\mathrm{i}}\,p(s ) , \end{array } \label{sw.3}\ ] ] where @xmath222 substitution into eq .
( [ sdreview.5.0 ] ) gives @xmath223^{r+1/2}}=\\ \\
= \displaystyle\frac{(r-1/2)!\,\,\,(5\,\mathrm{i})^{r+1/2}\,}{(2r)!}\ , \left\ { \displaystyle\frac{\mathrm{d}^{2r}}{\mathrm{d } s^{2r}}\ , \displaystyle\frac 1{[p(s)]^{r+1/2 } } \right\}_{s = s_n}. \end{array } \label{sw.5}\ ] ] equation  ( [ sw.5 ] ) can be further simplified by letting @xmath224 and then by changing @xmath225 into @xmath226 , which yields @xmath227_{t=0 } , \end{array } \label{generatingc.1.5.1}\ ] ] where @xmath228 note that eq .
( [ generatingc.1.5.1 ] ) can also be given the alternative following form : @xmath229 where the function @xmath230 is defined through @xmath231 which coincides with eq .
( [ generatingc.3 ] ) .
the starting point is the definition of @xmath161 through the generating function in eq .
( [ generatingc.3.1 ] ) which , once derived with respect to @xmath232 , gives @xmath233 and , by taking eq .
( [ generatingc.3.1 ] ) into account once again , leads to @xmath234 by operating the products term by term , by rearranging the series indices , after long but straightforward algebra eq .
( [ generatingc.4 ] ) follows .
00    r. borghi ,  joint use of the weniger transformation and hyperasymptotics for accurate asymptotic evaluations of a class of saddle - point integrals , " phys .
e * 78 , * 026703 ( 2008 ) .
m. born and e. wolf , _ principles of optics _
( cambridge university press , cambridge , 1999 ) .    m. v. berry ,  some quantum - to - classical asymptotics , " in les houches lecture series lii , eds .
giannoni , a. voros , and j. zinn - justin , chaos and quantum physics , elsevier north - holland , amsterdam , 251 - 304 ( 1991 ) .
f. j. wright ,  the stokes set of the cusp diffraction catastrophe , " j. phys . a * 13 , * 2913 - 2928 ( 1980 ) .
m. v. berry and c. j. howls ,  stokes surfaces of diffraction catastrophes with codimension three , " nonlinearity * 3 , * 281 - 291 ( 1990 ) .
e. j. weniger ,  nonlinear sequence transformations for the acceleration of convergence and the summation of divergent series , " comput .
phys . rep .
* 10 , * 189 - 371 ( 1989 ) .
e. j. weniger ,  mathematical properties of a new levin - type sequence transformation introduced by cizek , zamastil , and skala .
i. algebraic theory " , j. math .
* 45 , * 1209 - 1246 ( 2004 ) .    for an updated review about methods for decoding diverging series ,
see for instance e. caliceti , m. meyer - hermann , p. ribeca , a. surzhykov , and u. d. jentschura ,  from useful algorithms for slowly convergent series to physical predictions based on divergent perturbative expansions , " phys .
* 446 , * 1 - 96 ( 2007 ) .
arxiv:0707.1596v1 .    c. brezinski and m. redivo zaglia , _ extrapolation methods _ ( north - holland , amsterdam , 1991 ) .
u. d. jentschura , p. j. mohr , g. soff , and e. j. weniger ,  convergence acceleration via combined nonlinear - condensation transformations , " comput .
* 116 , * 28 - 54 ( 1999 ) .
d. jentschura ,  resummation of nonalternating divergent perturbative expansions , " phys .
d * 62 , * 076001 ( 2000 ) .
s. v. aksenov , m. a. savageau , u. d. jentschura , j. becher , g. soff , and p. j. mohr ,  application of the combined nonlinear - condensation transformation to problems in statistical analysis and theoretical physics , " comput .
* 150 , * 1 - 20 ( 2003 ) .
m. v. berry and c. j. howls ,  hyperasymptotics , " proc .
a * 430 , * 653 - 668 ( 1990 ) .
m. v. berry and c. j. howls ,
hyperasymptotics for integrals with saddles , " proc .
r. soc . a * 434 , * 657 - 675 ( 1991 ) .    m. v. berry and c. upstill ,  catastrophe optics : morphologies of caustics and their diffraction patterns , " prog
* xviii , * 257 - 346 ( 1980 ) .    i. s. gradshteyn and i. m. ryzhik , _ table of integrals , series , and products , _
, a. jeffrey and d. zwillinger , eds .
( academic , 2000 ) .
a. b. olde daalhuis ,
hyperterminants i , " j. comput .
appl . math .
* 76 , * 255 - 264 ( 1996 ) .
a. b. olde daalhuis ,  hyperterminants ii , " j. comput . appl . math . * 89 , * 87 - 95 ( 1998 ) .    a. p. prudnikov , yu
. a. brychkov , and o. i. marichev , _ integrals and series . vol .
( gordon breach publisher , new york , 1986 ) .
j. p. boyd ,  the devil s invention : asymptotic , superasymptotic and hyperasymptotic series , " acta applicandae mathematicae * 56 , * 1 - 98 ( 1999 ) .    m. v. berry and c. j. howls ,  infinity interpreted , " phys .
world * 6 , * 35 - 39 , ( 1993 ) .
it could be proved that the above  equivalence " between the airy and the instanton functions keeps its validity also for the 2nd - level h - wt .
j. n. l. connor , p. r. curtis , and d. farrelly , `` the uniform asymptotic swallowtail approximation : practical methods for oscillating integrals with four coalescing saddle points , '' j. phys . a : math
* 17 , * 281 - 310 ( 1984 ) .
j. f. nye , _ natural focusing and fine structures of light _ , ( institute of physics publishing , bristol , 1999 ) .
j. f. nye ,  dislocation lines in the swallowtail diffraction catastrophe , " proc .
a * 463 , * 343 - 355 ( 2007 ) .
r. borghi , `` on the numerical evaluation of cuspoid diffraction catastrophes , '' j. opt .
am . a * 25 , * 1682 - 1690 ( 2008 )
.    a. b. olde daalhuis ,
hyperasymptotics for nonlinear odes .
i. a riccati equation , " proc .
lond . a * 461 , * 2503 - 2520 ( 2005 ) .
a. b. olde daalhuis ,
hyperasymptotics for nonlinear odes ii . the first painlev equation and a second - order riccati equation , " proc .
a * 461 , * 3005 - 3021 ( 2005 ) .
m. v. berry and a. shelankov ,  the aharonov - bohm wave and the cornu spiral , " j. phys .
a * 32 , * l447-l455 ( 1999 ) .
a. sidi ,  a new method for deriving pad approximants for some hypergeometric functions , " j. comput .
appl . math .
* 7 , * 37 - 40 ( 1981 ) .
values , provided by the wt , of the contribution coming from the saddle @xmath112 in the evaluation of the airy function across its stokes set , when @xmath111 .
[ cols="^,^",options="header " , ]
fig . 1 - r. borghi


the cygnus loop is one of the brightest supernova remnants ( snrs ) in various wavelength domains due to its proximity of @xmath2540 pc @xcite and large angular size ( @xmath23@xmath3 in diameter ) .
it is considered to be a middle - aged snr ( @xmath21.4@xmath410@xmath5 yr ; levenson et al .
1998 ) , and to be interacting actively with the ambient medium .
these points make the remnant one of the most extensively observed and studied snrs .
many x - ray studies were performed to understand the spatial structures and spectral properties of the remnant @xcite . also , @xcite provided results from gamma - ray measurements . details of many optical filaments , including balmer filaments , were analyzed in various optical studies @xcite .
@xcite presented the optical h@xmath6 , [ ] , and [ ] images of the entire region .
the global ir images were presented by @xcite and @xcite , and the high resolution ir images of a non - radiative shock were presented by @xcite .
radio continuum studies at various frequencies provided results for spectral indices and polarization properties @xcite .
also , @xcite presented the results of the 21 cm observations over various velocity channels . in all wavelength ranges
, the global structure of the cygnus loop is well defined to be a shell - like feature with a south extension , and the outer rim of the shell is much brighter than the inner region .
this morphology can be explained by a scenario in which the cygnus loop exploded in a cavity , and its blast wave is now hitting the walls of the cavity @xcite .
many of the far - ultraviolet ( fuv ) observations for various optical filaments were made and several fuv emission lines originating from non - radiative or radiative shocks were detected and analyzed @xcite . however , the entire region of the cygnus loop was analyzed in only a few fuv studies .
@xcite reported maps of whole regions in the @xmath0977 and @xmath11032 , 1038 emission lines using the _ voyager 2 _ ultraviolet spectrometer , and @xcite presented those in the @xmath11032 , 1038 emission lines by a rocket - borne experiment , high resolution emission line spectrometer ( hires ) .
the spectroscopy of plasma evolution from astrophysical radiation ( spear ) , also known as far - ultraviolet imaging spectrograph ( fims ) , carried out the global fuv observation of the cygnus loop in two wavelength channels , 9001150  and 13401750  .
@xcite reported the early results using only the long wavelength channel data .
they showed the + @xmath11400 , 1403 ( unresolved ) , @xmath11548 , 1551 , @xmath01640.5 , and @xmath11661 , 1666 emission - line images , which were similar to the x - ray , optical , and radio - continuum images .
they also presented the spectra in the 13401750  range and estimated shock velocities from the line ratio of to .
in this paper , we analyze the short wavelength channel data of the spear / fims observations and obtain new spectral images in the long wavelength channel .
we create the @xmath0977 and @xmath11032 , 1038 emission - line images of the cygnus loop and compare them with the previous results .
we separate out the blended + image , originally presented in @xcite , into two individual images of and . a new image at @xmath01486 is also presented with the newly processed spear / fims data .
next , we compare the resulting fuv images with the x - ray , optical , and 21 cm images , and discuss the fuv properties of the cygnus loop on a global scale . from the newly processed data , we update the previous results for the fuv line intensities , ratios , and luminosities .
spear / fims is a dual - channel fuv imaging spectrograph for observing diffuse fuv emissions from the interstellar medium ( ism ) .
the short bandpass channel ( s - channel ) covers 9001150  with 4@xmath3.0@xmath44@xmath7.6 field of view , and the long bandpass channel ( l - channel ) covers 13401750  with 7@xmath3.4@xmath44@xmath7.3 field of view .
their spectral and spatial resolutions are @xmath8 and 5@xmath7 , respectively .
@xcite described the spear / fims mission and its science objectives in detail . the spear / fims instrument , its on - orbit performance , and the basic processing of the data are described in @xcite .    by combining the sky - survey data with the target - pointed data used in @xcite
, we utilized the data obtained from a total of 43 orbits .
we improved the accuracy of the attitude knowledge up to 5@xmath7 by using an updated version ( for sky - survey data ) of the software - correction code used in @xcite .
the procedure of correction was briefly described in @xcite .
we used the healpix pixelization scheme @xcite to make images and spectra .
a pixel size of @xmath23@xmath7.4 ( healpix resolution parameter ` nside ` = 1024 ) and a wavelength bin size of 0.5  were selected for the s - channel . on the other hand , ` nside ` = 2048 and 1  were used for the l - channel . to mask the pixels including fuv - bright stars
, we identified early - type stars in the area of present concern using the _ td-1 _ catalog , and removed a total of 15 stars in the l - channel and 3 stars in the s - channel .
for emission - line images , the spectrum extracted from each pixel was divided into small segments including individual emission lines . then each segment was fitted with gaussian line profiles plus a linear continuum . for better statistics ,
the pixel size was increased to @xmath26@xmath7.8 ( ` nside ` = 512 ) and the pixels with exposure time of @xmath91 second were excluded before the spectral fitting .
we fixed the central wavelengths and widths of the gaussian functions , to be the calibrated line centers and the spectral resolutions , respectively . the intrinsic line ratios of 2:1 for the @xmath11031.9 , 1037.6 and @xmath11548.2 , 1550.8 doublet lines approach 1:1 in an optically thick case due to resonant scattering @xcite .
therefore , we set the ratios for the doublet lines to vary between 1:1 and 2:1 . the lines at 1660.8 and 1666.1  were assumed to have an intensity ratio of 0.34:1 , according to their statistical weights .
we also attempted to fit the + line complexes as follows .
we assumed that the lines are a doublet line at 1393.8 and 1402.8  , while the lines are a quintet of lines at 1397.2 , 1399.8 , 1401.2 , 1404.8 , and 1407.4  .
however , since the @xmath01402.8  line was unable to be resolved from the quintet lines , only the @xmath01393.8  line was fitted along with the quintet lines .
the strength of the @xmath01402.8 line was assumed to be half that of the @xmath01393.8 line , although it would be valid only in an optically thin case .
then , the strength of the @xmath01402.8 line was subtracted from the value obtained for the quintet lines and added into that of the doublet
.                    * figure 1 .
* spear / fims ( a ) @xmath0977 , ( b ) @xmath11032 , 1038 , ( c ) @xmath11394 , 1403 , ( d ) @xmath01404 , ( e ) @xmath01486 , ( f ) @xmath11548 , 1551 , ( g ) @xmath01640.5 , and ( h ) @xmath11661 , 1666 emission - line images of the cygnus loop ( in units of 10@xmath10 erg s@xmath11 @xmath12 sr@xmath11 ) .
contour levels are stepped by seven equal intervals and range ( a ) from 15.0 to 36.0 , ( b ) from 25.0 to 70.0 , ( c ) from 1.3 to 6.1 , ( d ) from 1.5 to 9.5 , ( e ) from 1.2 to 4.2 , ( f ) from 3.5 to 20.0 , ( g ) from 3.0 to 12.0 , and ( h ) from 3.0 to 12.0 .
the solid lines indicate the outlines of the masked pixels .
figure 1 shows the resulting eight emission - line images .
the solid lines denote the masked bright pixels .
the images have been adaptively smoothed using a gaussian function with a position - variable kernel radius , which is the same as the adaptive kernel method adopted in @xcite .
the smoothing - kernel scales of @xmath23@xmath7 to @xmath230@xmath7 were used to obtain moderate signal - to - noise ratios ( s / n ) .
seven equally - spaced contour levels were overlaid on each image , and all pixels inside the lowest - level contours have s / n @xmath133 .
figures 1(a ) and ( b ) have lower spatial resolutions than the other images .
this is because the and lines belong to the s - channel ( with relatively lower sensitivity ) and consequently , the images have been smoothed with larger smoothing - kernel scales .
moreover , the eastern region has much lower exposure time for the s - channel observation , which may make the and local peaks in the eastern region weaker and more diffuse .
figures 1(c ) and ( d ) show a clear distinction between the deblended and images .
in particular , the emission line is found to be strongest around a local peak on the eastern edge , but relatively weak around the peaks in the western region .
even though the line was weakest in the l - channel , we could obtain an emission - line image in figure 1(e ) , which is compatible with the other emission - line images .
the strongest peak appears to be around the peak on the eastern edge . in figures 1(f)(h ) , the , , and emission - line images have been reproduced using the newly processed data , and do not show large differences from the images by @xcite .
* figure 2 . * ( a ) @xmath11548 , 1551 color - scale image overlaid on h@xmath6 grey - scale image .
( b ) , ( c ) @xmath01404 , and ( d ) ] @xmath11661 , 1666 color - scale images overlaid on [ ] @xmath05007 grey - scale images .
the color - scale images are the same as in figure 1 , and the grey - scale images are from @xcite .
the ellipses marked in ( a ) denote subregions discussed in the text .
( e ) @xmath11032 , 1038 contours overlaid on _ rosat _ pspc x - ray ( 0.12.4 kev ) intensity image and ( f ) 21 cm ( @xmath14 = -8.4 to 1.5 km s@xmath11 ) image with the 1420 mhz radio continuum contours @xcite .
the ellipse marked in ( f ) indicates the ring feature identified by @xcite .
the spear / fims fuv images are compared with the optical , x - ray , and radio images in figure 2 . in figures 2(a)(d ) , the , , and images are overlaid onto the h@xmath6 and [ ] @xmath05007 grey - scale images adapted from @xcite .
the figures show strong correlations between the fuv images and optical images , although their spatial resolutions are fairly different .
the ( the strongest - detected fuv emission line ) closely matches both optical images , not only in prominent features but also in faint features , as can be seen in figures 2(a ) and ( b ) .
ten subregions of the characteristic features , which will be discussed later , are denoted in figure 2(a ) . in figure 2(e ) , the contours of figure 1(b ) are overlaid onto the _ rosat _ pspc x - ray ( 0.12.4 kev ) intensity image plotted using the data from the skyview virtual observatory @xcite .
although the image has too low spatial resolution for detailed comparison , the clearly appears to peak around the x - ray peaks on the eastern , northern , and western edges .
additionally , the 21 cm image overlapped with the 1420 mhz radio continuum contours is shown in figure 2(f ) , which is adapted from @xcite .
the elliptical ring feature just inside the marked ellipse , is the identified cloud located on the front side of the cygnus loop @xcite .
we note that all of the spear / fims fuv images are faint in the direction of this cloud , as shown in figures 1 and 2(a)(d ) .
* figure 3 .
* spear / fims s - channel ( @xmath15 ) and l - channel ( @xmath16 ) spectra ( with 1 @xmath17 error bars ) for the whole region and several subregions indicated in figure 2(a ) .
the spectra are binned with 0.5 ( s - channel ) and 1 ( l - channel )  intervals , and smoothed with a boxcar average of 3 bins .
dashed lines denote the positions of the identified emission lines , and dotted lines indicate those of the hydrogen lyman series ( @xmath18 973  , @xmath19 1026  ) originating from the geocorona .    in figure 3 ,
the spear / fims fuv spectra for several subregions are plotted along with the spectra of the whole region .
the detector background estimated in @xcite was subtracted from the spectra .
the s - channel spectra ascertain the @xmath0977 and @xmath11031.9 , 1037.6 emission lines , together with the hydrogen lyman series at 973 and 1026  originating from the geocoronal airglow .
the l - channel spectra show the @xmath01393.8 , @xmath01404 , @xmath01486 , @xmath11548 , 1551 ( unresolved ) , @xmath01640.5 , and @xmath11661 , 1666 emission lines , all of which have already been reported in @xcite .
however , the and lines , unresolved in @xcite , are clearly distinguishable in our spectra ; particularly for subregions xa+ne and carrot+w+sw .
the spectra were fitted to calculate line intensities using the same method as for the emission - line maps .
the line intensities were reddening - corrected by adopting the extinction curve of @xcite with @xmath20 = 3.1 , and @xmath21 = 0.08 @xcite . although the extinction curve was provided only for wavelengths @xmath131000  , we extrapolated from the curve down to the wavelength of @xmath0977 . in this short wavelength domain ,
other known extinction curves are also highly uncertain . assuming a distance of 540 pc to the cygnus loop @xcite , the line luminosities were also calculated and are shown along with the values for the vela snr @xcite in table 1 .
the x - ray luminosity @xcite is also given for comparison .
the distance of 540 pc is slightly larger than that ( 440 pc ) adopted in @xcite , and is a more recent , revised value @xcite .
however , we note that there is some evidence that the distance to the cygnus loop may be larger than 540 pc @xcite .
the spear / fims , , and line luminosities of the cygnus loop have been newly obtained and the and line luminosities have been separately estimated in this study .
* figure 4 .
* line - ratio diagrams for a shock model and observed values .
the curves indicate the line - ratio values expected from mixture of shocks with velocity populations of 0.110 power - law indices .
the line ratios ( with 1 @xmath17 error bars ) observed for the whole region and each subregion indicated in figure 2(a ) are overplotted . the line ratios less than 1 @xmath17 error
are not displayed .
the values of the power - law indices to explain the observed / ratios in each subregion are denoted by numerical values on the left side in ( a ) , and the corresponding positions ( along with some integer values ) are marked on the curves in ( b)(d ) .
lccccccc @xmath0977 & @xmath91.64 & 3.44 @xmath22 1.00 & @xmath91.81 & 12.34 @xmath22 2.55 & 21.84 @xmath22 0.49 + @xmath11032 , 1038 & 2.95 @xmath22 1.35 & 4.68 @xmath22 1.35 & @xmath91.54 & 19.79 @xmath22 4.04 & 14.81 @xmath22 0.45 + @xmath11394 , 1403 & 0.27 @xmath22 0.05 & @xmath90.12 & @xmath90.07 & 0.41 @xmath22 0.14 & + @xmath01404 & 0.45 @xmath22 0.06 & 0.67 @xmath22 0.06 & 0.10 @xmath22 0.03 & 2.10 @xmath22 0.14 & 6.13 @xmath22 0.17 + @xmath01486 & 0.15 @xmath22 0.03 & 0.17 @xmath22 0.05 & 0.10 @xmath22 0.03 & 0.79 @xmath22 0.11 & 1.47 @xmath22 0.07 + @xmath11548 , 1551 & 1.23 @xmath22 0.07 & 1.90 @xmath22 0.08 & 0.36 @xmath22 0.04 & 6.73 @xmath22 0.17 & 20.28 @xmath22 0.12 + @xmath01640.5 & 0.31 @xmath22 0.05 & 0.48 @xmath22 0.05 & @xmath90.19 & 1.83 @xmath22 0.14 & 2.75 @xmath22 0.10 + @xmath11661 , 1666 & 0.51 @xmath22 0.06 & 0.58 @xmath22 0.07 & @xmath90.10 & 1.64 @xmath22 0.17 & 3.56 @xmath22 0.12 + x - ray & & & & 5.41 & 3.0 +    since the spear / fims data have limited spatial resolution , it would be impossible to estimate the velocity of an individual shock from the observed fuv line ratios . instead , we attempted to compare the ratios of fuv lines summed over each subregion with the values expected from mixture of radiative shocks with a power - law velocity population @xcite .
the line ratios predicted for each shock velocity were obtained using the shock+precursor model of the mappings iii code @xcite .
we used the models with a preshock density of 1.0 @xmath23 and a preshock magnetic field of 0.5 @xmath24 g , both of which do not significantly affect the results .
we also used the atomic abundance set `` dopita2005 '' ( with @xmath20.20.3 solar abundance for c , n , and o ; dopita et al .
2005 ) as defined in table 1 of @xcite to reflect the low abundance in the rim of the cygnus loop @xcite .
assuming a mixture of shocks with the velocity population p(v ) @xmath25 v@xmath26 , where shock velocities v are in the range of 100250 km s@xmath11 with an interval of 25 km s@xmath11 , we tried to determine the power - law indices ( @xmath27 ) , which match the observed line ratios . here ,
the maximum shock velocity was selected to exclude non - radiative shocks . with a higher power - law index ,
the mixture of shocks contains a relatively large fraction of slower shocks . in figure 4
, we plot the line - ratio curves obtained by varying the power - law indices in the 0.110 range .
the spectra extracted from the whole region and the subregions indicated in figure 2(a ) , were used to estimate each line ratio and the resulting line ratios with 1 @xmath17 error bars are overplotted .
we present only the ratios for the , , , , and lines .
the line was excluded because of high uncertainty regarding the reddening - correction , and the and lines were too weak to be reliable in most subregions . in figure 4(a ) , the to ratio for subregion ne is well matched with the highest power - law index ( @xmath137.4 regarding 1 @xmath17 error ) ; indicating a relatively larger fraction of slow shocks in the subregion than in the other subregions .
subregions w and carrot have the power - law indices of @xmath23.6 and @xmath22.7 , respectively . the lowest index ( @xmath91.6 regarding 1 @xmath17 error ) , meaning a large fraction of fast shocks , is found in subregion xa .
the whole region has the intermediate value of @xmath23.0 , which seems to be a reasonable mixture of the subregions .
however , the values of the power - law indices obtained from the to ratio do not coincide well with the results obtained with the , , or lines , as can be seen in figures 4(b)(d ) .
we will discuss the implications of these results in the next section .
using the s - channel data of spear / fims , we newly created the @xmath0977 and @xmath11032 , 1038 emission - line images of the cygnus loop in figures 1(a ) and ( b ) .
these images can be compared with the two previous results : the _ voyager 2 _ and images @xcite , and the hires contour map @xcite .
the spear / fims image has the brightest region along the northwestern side , and less bright regions on the eastern and northeastern sides .
the brightest region coincides well with the brightest one in the _ voyager 2 _ image . on the other hand
, there is no feature corresponding to the eastern bright patch of the _ voyager 2 _ image .
instead , a very diffuse feature is found in our image in the vicinity of the bright patch .
this difference could be due to a large smoothing scale in the eastern region of our image , as mentioned in section 3 , and/or the slit sampling effects of the _ voyager 2 _ data .
the spear / fims image shows bright regions along the eastern side and the rims from the north to the southwest . except for a faint region on the northeastern side and a somewhat diffuse , bright region on the eastern side ,
overall bright features agree well with those of the _ voyager 2 _ image .
the morphology of the image can also be compared with the hires contour map in more detail .
a local peak on the western edge , and the eastern and northwestern enhanced regions , coincide with the hires contour peaks .
however , another local peak seen on the northern edge has no hires counterpart , although it seems to be near the northern hot spot of the _ voyager 2 _ image noted in @xcite . with a slightly better resolution ,
the spear / fims and images are not affected by the slit sampling effects of the _ voyager 2 _ data .
our image shows better correlation with other wavelength images than the hires contour map does .
this will be further discussed later .
we note that our new , , and images in figures 1(f)(h ) are not exactly the same as those in @xcite . these differences are mainly due to differences in the data set and the line - fitting methods ( e.g. , a linear versus a constant continuum , and the intensity ratios of doublet lines ) .
most of our results for subregions indicated in figure 2(a ) are similar to those of @xcite , and thus we discuss only the distinctions and new results .
@xcite reported the detection of and originating from the northeastern , non - radiative region .
however , selecting a region with a slightly higher declination than the ne non - radiative region in @xcite , as shown in figure 2(a ) , results in the and detections with a significance @xmath92 @xmath17 : just 1 @xmath17 upper limits of @xmath90.42 and @xmath91.31 @xmath4 10@xmath10 erg s@xmath11 @xmath12 sr@xmath11 , respectively . according to figure 7 in @xcite , there are incomplete radiative shocks located just below the southern border of subregion non - radiative marked in figure 2(a ) . considering the spatial resolution of the spear / fims data , most of the and detections reported in @xcite may not originate from non - radiative shocks but from incomplete radiative shocks .
even the line , which is known to be emitted more strongly from non - radiative shocks , has not been detected in subregion non - radiative .
@xcite predicted the intensity of 1.212 @xmath4 10@xmath10 erg s@xmath11 @xmath12 sr@xmath11 from a sharp northeastern non - radiative shock .
if this intrinsic brightness is obscured by interstellar dust , then the intensity would become lower and it is likely impossible to detect by the spear / fims instrument with the limited spatial resolution .
subregions xa and ne are the brightest regions for all of the spear / fims fuv lines , except for the s - channel lines ( and ) , which have worse statistics there .
however , their peak positions show apparent disagreement , as evident in figure 1 .
the , , and lines are brightest in subregion xa , whereas the and lines are found to be stronger in subregion ne . the line separated from the line peaks around the middle point between the two subregions .
the newly observed and lines have been identified in subregions carrot , w , and breakout , as well .    for five 40@xmath7 diameter regions
, @xcite presented fuv images with the spatial resolution comparable to those of the optical h@xmath6 and [ ] @xmath05007 images .
their fuv images show very sharp features coinciding well with bright optical filaments , although they have wide bandpass including the , , , and lines . therefore , despite the coarse spatial resolutions of the fuv images , strong correlations between the fuv images and optical images in figures 2(a)(d ) , imply that most of fuv emissions seen in the figure likely originate from the shocks consisting of bright h@xmath6 and [ ] filaments . assuming this
, we will compare the fuv images with the optical images . since the resonance line is formed at roughly the same temperature as the [ ] @xmath05007 forbidden line @xcite , figure 2(b ) could directly provide clues to the regions where strong resonance - scattering occurs .
a vertical filamentary structure running from @xmath28 to @xmath29 at @xmath30 in the [ ] image seems to have no counterpart .
therefore , this region may be regarded as a candidate region of strong - resonance - scattering . however , the vertical feature is also absent in the other fuv line images in figure 1 , except for a very faint feature in the image .
therefore , the absence of the vertical feature in the image may not be due to the effect of resonance scattering , but instead be due to dust extinction by a foreground dust cloud . the optical , vertical filamentary feature
would not be extinguished even by the foreground dust cloud because of much weaker dust extinction in optical wavelengths .
this might be another case of the cloud displayed in figure 2(f ) , which has been considered to be on the front side of the cygnus loop , as noted in section 3 .
since the , , and [ ] lines are all free from resonance scattering , and their relative strengths are not affected by elemental abundances , the comparison between figures 2(c ) and ( d ) can be used to verify spatial variations of shock - velocity populations .
the obvious variations of the and intensities are found on two filamentary [ ] features : one running along subregions xa and ne , and the other running along subregions w and sw . on the former feature ,
the intensity is highest around the northern part of subregion xa , whereas the intensity becomes higher toward both ends of the filamentary feature . on the latter feature ,
the intensity becomes higher toward both ends of the filamentary feature , but the intensity is highest near the center of subregion w. in the shock models of figure 4(a ) , the line ratio / increases with the relative fraction of fast shocks . it should be noted that the positions with higher / on each filamentary feature are closer to the locations of bright x - ray peaks seen in figure 2(e ) .
this is compatible with the above argument about shock - velocity populations .
as mentioned above , figure 2(e ) shows close coincidence between the spear / fims and the x - ray peaks ; particularly the northern and the western peaks .
the bright region on the eastern side of the image is well correlated with the overlapped x - ray peaks there .
an extended , fainter region on the southwestern side also overlaps well with the v - shaped x - ray feature @xcite .
the similarity between the and x - ray images indicates that the observed emissions are mostly originating from the x - ray emitting gas .
the only enhanced region without bright x - ray counterparts is located on the northwestern side , which overlaps with the northern part of the optical features in subregion carrot . to estimate the amount of the contribution by x - ray emitting gas , we compared the observed intensity with the value expected from the observed x - ray emissions .
@xcite calculated the intensities predicted at the northern x - ray peak ( their field a ) using a few models .
their ei + nei model predicts a similar intensity ( 59 @xmath4 10@xmath10 erg s@xmath11 @xmath12 sr@xmath11 ) as the average intensity of @xmath250 @xmath4 10@xmath10 erg s@xmath11 @xmath12 sr@xmath11 observed at the northern peak in figure 1(b ) .
however , their two - ei model predicts only 8% of the observed intensity .
these indicate that most of the emissions near the northern peak can be interpreted to originate from the x - ray emitting gas of which the softer x - ray component is in the condition of equilibrium ionization , and the harder one is in the non - equilibrium ionization condition ( ei + nei model of @xcite ) .
the fuv line luminosities in table 1 are now compared with the previous results . for the and lines ,
we compare the observed intensities , which were not reddening - corrected , because of high uncertainty in the reddening - correction .
@xcite reported the average observed intensities of 4.0 @xmath4 10@xmath10 erg s@xmath11 @xmath12 sr@xmath11 for and 8.8 @xmath4 10@xmath10 erg s@xmath11 @xmath12 sr@xmath11 for . on the other hand , @xcite obtained the intensity of 3.0(@xmath221.2 ) @xmath4 10@xmath10 erg s@xmath11 @xmath12 sr@xmath11 .
fitting the spectrum for the whole region in figure 3 results in the observed intensities of 2.22(@xmath220.46 ) @xmath4 10@xmath10 erg s@xmath11 @xmath12 sr@xmath11 for and 4.22(@xmath220.86 ) @xmath4 10@xmath10 erg s@xmath11 @xmath12 sr@xmath11 for .
these values are roughly half the _ voyager 2 _ values @xcite .
however , @xcite have recently recalibrated the _ voyager 2 _ intensities and found the old values to be too high by a factor of 2 , meaning that our values are similar to the corrected _ voyager 2 _ values .
the hires value @xcite also agrees well with our intensity within 1 @xmath17 error .
we note that the spear / fims systematic errors for the sensitivity calibration were estimated to be @xmath225% @xcite .
if the same distance as in @xcite is adopted , the present luminosity for the whole region is almost exactly the same as that of @xcite .
the and luminosities , however , are @xmath27080% higher than those of @xcite . the difference in the line - fitting methods , as mentioned above , likely caused the larger gaps for the weaker emission lines .
@xcite assumed a fixed ratio of 1:5.5 for the blended + lines to calculate the to ratio .
this caused a somewhat higher line ratio in subregion xa than those in the other parts of the cygnus loop .
@xcite discussed that the to ratio in subregion xa would in fact , be lower than 5.5 , based on the result of @xcite .
the decomposed line luminosities of and in table 1 verify their argument .
the to ratio for the whole region is @xmath25.1 , whereas the ratio for subregion xa+ne is @xmath21.7 .
fitting the spectrum extracted from only subregion xa , results in the ratio of @xmath22.6 .
the relative increase of the line in subregion xa+ne can be recognized in the spectra of figure 3 .    in estimating shock velocities from comparisons between the observed fuv emission lines and radiative shock models
, we should consider the effects of resonance scattering , x - ray emitting gas , and non - radiative shocks @xcite .
the to ratio is not affected by the above effects .
therefore , the fact that the values of the power - law indices obtained from the to ratio in figure 4(a ) are not derived well from the other line ratios , implies that the above effects are indeed important in the line ratios containing the , , and lines .
we , therefore , can infer the importance of the above effects in figures 4(b)(d ) . in figure 4(b ) , the observed ratios / and / are found to be much lower than the values expected from the ratio / for most subregions .
this indicates that the intensities were significantly diminished by the resonance scattering , which is ubiquitous across the cygnus loop . the resonance scattering in subregions ne and w seems to be particularly strong ( reduced by a factor of @xmath133 and @xmath132 , respectively ) . within our wide subregions ,
overlapping projections of the less edge - on shocks may cause a significant contribution from diffuse emission , in which the resonance scattering effect is less important . the contribution would depend on the nature of the ripples in the shock front .
however , more numerous optical filaments in subregions ne and w than the other subregions , indicate the significance of the resonance scattering in subregions ne and w. a kinematic study of the whole cygnus loop region showed that the radial velocities around subregions ne and w are similar to the systemic velocity of the cygnus loop @xcite .
even the ratios for the whole region are reduced to @xmath265% of the expected values .
since the resonant scattering does not destroy photons , the decrease of the line ratios obtained from the whole region implies that the cygnus loop would neither be perfectly spherical nor be homogeneous .
another possibility would be the depletion onto dust grains .
a global ir study of the cygnus loop showed that there is plenty of dust correlated well with the x - ray and optical emission of the remnant @xcite .
also , some x - ray studies suggested that a possible explanation of the low abundance in the overall rim of the cygnus loop could be the dust depletion @xcite .
the difference in the degree of dust depletion for carbon and oxygen may affect the ratios / and /.
as @xcite presented clear evidence for resonance scattering of in a northeastern non - radiative shock , resonance scattering of can also be widespread throughout the cygnus loop .
however , figure 4(c ) does not show any decrease in the line ratio for .
instead , the line ratio seems to increase for some subregions , although their errors are somewhat large .
the ratios for the whole region are higher by a factor of @xmath132 than the values matching with the ratio /.
these enhancements can be explained by the above conclusion that large portions of the observed emissions would be due to the x - ray emitting gas . moreover , since the x - ray emitting gas fills more diffuse and larger volumes inside the radiative shocks @xcite , the emissions from the x - ray gas are less affected by resonance scattering .
figure 4(d ) for the line shows no clear tendency . only in subregion
n , ratios are found to be slightly higher than those expected .
this may be due to a relatively large number of non - radiative shocks within subregion n because some portions of the intensity can originate from non - radiative shocks @xcite .    in table 1 , the luminosities of the x - ray and the highest ionization fuv line ( ) for the cygnus loop are @xmath280% and @xmath230% higher than those for the vela snr , respectively . on the other hand ,
the luminosities of the other fuv lines for the cygnus loop are all @xmath23070% lower than those for the vela snr .
we note that the ratios of the x - ray and luminosities to those of the lower ionization lines are higher in the cygnus loop than in the vela snr . although the distances to the cygnus loop was recently found to be larger than the value adopted in the present paper @xcite , the luminosity ratios are independent of the distance .
therefore , our result suggests that the cygnus loop must possess a much larger fraction of shocks in the early evolutionary stages than the vela snr .
we conclude that the cygnus loop is globally younger than the vela snr from the point of view of shock evolution populations .
in fact , the vela snr exploded near the border of the @xmath31 velorum stellar wind bubble and large portions of the remnant have already passed through the dense wall of the bubble @xcite .
in contrast , the cygnus loop exploded near the center of a cavity formed by its progenitor star and it has recently encountered the cavity walls @xcite .
utilizing the short wavelength channel data of spear / fims , we have obtained new spectral images of the cygnus loop at @xmath0977 and @xmath11032 , 1038 .
these images supplement the previously reported images with the slit sampling effects or a little worse resolution .
the observed intensities agree with the hires value within 1 @xmath17 error .
they are also consistent with the recently recalibrated _ voyager 2 _ values .
in addition , the + lines unresolved in @xcite have been successfully decomposed .
a new image for the weak @xmath01486 line was also obtained by binning and fitting the spectra carefully with the newly processed data .
the deblended and line maps confirm a region with a relatively low to ratio , as predicted in @xcite .
we have compared the spear / fims fuv images with the x - ray , optical , and 21 cm images .
the very close coincidence between the and x - ray enhanced regions suggests that the observed emissions are most likely related to x - ray emitting gas . for the region around the northern peak , comparison of the observed intensity with the value predicted from the observed x - ray emissions indicates that most of the emissions observed there can be produced from the x - ray emitting gas .
comparisons with the optical h@xmath6 and [ ] @xmath05007 images have revealed clear variations of shock - velocity populations on the two filamentary features .
we also argue that the fuv lines may be strongly extinguished in the direction of the previously identified cloud , which would be located on the front side of the cygnus loop .
we have calculated the fuv line ratios for several subregions of the cygnus loop , and compared them with a radiative shock model . from this , we have diagnosed subregions containing relatively large fractions of slow or fast radiative shocks , and also obtained clues to the effects of resonance scattering , x - ray emitting gas , and non - radiative shocks . in particular , it has been confirmed that the resonance scattering in the regions viewed edge - on is much higher than in the other regions .
the x - ray and fuv luminosities of the cygnus loop calculated from a revised distance have been compared with those of the vela snr . in the cygnus loop ,
the x - ray and luminosities relative to those of the other lower ionization fuv lines are higher than the vela snr .
this suggests that the cygnus loop contains a much larger fraction of shocks in the early evolutionary stages than the vela snr .
spear / fims is a joint project of the korea astronomy and space science institute , the korea advanced institute of science and technology , and the university of california at berkeley , funded by the korea most and nasa grant nag5 - 5355 .
kim was supported by a national research foundation of korea grant funded by the korean government .
we thank the referee for comments that improved the manuscript .
allen , m. g. , groves , b. a. , dopita , m. a. , sutherland , r. s. , & kewley , l. j. 2008 , , 178 , 20 arendt , r. g. , dwek , e. , & leisawitz d. 1992 , , 400 , 562 aschenbach , b. , & leahy , d. a. 1999 , , 341 , 602 blair , w. p. , long , k. s. , vancura , o. , & holberg , j. b. 1991a , , 374 , 202 blair , w. p. , long , k. s. , vancura , o. , et al .
1991b , , 379 , l33 blair , w. p. , sankrit , r. , & raymond , j. c. 2005 , , 129 , 2268 blair , w. p. , sankrit , r. , raymond , j. c. , & long , k. s. 1999 , , 118 , 942 blair , w. p. , sankrit , r. , & tulin , s. 2002 , , 140 , 367 braun , r. , & strom , r. g. 1986 , , 164 , 208 cardelli , j. a. , clayton , g. c. , & mathis , j. s. 1989 , , 345 , 245 charles , p. a. , kahn , s. m. , & mckee , c. f. 1985 , , 295 , 456 cornett , r. h. , jenkins , e. b. , bohlin , r. c. , et al .
1992 , , 395 , l9 danforth , c. w. , blair , w. p. , & raymond , j. c. 2001 , , 122 , 938 danforth , c. w. , cornett , r. h. , levenson , n. a. , blair , w. p. , & stecher , t. p. 2000 , , 119
, 2319 decourchelle , a. , sauvageot , j. l. , ballet , j. , & aschenbach , b. 1997 , , 326 , 811 dopita , m. a. , groves , b. a. , fischera , j. , et al .
2005 , , 619 , 755 edelstein , j. , korpela , e. j. , adolfo , j. , et al .
2006a , , 644 , l159 edelstein , j. , min , k .- w .
, han , w. , et al . 2006b , , 644 , l153 grski , k. m. , hivon , e. , channelay , a. j. , et al .
2005 , , 622 , 759 green , d. a. 1990 , , 100 , 1927 hester , j. j. , raymond , j. c. , & blair , w. p. 1994 , , 420 , 721 katagiri , h. , tibaldo , l. , ballet , j. , et al .
2011 , , 741 , 44 katsuda , s. , tsunemi , h. , kimura , m. , & mori , k. 2008 , , 680 , 1198 katsuda , s. , tsunemi , h. , mori , k. , et al .
2011 , , 730 , 24 kim , i .- j . , seon , k .-
i . , min , k .- w . , han , w. , & edelstein , j. 2012 , , 761 , 135 ku , w. h .- m . ,
kahn , s. m. , & pisarski , r. 1984 , , 278 , 615 leahy , d. a. 2002 , , 123 , 2689 leahy , d. a. 2003 , , 586 , 224 leahy , d. a. 2004 , , 351 , 385 leahy , d. a. 2005 , , 130 , 165 leahy , d. a. , roger , r. s. , & ballantyne , d. 1997 , , 114 , 2081 levenson , n. a. , & graham , j. r. 2001 , , 559 , 948 levenson , n. a. , & graham , j. r. 2005 , , 622 , 366 levenson , n. a. , graham , j. r. , aschenbach , b. , et al .
1997 , , 484 , 304 levenson , n. a. , graham , j. r. , keller , l. d. , & richter , m. j. 1998 , , 118 , 541 levenson , n. a. , graham , j. r. , & snowden , s. l. 1999 , , 526 , 874 levenson , n. a. , graham , j. r. , & walters , j. l. 2002 , , 576 , 798 long , k. s. , blair , w. p. , vancura , o. , et al .
1992 , , 400 , 214 mcglynn , t. , scollick , k. , & white , n. 1998 , in iau symp .
179 , new horizons from multi - wavelength sky surveys , ed . b. j. mclean et al .
( dordrecht : kluwer ) , 465 miller , j. s. 1974 , , 189 , 239 miyata , e. , tsunemi , h. , kohmura , t. , suzuki , s. , & kumagai , s. 1998 , , 50 , 257 quemerais , e. , snow , m. , & bonnet , r .- m . 2013 , cross - calibration of far uv spectra of solar system objects and the heliosphere ( new york : springer ) rasmussen , a. , & martin , c. 1992 , , 396 , l103 raymond , j. c. , black , j. h. , dupree , a. k. , hartmann , l. & wolff , r. s. 1980a , , 238 , 881 raymond , j. c. , black , j. h. , dupree , a. k. , hartmann , l. & wolff , r. s. 1981 , , 246 , 100 raymond , j. c. , blair , w. p. , fesen , r. a. , & gull , t. r. 1983 , , 275 , 636 raymond , j. c. , blair , w. p. , long , k. s. , et al .
1997 , , 482 , 881 raymond , j. c. , davis , m. , gull , t. r. , & parker , r. a. r. 1980b , , 238 , l21 raymond , j. c. , ghavamian , p. , sankrit , r. , blair , w. p. , & curiel , s. 2003 , , 584 , 770 raymond , j. c. , hester , j. j. , cox , d. , et al . 1988 , , 324 , 869 salvesen , g. , raymond , j. c. , & edgar , r. j. 2009 , , 702 , 327 sankrit , r. , & blair , w. p. 2002
, , 565 , 297 sankrit , r. , blair , w. p. , cheng , j. y. , et al .
2007 , , 133 , 1383 sankrit , r. , blair , w. p. , raymond , j. c. , & long , k. s. 2000 , , 120 , 1925 sankrit , r. , williams , b. j. , borkowski , k. j. , et al .
2010 , , 712 , 1092 seon , k .- i . , edelstein , j. , korpela , e. , et al .
2011 , , 196 , 15 seon , k .-
han , w. , nam , u .- w . ,
2006 , , 644 , l175 shull , p. j. , & hippelein , h. 1991 , , 383 , 714 sun , x. h. , reich , w. , han , j. l. , reich , p. , & wielebinski , r. 2006 , , 447 , 937 sushch , i. , hnatyk , b. , & neronov a. 2011 , , 525 , a154 tsunemi , h. , katsuda , s. , nemes , n. & miller , e. d. 2007 , , 671 , 1717 uchida , h. , tsunemi , h. , katsuda , s. , & kimura , m. 2008 , , 688 , 1102 uchida , h. , tsunemi , h. , katsuda , s. , et al .
2009 , , 61 , 503 uyaniker , b. , reich , w. , yar , a. , kothes , r. , & frst , e. 2002 , , 389 , l61 vancura , o. , blair , w. p. , long , k. s. , & raymond , j. c. 1992 , , 394 , 158