in 1970 , efimov described a phenomenon occurring in a three - particle system of identical bosons , in which there is an accumulation of bound states near the point of zero energy and diverging scattering length of the interparticle interaction @xcite .
efimov was inspired by the physics of the nucleus of the tritium atom .
however , the experimental signature for the existence of such efimov states was not observed in nuclear physics , but in a gas of ultracold cesium atoms  @xcite by studying three - body recombination .
the reason that properties predicted for nuclear physics few - body systems can be derived from ultracold atoms lies in the fact that the two systems are connected through universal relations , that are irrespective of the specific nature of the interparticle interactions .
since the first observations in cesium , the efimov effect has received a lot of attention  @xcite .
several groups now have observed efimov physics by studying three - body recombination in a variety of atomic systems  @xcite .
evidence for the existence of trimers appears as some distinct features ( resonances ) in the recombination rate coefficient . in fig .
[ fi : efimovschematic ] the energies of the efimov trimers are shown , as a function of the inverse two - body scattering length @xmath0 .
a crucial advantage of ultracold atoms with respect to nuclear systems is that the interparticle interactions can be easily manipulated using feshbach resonances  @xcite , by varying the magnetic field .
the two - body scattering length , which is obtained from the scattering phase shift @xmath1 as @xmath2 , diverges on resonance to infinity , and is typically described by the dispersive formula    @xmath3    with @xmath4 the background value of the scattering length , @xmath5 the magnetic field , @xmath6 the resonance position , and @xmath7 the width of the resonance .
this possibility to change the scattering length allowed for the observation of distinct efimov features in the three - body recombination rate , connected to two different efimov states in a row @xcite , confirming the scaling law   given by    @xmath8    with @xmath9 a constant .
moreover , recent developments extended the universal few - body physics to the domain of four - body states which was theoretically predicted  @xcite and experimentally verified  @xcite .    in most of the experiments ,
efimov physics has been observed only indirectly via three - body recombination  . due to the short lifetime of these states , it is difficult to probe the efimov states directly .
these weakly - bound trimers are not only unstable in collisions with other particles , but they are also intrinsically unstable by itself , as trimers may decay to a deeper lying dimer state and an unbound atom . however , the group of jochim in heidelberg recently succeeded in associating efimov trimers directly via a new and promising approach using radio - frequency fields  @xcite
.    optical lattices might provide more experimental control , and potentially give rise to an enhanced lifetime of efimov trimers .
efimov physics can be studied in this way by loading three atoms in each lattice site .
when the lattice is deep enough , single sites act as harmonic oscillator potentials where individual trimers can be trapped .
a possible advantage is then that the interparticle separation of the trimers can be increased to the scale of the harmonic oscillator length , while their properties are still determined by efimov physics . in this way
the overlap of these states with deeply bound dimers is reduced , thereby suppressing the process of three - body recombination .
the influence of other limiting processes , where four or more particles involved , will be reduced as well .
, indicating dimer and efimov trimer states as a function of the inverse scattering length @xmath0 .
the horizontal line @xmath10 indicates the threshold to the three - body scattering continuum ( gray ) .
simultaneously , the blue region marks the atom - dimer continuum for @xmath11 , and the green zone indicates the trimer region .
the dimer binding energy scales as @xmath12 . only a few efimov states are shown .
the figure is not to scale : for clarity reasons we plotted instead of @xmath13 the quantity @xmath14 , and instead of @xmath0 the quantity @xmath15 ( cf .
eq .  ( [ eq : polarcord])).,width=302 ]    jonsell , heiselberg and pethick were the first to introduce an harmonic oscillator term in a three - body system  @xcite , by adapting the approach of ref .  @xcite . for diverging scattering length , they found exact expressions for the bound states of three identical bosons , while they obtained numerical results for finite scattering length . later , werner and castin @xcite obtained additional analytical results for diverging scattering length , and complimentary first - order corrections for large scattering length .
stoll and khler @xcite performed numerical simulations with a tunable separable potential to find the bound states of three particles in a trap .
et al . _  found numerical solutions of three particles in a harmonic trap  @xcite , mainly to study the non - universal corrections to the efimov effect that come from the effective range .
more recently , liu , hu and drummond derived an analytical description of three attractively interacting fermions in a harmonic trap  @xcite .
a convenient description of the efimov effect in free space is presented in ref .
, in which polar coordinates @xmath16 and @xmath17 are introduced , which are related to the wave number or energy parameter @xmath13 and the scattering length @xmath18 via    [ eq : polarcord ] @xmath19    using efimov s radial law , it is derived that @xmath20 for a function @xmath21 \to \mathbb{r}$ ] .
here @xmath22 is the three - body parameter that defines the trimer energies exactly on resonance as @xmath23 , with @xmath24 the mass of a single atom , and is a result of the short - range behavior of three particles being close together .
the combinations of @xmath13 and @xmath17 that yield bound states describe a family of curves that can be mapped on top of each other by multiplying the distance of the curves to the origin by a factor @xmath25 , see fig .
[ fi : efimovschematic ] .    in this paper
, we extend the description of ref .   and show how by introducing the harmonic oscillator length as a new parameter , and by using spherical coordinates in parameter space , the scaling behavior in free space very naturally generalizes to a situation with three particles in a harmonic trap . in particular , the curves that are shown in fig .
[ fi : efimovschematic ] will be replaced by surfaces that are mapped on top of each other . in the derivation ,
we use an approach slightly different from earlier work , that calculates the equivalent of the function @xmath26 by solving an ordinary differential equation . additionally , we also include the coupling terms , we investigate the stability of the trimers in the lattice , and we study the non - universal corrections coming from the effective range parameter .
we consider a system of three bosons with pair interactions trapped in a harmonic potential .
the corresponding hamiltonian is    @xmath27    where @xmath28 is the harmonic oscillator length , @xmath29 the trap frequency , and @xmath30 is the interaction between the particles @xmath31 , with @xmath32 a cyclic permutation of @xmath33 . for convenience ,
the hamiltonian is scaled by a factor @xmath34 . throughout this paper
, we will denote the eigenvalues of @xmath16 by @xmath35 , and the corresponding trimer energies by @xmath36 . to find the bound state energies , we closely follow the method described in ref .
@xcite and ref . .
we denote the permutation operator corresponding to one cyclic permutation of @xmath33 by @xmath37 and transform to the faddeev operator @xmath38 if @xmath39 is an eigenfunction of @xmath40 , @xmath41 is an eigenfunction of @xmath16 , provided it is nonzero . in case @xmath42 , @xmath39 is called a spurious solution .
we separate the center of mass motion , and consider the relative system in hyperspherical coordinates . for that
, we introduce jacobi - coordinates according to @xmath43 the hyperradius @xmath44 and the hyperangle @xmath45 are defined by @xmath46 the remaining hyperspherical coordinates are given by the unit vectors @xmath47 and @xmath48 .
we restrict ourselves to states with zero total angular momentum .
the advantage of considering the faddeev rather than the schrdinger operator is that subsystem angular momentum comes in only at second order in the short - range potential @xcite .
this motivates us to neglect the subsystem angular momentum , by averaging over the jacobi coordinates @xmath47 and @xmath48 . in this approximation ,
the schrdinger wavefunction @xmath49 is given in terms of reduced faddeev wavefunctions @xmath50 :    @xmath51    in the differential faddeev equation   we then make an adiabatic hyperspherical expansion    @xmath52    for every value @xmath44 , which is treated as a parameter , the functions @xmath53 of hyperangle @xmath45 form a complete set , and are solutions to the eigenvalue equation on @xmath54 @xmath55\phi_n(r,\alpha ) = - 2 r^2 v(\sqrt{2 } r \sin \alpha )   \left [ \phi_n(r,\alpha ) + \frac{4}{\sqrt{3 } } \int_{|\frac{1}{3}\pi-\alpha|}^{\frac{1}{2}\pi - |\frac{1}{6}\pi-\alpha| } \phi_n(r,\alpha ' ) d \alpha ' \right],\ ] ] with dirichlet boundary conditions @xmath56 .
the eigenvalues @xmath57 depend parametrically on @xmath44 .
the set of hyperangular functions @xmath53 are generally not orthogonal functions in @xmath45 .
therefore , we define @xmath58 which allows us to derive the following coupled set of differential equations for the functions @xmath59
@xmath60 f_n(r ) + \\+ \sum_m \left[2 p_{nm}(r ) \frac{\partial}{\partial r } + q_{nm}(r)\right]f_m(r ) = e f_n(r).\end{gathered}\ ] ]    the coupling terms in these equations are defined by    [ eq : couplingterms ] @xmath61    this framework is very similar to the work of refs .  , and the harmonic oscillator term is taken into account as in ref .  @xcite . for
now , we neglect the coupling terms conform the adiabatic hyperspherical approximation .
however , later on we wil study the coupling between different channels more carefully , as the harmonic oscillator gives rise to degeneracies between the efimov - like solutions in the zeroth channel and the continuum - like solutions of the first and higher channels .
+     + the effect of the interaction potential @xmath62 is now encoded in the eigenfunctions @xmath57 .
parametric solutions for @xmath63 can be found by fixing the parameter @xmath44 , while solving for the angular problem .
the functions @xmath57 will then serve as input for the radial equation ( [ eq : coupledfn ] ) , and when the coupling terms are ignored , solutions are obtained by solving an ordinary differential equation ( ode ) .    to find an approximation of @xmath63 , one can study solutions @xmath50 of eq .
( [ eq : eqhypwf ] ) for small values of @xmath45 and values of @xmath45 close enough to @xmath64 .
matching of the solutions at intermediate values of @xmath45 then yields an approximation of @xmath63 .
this procedure has been described in ref .
the equation found in this manner is @xmath65 where @xmath66 .
note that @xmath67 always solves this equation , but it gives rise to a spurious solution .
the right hand side of equation ( [ eq : approxlambda ] ) can be conveniently expanded using the effective range approximation for the s - wave scattering phase shift of the two - body problem @xmath68 with @xmath69 the effective range parameter . when this parameter is neglected , we find the approximating equation for @xmath70 that is at the basis of the derivation in ref .
@xcite : @xmath71 by taking into account the effective range as well , one obtains an approximation that is more accurate with respect to the two - body physics near narrow feshbach resonances .
this has been proposed in ref .
@xcite as a regularization procedure and has been readdressed in ref .
we will come back to this issue in section [ sec : effrange ] .
+   + with the approximation of the channel eigenvalues at hand , we may turn to the hyperradial problem . in the adiabatic hyperspherical approximation
, we ignore the coupling between the different hyperradial channels and only take the scattering length into account in the approximation of @xmath63 . in this setting ,
the efimov effect was derived in ref .
@xcite , and repeating it in this context allows for a comparison in the limit of @xmath72 .
the uncoupled system of equations is given by    @xmath73f_n(r)= 2 e f_n(r),\ ] ]    on @xmath74 , with a boundary condition @xmath75 and an appropriate boundary condition on the edge of the inner region @xmath76 where the short - range potentials become important , which is given in terms of the three - body parameter @xmath22 .
the exact boundary condition will be specified later .
we first concentrate on the zeroth channel , i.e.  @xmath77 , which is the only channel that has an attractive potential , and therefore the only channel in which the efimov effect is found .
we introduce spherical coordinates @xmath78 according to    [ eq : sphercor ] @xmath79    and write down the differential equation for @xmath80 as    [ subeq : problemforf ] @xmath81f(\sigma )   = 2\operatorname{sgn}{(\xi)}\sin^2 \vartheta \sin^2 \xi f(\sigma),\ ] ] with long - range boundary condition @xmath82 for small values of @xmath83 , the model does not provide us with a natural boundary condition .
solutions of eq .
( [ eq : diffeqforf ] ) have a log - periodic behavior as @xmath84 , given by @xmath85 for some value of @xmath86 . here
@xmath87 , and @xmath88    fortunately , the real short - range potentials @xmath30 provide us with a natural boundary condition at @xmath76 encoded by @xmath89 , where @xmath22 is the three - body parameter that accounts for the short - range physics , which was mentioned already in the introduction .
this gives rise to a short - range boundary condition at @xmath76 : @xmath90    , @xmath13 and @xmath91 ( blue ) .
only three surfaces are shown , but in fact there is an infinite number of them , which all can be mapped on top of each other by multiplying the distance to the origin by powers of @xmath25 .
these surfaces are sliced by a plane ( red transparant ) , for a given @xmath92 that corresponds to the height of the plane .
the intersection lines ( red ) represent bound state solutions for combinations of the scattering length @xmath18 and the wave vector @xmath13 , for this value of @xmath92 .
the picture is not too scale : we used similar scaling relations as used in fig .
[ fi : efimovschematic ] in order to show the surfaces more clearly.,width=302 ]    we may determine for which parameter combinations a solution exists by a shooting method from large values of @xmath83 .
a particular convenient method makes use of a certain phase variable @xmath93 .
one can set up an equation for @xmath93 that is much alike the variable phase equation , but is adapted to fit the boundary condition at short distances very well .
henceforth , it will play a role similar to the function @xmath26 as we mentioned in the introduction .
we introduce the transformation @xmath94 the functions @xmath95 , @xmath96 and @xmath97 may be chosen to our convenience .
continuity properties in the parameters are easier to derive considering @xmath93 than when one considers @xmath98 directly .
a particular choice @xmath99 is highly compatible with the boundary condition ( [ eq : bcsmallsigma ] ) .
the equation for @xmath93 becomes @xmath100 which is a nonlinear first order differential equation .
given values for @xmath17 and @xmath101 , we denote by @xmath102 the solution to ( [ eq : forphief ] ) that satisfies @xmath103 .
with this definition , the problem ( [ subeq : problemforf ] ) has a solution if and only if @xmath104 consequently , eigensolutions only exist if there exists an @xmath105 such that @xmath106 we would like to emphasize here the similarity with the description of the efimov effect in free space given by equation ( [ eq : descrefimovfreespace ] ) . using the asymptotic behavior of @xmath107 for small arguments
, one can easily show that @xmath108 converges to a finite value for @xmath109 .
the function @xmath110 can be shown to be continuous over the region . for @xmath111 and
@xmath112 $ ] , the efimov effect in free space is retained .
hence , @xmath113 is an approximation of @xmath114 .
+   + [ se : resultssimplemodel ]     s@xmath115 .
a deeper bound state , with much lower energy , can not be shown in this graph .
also indicated is the dimer binding energy ( solid blue line),width=302 ]    the trimer solutions are most easily represented by surfaces in three - dimensional @xmath116 space , as defined in eq .
( [ eq : sphercor ] ) . in fig .
[ fi : boundstatesurfaces ] such surfaces are shown for the lowest channel , which can be mapped on top of each other by multiplication of the distance to the origin by powers of @xmath25 . fig .
[ fi : boundstatesurfaces ] can be regarded as a generalization of the two - dimensional fig .
[ fi : efimovschematic ] .
when these surfaces are sliced with a plane corresponding to a fixed harmonic oscillator length @xmath92 , we obtain curves that express wave vectors @xmath13 in terms of the scattering length @xmath18 , and an example can be seen in fig .
[ fi : channelscombined ] .
when taking the limit of @xmath117 , the efimov effect in free space is obtained as depicted in fig .
[ fi : efimovschematic ] . in the region
@xmath112 $ ] , the intersections of the surfaces will tend to the curves that describe the efimov states . in the region @xmath118 ,
the energy levels come closer and closer together , expressing the fact that in the limit , this region is filled with continuum states .
the intersection of the surfaces with the @xmath119 plane precisely corresponds to the exact results obtained by jonsell _
et al . _ and werner  @xcite .
as an example , we show in fig .
[ fi : channelscombined ] the trimer binding energies of cesium as a function of scattering length with a fixed oscillator length , for the lowest and next to lowest hyperspherical channels .
the parameters are @xmath120  @xcite , with @xmath121 the bohr radius , @xmath122 , and @xmath123 , which corresponds to a trap frequency of approximately @xmath124 s@xmath115 .
it is interesting to note the change of energy scales in this figure . below the thresholds of the trimer and atom - dimer continuum ( cf .
[ fi : efimovschematic ] ) , which would exist for the free - space efimov effect , the bound states show the characteristic scaling law of eq .
[ efimovscaling ] . above these thresholds ,
the bound states show a spacing typical for the harmonic oscillator potential which is approximately @xmath125 .
this change of exponential scaling to linear scaling of the trimer energies makes that the accumulation of bound states in the point of zero energy and infinite scattering length disappears .
effectively , the efimov trimers are `` pushed '' out of the trimer region through the thresholds and aquire a strong harmonic oscillator state character .
are replaced by avoided crossings when the coupling between the zeroth and first hyperspherical channels is taken into account.,width=302 ]    there is also a clear difference between the zeroth and first hyperspherical channels . in this respect , it is interesting to note the different limiting behavior of the solutions for small positive and negative scattering length . for small and negative scattering length ,
there is no weakly - bound dimer state in the interatomic potential , and the solutions have the character of three weakly - interacting atoms in a harmonic oscillator potential . for small and positive scattering length ,
there is a weakly - bound dimer state , and for the zeroth hyperspherical channel the solutions have the character of a weakly - interacting dimer and an atom in a harmonic oscillator potential .
these regimes are connected through solutions which have a strong efimov character . however , for the first hyperspherical channel the situation is different . here
the solutions behave more as those of two particles undergoing a feshbach resonance  @xcite , which as a result interpolate between two different harmonic oscillator levels
.    it should be possible to investigate this change of character experimentally , by using similar techniques as have been done by stferle _
et al . _
@xcite . here
the creation of dimers was investigated in an optical lattice while sweeping the magnetic field through a feshbach resonance .
for the trimers , such adiabatic sweeps through trimer solutions should allow for a creation of states with a character changing between harmonic oscillator- , efimov- and dimer+atom- type solutions .
so far , we have applied the adiabatic hyperspherical approximation , in which we completely ignored the coupling between the various hyperspherical channels .
while this is generally a good approximation , we see from fig .  [ fi : channelscombined ] that there are degeneracies between the zeroth and first channel solutions , which should be lifted when the coupling is taking into account . in order to demonstrate this , we introduce an approximation of the coupling .
we calculate the coupling terms @xmath126 and @xmath127 of eq .
( [ eq : couplingterms ] ) using first order approximations of the functions @xmath128 , and find in this way combined solutions which involve both channels . in free - space , the coupling between these two hyperspherical channels is directly related to the process of three - body recombination , where three colliding particles relax into a weakly - bound dimer and a third particle , with a surplus of kinetic energy  .
interestingly , for our lattice solutions we observe that this coupling transforms the crossings in fig .
[ fi : channelscombined ] into avoided crossings , as shown in fig .
[ fi : avoidedcrossings ] . this could be experimentally interesting , since using magnetic field sweeps would allow for adiabatic transitions from one hyperspherical channel to an other .
of cesium efimov states in a trap , with @xmath129 s@xmath115 and @xmath130 , for the first ( red ) , second ( black ) , third ( green ) and fourth ( blue ) trimer states depicted in fig .
[ fi : channelscombined ] .
also , the width of a free efimov trimer corresponding to the first state is shown ( black dashed).,width=377 ]    the change of character of the trimer solutions , which was discussed above , will have a profound impact on the stability of these states .
so far , we implicitly assumed that the trimers in the harmonic oscillator potential are stable .
however , when the three particles are close together the trimer can dissociate into a deeply - bound dimer and a third atom , with large kinetic energy . here
, we investigate the stability of the trimers by applying the method of braaten and hammer  , that regards the trimers as scattering resonances in a dimer - atom scattering process . in this respect ,
the trimer energy @xmath131 acquires a width @xmath132 , where the latter can be used to calculate the lifetime @xmath133 of the trimers . in this method ,
an inelasticity parameter @xmath134 is introduced to account for the decay into various deeply bound states . which can be regarded as the imaginary part of a complex phase factor that describes the reflection of hyperradial waves in the inner region .
effectively , this comes down to changing the boundary condition at @xmath76 to @xmath135 for given values for the scattering length and the harmonic oscillator length , we are then able to find the associated complex eigenvalues    @xmath136    where the trimer energy @xmath137 and the width @xmath138 .    as an example
, we calculated the linewidths or cesium , with a trapping frequency @xmath139 s@xmath115 , @xmath140 , and @xmath130  @xcite .
the short - range boundary condition was implemented at a value @xmath141 .
for some of the efimov trimers shown in fig .
[ fi : channelscombined ] , the linewidths are shown in fig .
[ fi : stabplot ] .
the change of character of the trimers can be clearly seen from this graph : for values @xmath142 the states have a harmonic oscillator character , and the widths are relatively small .
however , the efimov effect shows up via a strong increase of the widths towards @xmath143 , which is the parameter region where in free - space a new efimov state would enter the tree - body system .
it is also the region where a local peak in the three - body recombination rate would occur ( trimer resonance ) . here
, we see that at this point the width of the lowest trimer keeps increasing , while for the higher trimer states the width saturates and displays an oscilating behavior .
this is not surprising , since the latter trimers are more similar to three - particle continuum states . for comparison ,
the linewidth of the free - space efimov trimer compared to the lowest lattice trimer state is shown in the region where it exists .
these lines overlap mostly , and can be more than two orders of magnitude above the widths for the higher trimer states . at the point where the trimer energies approach the atom - dimer threshold , around @xmath144 ,
the widths go down again to the level of the higher trimer states .
the width of the lattice trimer is there somewhat larger than the free - space trimer , probably due to the additional confinement of the weakly - bound dimer state caused by the lattice .    while the efimov effect shows up in the excited states by a clear deviation from the harmonic oscillator levels , these states benefit from a much smaller linewidth than the deeply bound efimov trimers , and they should therefore be easier to study experimentally . also the lowest trimer energy , which most closely follows a free - space efimov state , could be easier to study in experiment .
the trimer could first be prepared in the lowest harmonic oscillator level for @xmath145 , and then the scattering lenght can be changed via a magnetic field sweep to the region @xmath146 , where the decay is still manageable .
in order to improve the approximation of the channel eigenvalues @xmath63 arising in ( [ eq : eqhypwf ] ) , one might want to include the influence of the effective range @xmath69 .
this can be done by making the so - called effective range expansion , i.e.  using the expansion in eq .
( [ eq : effrangeexp ] ) up to second order . in literature
there have been various instances in which the effective range is taken into consideration in the determination of the three - particle spectrum @xcite .
having an additional interaction parameter at hand , raises the question whether in this regime efimov physics in ultracold gases can still be universal . on the experimental side ,
there are presently also some conflicting observations . in order to understand these issues ,
we first need to define what we mean by the universal efimov regime . in the above
, we have already seen that besides the scattering length , we need at least one additional parameter that fixes at least one of the efimov trimer energies .
the other ones can then be derived from this position .
this information is provided via the three - body parameter @xmath22 , which can be regarded as a boundary condition for the interaction potentials in the regime where the three particles are very close together .
it is related to the energy of the lowest efimov trimer via @xmath147 in the infinite scattering length limit , and with this only additional parameter the efimov physics is universal .     of @xmath148 depending on @xmath149 , as calculated in three different ways .
the blue lines represent the approximations of the eigenvalues by a truncated fourier expansion .
the black lines and red lines are the approximations obtained by ( [ eq : lambdaonlya ] ) and ( [ eq : lambdawitheffrange ] ) respectively .
we take for the square well parameters @xmath150 .
it can be clearly seen that the effective range gives a much better asymptotic behavior .
, title="fig:",width=264 ] +    however , some experiments on three - body recombination do not imply universality .
while the experiments at rice university did show the universal connection between efimov states in @xmath151li at one side of the resonance  @xcite , the universality did not extend to the other side of the resonance .
the rice group was not alone in this observation : also in florence  @xcite on @xmath152k , in innsbruck  @xcite on @xmath153cs , and in tokyo  @xcite on @xmath154li this non - universality was observed . on the other hand , in three other experiments universality across the resonance
was observed : in bar - ilan it was found that universality does extend across a feshbach resonance in @xmath151li , and even that the short range physics encapsulated in @xmath22 is independent of the nuclear spin configuration  @xcite .
furthermore universality was also observed in a three - component fermionic @xmath155li gas in heidelberg  @xcite and in state college  @xcite .
now the question arises whether the effective range paramer @xmath69 can explain , in addition to the three - body parameter @xmath22 , the non - universal behavior of some of the above - mentioned experiments .
however , it might also be that in a certain regime , the effective range will take over the role of the three - body parameter due to its regularising effect of the radial wavefunctions  @xcite . and
recently , it has been shown that small experimental deviations from universality  @xcite can be explained by effective range corrections using an effective field theory  @xcite .    as it has been shown by fedorov and jensen @xcite and more recently by platter _
@xcite , the introduction of the effective range in the approximate equation ( [ eq : approxlambda ] ) for the eigenvalues @xmath63 yields the relation    @xmath156    to get an idea on the validity of the approximation ( [ eq : lambdawitheffrange ] ) , we have compared the results with the eigenvalues obtained by a truncated fourier expansion for a square well potential with height @xmath157 and width @xmath76 , such that @xmath158 .
this expansion is very convenient when one considers the operator @xmath148 , which depends on the hyperradius @xmath44 and is defined via the eigenvalue equation ( [ eq : eqhypwf ] ) as follows @xmath159.\ ] ] the reason for this convenience is that in the basis of fourier expansion functions @xmath160 both the differential operator and the integral operator are diagonal .
the matrix elements are then given by @xmath161 where @xmath162 in fig .
[ fi : comparisoneffrange ] , we compare the numerically calculated eigenvalues @xmath57 of @xmath148 with results obtained from eq .
( [ eq : lambdaonlya ] ) , which depends on the scattering length only , and with the results from solving eq .
( [ eq : lambdawitheffrange ] ) , which depends both on scattering length and effective range . here
the values for @xmath18 and @xmath69 are taken from the effective range expansion of the square well potential .
it can be seen that asymptotically , that is for large hyperradius , the inclusion of the effective range in the model yields results that are considerably better than if one only takes the scattering length into account .
for relatively small values of the hyperradius , the predictions obtained by the fourier expansion and this simple model start to deviate .
hence , we expect that for short hyperradii , additional details of the interaction potential become important , i.e.  higher order expansion coefficients beyond the effective range approximation have to be taken into account .    for @xmath163 and a fixed value for @xmath18 , the lowest solution @xmath107 to eq .
( [ eq : lambdawitheffrange ] ) can be defined as a smooth function of @xmath44 for @xmath164 . moreover , @xmath165 as @xmath166 .
therefore , if we use this approximation in the hyperradial equation , we can find _ regular _ solutions instead of rapidly oscillating ones .
this could , at least in principle , remove the need for a renormalization boundary condition in terms of @xmath22 for @xmath166 . in the absence of a harmonic oscillator , and in the limits of small energy and @xmath167
, we will see that the solutions to the hyperradial equation are identical to those when the effective range is not taken into account , except that the value of @xmath22 is determined by @xmath69 by boundary layer behavior for small values of the hyperradius .
such a relation between @xmath22 and @xmath69 has been established before @xcite .
let us , therefore , study the equation ( [ eq : lambdawitheffrange ] ) more closely .
we have plotted the solution @xmath107 for the situation @xmath168 in fig .
[ fi : illlambda0 ] , to give a characteristic example . in the region
@xmath169 , the value of @xmath107 satisfying ( [ eq : lambdawitheffrange ] ) is close to that obtained by solving ( [ eq : lambdaonlya ] ) .
similarly , when @xmath170 , the solution @xmath107 to ( [ eq : lambdawitheffrange ] ) will be close to @xmath171 , where @xmath171 is the solution to    @xmath172    around @xmath173 , there is a transition from a boundary layer behavior similar to @xmath174 to behavior obtained by ignoring the effective range ( [ eq : lambdaonlya ] ) .
( purple , solid ) as a function of @xmath175 , for a fixed value of @xmath176 , with @xmath177 . for @xmath178 ,
the solution is very close to that obtained from ( [ eq : lambdaonlya ] ) ( red , dotted ) . for @xmath179 ,
there is a boundary layer , in which the solution resembles the solution @xmath174 of ( [ eq : mu ] ) ( blue , dashed).,width=302 ]    to find the low - energy behavior of the system , which is needed to find a proper relationship between @xmath69 and @xmath22 @xcite , we have to turn back to the hyperradial equation ( [ eq : hyperradial ] ) .
following from the above considerations , the solution @xmath180 that decays exponentially for large hyperradius is to leading order in @xmath181 , and for @xmath182 , equal to @xmath183 , the exponentially decaying solution to the hyperradial equation obtained when the effective range is ignored .
this outer solution @xmath183 will have log - periodic behavior as @xmath166 @xmath184 for some @xmath86 .
the solution @xmath180 will be similar for @xmath185 , and will therefore for low enough energies also show log - periodic behavior .
however , there is a boundary layer for @xmath170 , which makes that we can demand the solution @xmath186 to be regular for @xmath166 , that is , such that @xmath187 for @xmath166 .
this effectively determines the value of @xmath86 . to analyze the boundary layer
, we introduce the new variable @xmath188 and the function @xmath189 , @xmath190 . the hyperradial equation for @xmath189 reads @xmath191 \tilde{f}_0(\rho ) = 0,\ ] ] and we demand that @xmath192 behaves like @xmath193 for @xmath194 .
we consider now the case of a negligible harmonic oscillator potential , i.e. @xmath195 , and take the low - energy limit @xmath196 and a diverging scattering length .
the solution will then be equal to @xmath197 , which satisfies @xmath198 f_0^{\text{in}}(\rho ) = 0.\ ] ] for @xmath199 , @xmath200 we determine the phase @xmath201 numerically , by choosing a boundary condition @xmath202 for small enough values of @xmath203 , and by integrating the equation .
the boundary condition ensures that we pick up the right solution .
the result is @xmath204 . in order to find what value of @xmath22 ( @xmath205 in this limit )
is induced by @xmath69 , we have to match the inner and the outer solutions in the regime where @xmath185 .
this yields that @xmath206 .    remarkably , this numerical value for @xmath207 differs from the values obtained in ref .
@xcite and @xcite , which are @xmath208 and @xmath209 respectively .
it is unclear where the difference comes from , however , both papers use models different from the hyperspherical formalism that we used here , where we base our treatment of the effective range parameter on the descriptions in refs .
@xcite and @xcite .
more research is needed to resolve this discrepancy .
now we are in a good position to compare our model with earlier work on the efimov effect in a harmonic potential . in the work of thgersen , fedorov , and jensen ,
the influence of the effective range is studied on the universal efimov effect @xcite .
they considered a system with a gaussian interaction potential , given by @xmath210 where @xmath157 and @xmath76 are constants .
the range @xmath76 has been chosen to equal @xmath211 .    although the main purpose of ref .
@xcite is to show the influence of the effective range , their model contains a harmonic trap as well . in order to make a proper comparison with our model , we take @xmath212
the results are shown in figs .
[ fi : topth ] and fig .
[ fi : botth ] , and can be directly compared with fig .  1 in ref
.  @xcite .
we observe that the numerical results are reproduced quite closely and we can assign the numerically found bound states to the different hyperspherical channels . we also observed that if we include the above corrections that come from the effective range , the comparison is better than results obtained without these corrections .
we studied the efimov effect in a harmonic potential by using the hyperspherical formulation . in the first part of this paper
, we studied the universal regime by only taking a large scattering length into account , and by making the adiabatic hyperspherical approximation in that we ignored the coupling between different hyperspherical channels .
the model thus obtained is simple in the sense that the binding energies of the trimers result from an ordinary differential equation , with an appropriate boundary condition for the lowest channel encoded by the three - body parameter @xmath22 .
considering the lowest channel , we formulate an equation for a phase variable , the exponent of which is the radius in a spherical plot describing a surface of values in @xmath213-space ( inverse scattering length , wavenumber and inverse oscillator length ) for which bound states exists .
other such surfaces are obtained by multiplying the distance to the origin by a factor @xmath25 , and there is an accumulation point of surfaces near the origin , analogous to the efimov effect in free space . in the limit
@xmath214 , the efimov effect in free space is recovered .
interesting differences with the situation of three particles in free space can be observed when we focus at a fixed value of @xmath215 . for parameter combinations
where in free - space bound states would exist , the situation is very similar , except that the accumulation of energies only happens up to a typical energy of @xmath216 .
the rest of the efimov states are `` pushed out '' into the trimer continuum ( which is of course only a continuum in free space ) .
higher hyperspherical channels give rise to additional bound states , and by taking those into account as well we obtain results very close to numerical simulations of three - particle systems  @xcite .
better qualitative results are obtained by approximating the coupling between the hyperspherical channels , which lifts the degeneracies and creates avoided crossings .
the lifetime of the trimers is highest for states with a continuum character , i.e.  in the region which we would denote as trimer continuum in free space .
the reason for this is the small overlap of these states with deeply bound states .
experimental study of efimov - like states , however , could be feasible close to the transition from this continuum region to the trimer region , where the linewidth is still quite small .    for a more general efimov model that allows for deviations from the universal regime ,
a natural step is to improve on the approximation of the hyperangular problem by taking into account the effective range . by doing this
, we were able to find a relation between @xmath22 and the effective range @xmath177 if @xmath217 , where @xmath76 is a typical range of the potential .
we obtain the same qualitative behavior as ref .
@xcite and @xcite where different descriptions were used , however , the actual value of @xmath218 is also very different .
it appears that the actual value for this relationship depends strongly on the choice of the model and the corresponding approximations to describe the three - particle physics , and more research is needed to understand this discrepancy . in a different parameter regime for @xmath18 , @xmath69 and @xmath76 , such a relationship between @xmath22 and @xmath69 is not necessary .
then we have in addition to @xmath69 , the three - body parameter @xmath22 in the boundary condition for the hyperradial problem to our disposal . by comparing to numerical calculations
we find that this description indeed improves the predicted trimer energies .
99 v.  efimov , phys .
b * 33 * , 563 ( 1970 ) ; v.  efimov , sov .
. phys . * 12 * , 589 ( 1971 ) ; v.  efimov , _ ibid . _
* 29 * , 546 ( 1979 ) .
t.  kraemer , m.  mark , p.  waldburger , j.g .
danzl , c.  chin , b.  engeser , a.d .
lange , k.  pilch , a.  jaakkola , h .- c .
ngerl , and r.  grimm , nature * 440 * , 315 ( 2006 ) .
f. ferlaino and r. grimm , physics * 3 * , 9 ( 2010 ) . c.h .
greene , physics today * 63 * , 40 ( 2010 ) .
s.  knoop , f.  ferlaino , m.  berninger , m.  mark , h .- c .
ngerl , and r.  grimm , nature phys . * 5 * , 227 ( 2009 ) .
m.  zaccanti , b.  deissler , c.  derrico , m.  fattori , m.  jona - lasinio , s.  mller , g.  roati , m.  inguscio , and g.  modugno , nature phys .
* 5 * , 586 ( 2009 ) .
n.  gross , z.  shotan , s.j.j.m.f .
kokkelmans , and l.  khaykovich , phys .
103 * , 163202 ( 2009 ) .
pollack , d.  dries , and r.g .
hulet , science * 326 * , 1683 ( 2009 ) .
wenz , t.  lompe , t.b .
ottenstein , f.  serwane , g.  zrn , and s.  jochim , phys .
a * 80 * , 040702(r ) ( 2009 ) .
t. lompe , t. b. ottenstein , f. serwane , k. viering , a. n. wenz , g. zrn , and s. jochim , phys .
lett . * 105 * , 103201 ( 2010 ) .
huckans , j.r .
williams , e.l .
hazlett , r.w .
stites , and k.m .
ohara , phys .
102 * , 165302 ( 2009 ) .
williams , e.l .
hazlett , j.h .
huckans , r.w .
stites , y.  zhang , and k.m .
ohara , _ ibid . _
* 103 * , 130404 ( 2009 ) .
s.  nakajima , m.  horikoshi , t.  mukaiyama , p.  naidon , and m.  ueda , phys .
* 105 * , 023201 ( 2010 ) .
g.  barontini , c.  weber , f.  rabatti , j.  catani , g.  thalhammer , m.  inguscio , and f.  minardi , phys .
lett . * 103 * , 043201 ( 2009 ) .
h.  feshbach , ann . of phys .
( n.y . ) * 5 * , 357 ( 1958 ) ; * 19 * , 287 ( 1962 ) .
e.  tiesinga , a.  j.  moerdijk , b.  j.  verhaar , and h.  t.  c.  stoof , phys .
a * 46 * , r1167 ( 1992 ) ; e.  tiesinga , b.  j.  verhaar , and h.  t.  c.  stoof , phys . rev .
a * 47 * , 4114 ( 1993 ) .
e.  braaten and h .- w .
hammer , phys .
rep . * 428 * , 259 ( 2006 ) .
j.  von  stecher , j.p .
dincao , and c.h .
greene , nature phys .
* 5 * 417 ( 2009 ) .
f.  ferlaino , s.  knoop , m.  berninger , w.  harm , j.p .
dincao , h .- c .
ngerl , and r.  grimm , phys .
lett . * 102 * 140401 ( 2009 ) .
t.  lompe , t.b .
ottenstein , f.  serwane , a.n .
wenz , g.  zrn , and s.  jochim , science * 330 * , 940 ( 2010 ) .
s.  jonsell , h.  heiselberg , and c.j .
pethick , phys .
lett . * 89 * , 250401 ( 2002 ) .
fedorov and a.s .
jensen , phys .
lett . * 71 * , 4103 ( 1993 ) .
f.  werner and y.  castin , phys .
lett . * 97 * , 150401 ( 2006 ) .
m.   stoll and t.  khler , phys .
a * 72 * , 022714 ( 2005 ) .
m.  thgersen , d.v .
fedorov , and a.s .
jensen , phys .
a * 78 * 020501 ( 2008 ) .
x.  j.  liu , h.  hu , and p.
d.  drummond , phys .
a * 82 * , 023619 ( 2010 ) .
e.  nielsen and d.v .
fedorov , phys . rep . * 347 * 373 ( 2001 ) .
fedorov and a.s .
jensen , j. phys .
a , * 34 * 6003 ( 2001 ) . l.  platter , c.  ji , and d.r .
phillips , phys .
a * 79 * 022702 ( 2009 ) .
f.  werner , ph.d thesis _ trapped cold atoms with resonant interactions : unitary gas and three - body problem _ , ens , 2008 .
t.  busch , b .-
englert , k.  rzazewski , and m.  wilkens , found .
* 28 * , 549 ( 1998 ) .
t.  stferle , h.  moritz , k.  gnter , m.  khl , and t.  esslinger , phys .
lett . * 96 * , 030401 ( 2006 ) .
petrov , phys .
* 93 * , 143201 ( 2004 ) .
gogolin , c.  mora , and r.  egger , phys .
100 * 140404 ( 2008 ) .
n. gross , z. shotan , s. kokkelmans , and l. khaykovich , phys .
lett . * 105 * , 103203 ( 2010 ) .
n. gross , z. shotan , o. machtey , s. kokkelmans , and l. khaykovich , preprint arxiv:1009.0926 , accepted for comptes rendus physique . c.  ji , l.  platter , and d.  r.  phillips , epl * 92 * 13003 ( 2010 ) .
b.  marcelis , b. verhaar , and s.  kokkelmans , phys .
lett * 100 * , 153201 ( 2008 ) .


