recently , the proximity effect between a superconductor ( @xmath0 ) and a ferromagnet ( @xmath1 ) has attracted much attention because , due to advances in materials growth and fabrication techniques @xcite , well controlled structures , in which it can occur became available .
such @xmath1/@xmath0 hybrid structures are important from the point of view of their intrinsic scientific interest , as they allow the study of the interplay between ferromagnetism and superconductivity @xcite as well as of device applications in such areas of technology as magnetoelectronics @xcite or quantum computing @xcite .    in the present context the proximity , on one hand , means the leakage of superconductivity into the ferromagnet , and on the other , the spin polarization of the superconductor . in the normal metal - superconductor systems
the proximity effect has been studied over 30 years , and is now well understood @xcite to be governed by the andreev reflection @xcite processes , in which an electron at the interface is reflected as a hole while a cooper pair is transfered into superconductor .
this mechanism makes the movement of electrons and holes in the normal metal highly coherent .
if the normal metal is replaced by a ferromagnet a number of new phenomena may occur .
for instance one may expect a version of the @xmath2 state , predicted in sixties for a bulk superconductor in a strong spin exchange field , by fulde and ferrell @xcite and larkin and ovchinikov @xcite , to be realized . in the bulk , as it is well known , the exchange field tends to polarize the conduction electrons and holes . in particular
if these electrons are cooper paired , naively one would expect that either the exchange field is too week to break the pairs , or it leads to the first order phase transition to the normal state .
however , as was demonstrated in ref .
@xcite , for certain values of the exchange field a new superconducting depairing state is realized through first order phase transition from the @xmath3 state and it transforms continuously , by a second order phase transition , into the normal state as the strength of the exchange field is increased .
this @xmath2 state has a spatially dependent order parameter corresponding to the nonzero center of mass motion of the cooper pairs .
another novel feature of this state is a current flowing in the ground state .
the unpaired electrons tend to congregate at one portion of the fermi surface so a quasiparticle current is produced . in order to satisfy the bloch theorem : no current in the ground state , a supercurrent , generated by the nonzero value of the pairing momentum , flows in opposite direction , and the total current is zero .
interestingly , similar oscillations of the pairing amplitude ( cooper pair density ) have been predicted @xcite-@xcite in a ferromagnet proximity coupled to a superconductor .
it turns out that these oscillations are responsible for the oscillatory dependence of the critical temperature on the thickness of the ferromagnetic slab @xcite .
this effect as well as the corresponding oscillations of the density of states at the fermi level @xcite , have been observed experimentally @xcite . in this paper
we investigate further the ramifications of these interesting phenomena .
one of these is the formation of so called @xmath4 junction state @xcite which is also a fingerprint of the oscillatory behaviour of the pairing amplitude .
the @xmath4-state effect has been extensively studied in connection with the high-@xmath5 superconductors , where due to spatially inhomogeneity of the order parameter , the cooper pair wave function can change sign at the interface and this leads to the formation of a zero - energy mid - gap state . remarkably , this zero - energy state is unstable to the occurrence of a spontaneous current flowing parallel to and within a coherence length of the interface ( see eg .
@xcite and references therein ) . in what follows we report on our finding similar spontaneous currents but of different origin in a @xmath0/@xmath1/@xmath6 trilayer
.    we also find zero energy bound states due to the finite size of the @xmath1 slab .
in fact even for a @xmath7/@xmath0 heterostructure , in agreement with @xcite , we have found states within the @xmath0 gap .
although the energy of these states depend on the thickness of the normal metal , they never reach zero , namely @xmath8 .
interestingly , when the normal metal is replaced by a ferromagnet , these bound states split , and for certain values of the exchange field , they cross the fermi level @xcite . as we shall show , this circumstance leads to a current flow .
clearly , the spontaneous current is strictly related to the zero energy states , as in the case of the high - @xmath5 superconductors , but the origin of these states is completely different . as mentioned above ,
there are two competing parts to the total current in the @xmath2 state .
one is related to the unpaired electron movement and the other is a supercurrent .
so one may expect , that in the case of the @xmath1/@xmath0 heterostructure , these parts could be spatially unbalanced creating regions of net current flow .
indeed this is what we have discovered .    a brief report of our main results
has already been published @xcite . here
we wish to present a more detailed and systematic study of the zero energy andreev bound states and the corresponding spontaneous current in the @xmath1/@xmath0 heterostructures .
the paper is organized as follows : in sec .
[ sec2 ] the simple model which allows for self - consistent description of the @xmath1/@xmath0 heterostructure is introduced .
we also derive equations which have to be solved self - consistently .
in particular , these include the maxwell s equation ( ampere s law ) for the coupling of the current to the magnetic field .
then we describe some technical details concerning the principal layer technique , which allows us to describe a semi - infinite superconductor . in sec .
[ sec3 ] the nature of the andreev bound states in the ferromagnet is discussed .
the spontaneously generated current and corresponding magnetic field in the ground state are studied in the sec .
[ sec4 ] . in sec .
[ sec5 ] we provide some suggestions of how one might observe the spontaneous current and its polarization experimentally .
finally , we conclude in sec .
to study spontaneous currents and their polarization in the ferromagnet - superconductor heterostructure we have adopted a single orbital , nearest neighbour hopping hubbard model with negative @xmath9 on the @xmath0 side and zero @xmath9 otherwise .
additionally we have assumed the site energies @xmath10 to be exchange split in the ferromagnet and spin independent in the superconductor .
moreover , the simplest geometry allowing for a current flow is the @xmath11 system , shown in the fig.[fig1 ] , where the magnetic field in one direction , the vector potential and the current in another and the spatial modulation in a third orthogonal direction are explicitly indicated .    in short ,
our model hamiltonian is given in the form : @xmath12         c^+_{i\sigma } c_{j\sigma } +         \sum_{i\sigma } \frac{u_i}{2 } \hat n_{i\sigma } \hat n_{i-\sigma }   \label{hamiltonian}\end{aligned}\ ] ] where , in the presence of a vector potential @xmath13 , the hopping integral is given by @xmath14 for nearest neighbour lattice sites , whose positions are @xmath15 and @xmath16 , and zero otherwise .
the site energies @xmath10 are @xmath17 on the superconducting side and equal to @xmath18 on the ferromagnetic side , @xmath19 is the chemical potential , the on - site interaction @xmath20 is @xmath21 in the superconductor and zero elsewhere , @xmath22 , ( @xmath23 ) are the usual electron creation ( annihilation ) operators and @xmath24 is the electron number operator .
note that the above description of the electrons with charges @xmath25 includes a coupling to a magnetic field @xmath26 , which is necessary for calculating the effects of current on the electronic states .    in what follows
we shall work within spin - polarized - hartree - fock - gorkov ( @xmath27 ) approximation , which means that we have approximated the interaction term in hamiltonian ( [ hamiltonian ] ) by the usual mean field theory mapping : @xmath28 .
for the magnetic field we have chosen the landau gauge where @xmath29 and hence @xmath30 .
furthermore , we have assumed that the effective @xmath27 hamiltonian is periodic in the direction parallel to the interface and therefore we work in @xmath31 space in the @xmath32 direction but in real space in the @xmath33-direction ( see fig.[fig1 ] ) . labelling the  planes  ( lines ) in fig.[fig1 ] by integers @xmath34 and @xmath35 at each @xmath36 point of the brillouin zone we shall solve the following @xmath27 layer index ( @xmath34 ) matrix equation for the retarded green s function matrix @xmath37 : @xmath38 where @xmath39 is the pauli unit matrix and @xmath40 is of the form :    @xmath41    with @xmath42 .
our system consists of infinite number of layers ( @xmath34 ) , because of fact that we have a semi - infinite superconductor .
the @xmath1 region spreads from @xmath43 to @xmath17 , while superconductor is defined for @xmath44 ( see fig .
[ fig1 ] ) . in order to describe semi - infinite system
, we need infinite range matrix with block elements @xmath45 ( @xmath34 , @xmath35 going from @xmath46 through @xmath17 to @xmath47 ) . in short , strictly speaking , we have to invert the infinite range matrix :    @xmath48    to render the problem traceable , it is useful to define the surface green s function ( @xmath49 ) @xcite as @xmath50^{-1}\}_{n , n }   \label{sgfdef}\end{aligned}\ ] ] @xmath51 being semi - infinite submatrix for @xmath52 ( lower right corner of the matrix in eq.([hinf ] ) ) .
the physical meaning of the @xmath49 is that it represents a square diagonal subblock corresponding to the semi - infinite bulk shown with a shaded region in the fig .
[ fig1 ] .
now we can assume , that electronic properties in the @xmath0 differ from that in the bulk only over the finite distance from the interface ( @xmath53 ) .
@xmath54 can correspond to the distance of a few @xmath0 coherence lengths . for @xmath55
the electronic properties are that of the bulk superconductor . in other words , we assume that the system for @xmath55 is homogeneous .
namely , for @xmath55 we have : @xmath56 as a consequence the surface green s function ( [ sgfdef ] ) is also @xmath34-independent and can be obtained from the bulk @xmath57 @xcite : @xmath58 where @xmath59 is the fourier transform of the green s function of the homogeneous superconductor @xmath60 , namely @xmath61 .
clearly , the green s function @xmath60 for the case where the vector potential is equal to zero is given by the usual result :    @xmath62    where @xmath63 is the single electron energy spectrum .
the homogeneity of the system for @xmath55 leads to the conclusion , that the inversion of a semi - infinite block tridiagonal matrix ( [ hinf ] ) can be reduced to the inversion of the finite tridiagonal matrix :    @xmath64    in the above equation the effect of the semi - infinite homogeneous system is contained in the embedding potential @xmath65 , which is related to the surface green s function through : @xmath66 evidently , by increasing @xmath54 until @xmath67 in the surface region does not change would constitute an exact albeit numerical solution of the eq .
( [ hfg ] ) .
all calculations reported here were performed for @xmath68 .
however , we have checked that larger @xmath54 has no influence on the results .
the calculation of the various physical quantities , like particle concentration @xmath69 , spin polarization @xmath70 , order parameter @xmath71 , etc .
, needs the evaluation of integrals of the product of the green s function and fermi distribution function as in @xmath72 . while working on the real energy axis , we need a large number of integration points to get good accuracy , unless there are no singularities in the density of states
. however , in our calculations , there is a van hove singularity in the @xmath11 @xmath73 as well as @xmath3-like singularities .
in fact we need approximately @xmath74 - @xmath75 points to get good accuracy , which is very time - consuming in self - consistent calculations .
the situation is similar for summation over the matsubara energies @xmath76 because of the poles of the fermi function .
there is an infinite number of poles , and infinite sum does not always converge well , so numerical treatment needs a large number of data points .    to overcome these difficulties we follow ref .
@xcite and approximate the fermi distribution function by : @xmath77 to get the correct behaviour near the fermi energy ( @xmath78 ) we choose @xmath79 .
one can check , that the discrepancy between the approximate function and the true one , starts to play role for energies lower than @xmath80 .
if we choose @xmath81 , where @xmath82 is the half of the bandwidth , the approximation works very well . for example , for the temperature @xmath83
it is sufficient to consider only @xmath84 points to get very good accuracy .
the poles of the approximate function lie on a circle given by @xmath85 the solution of the above equation gives the poles in the form : @xmath86 in the limit @xmath87 , the poles ( [ polesenergy ] ) move to the matsubara energies , and the number of poles , @xmath88 , becomes infinite .
the calculation of any physical quantity becomes : @xmath89 where @xmath90 is given by eq.([polesenergy ] ) .      as usual in problems of interfaces , a large number of equations has to be solved self - consistently . in our case the parameters , to be determined self - consistently at each layer @xmath34 ,
are : the total electron density @xmath69 , the spin polarization @xmath70 , the @xmath0 order parameter @xmath71 , the current parallel to the interface @xmath91 and the vector potential @xmath92 .
most of these quantities can be expressed in terms of the corresponding layer - diagonal elements of the green s functions @xmath93 ( see eq.([hhfg ] ) ) .
the exception is the vector potential , for which we need an additional ( maxwell s ) equation .
these relations are ( see eq.([oaverage ] ) ) : @xmath94 @xmath95 @xmath96 where @xmath97 is the @xmath98-th element of the matrix green s function , which is the solution of the @xmath27 equation ( eq .
( [ hfg ] ) ) .
the current , for spin up electrons , in the @xmath32-direction can be calculated from the relation : @xmath99 which follows from the continuity equation for the charge ( @xmath100_- = - \sum_j j_{ij}$ ] ) .
this current will give rise to a vector potential @xmath101 which will have to be used to update @xmath102 in eq.([hfg ] ) at the end of each self - consistency cycle .
we calculated this new vector potential by solving numerically ampere s law , @xmath103 , which for the lattice problem at hand , is @xmath104 in all our calculations the eqs.([n])-([maxwell ] ) have been solved in each iteration step until self - consistency has been achieved .
the numerical calculations have been performed at temperature @xmath83 for @xmath84 energy points , @xmath105 - @xmath36 points and @xmath106 layers ( @xmath107 ferromagnetic and @xmath107 superconducting ) .
typically self - consistency , at the level of @xmath108 on all densities in eqs.([n])-([maxwell ] ) , has been achieved after @xmath109 - @xmath110 iterations .
since we have determined the ferromagnetic ( [ m ] ) and superconducting ( [ delta ] ) order parameters on both sides of the interface fully self - consistently , we were able to study both @xmath0 and @xmath1 proximity effects .
this means , that we were able to describe situations in which ferromagnetism and superconductivity coexist near the interface .
firstly , we want to discuss the magnetic proximity effect , namely the entering of the spin polarization into the superconductor .
the typical behaviour of the spin polarization ( @xmath111 ) is plotted in the fig .
[ fig2 ] as a function of the layer index @xmath34 for a number of exchange splittings @xmath112 .
we can see that the spin polarization exponentially decays over the distance of the @xmath0 coherence lengths @xmath113 , which in the present case is @xmath114 in units of the lattice constant @xmath115 .
we have checked this by explicit calculations of @xmath113 for various values of the superconducting energy gap .
so we can conclude that effect of proximity of the ferromagnet on the superconductor is very similar to case of the meissner effect , where an external magnetic field is excluded from the sample .
note however that while the magnetic field is excluded on the spatial scale of the penetration depth the effective exchange field drops to zero within the distance of the coherence length @xmath113 in the superconductor .
more interesting is the behaviour of the spin polarization on the ferromagnetic side of the heterostructure . as it is seen from fig .
[ fig2 ] , it oscillates around its bulk value with period depending on the exchange splitting @xmath112 but not on the superconducting gap .
however we were not able to fit the numerical points to any simple relation between parameters of the model .
the best fit we have is @xmath116 . whilst this means that the period of the oscillations is a linear function of the exchange splitting the above relation is valid only for the thickness of the @xmath1 slab @xmath117 layers . for different thicknesses we have different relations and different numbers of period within the slab .
for example , for @xmath118 , @xmath70 shows one period less for the corresponding exchange splittings @xmath112 .
we now turn to the superconducting proximity effect , namely the leakage of the @xmath0 properties into the non - superconducting region . in our case , although @xmath119 on the ferromagnetic side , due to the fact that for @xmath120 , @xmath121 , the pairing amplitude @xmath122 can , in general , be non - zero and it usually is .
the proximity effect is well understood in the case of non - magnetic metal - superconductor interface ( @xmath123/@xmath0 ) @xcite in terms of the andreev reflections @xcite .
the andreev process describes a situation in which an electron impinging onto the @xmath123/@xmath0 interface is reflected back as a hole ( with opposite spin ) and a cooper pair is created in the @xmath0 .
moreover , the movement of the electrons and holes in the @xmath123 close to the interface ( over the distance of the @xmath0 coherence length @xmath113 ) is highly correlated .
the pairing amplitude , in this case , enters the normal metal and decays exponentially over the @xmath113 .
when the normal metal is replaced by a ferromagnet one might naively expect that the proximity effect will be suppressed due to the pair breaking effect of the exchange field @xcite .
this suggests that the pairing amplitude @xmath124 will decay over a distance much shorter than the superconducting coherence length .
remarkably , as we can see in the fig . [ fig3 ] , this is not the case .    instead , we can observe a very long range proximity effect with an oscillating pairing amplitude .
this effect has been first noted by a. i. buzdin _
_ @xcite and have been studied in several recent publications @xcite . according to them ,
it can be attributed to a @xmath2 - like phenomena @xcite in @xmath1/@xmath0 heterostructures .
the period of the oscillations of the @xmath124 depends on the exchange splitting @xmath112 ( see fig .
[ fig3 ] ) and , similarly to the oscillations of @xmath70 , does not depend on the superconducting energy gap .
clearly this behaviour is different from @xmath123/@xmath0 case , where the proximity effect depends on the @xmath0 coherence length @xmath113 . here
, it depends only on the properties of the ferromagnet . moreover ,
the numerically obtained behaviour of the pairing amplitude @xmath124 is consistent with the analytical formula @xmath125 with @xmath126 being ferromagnetic coherence length . reassuringly
this result is fully consistent with those of ref .
the comparison our numerical results with the formula ( [ chi_analit ] ) for the exchange splitting @xmath127 , which gives @xmath128 is depicted in the inset of the fig .
in fact we had to change @xmath129 by @xmath130 , in order to get better fit .
so the solid line in the inset of the fig .
[ fig3 ] corresponds to the formula @xmath131 with @xmath132 .
it turns out that as exchange splitting becomes smaller the fitting to the analytical formula ( [ chi_analit ] ) becomes better .
the small discrepancy in the coherence length may come from the fact , that we have used bulk value of the exchange field while calculating the coherence length .
however in our system , the spin polarization oscillates around its bulk value , and if we calculate the average value of the spin polarization @xmath133 , we get a little bit smaller value than the bulk one .
so , we conclude , that the effective exchange splitting @xmath134 also has to be smaller than its bulk value .
indeed , when we used @xmath134 to calculate the coherence length , the agreement is improved .
another reason may be the approximate form of the formula for the coherence length in the lattice case @xmath135 .
the proximity effect manifests itself not only in the pairing amplitude and the spin polarization characteristics but also in the density of states ( @xmath73 ) .
the layer resolved @xmath73 is defined as @xmath136 with @xmath137 being the @xmath98-th element of the matrix green s function - solution of the eq .
( [ hfg ] ) .
we shall also use the density of states integrated over the total sample : @xmath138 .
at this point we want to stress that , when we calculate thermodynamic quantities ( eq.([n])-([current ] ) ) we work on the complex energy plane , but calculating the @xmath73 we solve eq .
[ hfg ] on the real energy axis only once .
in fact we have added a small imaginary part ( @xmath139 ) for numerical purposes .
the quantities : @xmath69 , @xmath70 , @xmath71 , @xmath92 as determined by the previous self - consistent procedure allow us to find the green function on the real energy axis in only one step .
we also check if the new quantities like @xmath69 , @xmath70 , etc . agree with the old ones . for all such calculations
the agreement has been found to be excellent .
while the general solution of eq.([hfg ] ) allows for a current flow , in this section we show the results for a solution for which the current has been constrained to be zero .    to appreciate the effect of spin polarization in the @xmath1 slab ,
let us recall again the case of no exchange splitting .
if the non - magnetic metal ( @xmath123 ) is in the contact with superconductor ( @xmath0 ) , the density of states on the @xmath123 side shows features of the @xmath0 @xmath73 , i.e. the pseudo - gap opens up around the fermi energy in the @xmath73 .
however it never turns into a true gap , namely the @xmath123 @xmath73 always possess low energy excitations .
moreover , the gap in the @xmath0 @xmath73 also starts to fill up , due to the proximity of the normal metal .
so in this case we deal with gapless superconductivity close to interface of the @xmath123/@xmath0 heterostructure .
compared to the above summary , the ferromagnet - superconductor interface is enriched with a number of new features .
for instance , the @xmath0 @xmath73 splits , if the exchange field is small . for strong enough exchange field , @xmath0 gap
can be suppressed , indicating transition to the non-@xmath0 ( normal ) phase .
an example of the layer resolved density of states for spin up ( solid line ) and down ( dashed ) is shown in the fig .
[ fig4 ] .
evidently there is a pseudogap in the @xmath1 @xmath73 , indicating that superconducting correlations are present in the @xmath1 layers .
this effect is particularly strong for the layer index @xmath140 and @xmath141 . on the other hand ,
the proximity of the ferromagnet leads to existence of the low energy states in the @xmath0 @xmath73 .
this effect is also observed in @xmath123/@xmath0 heterostructure but here we also notice differences in the spin up and spin down @xmath73 on the @xmath0 layers close to the interface .
these are clear manifestations of the @xmath1 proximity effect .
the most interesting physics of the heterostructures we consider , is the formation of andreev bound states in the @xmath123 or @xmath1 layers .
they are  particle in a box  like states which can be associated with the semiclassical orbits bouncing back and forth between the @xmath0 and @xmath6 regions , as depicted in the fig .
[ fig5 ] . as is well known in the case of @xmath123 layers the reflections at the @xmath0 are andreev ones , while those at the @xmath6 are normal reflections @xcite .
the arrows in the fig .
[ fig5 ] point along the momenta of the particle segment of the orbit .
its @xmath32 component is @xmath36 which labels the state .
clearly , the collection of states labelled in this way forms an andreev band .    for the normal metal - superconductor heterostructure , andreev bound states ( bands )
are symmetrically placed with respect to the fermi energy .
of course the position of these states changes with the thickness of the sample . as we increase the size of the sample , the andreev bound states approach the fermi energy , but they will never reach it . in other words , there is no possibility of crossing the fermi level . somewhat surprisingly , the situation is quite different for the @xmath1/@xmath0 structure @xcite . in this case the bound state energies can cross the fermi level . in short
they can become the zero - energy mid - gap states @xcite .
this situation is illustrated , for @xmath142 , in the fig .
[ fig6 ] , where the bound state energies are plotted against the exchange splitting .
if we take into account states for different @xmath36 ( andreev band ) , the situation is more complicated as in general the states with different symmetry of the wave function can mix , and so we are not able to separate andreev bands so easily .
the example of such andreev bands rather than states is shown in the fig .
[ fig7 ] , where the density of states of the spin up electrons is plotted for a few values of the exchange splitting .
as we shall report presently the fact that the andreev bound states can cross the fermi level has remarkable consequences : it leads to the generation of the spontaneous current in the ground state .
in fact the situation is rather similar to that of exotic superconductors ( with non-@xmath143-wave order parameter ) , where the surface @xmath73 shows the zero - energy state , leading to the spontaneously generated current .
however the origin of this state is quite different , as there it comes from the symmetry related sign change of the order parameter at the surface @xcite . here , in the @xmath1/@xmath0 heterostructure
, the zero - energy state comes from the oscillatory behaviour of the pairing amplitude ( @xmath144 ) in @xmath1 , which changes its phase by @xmath4 each time crosses zero .
nevertheless it leads to the same effect : a spontaneously generated current .
the most remarkable feature of our calculations is that the solution of the @xmath27 equations frequently converges to a solution with the finite current @xmath145 even though the external vector potential is zero .
the generation of the spontaneous current is strictly related to the crossing of the andreev bound states through the fermi energy . in the fig .
[ fig7 ] , we have plotted the @xmath73 for different values of the exchange splitting for calculations where the current was allowed to develop in the @xmath32 - direction . for the @xmath146 we observe that spontaneous current is generated in the ground state .
we have checked that the solution with the current flowing is the true ground state , as it has energy lower than solution without the current .
physically , in the presence of the current , the zero energy states split because of the momentum dependence of a @xmath147 contribution to the quasiparticle energies where the velocity vector is defined by the current @xmath148 .
the example of such a doppler splitting in the @xmath73 is shown in fig .
[ fig8 ] for the exchange field @xmath146 .
we can see that , besides the zero energy state , the other andreev states also are split .
it turns out that the splitting of these states is only weakly sensitive to the exchange field @xmath112 , but it strongly depends on the thickness of the ferromagnetic slab and it decreases as the thickness increases .
interestingly , we have found a simple relation between the splitting of the zero energy state @xmath149 and the vector potential @xmath92 .
this is given by : @xmath150 where @xmath151 is the value of the vector potential @xmath92 averaged over the @xmath1 side of the system only , i.e. @xmath152 .
corresponding to the splitting of the zero energy states , when the current flows , we find differences in the @xmath73 at the fermi level .
in fact , the difference is quite small .
the main behaviour remains the same .
the comparison of the density of states at the fermi energy for both cases ( @xmath153 and @xmath154 ) is shown in the fig .
[ fig9 ] .
the @xmath73 for the @xmath123/@xmath0 structure ( dashed line with circles ) is also drawn .
as we can see , the main features remain more or less the same , in particular the period of the oscillations . similarly like the spin polarization ( fig .
[ fig2 ] ) and pairing amplitude ( fig .
[ fig3 ] ) , the density of states at the fermi level shows decreasing period of the oscillations as the exchange splitting increases .
as we already discussed , the zero energy state is responsible for the generation of the spontaneous current in the system .
the typical example of such a current , flowing parallel to the @xmath1/@xmath0 interface , ( @xmath155 ) is shown in the fig .
[ fig10 ] for a few values of the exchange splitting .
evidently , the behaviour of the current , as a function of the layer index , is very similar to the density of states at the fermi level .
for example , if we compare the dashed curve with the triangles in fig .
[ fig10 ] to the dashed or solid line with triangles in fig .
[ fig9 ] , we readily notice the similarity in the oscillating nature of the current and @xmath73 respectively .
another important issue is the distribution of the current through the whole trilayer structure .
we find that it flows mostly in the positive @xmath32 direction on ferromagnetic side and in the negative direction in the superconductor .
notably the total current , integrated over the whole sample , is equal to zero within numerical accuracy .
this is as it should be for the true ground state and found to be in the @xmath2 state , where the current associated with the unpaired electrons is balanced by the supercurrent flowing in the opposite direction .
while the total current in the system is equal to zero , the current integrated over @xmath1 side only ( @xmath156 ) has a finite value , almost independent of the exchange splitting . but that current is very sensitive to the thickness of the @xmath1 slab , and for example for @xmath157 @xmath158 ( @xmath159 per layer ) , while for @xmath160 we have @xmath161 ( @xmath162 per layer ) .
this suggests that effect associated with the spontaneous current is very important in the samples of the small thicknesses and it seems to play minor role if the size of @xmath1 part of the system is large .
obviously , the above current distribution should lead to the generation of the magnetic flux through the sample .
in fact we have found that again , as the current , the magnetic flux bounded by the layers only weakly depends on the exchange splitting but it does change with the distance from the interface . for @xmath157
the magnetic flux @xmath163 is nearly equal to the half of the flux quantum @xmath164 , @xmath165 .
this gives a flux per plaquette , associated with the layers @xmath34 and @xmath166 , @xmath167 .
on the other hand , for @xmath160 we have @xmath168 and @xmath169 respectively .
the typical layer dependence of the spontaneous magnetic flux @xmath170 ( magnetic flux through a plaquette ) is shown in the fig .
[ fig11 ] .
we can see that while there are small differences in the behaviour on the @xmath1 and @xmath0 sides of the sample , related to the different current distributions ( see the fig .
[ fig10 ] ) , the total magnetic flux @xmath171 is roughly the same and more or less independent of the exchange splitting .
as we have recalled above , in the bulk @xmath2 state the total current vanishes . moreover , the quasiparticle current @xmath172 is exactly equal to the supercurrent @xmath173 ( related to the @xmath174 in eq .
( [ hfg ] ) ) with opposite sign .
this means that in a bulk layered superconductor , these two currents cancel each other layer by layer , i.e. @xmath175 .
evidently this is the case because the exchange field and superconducting pairing potential are present in the whole system .
by contrast , in the case of the @xmath1/@xmath0 heterostructure , the exchange field and the pairing potential are spatially separated . in this case
@xmath176 is not equal to @xmath177 within each layer , but it is when integrated over the whole sample , i. e. @xmath178 thus the lack of exact cancellation of the current layer by layer leads to a finite current on both sides of the interface , but zero in the whole sample .
an example of such behaviour is shown in fig .
[ fig12 ] .
we see , that the quasiparticle part of the total current ( solid line ) is almost in anti - phase to the supercurrent ( dashed line ) .
moreover , on the superconducting side , the current is carried mainly by quasiparticles .
this is rather surprising since usually one expects current to be related to the cooper pairs .
however due to the proximity of the ferromagnet , there are also  normal  particles in the superconductor near the interface .
these certainly contribute to the total current , like in the bulk @xmath2 state . in our system , only the cooper pairs close to the interface can feel the exchange field of the ferromagnet and have moving center of mass . note that magnetization goes to zero very quickly as we go away from the interface ( see fig .
[ fig2 ] ) .
clearly this makes the supercurrent equal to zero except within a few @xmath0 layers .
interestingly we also found that whenever the pairing amplitude @xmath124 changes sign , the local quasiparticle current as well as the local supercurrent becomes large .
in other words , the maxima of the @xmath176 and @xmath179 correspond to the zeros of the pairing amplitude @xmath124 .
as we have mentioned above , the ground state does not carry current for all values of the exchange splitting .
we have shown , that whenever there is a current flowing , there are the andreev bound states crossing the fermi level .
also we have found another quantity which is related to the current flowing .
namely , the pairing amplitude at the @xmath1/@xmath6 interface ( i.e. @xmath180 - see fig .
[ fig1 ] ) .
if @xmath181 changes its sign , the spontaneous current is generated .
this happens exactly for the same values of the exchange splitting , for which a band of the andreev bound states crosses the fermi energy .
moreover , the presence or absence of the current has a dramatic effect on the pairing amplitude @xmath181 .
this can be seen in the fig .
[ fig13 ] , where @xmath181 is plotted as a function of the dimensionless parameter @xmath182 .    as implied by eq.[delta ]
, @xmath124 is the real part of a complex order parameter which has an amplitude @xmath183 and a phase @xmath184 . clearly , fig .
[ fig13 ] can be red as @xmath185 flipping from @xmath186 where @xmath181 has the same sign as @xmath144 in @xmath0 , to @xmath187 , where @xmath181 has the opposite sign .
evidently , the pairing amplitude is pinned to zero by the spontaneous current in the system . introducing the dimensionless quantity @xmath188 we can say , that system is in a @xmath4-like phase ( for odd @xmath188 ) or in the @xmath17-phase ( even @xmath188 ) , in close analogy to studies on the @xmath4-effect in the @xmath0/@xmath1/@xmath0 structure by chtchelkatchev et al .
@xcite . in ref .
@xcite , we have shown a similar picture for the thickness of the @xmath1 slab equal to @xmath189 layers . here we show the corresponding results for @xmath190 layers . to be precise in the present case
@xmath188 has slightly different value with prefactor @xmath191 ( not @xmath192 as it was for @xmath189 layers case ) .
this discrepancy may come from the fact , that effective exchange splitting is little bit smaller than its bulk value ( see discussion in the paragraph [ subsec3.a ] ) .      as our system consists of ferromagnet and superconductor
, we can expect the spontaneous current to be polarized .
however , up to now , we have presented results for the special case of the particle - hole symmetry ( @xmath193 or @xmath194 ) and we found that the spontaneous current was , in fact , not spin - polarized .
the reason is as follows .
within linear response theory , the total current can be divided into two parts : a diamagnetic one - giving the response of the bulk density , and a paramagnetic one , which is proportional to the density of states at the fermi energy , as it comes from the deformation of the wave function at the fermi surface @xcite .
now if we note , that for the particle - hole symmetric case ( @xmath193 ) , the spin up @xmath195 and spin down @xmath196 @xmath73 at the fermi energy are equal , we are not surprised that the polarization of the current is equal to zero . as we go away from the half filling , the difference between spin up and spin down @xmath73 ( @xmath197 starts to play a role as it leads to the polarization of the current .
a typical example of the polarization of the current @xmath198 is shown in the fig .
[ fig14 ] for concentration of electrons @xmath199 and @xmath200 .    to make the two calculations comparable
we have adjusted the parameters so the value of the @xmath0 order parameter would be the same .
this was achieved by shifting of the coupling constant @xmath201 to @xmath202 for @xmath203 , @xmath204 for @xmath205 as well as chemical potential to @xmath206 and @xmath207 respectively .
the total polarization of the current @xmath208 is of order @xmath209 and does not depend much on the concentration of electrons @xmath210 .
moreover , we have found a very interesting property of this current , namely , the sign of the polarization of the current depends on whether @xmath211 or @xmath212 .
this means , that for @xmath213 we get the same results as for @xmath214 , except of the sign of the polarization of the current .
this effect can be explained very easily if we recall that polarization of the current is proportional to the difference in the @xmath73 at the fermi energy , and one can check , that @xmath215 is an antisymmetric function under the shift @xmath216 , i.e. @xmath217 .
clearly , the important issue raised by the above results is `` how to observe experimentally the predicted current ? '' but first of all let s estimate how big this current is . if we assume a bandwidth equal to @xmath218 , the average current per layer is of order of tenths of @xmath219 . in general the polarization of the ferromagnet
can be large , however it can not be equal to @xmath220 .
there must be a small number of the minority electrons .
in other words it does nt matter how big the ratio @xmath221 is , we should always observe the current .
the next point is the thickness of the @xmath1 region .
namely , the current is bigger for a thinner sample , but even for thickness of order of several hundreds of @xmath129 , we expect the current to be easily detectable
.    clearly , the presence of the current should manifest itself in a number of experiments .
for example , since there is a net magnetic flux associated with the current , it could be detected by @xmath222 experiments .
another possibility is to measure the conductance across the @xmath1/@xmath0 junction @xcite . in such experiments
the splitting of the zero energy state in the @xmath73 is expected to be observed and , as indicated by eq.[splitting ] the energy distance between the split peaks measures the current flowing .
finally one can also use some local probes like @xmath223 techniques to measure the layer resolved @xmath73 .    confirming that the current is polarized is more troublesome .
as we have shown , to get polarized current the system has to be away from the @xmath193 point . in our case , to avoid a fermi energy mismatch , both the @xmath1 and @xmath0 are moved from @xmath193 .
fortunately , this is only necessary for the ferromagnet .
the superconductor can be in a particle - hole symmetric state . in short ,
realistic estimations of the spin polarization of the spontaneous current require more material specific calculations .
in conclusion , we have studied properties of the insulator ( or vacuum ) - ferromagnet - superconductor trilayer .
we have shown , that such a structure supports andreev bound states forming andreev bands , whose position can be tuned by thickness of the sample or exchange splitting . moreover ,
when a band crosses the fermi energy , a spontaneous current ( and magnetic flux ) is generated in the ground state . we have found a relation between the pairing amplitude in the ferromagnet , which has oscillatory behaviour , and the andreev bound states .
namely , when the andreev band crosses the fermi energy the pairing amplitude at the surface opposite to the @xmath1/@xmath0 interface changes sign , and as long as current flows it is pinned to zero .
the polarization of the current strongly depends on the band filling and is related to the difference in the spin up and spin down density of states at the fermi level .
we also have discussed the question of the supercurrent vs quasiparticle current and gave some experimental clues on how to observe that current .
the presented results can be , and we have argued that they are , attributed to the @xmath2 effect influenced andreev bound states and can be verified experimentally .
p. m. tedrow , r. meservey , phys . rep . * 238 * , 173 ( 1994 ) ; v. t. petrashov , v. n. antonov , s. v. maximov , r. s. shaikhaidarov , jetp lett . * 59 * , 551 ( 1994 ) ; v. t. petrashov , i. a. sosnin , i. cox , a. parsons , c. troadec , j. low temp
* 118 * , 689 ( 2000 ) ; h. -u .
habermeier , g. cristiani , r. k. kremer , o. lebedev , g. van tendeloo , physica * c364 - 365 * , 298 ( 2001 ) .
n. f. berk , j. r. schrieffer , phys .
* 17 * , 433 ( 1966 ) ; c. pfleiderer , m .
uhlarz , s. m. hayden , r. vollmer , h. v. lohneysen , n. r. bernhoeft , g. g. lonzarich , nature * 412 * , 58 ( 2001 ) ; d. aoki , a. huxley , e. ressouche , d. braithwaite , j. flouquet , j. -p .
brison , e. lhotel , c. paulsen , nature * 413 * , 613 ( 2001 ) .
g. e. w. bauer , yu .
v. nazarov , d. huertas - hernando , a. brataas , k. xia , p. j. kelly , materials sci . eng . *
b84 * , 31 ( 2001 ) ; s. oh , d. youm , m. r. beasley , appl .
. lett . * 71 * , 2376 ( 1997 ) ; l. r. tagirov , phys . rev . lett . *
83 * , 2058 ( 1999 ) .
h. k. wong , b. y. jin , h. q. yang , j. b. ketterson , j. e. hilliard , j. low temp . phys . * 63 * , 307 ( 1986 ) ; j. s. jiang , d davidovi , d. h. reich , c. l. chien , phys .
lett . * 74 * , 314 ( 1995 ) ; v. mercaldo , c. attanasio , c. coccorese , l. maritato , s. l. prischepa , m. salvato , phys . rev . *
b53 * , 14 040 ( 1996 ) ; th .
mhge , n. n. garifyanov , yu . v. goryunov , g. g. khaliullin , l. r. tagirov , k. westerholt , i. a. garifullin , h. zabel , phys .
* 77 * , 1857 ( 1996 ) ; c. l. chien , d. h. reich , j. magn .
magn . mater . *
200 * , 83 ( 1999 ) .
v. proki , a. i. buzdin , l. dobrosavljevi - gruji , phys .
* b59 * , 587 ( 1999 ) ; a. kadigrobov , r. i. shekhter , m. jonson , z. g. ivanov , phys . rev . *
b60 * , 14 593 ( 1999 ) ; l. dobrosavljevi - gruji , r. ziki , z. radovi , physica * c331 * , 254 ( 2000 ) ; m. zareyan , w. belzig , yu . v. nazarov , phys .
lett . * 86 * , 308 ( 2001 ) .


galaxy clusters are the largest known gravitationally bound systems in the universe .
the detailed analysis of the mass distribution of clusters is thus important in the process of understanding the large scale structure , and the nature of dark matter .
the three main methods of measuring galaxy cluster masses : virial masses from velocity dispersions of cluster galaxies , x - ray imaging and spectroscopy of the intra - cluster medium ( icm ) emission , and the gravitational lensing of background galaxies , have been found in recent years to be in disagreement for some clusters . generally , the x - ray estimates are in good agreement with gravitational lensing for clusters with a high concentration of central x - ray emission ( the so - called `` cooling flow '' clusters ) but seemingly in disagreement for other , less centrally peaked objects @xcite . to obtain the estimate of the total mass ( including that due to dark matter ) of a galaxy cluster from its x - ray emission  commonly assumed to be from optically thin , hot plasma that subtends the space between galaxies  it is necessary to make the assumption of hydrostatic equilibrium .
this is appropriate of course only for clusters that have had time to relax into equilibrium and have not experienced any recent merger events .
generally , it is assumed that clusters with circular isophotes meet this criterion .
cold dark matter ( cdm ) hierarchal clustering is the leading theory describing the formation of large scale structure quite well .
in particular , the numerical simulations such as @xcite ( nfw ) successfully reproduce the observed dark matter halo profiles , which appear be independent of the halo mass , initial power spectrum of fluctuations , and cosmological model .
however , observations often disagree with the numerical models .
for instance , one disagreement regards the rotation curves of dwarf elliptical galaxies which appear to be the result of a constant - density core whereas numerical simulations predict cuspy dark matter halo profiles @xcite .
in addition , observations show fewer milky way satellites than predicted by the models @xcite . clearly , to understand the nature of galaxy clusters and the dark matter they consist of
, it is important to measure the matter distribution in clusters via all available means .
fortunately , there are two superb x - ray observatories , chandra and xmm - newton , featuring excellent angular resolution and exceptional effective area coupled with good spectral resolution , and those are well suited for detailed analysis of the x - ray emitting gas of galaxy clusters    abell 1689 is a cluster showing a large discrepancy among various mass determinations , and we chose it for a detailed study .
it is a rich cluster , @xmath5 , without a pronounced cooling flow but with an approximately circular surface brightness distribution suggesting a relatively relaxed structure .
the large mass , @xmath6 @xcite , and apparent symmetric distribution of abell 1689 make it a suitable cluster for gravitational lensing measurements as well as x - ray measurements .
however , the type of clusters that are believed to be the most relaxed have a cool central component of enhanced surface brightness .
the absence of such a component in abell 1689 suggests that the cluster is not fully relaxed .
also , the galaxy content of the cluster is unusually spiral rich for a cluster with high spherical symmetry and richness , with a galaxy type ratio e : s0 : sp of approximately 22 : 22 : 28 plus 25 unidentified galaxies @xcite .
@xcite measure a velocity dispersion of @xmath7 for 66 cluster members , unusually high for a cluster of this temperature . positions and redshifts from @xcite for all cluster members with red magnitude @xmath8 and within @xmath9 of the brightest central galaxy are shown in fig .
[ figoptical ] together with logarithmic x - ray intensity contours from xmm mos 1 .
the gravitational lensing estimate from 6000 blue arcs and arclets calibrated by giant arcs at the einstein radius of abell 1689 gives a best fit power - law exponent of @xmath10 for the projected density profile from @xmath11 to @xmath12 @xcite .
( unless otherwise stated , we assume an einstein - desitter ( eds ) cosmology with @xmath13 and @xmath14 . ) this is steeper than the profile of an isothermal sphere @xmath15 .
the strong lensing analysis of two giant arcs directly gives @xmath16 @xcite .
the mass profile derived from the deficit of lensed red galaxies behind the cluster due to magnification and deflection of background galaxies suggests a projected mass profile of @xmath17 for @xmath18 , close to that of an isothermal sphere @xmath19 .
the mass interior to @xmath20 from this method is @xmath21 @xcite .
measurement of the distortion of the luminosity function due to gravitational lens magnification of background galaxies gives @xmath22 @xcite .
finally , the weak gravitational shear of galaxies in a eso / mpg wide field imager @xmath23 image gives a mass profile with best fit nfw profile with @xmath24 and @xmath25 or a best fit sis with @xmath26 @xcite .
there is a good indication from the optical data that the cluster consists of substructures .
@xcite suggest a strong lensing model with two clumps in order to reproduce the positions of the brightest arcs .
a larger mass clump @xmath27 is centered on the brightest cluster galaxy while a smaller clump @xmath28 is located @xmath29 northeast of the main clump .
they arrive at a mass a factor 2 - 2.5 lower for their x - ray estimate compared to their gravitational lensing model .
@xcite identify two distinct substructures centered on redshifts @xmath30 and @xmath31 using positional and redshift data of cluster galaxies from @xcite , providing further evidence that the cluster is not relaxed .
these clumps are also aligned in the southwest  northeast direction but the locations do not agree with the ones of @xcite .
both structures are found to have @xmath32 yielding virial masses several times smaller than those derived from lensing and x - ray estimates .
can the x - ray observations provide any evidence for a substructure in abell 1689 , and what are the implications on the mass inferred from the x - ray data ? in an attempt to answer this ,
we analyze xmm - newton epic pn and epic mos data to measure the mass profile of a1689 and to investigate the spatial structure of the cluster .
2 contains the details of the observations and data reduction with the xmm - newton as well as with summary of the rosat , asca and chandra data ;  3 covers the methods of spectral fitting of the xmm data and presents the analysis of cluster asymmetry ;  4 derives the mass and the slope of the mass distribution in the core ; and  5 presents the inferences about the structure of the cluster inferred from the spatial analysis .
the paper concludes with the summary in  6 .
abell 1689 was observed with xmm - newton for 40 ks on december 24th 2001 during revolution 374 .
for imaging spectroscopy we use data from the european photon imaging camera ( epic ) detectors mos1 , mos2 and pn .
both mos cameras were operating in the full frame mode whereas pn was using the extended full frame mode .
the extended full frame mode for pn is appropriate for studying diffuse sources since it has lower time resolution and so it is less sensitive to contamination from photons being detected during readout of the ccds .
these events ( so called out - of - time events ) show up as streaks across the x - ray image and are especially bothersome when the goal of an observation is spatially resolved spectroscopy of diffuse sources .
all cameras used the thin filter during the observation .
epic background is comprised mainly of three components .
the external particle background consists primarily of soft protons @xmath33 being funneled through the mirrors and causing a time variable flaring signal in the detector .
the internal particle background is mainly due to high energy particles interacting with the detector material and causing a roughly flat spectrum with flourescent emission - lines characteristic to the detector material .
this component varies over the detector surface .
the third source of background is the cosmic x - ray background ( cxb ) which is roughly constant in time but varies over the sky .
for all data reduction we use the software and calibration data implemented in xmm science analysis software ( sas ) 5.4.1 . to exclude the events contaminated by proton flares , we produce light curves in the @xmath34kev band where the true x - ray signal is low .
we screen the data using a constraint on the total count rate of less than 1.5 ct s@xmath35 for mos and 1.1 ct s@xmath35 for pn in this band , leaving an effective exposure time of 37 ks for mos and 29 ks for pn .
this screening corresponds to a limit on the count rate of approximately @xmath36 above the quiescent period in the @xmath37kev band .
the effective mirror collecting area of xmm - newton is not constant across the field of view : it decreases with increasing off - axis angle and this decrease is energy dependent .
this results in a position dependent decrease in the fraction of detected events and when doing imaging spectroscopy for extended sources , we need to account for this effect . by generating an ancillary response
file ( arf ) for each source spectrum region , using xmm sas 5.4.1 command ` arfgen ` we calculate an average effective area for each region considered by us ( see below ) .      in order to correctly account for particle - induced and cosmic x - ray background
, it would be optimal to extract a background spectrum from the same detector region collected at the same time as the source spectrum .
this is of course impossible , and the background can be taken from a source - free region of the detector ( other than the target , but near it ) , or can be estimated using blank sky data .
we adopt the former method , noting that the background is not entirely constant over the field of view ; however , it can be assumed to be approximately constant with the exception of the fluorescent cu line at @xmath38 in pn which is the strongest contaminant emission line . using this assumption
, we can effectively subtract the particle background by extracting a spectrum in a source free region in the same exposure .
for pn data in the @xmath39 range is excluded due to the strong spatial dependence of the cu internal emission .
using an arf generated for a source region on a background subtracted spectrum will not take into account that the spectrum used as background was extracted from a different region where the vignetting was higher , due to the larger off - axis distance of the background region .
this will have the effect that the cxb component in our source spectrum will be under - subtracted and the net spectrum will contain some remnant of the cxb .
the particle induced background however is not vignetted and therefore should leave no remnant .    to estimate the cxb component in our exposure we take the events outside of the field - of - view to be the particle background .
another spectrum is then extracted from a large source free area , located away from the cluster .
the particle background is normalized and removed from this spectrum .
the resulting spectrum is fitted to a broken power law model to determine the cxb flux .
the incorrect vignetting correction of the cxb is found to cause an at most 2 % over - estimation of the flux which is the case in the outermost annulus in our analysis ( see below ) .
however since the vignetting is energy - dependent , the incorrect vignetting may cause a small ( @xmath40 kev ) shift in temperature for the outermost region .
the overall uncertainty in the background is estimated to be at most a few percent .
we compare our background subtraction method with the method of using xmm standard blank sky data compiled by @xcite ; we find that both methods give similar fit parameters .
however , given the difficulty of normalizing the cxb , the different high - energy leftovers from proton - flare subtraction and the different internal particle background which occurs when using background data from a different exposure , we chose to use in - field ( rather than blank - sky ) background .
we find this is more robust in keeping the overall shape of the source spectrum uncontaminated . in all our analysis
we use a background region of a circular annulus with inner and outer radii of @xmath41 and @xmath42 respectively .
in fact , recent work by @xcite suggests that the in - field background method is probably more accurate .
the clusters in the above paper , however , have smaller apparent angular size which makes the background method more reliable .      in analyzing the vignetting corrected , background subtracted radial count - rate profile , we fit it to a conventional beta model , @xmath43 , where @xmath44 is the source surface brightness at radius @xmath45 .
the fit gives @xmath46 and @xmath47 with @xmath48 using data out to @xmath49 .
clearly this model is not a very good fit ; we show it here only for completeness and comparison with previous work , and note that it is not used in the subsequent analysis .
the bad fit above results from the fact that the cluster emission is more peaked in the core than the best fit beta model . to obtain a general idea of the properties of the cluster and compare this with previous results
, we extracted spectra for the central @xmath50 @xmath51 region , centered on the x - ray centroid at @xmath52 @xmath53 .
this radius corresponds to @xmath54 or @xmath55 .
background spectra were extracted from source - free regions from the same exposure .
for xmm pn we use single and double pixel events only whereas for xmm mos we also use triple and quadruple pixel events .    for spectral fitting
we use the xspec @xcite software package .
we fit the data in the @xmath56 range using the mekal @xcite model for the optically thin plasma and galactic absorption . with the absorption fixed at the galactic value , @xmath57 @xmath58 @xcite , and assuming a redshift @xmath59 we arrive at a temperature of @xmath60 and a metal abundance of @xmath61 solar for both mos cameras with @xmath62 . from the pn camera
we get @xmath63 and @xmath64 solar with @xmath65 .
the best fit models with residuals can be seen in fig .
[ mospn ] for mos ( left ) and pn ( right ) .
the pn temperature is in disagreement with mos data and the reason for this effect can be seen from the residuals below 1 kev for pn ( fig .
[ mospn ] ( right ) ) . fitting the pn data above 1 kev we get @xmath66 and @xmath67 solar with @xmath68 , in better agreement with mos .
the temperature from mos is in agreement with that found by asca , @xmath69 @xcite .
this observation does not suffer from pile - up , nor is the low energy spectrum sensitive to background subtraction .
the background uncertainties below 1 kev for this pn spectrum are less than 1 % whereas the pn soft excess is sometimes higher than 10 % .
the excess is certainly not background related .
it is possible that the pn low - energy discrepancy can be due to incorrect treatment of charge collection at lower energies ( s. snowden , priv .
comm . ) . to resolve the discrepancy between mos and pn we decided to compare with asca gis / sis , rosat pspc and chandra acis - i data .
besides the discrepancy regarding the softest x - ray band for the abell 1689 data between the mos and pn detectors aboard xmm - newton , the spectral fits to the chandra data presented by @xcite imply a higher absorbing column , @xmath70 @xmath58 than the galactic value of @xmath71 @xmath58 @xcite . since such excess absorption is not commonly detected in x - ray data for clusters , this requires further investigation . to determine if there is indeed any additional component of absorption beyond that attributable to the galactic column  and assess the reliability of the softest energy band of the pn vs. mos data
we used the most sensitive soft ( @xmath72 kev ) x - ray data for this cluster obtained prior to the xmm - newton observations , collected by the rosat pspc . the rosat observation conducted during july 18 - 24 , 1992 ( available from heasarc ) yielded 13.5 ks of good data .
we extracted the rosat pspc counts from a region @xmath73 in radius , centered on the nominal center of x - ray emission . for background
, we selected a source - free region of the same pspc image . using these data over the nominal energy range 0.15 - 2.1 kev with the standard pspc response matrix applicable to the observation epoch
, we performed a spectral fit to a simple , single - temperature mekal model with soft x - ray absorption due to gas with solar abundances , using the xspec package as above . in the fit we use metal abundances of 0.27 solar
obtained in the xmm mos fit above .
while the limited bandpass of rosat pspc precludes an accurate determination of the temperature ( the best fit value is @xmath74 kev , 90% confidence regions quoted ) , the pspc data provide a good measure of the absorbing column : the best value is @xmath75 @xmath58 , certainly consistent with the galactic value .
we note that the difference between the temperature inferred from the pspc fit and that obtained from the xmm - newton data as above is a result of the limited bandpass of the pspc , located much below the peak of the energy distribution of the cluster photons .
the measurement of the absorbing column , however , clearly indicates that the column inferred from the chandra observation by @xcite is not correct , and might be due to instrumental effects .
we conclude that the absorbing column is consistent with the galactic value .    to obtain further constraints on the absorbing column
, we also used the asca gis and sis data together with the pspc data for an independent constraint on the continuum radiation in the fitting process .
we performed standard extraction of data from all asca detectors , also from a region from a region @xmath73 in radius for the source , and a source - free region of the same image for background .
we performed a spectral fit simultaneously to data from the pspc and four asca detectors . to account for possible flux calibration differences ,
we let the normalization among all the different detectors run free .
we used the energy range of @xmath76 and @xmath77 for the asca gis and sis cameras respectively .
since asca detectors ( and in particular , the siss ) often return spectral fits with excess intrinsic absorption @xcite we also let the absorbing column be fitted independently for the gis , sis and pspc detectors .
temperature and metal abundances were tied together for all datasets in the fit .
the optical redshift @xmath0 was used .
the joint fit of rosat and asca data gives us the best fit temperature of @xmath78 kev and abundances of @xmath79 solar , in agreement with the values quoted by @xcite , and the absorption ( for the pspc data ) @xmath80 @xmath58 , in agreement with the galactic hi 21 cm data .
we shall use this value for the absorption in the subsequent analysis .
this best fit rosat / asca model was compared with the data from the same region in the xmm - newton cameras giving an unfitted reduced @xmath81 of 1.39 for pn and 1.06 for both mos cameras combined .
the ratio of these spectra to the rosat / asca model can be seen in fig .
[ mospndisc ] . from this result
we conclude that mos low energy response is more consistent with previous data and subsequently , we choose to ignore all pn data below @xmath82 .
re - fitting the xmm data from the above region using @xmath83 for mos data and @xmath84 for pn leaving the absorbing column as a free parameter yields @xmath85 @xmath58 , a temperature of @xmath86 and metallicity of @xmath87 solar . while this fitted value of absorption is formally inconsistent with the galactic and rosat - inferred values , this is a relatively small difference , which might be due to the slightly imperfect calibration of the mos detectors or the assumption of isothermality made by us for this fit ( since @xmath88 and @xmath89 are correlated in the fitting procedure ) .
we note that using the rosat value for absorption will give us a somewhat lower measure on the temperature ( see below ) .
finally , we reduced the chandra data for abell 1689 using the most recent release of data reduction software ; specifically , due to the degradation of the chandra acis - i low - energy response correction for the charge transfer inefficiency ( cti ) is necessary , and we applied this to the chandra data .
as of chandra data analysis software package ciao ver .
2.3 this correction can be applied in the standard event processing .
it is also necessary to account for the acis excess low - energy absorption due to hydrocarbon contamination .
we used the ` acisabs ` code for the correction to the auxiliary response function .
the event grades used in the acis analysis were grade=0,2,3,4 and 6 .
the cluster was observed with the chandra acis - i detector array at two separate occasions for @xmath90ks each on 2000 - 04 - 15 and 2001 - 01 - 07 .
spectra were extracted for the central @xmath50 region centered on the x - ray centroid at @xmath52 @xmath53 following the ciao 2.3 science threads for extended sources .
[ chandra ] shows the ratio of the combined chandra data to the best fit model determined from rosat and asca above .
the left spectrum is before and the right after the ` acisabs ` correction .
fitting the corrected spectrum to an absorbed mekal model in the energy range 0.3 - 7.0 kev yields absorption of @xmath91 @xmath58 , a temperature of @xmath92 kev and abundances of @xmath93 solar using the optical redshift z=0.183 .
the discrepant absorption found by @xcite ( using data 0.7 - 9.0 kev ) is apparently corrected for by ` acisabs ` and the value is in agreement galactic absorption .
however there is still a large discrepancy in temperature , which also can be seen from the high energy ends of fig .
[ chandra ] .
part of this effect can possibly be attributed to the high energy particle background but more likely to uncorrected instrumental effects .
repeating all above steps for single pixel events only ( grade=0 ) in order to achieve higher spectral accuracy gives us a best fit temperature of @xmath94 kev .
we can not account for the differences between the two chandra data sets ( using grade=0 vs. grade=0,2,3,4 and 6 events ) nor between the results of the chandra and xmm spectral fits .
we note here that the photon statistics resulting from the xmm observation is superior to that in the chandra data , and since our analysis does not require the superior angular resolution of the chandra mirror , we limit the analysis below to the xmm data .
to obtain a radial profile of cluster gas temperature , abundance , and density , we first make the assumption that the cluster is spherical and that above properties are only functions of radius . for this
, we divide the image of the cluster into 11 concentric annuli out to @xmath95 @xmath96 centered on the x - ray centroid . for each annulus , we extract spectra from all epic cameras , and we set the inner and outer radii of each region by requiring that each annulus contains at least 9000 counts per each mos camera and 13000 for pn .
this allows us to derive a reliable estimate of temperature in each region .
point sources with intensities greater than @xmath97 over the average are excluded .
the outer radii of the annuli are as follows : @xmath98 , @xmath99 , @xmath100 , @xmath101 , @xmath102 , @xmath103 , @xmath104 , @xmath105 , @xmath106 , @xmath107 and @xmath108 .
average cluster properties were determined using all annular spectra simultaneously out to @xmath109 ( @xmath110 , @xmath111 ) using the same spectral fitting procedures as in section [ fitting ] .
we use the optical redshift of @xmath59 @xcite and the line of sight absorption of @xmath112 as measured from rosat data , and also consistent with the galactic column density ( section [ rosasc ] ) .
the mean temperature of the cluster is found to be @xmath113 and the mean abundance @xmath114 solar .
we note that leaving the redshift as a free parameter gives a best fit redshift @xmath115 ( 90% confidence range ) which is considerably less than measured via optical observations . considering mos and pn data separately gives @xmath116 and @xmath117 respectively .    in order to take into account the three - dimensional nature of the cluster we consider the spectrum from each of the 11 annuli to be a superposition of spectra from a number of concentric spherical shells intersected by the same annulus .
the spherical shells have the same spherical radii as the projected radii of the annuli .
we assume that each shell has a constant temperature , gas density and abundance .
the volume for each annulus / shell intersection is calculated to determine how large a fraction of emission from each shell should be attributed to each annulus .
assigning a spectral model to each spherical shell we can simultaneously fit the properties of all spherical shells .    in practice
we use the method of @xcite where we have a matrix of @xmath118 mekal models with absorption in xspec . in the process of fitting the data in xspec , each datagroup ( consisting of three data files : pn , mos1 and mos2 spectra for each annulus )
is fitted using the same set of models . for the central annulus ( the one with zero inner radius
so it s actually a circle ) , which will intersect all 11 spherical shells , we will need to apply 11 mekal models to represent these .
this means that we have to apply 11 mekal models to each datagroup ( annulus ) with absorption where each model represents the properties of one spherical shell .
the normalization of each shell model is set to be the ratio of the volume that shell occupies in the cylinder that is the annulus / shell intersection to the volume it occupies in the central annulus .
abundance and temperature are tied together for the models representing the same spherical shell .
of course not all annuli intersect every shell , and for those shells not intersected by the annulus to which they are attributed , the model normalization will be zero .
the matrix of mekal models is thus triangular and can be fitted directly to the spectra we have extracted from annular regions in the data .
this approach allows for all data to be fitted simultaneously , and we do not have the problem with error propagation which occurs when starting to fit the outermost annulus and propagating inward subtracting contributions from each previous shell .    the metallicity profile ( fig .
[ abun ] ) from the deprojection shows signs of increasing abundance toward the center of the cluster . in temperature
[ temps ] ) we find an apparent decrease for large radii ( @xmath119 ) , an effect that has been seen in analysis of other clusters with xmm ( see e.g. @xcite ) .
gas dynamic simulations of the formation of galaxy clusters also show a decline of temperature at large radii @xcite .
we do not find a significant cooling in the cluster core with the highest temperature ( @xmath120 ) near the core radius : in fact , we will show in section [ xyana ] that the temperature profile is not symmetric around the cluster center .
we note that for completeness , we also performed the above analysis with the best - fit value of the redshift inferred from the x - ray data alone , and while the exact values of temperature and elemental abundances are slightly different , about 0.2 kev lower for temperature , the general trends in the radial runs of the parameters are the same .
the limited point spread function ( psf ) of the xmm mirrors is a potential problem especially for the annuli located close to the center since those are not much larger than the psf fwhm of @xmath121 .
some of the flux incident on the central ( circular ) region will be distributed over the outer annuli and vice versa .
this flux redistribution will have the effect of smoothing out the measured temperature profile since all annuli will have some flux that originally belong in other annuli .
this effect has been studied in detail by e.g. @xcite who find that correcting for the psf redistribution gives a profile that is consistent with an uncorrected profile .
abell 1689 has a temperature profile without large temperature variations and no large central flux concentration .
we conclude that the effect of flux redistribution in our case will be small and we do not attempt to correct for this . the psf is also weakly energy dependent , and to quantify its possible effect on the observed temperature profile , we calculate the energy dependent flux loss from the central annulus and the effect on the central temperature .
the difference in flux loss between various energy bands ( ranging from 1.5 to 7.5 kev ) for the on - axis psf is @xmath122 .
we find that for a cluster with an assumed temperature of 9 kev , this could give an error of the central temperature by at the most @xmath123 .
we note that this is a maximum difference since in practice , the flux gained from outer annuli could somewhat reduce this effect by working in the opposite manner .
the luminosity of the cluster in the @xmath124 band , calculated from the best fit model above ( with @xmath59 ) gives @xmath125 for an eds ( @xmath126 , @xmath127 ) universe or @xmath128 for a cdm ( @xmath129 , @xmath130 ) universe .
this corresponds to a bolometric luminosity of @xmath131 or @xmath132 .
all above values should include 10% as the absolute calibration error of xmm . from chandra analysis ,
@xcite find @xmath133 whereas @xcite find @xmath134 from asca data , both in agreement with our results .
this does not provide any new information regarding the location of abell 1689 in the luminosity - temperature @xcite relation and it is still in a close agreement with the trend suggested by other clusters .
llll x - ray ( rosat ) & @xmath135 & @xmath136 & 1 + x - ray ( chandra ) & @xmath137 & @xmath138 & 2 + lensing & @xmath139 & @xmath140 & 1 + optical & @xmath141 & @xmath142 & 3 + x - ray ( xmm ) & @xmath52 & @xmath53 & 4 +      with the good quality xmm data , it is possible to verify the result of @xcite that there is no discrepancy between the optical and x - ray centers of the cluster .
we determine the center of x - ray emission for abell 1689 using ximage command ` centroid ` .
we also include a measurement of the lensing center from @xcite and the x - ray measurement from rosat by @xcite ( table [ centers ] ) .
we find that all values agree within @xmath143 except the rosat estimate , the offset of which we attribute to uncertainties in rosat hri astrometry .
this apparently perfect agreement among x - ray , lensing , and optical centers leads us to conclude that the icm density peak and the central dominant galaxy is probably located at the bottom of the dark matter potential well . still , the apparent non - uniformity of the icm radial temperature distribution as well as the offset of optical and x - ray redshifts prompted us to analyze the spatial structure of the cluster . in fig . [ btemps ]
we show the spatial temperature distribution .
spectra were extracted in rectangular regions and fitted using the same method as in section [ fitting ] .
the temperature in the figure scales linearly from 6 kev ( black ) to 10 kev ( white ) .
the errors on the temperature of the inner 16 regions are @xmath144 , while the errors on the outer 8 regions are @xmath145 .
we clearly see an asymmetry in the temperature around the cluster center with an overall increase toward the northeast .    to check the consistency of these results and to increase our accuracy , we re - group the data in larger spatial regions and perform a fit using the same above procedure .
we first fit the data keeping the redshift frozen at the optical value .
temperatures derived from this fit are shown in fig .
[ urdl ] ( left ) .
thereafter , we leave the redshift as a free parameter and re - fit the data .
the fitted redshifts for these regions are shown in fig .
[ urdl ] ( right ) .
all errors in fig .
[ urdl ] are 90% confidence limits . in the temperature map ,
we see a clear discrepancy between the northern and southern part of the cluster with a hint of a temperature gradient in the southwest  northeast direction .
the redshift map reveals a high redshift structure to the east at @xmath146 separated from the rest of the cluster at @xmath147 . analyzing the mos and
pn data separately for this high redshift region gives a broad minimum in @xmath81 at @xmath148 for mos whereas pn shows several minima in the @xmath149 range .
this region may indicate a subcluster falling inward away from the observer .
this is further supported by the optical data , indicating that there are also high - redshift giant elliptical galaxies in this region ( fig .
[ figoptical ] ) : it is interesting to speculate if this is actually the remains of a cluster core ? especially intriguing is the apparent coincidence between the smaller subcluster as suggested by strong lensing @xcite and our high redshift gas region approximately @xmath150 northeast of the main cluster .
another possibility is that the redshift variation is due to large bulk motions of the intra - cluster gas .
it has been shown in gas dynamic simulations that clusters with apparently relaxed x - ray profiles can have complex gas - velocity fields and be far from relaxed @xcite .
this kind of motion could give rise to non - thermal emission from shocks etc .
in our analysis we can not distinguish between the two possibilities of bulk motion and subclustering .
this measurement of non - uniform redshift distribution of the x - ray emitting gas is important , and to verify if it could be due to instrumental effects , we investigated this in more detail .
the data in the regions above are from different ccd chips in the pn data whereas for mos all data are from the same chip .
since the pn camera provides about half of the data , we want to verify that there are no gain shifts between the ccds , which , if present , could easily cause such an effect .
most of the cluster emission is on the ccd chips 4 and 7 . to test for any possible gain offset
, we extracted the data from each of the chips individually to verify the position of the internal fluorescent @xmath151 line , with the dominant component at @xmath152 .
the spectra in the range @xmath39 are fitted to a gaussian profile : both datasets yield essentially the same peak position of @xmath153 .
this corresponds to a possible artificial redshift offset of maximum @xmath154 .
to explain the difference in the offset measured by us as an instrumental effect , we would have to have an offset ( say at the fe k line at 6.7 kev ) of 62 ev , and not on the order of 1 ev , as inferred from the cu k instrumental line .
hence we conclude that there is no gain shift that could alter our redshift measurements between pn ccds 4 and 7 .
according to xmm calibration documentation @xcite the magnitude of calibration errors for pn & mos should be no larger than 10 ev .
if we assume that the cluster is spherical with a smooth static gravitational potential and that the intra - cluster medium is a pressure - supported plasma , we can employ the hydrostatic equilibrium equation .
the circular x - ray isophotes ( fig .
[ figoptical ] ) generally indicate that a cluster is in dynamical equilibrium .
the hydrostatic equation can be written as @xcite : @xmath155 where @xmath156 is the enclosed total gravitating mass enclosed within a sphere of a radius @xmath45 , @xmath157 and @xmath158 are temperature and density of the icm at @xmath45 , @xmath159 is the mean particle weight and @xmath160 is the proton mass .
using the temperature and normalizations from the spectral deprojection fitting we calculate the total gravitating mass .
errors are treated by error propagation .
the mass distribution is fitted to a singular isothermal sphere ( sis ) @xmath161 where @xmath162 is the 1-dimensional velocity dispersion , which is used here for a comparison with previous mass estimates .
the predicted density profile from cdm hierarchal clustering according to @xcite for dark matter halos is @xmath163 where @xmath164 is the critical energy density at halo redshift @xmath165 and @xmath166 is characteristic density defined by @xmath167 } \label{eqnfw2}\ ] ] where @xmath168 is the concentration of the halo defined as the ratio of the virial radius @xmath169 to @xmath170 , which in turn is a characteristic radius in the nfw model .
the critical density at redshift @xmath165 for a flat ( @xmath171 ) universe is @xmath172 \label{rhocrit}\ ] ] where @xmath173 is the hubble constant , @xmath174 and @xmath175 are the current contributions of matter and vacuum energy respectively to the energy density of the universe .
more recent numerical studies suggest a steeper core slope and a sharper turn - over from small to large radii @xcite .
both models can be generalized as @xmath176^{(\beta - \gamma)/\alpha } } \label{eqgeneral}\ ] ] @xcite , where @xmath177 and @xmath178 characterize the density slope at small and large radii respectively whereas @xmath179 determines the sharpness in the turn - over .
most studies agree on @xmath180 but the value of @xmath177 is still being debated .
the nfw and moore profiles fit into the parameter space @xmath181 as @xmath182 and @xmath183 respectively .
we choose to fit our data to the nfw model ( eq . [ eqnfw1 ] ) which , when integrated over @xmath45 , yields @xmath184 \label{massnfw}\ ] ] where @xmath185 .
we find that the data give the best fit for the nfw model with @xmath1 and @xmath186 whereas the sis fit gives @xmath187 . the total mass data and models are shown in fig .
[ mass3d ] together with the mass of the x - ray emitting gas @xmath188 .
model parameters are summarized in table [ bestfit ] . for a cosmology with @xmath129 and @xmath130 , the best fit nfw model changes to @xmath189 and @xmath190 .
another cosmologically important quantity for clusters is the fraction of the total mass that is in the x - ray emitting gas , @xmath191 .
@xcite analyze data from the chandra observatory for 10 dynamically relaxed clusters between @xmath192 and @xmath193 and measure an average redshift independent @xmath194 at the @xmath195 radius ( where the total mass density is 2500 times the critical density at the redshift of the cluster ) .
the cosmology where this is valid is a flat @xmath196 cosmology with @xmath197 assuming @xmath198 and @xmath199 .
the error on @xmath200 above is the rms dispersion of the @xcite sample which is comparable to the individual errors on @xmath200 . with our best fit nfw model in the above cosmology
we find @xmath201 for abell 1689 at the @xmath195 radius .
this is lower than all the clusters in the @xcite sample and significantly lower than the mean .
we find that in our estimate @xmath200 has not converged to a constant and this might help explain part of the discrepancy .
however , it may be the case that for many clusters @xmath200 does not converge until well beyond @xmath195 .
cccc @xmath202 & @xmath203 & @xmath204 & 7.64 ( 8) + @xmath205 & @xmath203 & @xmath206 & 7.64 ( 8) + sis & @xmath203 & @xmath207 & 61.4 ( 9 ) + pow & @xmath208 & @xmath209 & + pow & @xmath210 & @xmath211 & +    comparing the mass and temperature at @xmath195 of abell 1689 to the m - t relation derived for a set of relaxed clusters @xcite shows a low mass for the temperature of abell 1689 .
the m - t relation predicts @xmath212 for a @xmath213 cluster where @xmath214 is the hubble constant at the redshift of the cluster . for abell 1689
we find @xmath215 , significantly lower .
the above values were derived assuming a flat @xmath196 cosmology with @xmath129 , @xmath130 and @xmath216 .
the unusually low mass may be due to the fact that the mass of abell 1689 seems to increase steadily beyond @xmath195 .
however , we note that a lower mass ( than would be predicted by the m - t relation ) is not an uncommon feature for unrelaxed clusters ( cf .
@xcite ) .    for completeness
, we note that the calculated total mass includes the icm and galaxy mass contributions as well as the dark matter .
the proper way of fitting the nfw model would be to subtract these contributions prior to performing the fit .
the nfw model for the actual dark matter halo is not used in this paper ; we only calculate the total mass profile .
the mass data were fitted to a simple power law ( @xmath217 ) in the ranges @xmath218 and @xmath219 .
we find the best fit of the slope of the matter profile to to be @xmath220 and @xmath221 for small and large radii , respectively .
this corresponds to total mass density slopes ( @xmath222 ) of @xmath223 and @xmath224 , in good agreement with what is expected from numerical simulations of cdm hierarchal clustering .
we note here that @xcite have measured the density profile of abell 1689 for @xmath225 using chandra data and found @xmath226 .
we do not observe a flattening of the core density profile but find a slope close to the core to be somewhere in between the preferred moore and nfw profiles .
thus we do not claim to be able to detect nor dismiss any kind of modification to standard cold dark matter . nonetheless , the x - ray data imply an upper limit on the self - interaction of dark matter ; see , e.g. , the discussion in @xcite .
such comparisons are more meaningful for cooling - flow clusters which are presumably more relaxed objects .
to compare our results with those obtained from gravitational lensing , we reprojected our derived mass distribution into a two dimensional distribution by summing up the contributions from each shell along the line - of - sight .
of course this method assumes that the outermost shell is the absolute limit of the cluster , and we recognize that this will not be entirely accurate .
therefore we also include our projected best - fit nfw model which has an analytic expression . reprojected mass and nfw model are shown in fig .
[ mass2d ] where pluses with error bars are the reprojected mass , the solid line is our nfw model , the triangle is the strong lensing result @xcite , the asterisks are lensing magnification results measured by the distortion of the background galaxy luminosity function @xcite , the dot - dashed lines are lensing magnification results measured by the deficit of red background galaxies @xcite , and the dashed line is the best fit nfw model from weak gravitational shear analysis @xcite .
the comparison of our results with measurements from gravitational lensing is summarized in table [ oldres ] .
the x - ray mass appears to be in good agreement with that derived from weak gravitational shear but can not be reconciled with the strong lensing data nor with the data from gravitational lens magnification .
this discrepancy in lensing was noted earlier by @xcite as well as @xcite who can not find an agreement using any realistic corrections . the high velocity dispersion measured by @xcite and the apparent grouping of galaxies along the line of sight @xcite may indicate that this cluster is not as regular as we expect from its smooth circular surface intensity .
the two component model of @xcite also strongly suggests non - uniformity .
this might explain the discrepancies between mass estimates from gravitational lensing and from x - rays .
a possibility could be that this is a cluster undergoing a major merger , where a sub - clump close to or along the line - of - sight is being stripped from gas or just has very low x - ray luminosity .
cccc @xmath227 & @xmath228 & 1 & strong lensing + @xmath229 & @xmath228 & 4 & x - ray + @xmath230 & @xmath231 & 2 & weak lensing magnification + @xmath232 & @xmath233 & 3 & weak lensing magnification + @xmath234 & @xmath233 & 4 & x - ray +      to illustrate the possible explanation that a merger might be taking place in this cluster and motivate why the x - ray derived mass should be lower than expected we employ a simple model .
we assume that we have two perfectly spherical clusters aligned exactly along the line of sight having identical mass ( @xmath235 ) and gas density ( @xmath236 ) . since the thermal bremsstrahlung emissivity scales as @xmath237 we will measure a surface intensity @xmath238 and from this infer a gas density @xmath239 .
in this example however we have @xmath240 and hence we will conclude that @xmath241 .
this means that we underestimate the total mass of the x - ray emitting gas by a factor @xmath242 .    from the hydrostatic equation ( eq . [ mass ] ) we can write @xmath243 . from the gradient of the measured surface brightness @xmath244
we infer a gradient on the gas density . in our scenario
it is true that @xmath245 which yields @xmath246 . for the total mass we will measure @xmath247 .
the actual mass of the cluster pair is of course @xmath248 which is twice what we measure .
we have underestimated the total gravitating mass by a factor @xmath249 .
this is in agreement what we find in our comparison with gravitational lensing derived masses which would be able to measure the total mass accurately in this example .    while this exact scenario is not very probable it shows that a close configuration of clusters will underestimate the x - ray mass , maybe by a factor as large as 2 . in the above example the gas mass fraction
would actually be overestimated which is the opposite of what we find .
however if a merger was in fact taking place the hydrostatic equation would probably be quite inaccurate estimation of the mass .
there would be other sources of pressure and probably the cluster would not be in equilibrium : examples here might be magnetic fields , or additional pressure support from non - thermal particles , likely to be accelerated in shocks that arise during a merger . also since we do not detect two separate peaks in the surface intensity map , it is much more likely we have a lower density companion cluster or one merging irregular system .
there is a good evidence suggesting that abell 1689 is undergoing a merger .
we find that the x - ray measurement yields the value of mass that is low in comparison to gravitational lensing and the m - t relation for relaxed clusters , both by about a factor 2 .
we have shown that this may be the result of grouping or elongation along the line of sight .
the cluster has a low gas mass fraction compared to other clusters , which is perhaps the result of large scale gas motion .
the velocity dispersion of the member galaxies is very high and there are hints of subgroupings in the redshift space .
the high number of spiral galaxies is unusual for a rich cluster and suggests a lower density region , where the spiral fraction would be higher , perhaps in front of or behind the main cluster .
a state of complex dynamics is supported by the non - uniform temperature distribution and maybe most importantly by the variation in redshift of the x - ray emitting gas across the cluster .
finally , recent measurements of the sunyaev - zeldovich effect by @xcite indicate that the inferred optical depth of the comptonizing gas might actually be higher than one would infer from the simple spherically symmetric model obeying hydrostatic equilibrium , and a clumpy or elongated structure ( aligned along the line of sight ) would alleviate the discrepancy .
clearly , this cluster deserves further detailed studies ; it is a very interesting potential target for the astro - e2 mission , where the x - ray calorimeter will provide an unprecedented spectral resolution , capable of clearly confirming the complex redshift structure of the intra - cluster gas , hinted by the xmm data .
if the cluster is indeed undergoing a merger , this might provide an environment where at least some fraction of the gas is accelerated ( for instance , via shocks ) to form a non - thermal distribution of particles , which might produce radiation detectable in radio ( via synchrotron emission ) , in hard x - rays ( by comptonizing the cosmic microwave background ) , and in gamma rays , potentially detectable by the future mission glast .
the superior effective area of the xmm - newton telescope has allowed us to perform a detailed analysis of abell 1689 . comparing with the data from rosat , asca and chandra we verified that the data are consistent with earlier observations .
importantly , there is no indication of any additional absorbing component in the xmm data , attributing the low - energy excess absorption found by @xcite in the chandra data to uncorrected instrumental effects .
the now available ` acisabs ` code successfully corrects for these effects .
we confirm earlier findings that there is a large discrepancy , of a factor 2 or more , between mass estimates from gravitational lensing and x - ray derived mass using an unparameterized deprojection technique .
our analysis indicates that this discrepancy is true for the central part of the cluster , but also might be the case for the entire observable cluster ; our finding is in contrast to @xcite , who conclude that at large radii , there is no disagreement between the x - ray and lensing masses .
although the x - ray determined mass appears to be discrepant from the values determined from most lensing techniques , it seems to be in good agreement with that derived from weak gravitational shear .
we compare the gas mass of the cluster with the total mass , and find that for abell 1689 , @xmath191 is @xmath250 , significantly less than @xmath251 , the value derived by @xcite for 10 dynamically relaxed clusters .
our calculation of the asymmetric temperature distribution of the cluster provides further evidence that this cluster is not in a relaxed state : the lower than expected gas mass fraction is yet another piece of evidence .
we also present the first measurement of asymmetries in redshift for different regions of the icm in the cluster determined from the x - ray data alone .
a similar analysis was done by @xcite who claim a 3-@xmath252 detection of bulk motion from asca observations of the centaurus cluster .
we argue that the redshift variation detected here might be either due to line - of - sight clustering or possibly due to large bulk motions of the gas . even though this cluster is clearly not as relaxed as might be expected from its apparent spherical form , we find the slope of the total mass in the central region to be in good agreement with what is expected from numerical simulations of structure formation .
the density slope for @xmath253 is @xmath254 . while our current understanding of the structure and dynamics of galaxy clusters is insufficient to put limits on the self - interaction of dark matter , better data and more accurate simulations appear promising for the future .
we would like to thank yasushi ikebe and steve snowden for kindly helping in answering questions regarding xmm - newton data reduction and calibration issues , and james chiang for the suggestions regarding spectral deprojection .
we are grateful to the xmm sas software team for supplying the data reduction software .
k. andersson would also like to thank per carlson and tsuneyoshi kamae for making this research possible .
x - ray data used in this work was extracted from nasa s heasarc archives .
we acknowledge the use of the digitized sky survey , produced at the space telescope science institute under us government grant nag w-2166 , based on the original photographic data obtained by the uk schmidt telescope .
this work was partially supported by the alice & knut wallenberg foundation via the wallenberg research link , by nasa chandra funds administered via sao as go1 - 2113x , and by department of energy contract de - ac03 - 76sf00515 to the stanford linear accelerator center .
we are grateful to the anonymous referee for many constructive suggestions and comments .