the orion nebula ( m42 ) is the most active site of star formation of the orion molecular cloud and contains many young stars , among them the trapezium cluster .
its high mass young stellar members generate an intense ultraviolet radiation field that photodissociates and photoionizes the nearby material .
these mechanisms promote the appearance of some complex structures as the so - called proplyds @xcite .
proplyds are low mass ysos that are being exposed to an intense ultraviolet radiation field which renders them visible . in orion , there are approximately 160 proplyds @xcite which are being photoionized mainly by @xmath6 ori c , an o6 spectral type star .
there are several objects that have been identified as being proplyds not only near the trapezium cluster @xcite but also in other star forming regions ( like ngc 3372 , * ? ? ?
most of them share the same features : a bow - shaped head that faces the ionization source , a tail , that is primordially directed away from the source , a young star , that may ( or not ) be visible , and a disk , sometimes seen in silhouette against the hii region ( e.g. , * ? ? ?
* ; * ? ? ?
* ; * ? ? ?
* ; * ? ? ?
early studies of these objects , however , were only able to determine their apparent ubiquity close to the @xmath6 ori c star , their high ionization level and , in later studies , the presence of high velocity structures ( outflows ) associated with these condensations ( e.g. , * ? ? ?
* ; * ? ? ?
* ; * ? ? ?
the proplyds are presently explained by a set of models which include a photoevaporated wind , an ionization front and a photoionized wind @xcite .
@xcite proposed models in which far ultraviolet radiation ( fuv ) and extreme ultraviolet radiation ( euv ) are responsible for the mass loss rates of the proplyds depending on the distance of the proplyd to the ionization source .
for the proplyds situated at intermediate distances from the ionizing star , fuv photons penetrate the ionization front and photodissociate and photoevaporate material from the accretion disk surface , generating a supersonic neutral wind that passes through a shock front before reaching the ionization front as a neutral , subsonic wind . at the ionization front ,
euv photons ionize the wind , and the material is then reacelerated to supersonic velocities .
the interaction of this supersonic ionized wind with the star wind generates bow shocks seen in h@xmath2 and in [ o iii]@xmath75007 in some proplyds @xcite . for proplyds
closer to the ionizing star , the ionization front reaches the disk surface , and the resulting wind is initially subsonic , becoming supersonic further out .    in general
, the models cited above are able to explain the main observational features of the orion proplyds .
the region close to the ionization front is responsible for most of the emission of these objects . @xcite
have shown that heating and dissociation of h@xmath8 can explain the observed [ o i]@xmath76300 emission at the disk surface @xcite .
@xcite proposed that neutral gas compressed between the shock front and the ionization front can explain the [ o iii]@xmath75007 emission seen near the ionization front in some objects . @xcite were able to reproduce the observed h@xmath2 intensity profile of the proplyds in orion in terms of accelerating photoevaporated flows .
numerical 2d hd simulations with a treatment of the radiative transfer @xcite also reproduce the proplyd morphology and emission , although in this work the treatment of the diffuse radiation is somewhat simplified .
there are several open issues about these objects .
one of them is related with the calculation of the mass loss rate of the proplyds , that is strongly model dependent and which may pose severe constraints on the age of the orion proplyds ( as well as for @xmath6 ori c ) . as the flux of fuv photons dissociates the molecules in the disk of the proplyds , the sub- or trans - sonic _ wind _
is accelerated at the proplyd s ionization front ( if ) ( see henney et al .
2002 for a discussion ) , leading to derived mass - loss rates of @xmath9 yr@xmath0 which imply a short lifetime for these systems @xcite .
it is hard to reconcile such short lifetimes with the fact that the region of the trapezium cluster is populated by several proplyds .
improvements in both models and observational techniques have slightly reduced the calculated mass - loss rates for these systems , and @xcite find @xmath10 yr@xmath0 for the 167 - 317 ( lv2 ) proplyd , implying an age of less than @xmath11 yr for @xmath6 ori c , and @xmath12 yr@xmath0 for the 170 - 377 , 177 - 341 , 182 - 413 , 244 - 440 proplyds ( e.g. , * ? ? ? * ) .
it is presently known that many proplyds show , besides protostellar features such as accretion disks and a young low mass protostar , the presence of jets .
@xcite , in a survey carried out using the wfpc2 camera on the hst , find 23 objects which appear to have collimated outflows seen as one - sided jets or bipolar chains of bow shocks .
@xcite also found kinematical traces of jets associated with lv5 ( 158 - 323 ) and gmr 15 ( 161 - 307 ) .
the proplyd 167 - 317 shows evidence of the existence of a collimated outflow , detected spectroscopically @xcite . until now
, this outflow was detected as a one - sided jet , with a p.a .
of @xmath13 and a propagation velocity of about 100 km s@xmath0 .
recently , @xcite and @xcite , using the wide field camera of the advanced camera for surveys ( acs / wfc ) on hst , found jets associated with silhouette disks in the outer regions of orion .    in this work , we present the first integral field unity ( ifu ) observation of a proplyd
the observed object is 167 - 317 , one of the brightest proplyds of the trapezium .
we discuss one of the very first results from the gemini multi - object spectrograph , in its ifu mode ( hereafter , gmos - ifu ) , at the gemini south telescope . in  2 , we present the observations and the steps to obtain a calibrated , clean datacube . in
3 we present the spectral analysis for the observed region of the 167 - 317 proplyd , as well as the observed lines and intensity maps .
we use the observed line profiles to identify the outflows in the system . in  4 , we present the discussion and the conclusions .
the data were taken during the system verification run of the gmos - ifu at the gemini south telescope ( gst ) , under the gs-2003b - sv-212 program , on 2004 february 26@xmath14 and 27@xmath14 .
the science field of view ( fov ) is @xmath15 with an array of 1000 lenslets of @xmath16 .
the sky is sampled with 500 lenses which are located @xmath17 from the science fov .
we used the r831_g5322 grating in single slit mode , giving a sampling of 0.34  per pixel ( @xmath18 km s@xmath0 per pixel at h@xmath2 ) and a spectral coverage between @xmath19   to @xmath20 . the spectral resolution ( i.e. , the instrumental profile ) is 47 km s@xmath0 @xmath21 fwhm @xmath22 km s@xmath0 ( the fwhm decreasing with increasing wavelength ) . a field centered on the 167 - 317 proplyd
has been observed with an exposure time of 300s . a 60 s exposure of the standard star
the hiltner 600 has also been taken in order to derive the sensitivity function .
the data were reduced using the standard gemini iraf v1.6 routines .
an average bias image has been prepared using the gbias task .
the extraction has been worked out using the gcal flat and the response using the twilight flat .
spectra extraction has been performed and wavelength calibration done using an arc taken during the run . the 167 - 317 field and hiltner 600 star have been processed in the same way .
the sensitivity function derived from the standard star has been applied to the observed field .
finally , data cubes have been built using the gfcube routine .
we maintain the ifu original resolution of @xmath23 px@xmath0 although an interpolation was performed in order to turn the ifu hexagonal lens shape into squared pixels .
because some of the emission lines are very intense , it was not possible to eliminate cosmic rays using the standard gscrrej routine of the iraf gemini v1.6 package . in order to remove cosmic rays
we have then programmed an idl routine .
this routine works directly on the data cube by taking out impacts that are above a certain convenient level calculated based on the mean intensity of nearby pixels .
the cosmic ray is then eliminated doing a linear interpolation in the wavelength direction . with this procedure , it is possible to keep very intense lines and eliminate lower level cosmic ray impacts that could be misinterpreted as low emission lines .
further manipulation of the cube , such as the production of channel maps , emission line maps , dispersion and velocity maps , have been done using starlink , idl routines and fortran programs developed specially for this purpose .
the subtraction of the background emission from the proplyd spectrum is a very challenging task @xcite .
this problem arises since the background nebular emission is strongly inhomogeneous on an arc - second scale . in this work ,
we have used standard @xmath24-method in order to fit an intensity versus position plane in a semi - rectangular field located @xmath25 away from the peak proplyd emission in each of the individual velocity channel maps .
we assure that we only take samples of the background avoiding the region of the proplyd .
these planar fits to the nebular emission are then subtracted from the corresponding velocity channel maps . in figure
[ f1 ] , we show the observed field , the object orientation in the plane of the sky , and also the region that defines our background ( region 4 , labeled as r4 in figure [ f1 ] ) .
region 4 ( see figure [ f1 ] ) contains 76 spectra which are used to define the background .
this figure also shows an internal box ( labeled as r1 , or region 1 ) , for which we have defined an arbitrary @xmath26-coordinate system limited by @xmath27 ( in the east direction ) and @xmath28 ( in the north direction ; see figure [ f1 ] ) . in the following sections ,
we show maps of spatially distributed physical variables that will be constrained to the limits defined by this box .
the spatial positions inside this domain will be defined using the @xmath26-coordinate system described above .
regions 2 and 3 ( labelled as r2 and r3 in figure [ f1 ] ) define , respectively , a region near the center of the 167 - 317 proplyd and a region which has a high - velocity , redshifted feature ( see below ) .
the lines that we identify in the spectra have already been reported in previous papers of the orion nebula ( e.g , . * ? ? ? * ) .
we select 4 different regions to measure the line intensities , namely , regions r1 , r2 , r3 and r4 ( see fig .
[ f1 ] and  2 for the coordinate system definition ) .
four spectra are then obtained by co - adding the spectra of the pixels included in each region .
the spectrum of region 1 is the result of the sum of 187 pixels , for region 2 the sum of 6 pixels , region 3 represents the spectrum of a single pixel and region 4 the sum of 76 pixels .
a set of 38 lines has been selected from a manual search for emission features in the spectrum integrated over the selected fields .
the intensities of the emission lines are determined by fitting a flat continuum and integrating over the whole emission feature .
then , a mean flux was obtained dividing the integrated flux by the total number of pixels of each region .
table  [ lines ] gives the list of the observed lines together with the mean flux of each region ( after subtraction of the emission from the background , region r4 ) and the ratio to h@xmath2 ( normalized to h@xmath29 ) .
the fourth and eigth columns show the background emission and ratios to h@xmath2 , respectively , for each line .
the mean fluxes are given in units of @xmath30 erg @xmath31 s@xmath0 px@xmath0 , and are not corrected for reddening .
the absolute error of the mean flux is given in parentheses .
this error was calculated taking into account the variations of the local continuum .
the errors can be very large for the weaker lines .
also , the detection of spectral lines in region r3 is affected by the low signal to noise ratio .
there are some lines which show stronger mean fluxes in the background , namely , the lines of si iii @xmath32 , si ii @xmath33 , [ ni ii]@xmath34 and [ fe ii]@xmath35  ( see table [ lines ] ) .
these lines have low ionization potentials and are expected to appear mainly in regions of lower ionization degree .
the lower line intensities within the proplyd indicate that these lines are mainly emitted by the background nebula , and that they are absorbed at least partially by the dust in the proplyd .    in figure [ f2 ] , we show the full extracted spectra from the data cube , for regions r1 to r4 as defined in figure [ f1 ] . in this figure ,
the spectra for regions r1 , r2 , r3 and r4 are depicted from bottom to top in each diagram .
the data shown here are not background subtracted .
most of the detected 38 lines can be clearly seen in the spectra .
it is also evident that the s / n is higher in regions r1 and r4 ( which is due to the increase of the pixel number of these regions ) . here
, a third order cubic spline was fitted to the continuum and then subtracted from the data . because of the high order spline polynome subtraction , some minor variations are still present near the h@xmath2 line .
figure [ f3 ] depicts ( for region r1 ) , 2d - intensity maps of the lines [ n ii]@xmath75755 , [ n ii]@xmath76548 , h@xmath2 , [ n ii]@xmath76583 , [ he i]@xmath77065 and [ ar iii]@xmath77135 , superimposed by line profiles .
the data shown here are background subtracted , as described in  2 .
the line profiles can be better seen in figure [ f4 ] , in which we show the mean line profile for each line and for each of the regions r1 , r2 , r3 and r4 , in the same diagram . as could be anticipated , the flux is more intense in region r2 ( dotted line ) , where the proplyd is located , except for the [ n ii]@xmath36 line , which shows a stronger flux in region r3 . for the h@xmath2 and [ n ii]@xmath37 lines the emission of region r3 is more intense than the emission of region r1 . for lines with good
s / n , the presence of a redshifted feature starts to become clear .
an analisys of the profile components will be given in more detail below .    in order to see how the spectra of the different lines change spatially , in figure [ f5 ] we present an intensity _ versus _ position plot for the same lines shown in figures [ f3 ] and [ f4 ] .
we show how the intensity of these lines changes along a diagonal line , that crosses region r1 passing through the pixel ( 3,4 ) ( region r3 ) and through the centre of region r2 .
we note that there are 2 intensity maxima for the h@xmath2 , [ n ii]@xmath38 and [ n ii]@xmath39 lines , clearly showing in which lines region r3 is clearly visible ( also see figure [ f3 ] ) .
we associate region r3 with the redshifted jet of the proplyd , although in figure [ f9 ] we will show that the emission of the jet is not restricted only to this pixel .
a redshifted jet associated with the 167 - 317 proplyd was previously detected , with a propagation velocity of @xmath40 100 km s@xmath0
( e.g. , * ? ? ? * ) .
we have also detected this high velocity feature in our data . in order to see this , we examine the behaviour of the moments of the radial velocity distribution .
these are the flux - weighted mean radial velocity @xmath41 , the flux - weighted rms width of the line @xmath42 and the skewness @xmath43 , that are given , respectively , by :    @xmath44    @xmath45    @xmath46    where    @xmath47    the integrated intensity ( see eq .
[ eq1 ] and figure [ f3 ] ) shows that if we move away from the proplyd position ( r2 ) , the emission for several lines ( for example , h@xmath2 and [ ar iii ] ) rapidly drops to low values ( close to the background value ) .
we also note that the maps for the flux - weighted mean radial velocity km s@xmath0 ( e.g. , * ? ? ?
the error in the radial velocity measurements is of the order of 2.5 km s@xmath0 . ] ( not shown here ; see eq .
[ eq2 ] ) show that @xmath48 km s@xmath0 .
the flux - weighted rms width of the line @xmath49 and the skewness @xmath50 moment of the radial velocity distribution ( see equations [ eq3 ] and [ eq4 ] , respectively ) are depicted in figure [ f6 ] , for the h@xmath2 ( left ) and [ ar iii ] ( right ) lines .
both of them were computed after carrying out both background and continuum subtractions . from figure [ f6 ]
, we see that there is a clear enhancement of both of these moments from the proplyd position ( r2 region ) towards the r3 region ( i.e. , towards the se direction , at a pa @xmath51 ) .
this behaviour of the moments leads us to infer that : 1 ) there is an increase in the line width as we go from the proplyd position to the se direction ; and 2 ) that there are redshifted wings in the line profiles .
these results indicate the presence of a redshifted outflow in this region .
we note that the same behaviour is seen in the he i @xmath76678.15 , [ n ii]@xmath76583.46 lines ( not shown here ) .    in order to confirm the presence of high velocity components in the observed field , as well as to see the behaviour of the photoevaporated proplyd flow and the background emission , we have computed three - component gaussian minimum @xmath24 fits for each position - dependent line profile .
the profiles of all emission lines have a major peak at @xmath52 km s@xmath0 , the exact radial velocity of this peak changing with spatial position . in some spectral lines
, there is an evident second peak at redshifted velocities ranging from 100 to 150 km s@xmath0 , corresponding to the jet associated with the proplyd .
finally , we have been able to detect a blueshifted component , which is fainter than the redshifted components . as an example of our three - component gaussian fit , in figure [ f7 ]
we show the data ( full line ) and the fit ( crosses ) , for the h@xmath2 and [ ar iii]@xmath77135 lines in the ( @xmath53,@xmath54)=(6,7 ) position . in this figure , for each emission line , the top - left panel represents the data and the three gaussian fits , the top - right panel shows the main , low velocity component ( and the gaussian fit for this component ) , the bottom - left panel depicts the data minus the main component together with the fits for the blue- and red - shifted components and , finally , the bottom - right panel shows the residual , obtained by subtracting the three - component fit from the observed line profile .
the fits depicted in this figure show the presence of redshifted and blueshifted components .
the redshifted component is present in several lines in pixels around the se direction , as already mentioned before .
on the other hand , blueshifted emission can be found in several pixels to the nw of the proplyd peak emission as can be seen in the left bottom panel of figure [ f7 ] and this could represent the first spectroscopic determination of the presence of a blueshifted counter - jet , with systemic velocities of -60 km s@xmath0 @xmath55 km s@xmath0 .
this blueshifted component is also suggested by the fitted profiles of hei @xmath77065 ( not shown here ) .
figure [ f8 ] depicts the spatially distributed intensity over the r1 region ( see figure [ f1 ] ) , of the main peak for the [ n ii]@xmath75754 ( top - left ) , [ n ii]@xmath76548 ( top - middle ) , h@xmath2 ( top - right ) , [ n ii]@xmath76583 ( bottom - left ) , hei @xmath77035 ( bottom - middle ) and [ ar iii]@xmath77135 ( bottom - right ) lines .
the values shown in this figure were obtained from three - component gaussian fits to the observed line profile .
all of the lines peak at the proplyd position , at ( @xmath53,@xmath54 ) = ( 6,7 ) , inside the r2 region , and three of them ( namely [ n ii]@xmath76548 , top - middle ; h@xmath2 , top - right ; and [ n ii]@xmath76583 , bottom - left ) , have a secondary , less intense peak at ( @xmath53,@xmath54)=(3,4 ) ( i.e. , in the r3 region ; see figure [ f1 ] for the definition of the coordinate system ) .
this secondary peak is probably related to the jet , since , as we have seen before , the jet propagates from the r2 region towards the se direction .
it is also interesting to note that there is a tail of faint emission ( compared with the maximum in a given map ) that extends in the ne direction , and that can be seen in both of the [ n ii ] lines that bracket h@xmath2 ( top - middle and bottom - left maps in figure [ f8 ] ) .    in figure
[ f9 ] , we show the spatially dependent intensity ( fig . [ f9]a ) and the central velocity ( fig .
[ f9]b ) of the redshifted component of the fitted h@xmath2 profile .
we limit the maps to the spatial pixels which have a skewness @xmath56 km s@xmath0 .
the pixels which satisfy this criterium have a well defined high velocity , redshifted component .
these figures show that the high velocity component has an intensity maximum near the center of the proplyd , with an extension towards the se direction ( along the jet axis ) , surrounded by a region of decreasing fluxes .
the high intensity spike seen in h@xmath2 , [ n ii]@xmath57 and [ n ii]@xmath58 lines ( figure [ f3 ] ) can be seen here .
the jet shows a spread in velocity values , although most of the pixels present velocities around 120 - 140 km s@xmath0 .
this figure also shows that the jet velocity decreases with increasing distance from the proplyd .
this trend can be seen in the figure [ f9]c , where we show the radial velocity of the redshifted component as a function of distance from the proplyd position , ( @xmath53,@xmath54 ) = ( 6,7 ) in region r1 .
to build this figure , we assume that the jet is propagating in a pa @xmath51 position angle .
we define a box with one axis aligned with the jet axis , and the second , perpendicular dimension extending two pixels to each side of the jet axis .
we then take averages ( of the mean velocity of the redshifted component ) perpendicular to the jet direction and plot the resulting velocity as a function of distance from the proplyd ( see the bottom frame of figure [ f8 ] ) .
there is an indication that the jet velocity is slowly diminishing as a function of distance from the source .      in figure [ f10]a
, we show a map for the r1 region of the line ratio :    @xmath59\lambda 6548 + [ { \rm n~ii}]\lambda 6583}{[{\rm n~ii}]\lambda 5754 } = \frac{a_{6548 } \sigma_{6548 } + a_{6583 } \sigma_{6583}}{a_{5754 } \sigma_{5754}},\ ] ]    where @xmath60 and @xmath61 are the height and the dispersion of the fitted gaussian profile ( of the main , low velocity component ) .
the [ nii](6548 + 6583)/5754 ratio is a classical electron temperature diagnostic of low / medium density nebulae ( @xmath62 @xmath63 ) . at high densities , due to the collisional desexcitations
, this ratio can not be used for the determination of the electron temperature , but it turns out to be a quite good electron density diagnostic .
this can be seen , for example , from the @xcite equation :    @xmath64(6548 + 6583)/5754 = \frac{6.91 \exp ( 2.5 \times 10^{4 } / t_{e})}{1 + 2.5 \times 10^{-3 } ( n_{e } / t_{e}^{1/2})}\ ] ]    figure [ f10]b shows the lines corresponding to [ nii](6548 + 6583)/5754 = 9 , 10 , 11 , 80 , 90 , 100 in a @xmath65 plane .
we can see that for the values obtained in our case ( line ratios close to 10 ) , the [ nii](6548 + 6583)/5754 ratio is strongly dependent on the electron density , and that the dependence on the electron temperature is small .    in figures [ f10]c and [ f10]d
, we show the densities computed from the [ nii](6548 + 6583)/5754 ratio for @xmath66 k and @xmath67 k , respectively .
although presenting a very complicated pattern , these figures show that the region coinciding with the proplyd center , and extending towards the se direction ( the jet propagation direction ) , shows the lowest line ratio , indicating electron densities ( @xmath68 ) of the order of @xmath69 @xmath63 .
this value is at least in qualitatively agreement with the @xmath70 @xmath63 obtained by @xcite considering a @xmath71 k temperature for the ionization front ( that is inside the r2 region ) .
we can also note that both figures are similar , presenting few differences in the density values .
the extended structure from the r2 region towards the se direction is surrounded by a region in which the nitrogen line ratio reaches values of up to @xmath72 . in this limit , @xmath73 @xmath63 .
it is interesting to note that in the clump associated with the redshifted jet , the nitrogen line ratio increases substantially , indicating a decrease in the electron density when compared with the same values near the center of the proplyd .
we have also obtained the [ s ii]@xmath74 emission maps .
unfortunately , we are not able to use these lines as a diagnostic because neither the ratio of this doublet is a good electron density indicator for this case ( since it is constant for densities higher than 10@xmath75 @xmath63 ) nor do we have a good enough signal to noise ratio for these lines . furthermore , the background subtracted spectra for these two lines give us negative fluxes in the proplyd region ( see table  [ lines ] ) .
the presence of these negative features in the line profiles is probably related with the presence of dust , as discussed by @xcite .
we have presented in this work the first integral field unit ( ifu ) spectroscopic observations of a photoevaporating disk immersed in a hii region . in particular ,
we have taken advantage of the system verification run of the gemini south telescope multi object spectrograph ( gmos ) to obtain spectra of the 167 - 317 proplyd in the orion nebula .
the 167 - 317 proplyd , also known as lv2 ( from the pioneering work from * ? ? ?
* ) , is one of the brightest and best studied proplyds . in a single exposure , we took 400 spectra with a spatial resolution of @xmath16 .
these spectra have been combined in order to optimize this instrumental feature , and , before discussing of our results , we want to make a few comments concerning the potential use of such ifu data - cubes for these objects . as previously discussed in the literature @xcite , the subtraction of the background emission is a very challenging task .
this problem arises since the background near these proplyds ( which are immersed in a highly non - homogeneous , photoionized ambient medium ) is very complex and changes in an arc - second spatial scale .
a careful computation of the background contribution to the spectra of these objects is one of the most important tasks in order to have reliable radial velocity measurements for both the proplyd photoevaporated flow and the high velocity features ( jets ) present in these systems .
here , we have used the 76 spectra , defined in figure [ f1 ] as region r4 , in order to obtain a planar fit for the background emission .
we think that such a definition of the background emission is more precise , since we use several regions to calculate a better function to describe the background .
one more thing to point out is that , with spatially distributed spectra , we are able to confidently separate the contribution of the background from the contribution of the object itself
.    the data from the ifu observation of the 167 - 317 proplyd also allow us to investigate the spatially dependent properties of the outflows associated with this object .
in particular , the 167 - 317 proplyd is known to have a redshifted , collimated jet that propagates with heliocentric velocities of @xmath76 km s@xmath0 towards the se direction , and with a spatial extension of @xmath77 ( e.g , * ? ? ? * ; * ? ? ? * ; * ? ? ? * ; * ? ? ? * ; * ? ? ? * ) .
we find that a prominent , high velocity redshifted component can be detected in some emission lines , particularly in h@xmath2 , in this se direction .
the redshifted jet has a trend in its radial velocity ( as a function of distance from the proplyd ) with higher velocities close to the proplyd and lower velocities at increasing distances .
we note that @xcite had previously suggested a variation in the jet velocity as a function of distance from the source .
there is a subtle peak ( in several emission lines , particularly in [ n ii]@xmath75754 and h@xmath2 ) that is associated with the se jet emission , and located at @xmath78 se of the center of the proplyd [ at ( x , y)=(3,4 ) ; the r3 region in figure [ f1 ] ] . for this intensity peak , and also using the h@xmath2 profiles of the neighbouring pixels , we find a mean heliocentric velocity @xmath79 km s@xmath0 .
it is interesting to note that previous hst images and spectroscopic analyses @xcite also reveal a spike in the jet emission at @xmath80 from the proplyd cusp .
if we associate these two emission regions as belonging to the same condensation , we can estimate a proper motion of @xmath81 km s@xmath0 , which combined with the @xmath82 km s@xmath0 radial velocity give a full jet velocity of @xmath83 km s@xmath0 ( where the error was inferred from the spatial resolution of our data ) .
this velocity is of the same order as the values inferred for herbig - haro jets associated with t tauri stars .
however , the exact nature of this spike is unknown .
we could relate this spike with a bow shock of the jet , since the line profiles for both are similar @xcite , showing a low velocity , more intense peak together with a high velocity , less intense peak ( similar to the line profiles seen in figures [ f3 ] , [ f4 ] and [ f7 ] ) .
we have found evidence for a blueshifted component in several emission lines around the northwest tip of the proplyd position . in particular , the h@xmath2 profiles reveal the presence of a blue - shifted component with systemic velocity @xmath84 km s@xmath0 , which indicates that its velocity is similar though smaller than the one of the redshifted jet . however , the intensities of the blue - shifted features are much lower than the intensity of the redshifted components associated with the se jet : the counter - jet is at least four times less intense than the redshifted jet ( in h@xmath2 ) .
the presence of such a blue - shifted component very close to the lv2 proplyd , together with previously reported evidence that ( at several arcseconds to the nw of this proplyd ; see massey & meaburn 1995 ) there is a blueshifted component ( in the [ o iii]@xmath75007 line ) strongly suggests the existence of a faint counterjet .
it is interesting to estimate the mass loss rate for the proplyd and the associated redshifted jet . in order to do this
, we consider a simple model that assumes that the proplyd flow arises from a hemispherical , constant velocity wind that originates at the proplyd ionization front , at a radius @xmath85 .
using the h@xmath2 luminosity @xmath86 , we can then obtain the particle density @xmath87 at this point and the mass loss rate , from :        where @xmath90 accounts for the extintion correction , @xmath91 @xmath92 is the effective recombination coefficient for h@xmath2 , @xmath93 is the energy of the h@xmath2 transition , @xmath94 is the mean molecular weight and @xmath95 km s@xmath0 is a typical value for the sound velocity in an hii region .
the extintion correction is obtained from the base 10 logarithm of the extintion at @xmath96 ( @xmath97 ) which for the proplyd 167 - 317 is equal to 0.83 @xcite through the relation , @xmath98 @xcite , where here we take @xmath99 = 0.56 .
assuming for the ionization front radius a value @xmath100 cm @xcite we obtain @xmath101 @xmath63 and @xmath102 m@xmath4 year@xmath0 . @xcite have obtained a @xmath103 m@xmath4 year@xmath0 @xmath104 , and the small differences between both results may arise from the different observational techniques employed in both cases .            where @xmath108 cm is the radius of the blob , @xmath109 km s@xmath0 and @xmath110 cm are the propagation velocity and length of the jet , respectively , derived above .
with these values , we obtain for the jet a mass loss rate equal to @xmath111 m@xmath4 year@xmath0 , which is similar to the mass loss rate of typical hh jets from t tauri stars .
we have been able to construct a ( [ nii]@xmath36 + [ nii]@xmath37)/[nii]@xmath112 line ratio map .
this ratio indicates densities higher than @xmath11 @xmath63 for the proplyd emission region , consistent with values previously reported in the literature .
we find a subtle enhancement of this ratio in the region of the redshifted jet .
further observations of electron temperature and density diagnostic lines would allow a more accurate determination of the spatial dependence of these important parameters .
we thank the anonymous referee for his / her comments and suggestions , which have improved significatively the presentation of the paper .
we are greatful to w. henney for enlightning discussions .
m.j.v . and a.h.c .
would like to thank the staff of the instituto de ciencias nucleares , unam , mxico , for their kind hospitality , as well as for partial financial support during our visit .
we would like to thank bryan miller , tracy beck and rodrigo carrasco from the gemini observatory , and bruno castilho from lna / cnpq , for help us in the early stages of ifu - data processing .
the work of m.j.v . ,
a.h.c . and h.p .
was partially supported by the bahia state funding agency ( fapesb ) , by the universidade estadual de santa cruz ( uesc ) and by the milenium institute ( cnpq project n. 62.0053/01 - 1-padct
iii / milnio ) .
this work was partially funded by the conacyt grants 41320-f and 43103-f and the dgapa ( unam ) grants in112602 and in111803 .
integrated emission map ( shown with a linear grey scale in arbitrary units ) showing the location of the 167 - 317 proplyd with respect to the scientific field of view ( sfov ) , and its orientation : north is up and east is left , as indicated by the arrows on the top - left corner of the figure .
the field has 16 pixels in the east direction and 25 in the north direction . the r4 frame located at the left and top parts of the figure represents the region choosen to define the background , which covers two columns of pixels on both the north and east directions .
an internal box ( r1 ) , approximately centered on the object emission defines a coordinate system , and the numbers on the bottom axis and on the left axis of this internal box are referred to in the text as the @xmath53- and @xmath54-coordinates , respectively .
this spatial sub - sample of spectra will be used throughout the paper .
the @xmath6 ori c star is located at @xmath113 from the proplyd , at a p.a .
= 236@xmath114 .
four regions are defined in this figure , namely r1 ( the internal box ) , r2 ( the proplyd center ) , r3 ( a clump associated with a redshifted jet ) and r4 ( the background ) .
we will refer to these regions throughout the paper .
, width=340 ]     has been added to the spectra of regions r2
( @xmath115 ) , r3 ( @xmath116 ) and r4 ( @xmath117 ) in order to produce an intensity offset between the successive spectra .
wavelengths ( @xmath7 ) are in  , and the fluxes are in units of @xmath30 erg @xmath31 s@xmath0 @xmath0 . , width=566 ]    -3truecm ) , showing the integrated flux ( the grey - scale map ) for the following lines : [ n ii]@xmath75755 ( top - left ) , [ n ii]@xmath76548 ( top - center ) , h@xmath2 ( top - right ) , [ n ii]@xmath76583 ( bottom - left ) , [ he i]@xmath77065 ( bottom - center ) and [ ar iii]@xmath77135 ( bottom - right ) .
the linear bars give the grey - scale in unit of @xmath118 erg @xmath31 s@xmath0 px@xmath0 .
the spectra in each position was background subtracted .
superimposed on the intensity maps are the line profiles for each box .
see the text for a discussion .
the orientation is the same as in figure [ f1 ] ( i.e. , north is up and east is left ) .
the original pixels ( see figure [ f1 ] ) have been resampled in order to smooth the emission features , and each box in these panels are slightly larger than the original pixels in fig .
title="fig:",width=566 ]    -3truecm .
each diagram shows the line profile for regions r1 ( continuous line ) , r2 ( dotted line ) , r3 ( dashed line ) and r4 ( dash - dotted line ) .
the emission line is indicated on the top - left corner of each plot .
the fluxes are in units of @xmath118 erg @xmath31 s@xmath0 per    and tey are plotted against the radial velocity ( in units of km s@xmath0 ; not corrected for the heliocentric velocity ) .
, title="fig:",width=566 ] -5truecm     as a function of the pixel position along a diagonal line that crosses region r1 , beginning at pixel ( 1,2 ) , passing through region r3 ( pixel 3,4 ) and finishing at pixel ( 11,12 ) .
each type of line shows the intensity of a different spectral line : the solid line shows [ n ii ] @xmath119 ; the short dashed line , the intensity of [ n ii ] @xmath39 ; the dot - dashed line , h@xmath2 ; the triple dotted - dashed line , [ n
ii ] @xmath38 ; the long dashed line , he i @xmath120 and , finally , the dotted line shows [ ar iii ] @xmath121 . in order to show all of the intensities in the same plot ,
we have divided the h@xmath2 intensity by a factor of 10 .
, width=566 ]    -2truecm ) and ( [ eq4 ] ) , respectively , for the h@xmath2 ( left ) and [ ar iii]@xmath77135.8 ( right ) lines .
the grey scale on the right side of each map is in unit of km s@xmath0 .
we note that both integrals show positive values and peak in a se region , indicating the presence of a redshifted component towards the se direction .
the velocities are not corrected by the systemic radial velocity , and the integration limits in all maps are from @xmath122 km s@xmath0 .
the white pixels have values close to zero . , title="fig:",width=377 ] -4truecm     line at position ( @xmath53,@xmath54)=(6,7 ) ( see figure [ f1 ] ) . @xmath123 line profile and the three gaussian fit .
@xmath124 main , low velocity line profile and its gaussian fit .
@xmath125 the data profile subtracted from the main fitted profile , showing evidence for a blue- and a red - shifted component
. @xmath126 the data profile subtracted from the three - component gaussian fit .
the residual profile is a random fluctuation around a flux of 0.5@xmath127 erg @xmath31 s@xmath0 .
these fluctuations corresponds to 10% of the flux for the blueshifted component , 0.5% of the flux for the main component and 10% of the flux for the redshifted component .
b ) the same but for the [ ar iii]@xmath77135 line . , title="fig:",width=453 ] -5truecm     erg @xmath31 s@xmath0 px@xmath0 ) as a function of position for the [ n ii]@xmath1285754 ( top - left ) , [ n ii ] @xmath76548 ( top - center ) , h@xmath2 ( top - right ) , [ n ii]@xmath76583 ( bottom - left ) , hei @xmath76678 ( bottom - center ) and [ ar iii]@xmath77135 ( bottom - right ) lines .
almost all of the lines peak at the proplyd position , and three of them , namely the [ n ii]@xmath76548 ( top - centre ) , h@xmath2 ( top - right ) and [ n ii]@xmath76583 ( bottom - left lines ) , have a secondary , less intense peak in region r3 ( see figure [ f1 ] ) .
, title="fig:",width=566 ] -9truecm    -2truecm   erg @xmath31 s@xmath0 px@xmath0 , as a function of position in region r1 .
b ) the velocity , in km s@xmath0 for the redshifted gaussian fit to the h@xmath2 profiles .
c ) the jet velocity variation as a function of distance from the proplyd ( see the text for a discussion ) . the velocities take into account the systemic radial velocity of @xmath129 km s@xmath0 .
the error associated with the radial velocity measurement is of the order of 2.5 km s@xmath0 ( see the text for a discussion ) . , title="fig:",width=566 ] -6truecm     @xmath130([nii]@xmath36 + [ nii]@xmath37)/[nii]@xmath112 ratio map of region r1 ( see figure [ f1 ] ) of the observed field .
@xmath131 a [ @xmath132 diagnostic diagram showing the locci corresponding to different values of the nitrogen line ratio .
@xmath133 the density for an assumed temperature of @xmath66 k. @xmath134 the density for an assumed temperature of @xmath135 k. , title="fig:",width=604 ] -8truecm