this paper describes a geometrical framework for the study of conformally invariant fields in conformally flat spaces in four dimensions .
applications to scalar and electromagnetic fields are made .
in particular , we show that the equation @xmath0 in @xmath1 , @xmath2 , where @xmath3 is the laplace - beltrami operator , leads under very moderate assumptions to both the maxwell equations and the conformal eastwood - singer gauge conditions @xcite on @xmath4d conformally flat spaces .
in addition , our proof of this result allows us to propose a new fiber bundle in which , broadly speaking , the maxwell equations on a conformally flat space are converted into constrained scalar equations on minkowski space .
this drastically simplifies the practical calculations .
the basic geometrical idea is to build four - dimensional spaces as the intersection in @xmath1 , of a surface and the five - dimensional null cone ( invariant under the linear conformal group so@xmath5 ) .
the metric on such a space is induced from that of @xmath1 . in particular , for any given robertson - walker ( rw ) metric , one can always find a surface such that the induced metric is the rw metric .
now , to each point of a space obtained in this way corresponds a half line of the cone , hence all spaces can be realized as subsets of the set of the half lines : the cone modulo the dilations .
in addition , in the intersection of two such subsets the spaces are related through a weyl rescaling . as a special case
, the minkowski space can be obtained in that scheme .
this allows us to build a particular atlas of the cone modulo the dilations . in effect
, copies of that minkowski space can be obtained by displacing the surface intersecting the cone thanks to the action of o@xmath6 .
this generates a covering of the set of the half lines .
one can go a step further by introducing coordinates of @xmath7 with the specific property that their restriction to a particular minkowski space of the covering are the usual minkowskian cartesian coordinates .
each space of the covering together with its corresponding set of minkowskian cartesian coordinates can be turned into a local chart of the cone modulo the dilations , which is thus endowed with a minkowskian atlas .
this atlas is one of the key ingredients to perform simple calculations . in particular
, it allows us to use the local weyl rescaling to handle global problems .
the study of conformal fields , as the electromagnetic field , viewed as restrictions to four dimensional spaces of fields on @xmath7 can be traced back to the seminal paper of dirac @xcite .
there , he introduces the `` six cone formalism '' which sets the conditions , mainly the homogeneity of the fields on @xmath7 , to obtain conformal minkowskian fields .
this approach takes into account the so@xmath5 symmetry of the minkowskian equations under consideration from the beginning .
this formalism has been extended in group theoretical context by mack and salam @xcite almost @xmath8 years later and in study of conformal generalizations of qed on minkowski and anti - de  sitter spaces in the mid 80 s by many authors ( see @xcite and references herein ) .
the fact that the conformal group so@xmath5 includes as subgroups , besides the poincar and the anti - de  sitter groups , the de  sitter group was our starting point in the study of the relation between conformal scalar field in the de  sitter and minkowski space @xcite,@xcite .
the generalization to the electromagnetic field in a conformal gauge on the de  sitter space was tackled in @xcite . there following @xcite we used a method of auxiliary fields to obtain a two - point function in a conformal gauge which reduces to the eastwood - singer gauge @xcite .
these works was concerned by free quantum fields , but part of the methods used also applied in the case of classical fields with sources . in that context , the tools previously build , especially the use of minkowskian charts on the set of half lines , allowed us to reproduce , with simplified calculations , the result obtained by higuchi and cheong @xcite for the problem of two charges in de  sitter space .
the geometrical framework we build up in the present work encompass and generalizes the method we developed in these previous studies .
it extends the formalism to a class of conformally flat spaces which contains in particular the spaces endowed with a rw metric .
it provides an atlas which permits a global use of the weyl relation to minkowski space .
last but not least , it also provides a deeper view on the geometrical nature of the objects involved .
this article is organized as follows . the sec .
[ sec - geometry ] is devoted to the geometrical framework .
we first gives a global view and motivate the main definitions .
then we turn to the geometrical formulation and consequences on tensors fields of the assumptions of transversality and homogeneity .
this allows us to see how the weyl rescaling emerges in this context .
the action of the so@xmath5 group on tensor fields is then discussed . the construction of the minkowskian atlas , and the expression of some properties in it follows .
this section ends by the proof that all the rw metrics can be obtained in the present formalism . in sec .
[ sec - scalarfield ] , we discuss the conformal scalar fields , we show how both equations and fields on @xmath7 and on a conformally flat space are related . the sec .
[ sec - maxwell ] contains the proof of the proposition that the maxwell equations and the eastwood - singer gauge condition are obtained from homogeneous one form fields satisfying the laplace - beltrami in @xmath7 and some transversality requirement .
although the result is purely geometrical , the proof makes use of calculations performed in the minkowskian atlas .
this make apparent the practical calculation method inherited from this geometrical framework .
some generalizations are discussed in sec .
[ sec - somecom ] . some properties of homogeneous tensors
are reminded in appendix [ sec - homogtensor ] .
appendix [ sec - appcoordmink ] resembles a few additional comments on some particular coordinate systems for the sake of completeness .
definitions and conventions are collected in appendix [ sec - appendicedef ] . for convenience , the conventions for indices are repeated here : @xmath9 the space @xmath7 is provided with the metric @xmath10 .
a point of @xmath7 is denoted by @xmath11 or @xmath12 , the letter @xmath12 refers most often to the linear structure of @xmath7 .
this section is intended to give the reader a global view of the formalism .
we motivate the introduction of key structures as the cone modulo the dilations or the bundle @xmath13 .
we give their definitions and their main properties , leaving the proofs to the forthcoming sections . for convenience a reminder of definitions and conventions
is given in appendix [ sec - appendicedef ] .
the main geometrical construction is pictured in fig . [ fig - globalpicture ] .    a four dimensional manifold @xmath14
is obtained as the intersection of the five dimensional null cone @xmath15 of @xmath7 and a surface @xmath16 , in which @xmath17 is homogeneous of degree one . on a manifold @xmath14 each point is intercepted by a single half line of @xmath15 .
conversely , thanks to the homogeneity of @xmath17 , each half line intercepts at most once the manifold @xmath14 . as a consequence each manifold @xmath14 can be realized as a subset @xmath18 of a single set : the set @xmath19 of the half lines of @xmath15 , namely the cone modulo the dilations .
this construction can be translated in relations between spaces which are collected in the following diagram @xmath20 }   \xymatrix { \setr^6 & \\ x_f \ar[u]^{l_f } & x'_f\subset \mathcal{c ' } \ar[ul]_{l'_f } \ar[l]^{\lambda_f } } \ ] ] in which @xmath21 & \mapsto x_f : = [ x ] \cap x_f , \nonumber\end{aligned}\ ] ] where @xmath22 $ ] is the half line of @xmath15 ( that is a point of @xmath19 ) which contains the point @xmath11 , @xmath23 is the canonical injection , and @xmath24&\mapsto x_f : = [ x ] \cap x_f,\nonumber\end{aligned}\ ] ] which is a diffeomorphism .
note that , thanks to the homogeneity of @xmath17 one has @xmath25 , and thus @xmath26    now , a set of half lines may cross several manifolds @xmath14 , in other words the realizations @xmath18 on @xmath19 of these manifolds have a non - empty intersection , in that case one can define some common coordinate systems ( at least locally ) .
this property allows us to drastically simplify many practical calculations , this can be explained as follows .
firstly , as we will prove in sec .
[ sec - minkatlas ] , the cone modulo the dilations as a manifold can be covered by a collection of realizations @xmath18 of minkowski spaces , each endowed with a cartesian system of coordinates .
we will call minkowskian atlas this covering together with these coordinates , and minkowskian charts the elements of this atlas . secondly
, we will prove in sec .
[ sec - weylrelation ] that in their intersection two realizations @xmath18 are related through a weyl rescaling . finally , due to the existence of the above covering all the manifolds @xmath14
are weyl related to a minkowski space , that is are conformally flat .
in addition , and ( this is the point for the simplifications ) a conformal equation reads under its minkowskian form in a minkowskian chart .    amongst the spaces obtained within this scheme
the space @xmath27 , which is obtained as the intersection of the cone with the @xmath28-sphere @xmath29 , provides a realization of the abstract manifold @xmath19 itself .
indeed , @xmath29 is intercepted only once by each half line of @xmath15 ( fig .
[ fig - globalpicture ] ) .
we will use this realization of @xmath19 , which is included in the einstein space , for drawing conformal diagrams in sec .
[ sec - minkatlas ] and some practical calculations ( sec .
[ sec - scalarfield ] ) .
other important structures , more specifically related to the maxwell field , have to be introduced . in the scheme we use , which is inspired by the dirac s `` six cone formalism '' @xcite , the maxwell field @xmath30 , a one form field on @xmath14 ,
is obtained from a one form field @xmath31 on @xmath7 homogeneous of degree zero .
we remind ( see appendix [ sec - homogtensor ] ) that the cartesian components of @xmath31 are then homogeneous of degree -1 . in sec .
[ sec - maxwell ] we will prove the following property  : for @xmath31 homogeneous of degree zero , @xmath32 homogeneous of degree @xmath33 satisfying @xmath34 the field @xmath35 and the current @xmath36 defined on @xmath14 satisfy the maxwell equations and the eastwood - singer gauge condition @xcite  : @xmath37 where @xmath38 is the laplace - beltrami operator on @xmath14 . note that for maximally symmetric spaces the eastwood - singer gauge condition reduces to the more familiar form @xmath39 with a constant ricci scalar @xmath40 .
note also that the two conditions appearing in ( [ eq - box6aetdieze - overviewformalism ] ) mean that the vector fields @xmath41 and @xmath42 are tangent respectively to @xmath15 and to @xmath14 .    beside the property by itself
its proof makes use of a new geometrical structure : the fiber bundle @xmath43 it is obtained in restricting the base space of the cotangent bundle @xmath44 to @xmath14 .
this object allows us to solve some problems in a more simple way than by tackling them directly on @xmath14
. it also describes more accurately the transition between the fields on @xmath7 and those on @xmath14 .
more precisely , although under the assumptions made in ( [ eq - box6aetdieze - overviewformalism ] ) the maxwell field @xmath30 is directly obtained as the pullback of @xmath31 .
the proof of this result involves an intermediate step in which the equations @xmath45 are , in a sense which will be made precise in sec .
[ sec - maxwell ] , restricted to the six conformal scalar equations @xmath46 where @xmath47 is the scalar laplace - beltrami operator on @xmath14 and in which both the @xmath48 s and the @xmath49 s belongs to the set @xmath50 of the sections of @xmath51 .
these equations reduce ( when the two constraints of ( [ eq - box6aetdieze - overviewformalism ] ) are applied ) to the maxwell equations and the eastwood - singer condition ( [ eq - maes - overviewformalism ] ) .
however they are more easier to handle .
even more , when quantizing the electromagnetic field , these equations account directly for the gauge fixing : part of the @xmath52 s components leads to the maxwell field , the others carry constraints @xcite .    the field @xmath48 is obtained from the field @xmath31 by restricting the base of @xmath44 to @xmath14 and the field @xmath30 from the field @xmath48 by further restricting the fiber @xmath53 to @xmath54 .
the relations between @xmath55 and @xmath30 ( also true for @xmath56 and @xmath57 ) are collected in the following diagram @xmath58_{\widetilde{l_f } } \ar[dr]^{l^*_f } & \\
a^f \in \mathcal{b}_f^ * \ar[r]_{{r_f } } & \omega_1(x_f ) \ni
a^f   } \ ] ] in which @xmath59 such that : @xmath60 , @xmath61 = a(x ) [ v]$ ] ,    @xmath62    such that : @xmath63 , @xmath64 = a(x)[u ] $ ] .
this is a condition on tensor fields which will appear very often . roughly speaking , for a tensor @xmath65 fulfilling this relation , the map @xmath66 has a good behavior with respect to group invariance , weyl relation , equations
... this condition reads as follows : a tensor @xmath67 fulfills the transversality condition if and only if , for any @xmath68 and for any @xmath69 , @xmath70as soon as one of the arguments @xmath71 is equal to the dilation field @xmath72 ,  @xmath73 .
note that , when @xmath74 , one often finds the little bit stronger condition @xmath75 which , for @xmath76 , reduces to @xmath77 ( from which the name `` transversality condition '' originates after dirac s paper @xcite , altough in the context of mechanics this kind of condition is usualy termed `` horizontality '' ) , equivalently in an index - free notation : @xmath78 .
this condition will eventually be required for any field considered here . for quantum fields , as usual
when a gauge condition is present , the implementation of the constraint at the quantum level ( in a gupta - bleuler quantization scheme ) is done after the quantization .      in this section ,
we consider two manifolds , @xmath14 and @xmath79 defined , as in sec .
[ sec - overviewformalism ] , as the intersections of @xmath15 and the manifolds @xmath80 .
let @xmath67 be an homogeneous tensor field of degree @xmath81 , we want to compare the fields @xmath82 and @xmath83 .
for this purpose we realize both fields on @xmath19 and obtain a weyl relation between @xmath84 and @xmath85 .
let us show that if @xmath65 fulfills the transversality condition then , on @xmath86 , the following property holds : @xmath87\right ) = \left(k^f_h([x])\right)^{d(t ) }   \left(t'^h\right)\left([x]\right).\ ] ] where we defined @xmath88 ) : = \frac{h(x)}{f(x)}.\ ] ] thus the function @xmath89 appears as the conformal factor . for further references
we also define the related one form on @xmath19 @xmath90 note that , no reference to a metric structure is made in these two definitions .    in order to prove ( [ eq - weyltenseurshomogenes ] ) , we begin with showing that , for @xmath91 and @xmath92}(\mathcal{c'})$ ] one has @xmath93)\left [ v'\right ] = h(x)l'_{h*}([x])\left[v'\right ] + n,\ ] ] where @xmath94 is a null vector which belongs to a half - line @xmath22 $ ] .
we first remark that @xmath95)= \frac{f(x)}{h(x)}l'_f([x])=k_f^h([x])l'_f([x]).\ ] ] then , differentiating the rightmost term , one obtains , for any @xmath68 and @xmath96}(\mathcal{c'})$ ] , @xmath97)[v']= l'_{k_f^h*}([x])[v']l'_f([x])+k_f^h([x])l'_{f*}([x])[v'].\ ] ] the first term of the r.h.s . of this expression is proportional to @xmath11 , and as a consequence , belongs to the half - line @xmath22 $ ] .
the expression ( [ eq - lf*etlh * ] ) follows at once .
now , let us consider @xmath98}(\mathcal{c'})$ ] .
then using in succession the definitions ( [ eq - deflpf])-([eq - deflamdaf ] ) , the homogeneity of @xmath65 , the property ( [ eq - lf*etlh * ] ) and finally the transversality condition for @xmath65 ( sec . [ sec - transvcond ] ) , one has @xmath99)\left[v'_1,\ldots , v'_p\right]\\ & = \left(l'^*_f t\right)\left([x]\right)\left[v'_1,\ldots , v'_p\right]\\   & = t\left(\frac{x}{f(x)}\right)\left[l'_{f*}v'_1,\ldots , l'_{f*}v'_p\right]\\ & = \left(\frac{1}{f(x)}\right)^{d(t ) - p}\!t(x)\left[l'_{f*}v'_1,\ldots , l'_{f*}v'_p\right]\\ & = \left(\frac{1}{f(x)}\right)^{d(t)}\!\!\!t(x)\left[f(x)l'_{f*}v'_1,\ldots , f(x)l'_{f*}v'_p\right]\\ & = \left(\frac{1}{f(x)}\right)^{d(t)}\!\!\!t(x)\left[h(x)l'_{h*}v'_1 + n_1,\ldots , h(x)l'_{h*}v'_p + n_p\right]\\ & = \left(\frac{1}{f(x)}\right)^{d(t)}\!\!\!t(x)\left[h(x)l'_{h*}v'_1,\ldots , h(x)l'_{h*}v'_p \right]\\ & = \left(\frac{h(x)}{f(x)}\right)^{d(t)}\!\!\!t'^h([x])\left[v'_1,\ldots , v'_p\right],\\\end{aligned}\ ] ] which is the announced result .
finally , thanks to the isomorphism @xmath100 between the @xmath14 s and the @xmath18 and using the notation @xmath101 for the tensors on the @xmath18 s instead of @xmath102 s , one can recast the above relation ( [ eq - weyltenseurshomogenes ] ) under the more familiar form @xmath103 in which @xmath104 stands for @xmath105)=t^f(x_f)$ ] . in particular , between metrics
this relation specializes to @xmath106 which makes apparent that the eq .
( [ eq - weyltenseurshomogenes ] ) is in fact a weyl relation in the usual sense .      in this section
we specify the action of so@xmath5 on the various objects defined above .
in particular , we define the action of so@xmath5 on the sections @xmath48 of the fiber bundle @xmath107 and show that the transversality condition ensures the so@xmath5 invariance of the construction .
the natural action @xmath108 of so@xmath5 on @xmath7 yields an action @xmath109 on @xmath14 defined through : @xmath110 setting @xmath111 for any @xmath112 , one obtains the action of the conformal group on the space time : @xmath113      let us recall that , for a scalar field @xmath114 homogeneous of degree @xmath115 the operator @xmath116 is defined through : @xmath117 also , the group acts on @xmath114 through the natural representation @xmath118 .
the field @xmath119 is defined on @xmath14 and we now define the representation @xmath120 of the conformal group on it .
we just impose that the operator @xmath116 intertwines @xmath120 and the natural representation @xmath121 on @xmath114:@xmath122 that is to say , using that @xmath123 for any @xmath124 , @xmath125 which is the well - known action of the conformal group on conformal pseudo scalar fields of weight @xmath115 .
the group acts on a 1-form @xmath31 of @xmath7 through @xmath126 where @xmath127 stands for the spinorial action on the 6 components of @xmath31 .
we are going to define a representation @xmath128 of so@xmath5 on the 1-forms of @xmath14 as well as a representation @xmath129 on the sections of @xmath51 .
moreover we will prove that the map @xmath130 intertwines these representations as soon as the form @xmath31 fulfills the transversality condition : @xmath131 which reads @xmath132= ( l_g^f)^ca^f(x)[u],\forall x\in x_f,\ \forall u\in t_x(x_f),\ ] ] as soon as @xmath133 .
we first build the representation @xmath129 on the field @xmath48 in a very similar way as for the scalar field .
we just impose that @xmath134 intertwines the representations @xmath129 and @xmath135 and obtain straightforwardly : @xmath136=(\omega_g^f(g^{-1}x))^{r-1}g\cdot a^f(l_{g^{-1}}^fx)[v],\ ] ] for any @xmath137 , which is _ a fortiori _
true for any @xmath138 .
we now define the representation @xmath128 of the group on the forms @xmath30 . as a result
, we can not define a representation on the @xmath30 in the same way as above . as we will see
, one must impose a condition on @xmath31 . in place of this construction
, we define directly the representation @xmath128 through @xmath139= ( \omega_g^f(g^{-1}x))^ra^f(l_{g^-1}^fx)[(l_{g^{-1}}^f)_*u],\ ] ] for @xmath124 and @xmath138 .
we begin with calculating @xmath140 on any @xmath138 using the leibniz rule : @xmath141=\frac{g\cdot u}{f(gx)}-g x\left(\frac{\partial f(g x)[u]}{f(gx)^2}\right).\ ] ] the crucial remark is that the second term of the r.h.s . belongs to a line of the cone .
this term will be denoted as @xmath142 in the following . at this time
we can set down the representation @xmath128 , putting @xmath143 for the readability : @xmath144&=\omega^ra^f(l^f_{g^{-1}}x)\left[\frac{1}{\omega } g^{-1}\cdot u+f\right]\\&=\omega^ra(l^f_{g^{-1}}x)\left[\frac{1}{\omega } g^{-1}\cdot u+f\right]\\ & = \omega^ra(l^f_{g^{-1}}x)\left[\frac{1}{\omega } g^{-1}\cdot u\right]\\ & = ( \omega_g^f(g^{-1}x))^{r-1}g\cdot a^f(l^f_{g^{-1}}x)[u ] , \end{aligned}\ ] ] the before last equality being due to the transversality condition .
the result ( [ eq - intertwine ] ) follows immediately from ( [ eq - rho_tilde ] ) and ( [ eq - rho ] ) . as a consequence ,
the map @xmath145 intertwines the representations @xmath128 and @xmath146 as soon as @xmath31 fulfills the transversality condition .
note that the identification @xmath147 allows to realize the group action on @xmath18 as well .      as explained in sec .
[ sec - overviewformalism ] all the spaces @xmath14 are realized as subsets @xmath18 of @xmath19 through the one to one map @xmath100 .
then , one can use the weyl relation ( [ eq - weylmetricxfxh - short ] ) in order to simplify practical calculations in problems involving conformally invariant equations . as well known
, these calculations can be much simpler than those performed on @xmath14 especially when the second space @xmath79 is a minkowskian one .
unfortunately , one can not in general recover the space @xmath14 with only one minkowskian space .
however , we will prove that the covering of the whole @xmath19 is possible using four minkowskian spaces .
these can be endowed with coordinates systems in order to form an atlas .
in addition , these systems can be chosen to be cartesian minkowskian coordinates to makes calculations simpler .
now , the formalism we develop use @xmath7 as framework , in practical calculations the coordinates basis on which fields can be expanded are chosen on @xmath7 . as a consequence ,
the cartesian minkowskian coordinates of the atlas have to be related to coordinates in @xmath7 . in the sequel , we will call minkowskian systems both the systems in @xmath7 and those deduced from them in the atlas of @xmath19 .
this atlas will also be named minkowskian as well as charts that compose it .
we will prove that one minkowskian plane @xmath148 together with one minkowskian coordinate system , the @xmath94-system , are sufficient to generate the whole minkowskian atlas .
note that , some general comments in relation with such kind of coordinates are made in appendix [ sec - appcoordmink ] .
first let us introduce a useful graphical representation of @xmath19 as conformal diagram .
remind that @xmath19 can be realized as the intersection @xmath27 of the @xmath28-sphere , obtained through @xmath149 , @xmath150 , with the cone @xmath15 ( see fig .
[ fig - globalpicture ] ) .
points of @xmath27 and @xmath19 are then identified through the map @xmath151 .
this space @xmath27 is naturally endowed with a very convenient global coordinate system @xmath152 obtained by setting @xmath153 in the hyper - spherical coordinate system @xmath154 defined through @xmath155 in which @xmath156 , \varphi\in[0,2\pi]$ ] and @xmath157 correspond to the usual spherical coordinates on @xmath158 .
the condition @xmath159 gives a system on the cone , @xmath160 gives a system on @xmath27 .
now , since as a result of the considerations of the sec .
[ sec - weylrelation ] a space @xmath14 is related to @xmath27 through a weyl rescaling , the subsets @xmath18 can be conformally pictured in the plane @xmath161 . in the sequel
we will call `` @xmath161-diagram '' this conformal mapping ( see fig . ( 2 ) for an example . ) .
one may note that the metric element on @xmath27 in the @xmath152 system is @xmath162 this form shows that @xmath27 is indeed included ( @xmath163 belongs to a subset of @xmath164 ) in the static einstein space of positive curvature .
now , let us consider the space @xmath165 , @xmath166 . as shown in @xcite
this space is a minkowski space .
the coordinate system @xmath167 defined by @xmath168 ( whose inverse is given for convenience in appendix [ sec - appcoordmink ] ) provides on @xmath169 a chart , called @xmath94-chart in the sequel .
the restriction to the cone is obtained for @xmath170 and that for @xmath148 for @xmath171 .
the metric induced from @xmath7 on @xmath172 is , in the @xmath173-coordinate basis , @xmath174 , the usual form of the minkowski metric .
the points where the above system becomes singular , that is on the subset @xmath175 , correspond to points at infinity in the minkowski space @xmath165 .
they are mapped to the boundaries of the @xmath94-chart on the conformal @xmath176-diagram in fig .
( [ fig - northandsouthchart ] ) .
similar considerations apply to the space @xmath177 defined by @xmath178 and the coordinate system @xmath179 is obtained by replacing @xmath94 by @xmath29 in ( [ eq - coordnord ] ) .    a set of minkowskian charts covering @xmath19 is finally obtained by moving the two surfaces @xmath148 , @xmath180 with some elements of so@xmath5 .
more precisely , we obtain new minkowski spaces through the action of the one parameter subgroup of so@xmath5 generated by @xmath181 .
the manifolds @xmath182 ( resp .
@xmath183 ) defined through @xmath184 are minkowski spaces . corresponding to each case , one can define a coordinate system @xmath185 analogous to @xmath186 , replacing @xmath187 by @xmath188 in ( [ eq - coordnord ] ) and similarly for @xmath29 in place of @xmath94 . as a consequence
, we obtain an atlas of @xmath19 , in the usual sense , in which the systems of coordinates make the weyl relations very transparent .
now every @xmath18 space is an open subset of @xmath19 and can be endowed with the above atlas .    as a final remark
, the space @xmath180 can be obtained from the space @xmath148 through the transformation @xmath189 which is also an element of o@xmath6 . as a consequence
, all the elements of our minkowskian atlas can be obtained from @xmath148 through a transformation of o@xmath6 that is to say through an ( inoffensive ) cartesian isometric change of variables in @xmath7 .
thus , the whole atlas is generated from @xmath148 and the @xmath94-coordinates as announced in the beginning of this section .
we note that , for @xmath124 , @xmath190    also , in @xmath94-coordinates the dilation field reads @xmath191 .
now , a one - form @xmath74 homogeneous of degree zero , is expanded in the coordinate basis @xmath192 as @xmath193 the last term is for convenience divided by @xmath194 in order that @xmath195 be of the same degree of homogeneity ( zero ) as the other components .
one has @xmath196 in particular , the transversality condition reads @xmath197 .
the same result applies in others minkowskian coordinates .
here we prove that for any robertson - walker ( rw ) metric @xmath198 one can find a function @xmath199 such that the metric induced from @xmath7 on the conformally flat space @xmath200 be the robertson - walker metric @xmath198 .
in other words , each robertson walker space can be realized at least locally as a space @xmath201 .
following ibison @xcite the set called @xmath202 of conformally flat metrics defined through  : @xmath203 @xmath204 being a conformal factor and @xmath205 the cartesian minkowskian coordinates , includes the rw metrics .
the notation @xmath206 for the arguments of the function @xmath207 indicates that it depends separately of the time coordinates and of the space coordinates and that the space coordinates appears only through the radius @xmath208 .
let us assume now that @xmath207 for some rw metric is given in a minkowskian coordinate system , say the system @xmath209 ( [ eq - coordnord ] ) .
the function @xmath210 is invariant under the subgroup so@xmath211 of so@xmath5 generated by @xmath212 .
the equation @xmath213 defines a so@xmath211 invariant space endowed with a metric which satisfies ( [ eq - rtmetric ] ) , thus a rw space .
in this section , we consider scalar fields on @xmath7 , homogeneous of degree @xmath115 .
these fields yield scalar fields on @xmath14 .
we show that these fields fulfill the scalar conformal equation as soon as the original field fulfill the laplace - beltrami equation on @xmath7 .
let @xmath114 and @xmath32 two scalar fields on @xmath7 homogeneous of degree -1 and -3 respectively , and suppose that @xmath214 we claim that @xmath215 where @xmath38 is the laplace - beltrami operator and @xmath40 the scalar curvature for the space @xmath14 .
we consider @xmath27 the realization of the cone up to the dilation defined sec .
[ sec - minkatlas ] , on which all @xmath14 spaces realize as subsets . setting @xmath216 and @xmath217 in the hyperspherical coordinates ( [ eq - coordalpbetsurc ] ) and using the homogeneity of @xmath114 ,
the equation ( [ eq - equationscalairer6 ] ) reads @xmath218 the restriction to @xmath27 ( @xmath219 ) then gives @xmath220 which is the desired result .
we obtain the similar result for the other hyper surfaces @xmath14 by using the weyl correspondence between the spaces @xmath14 .
note that although straightforward the change of variables leading to the equation ( [ eq - equationscalairecoordr1r2 ] ) is rather cumbersome .
in this section we deal with one - forms @xmath31 and @xmath32 of @xmath7 satisfying @xmath221 and such that @xmath31 is homogeneous of degree 0 and @xmath32 homogeneous of degree @xmath33 ( remind that this implies that @xmath222 and @xmath223 are homogeneous of degree @xmath224 and @xmath225 respectively ) .
we first show how in the present formalism the equation @xmath221 leads to a set of six copies of the equation of the conformal scalar field .
then , we give the proof that @xmath226 satisfies the maxwell equations as soon as @xmath31 fulfills the transversality condition , and that an additional condition on @xmath32 allows us to recover a conformal gauge condition .
suppose that @xmath51 is endowed with a coordinate system @xmath228 where @xmath229 is any coordinate system of @xmath14 and @xmath230 a cartesian system of coordinates on @xmath231 .
in such a coordinate system the equation on @xmath48 reads @xmath232 where @xmath47 is the scalar laplace - beltrami on @xmath14 and @xmath40 the ricci scalar .
this is a straightforward consequence of the fact that , on @xmath7 in cartesian coordinates one has @xmath233 : since in a cartesian coordinate system @xmath234 of @xmath7 the equation @xmath235 reduces to six copies ( one per components of @xmath31 and @xmath32 ) of the scalar equation ( [ eq - equationscalairer6 ] ) , each components @xmath222 and @xmath223 being homogeneous of degree @xmath236 and @xmath237 , one can thus apply the result ( [ eq - equationscalairexpf ] ) of the sec .
[ sec - scalarfield ] .
two comments are in order concerning equations ( [ eq - equationsurfibre ] ) . in first
, they do not present any gauge ambiguity and we proved in @xcite that they allow a quantization of the maxwell field on de sitter space using a gupta - bleuler scheme where the condition of transversality is translated to a condition on states after quantization . in second
, they can be used to solve classical propagation problems for the maxwell field : for instance , we considered in @xcite the two - charges problem for the de sitter space .      in this section , we implicitly identify the spaces @xmath14 and @xmath18 , writing all the objects on the common manifold @xmath238 .
we now prove the statement of sec .
[ sec - overviewformalism ] : for @xmath31 homogeneous of degree zero , @xmath32 homogeneous of degree @xmath33 satisfying @xmath239 the field @xmath35 and the current @xmath240 defined on @xmath14 satisfy the maxwell equations and the eastwood - singer gauge condition @xcite  : @xmath241 we first prove the result on the minkowskian chart @xmath165 , after what , using the o@xmath6 invariance of the hypothesis and of the conclusion , and the properties of the minkowskian charts , we obtain this result on all the minkowskian charts . in the following , for readability , we note @xmath242 , for @xmath243 and @xmath244 .
we first consider the equation @xmath45 , and apply ( [ eq - equationsurfibre ] ) on the minkowski space @xmath165 defined by @xmath187 ( used to build the n - chart ) , this yields the system @xmath245 where @xmath246 is the scalar laplace - beltrami operator on @xmath247 , @xmath248 and @xmath249 . as a second step we express in the @xmath94-chart of @xmath19 ( [ eq - coordnord ] ) the above system of equations ( [ eq - equationscalaire - an ] ) .
this leads to the system already obtained in @xcite ( with slightly different notations ) @xmath250    then , we express the two constraints appearing in ( [ eq - box6aetdieze - maxwell ] ) in the @xmath94-chart . the first one is the transversality condition applied to the field @xmath31 , using the formula of ( [ eq - transvcondcn ] ) one has @xmath251 the second condition @xmath252 implies the transversality of @xmath32 since @xmath253 , thus @xmath254 in addition , @xmath252 rewrites @xmath255 = 0 $ ] , which in the @xmath94-coordinates on the cone ( @xmath256 ) reads @xmath257 thanks to the homogeneity of @xmath17 this equation rewrites @xmath258 then , using the definition ( [ eq - defw ] ) , it becomes @xmath259 where the dot refer to the minkowskian metric @xmath260 .    using the constraint ( [ eq - aplus=0]-[eq - jcw ] ) ,
the system ( [ eq - systauxmink ] ) becomes @xmath261 owning to the map @xmath262 the above system leads to @xmath263 on the minkowski space @xmath165 .
note that the map @xmath262 becomes obvious in the n - coordinates ( [ eq - coordnord ] ) : for all @xmath264 one has @xmath265 .    for the sake of argument , let us introduce a self explanatory symbolical notation for the two operators appearing in the l.h.s . of the maxwell equations and the eastwood - singer gauge condition .
with them the above system reads : @xmath266\right)_\mu = j^{{\scriptscriptstyle n}}_\mu\\ & es_{{\scriptscriptstyle n}}[a^{{\scriptscriptstyle n } } ] = j^{{\scriptscriptstyle
n } } \cdot w^{(f , f_{{\scriptscriptstyle n } } ) } , \end{aligned } \right .\ ] ] this result is available for the @xmath94-chart .
nevertheless , the o@xmath6 invariance of hypothesis and conclusion and the fact that any minkowskian chart can be deduced from the @xmath94-chart by mean of an o@xmath6 transformation , proves that this result is true for any minkowskian chart of our atlas .
now , using the same symbolical notation for the two operators appearing in the l.h.s of ( [ eq - maes - maxwell ] ) , and following @xcite to apply a usual ( local ) weyl transformation between @xmath165 and @xmath14 to each eqs .
of the system ( [ eq - systmink - w ] ) , one obtains @xmath267\right)_\mu = j^f_\mu\\ & es_f[a^f ] + w^{(f , f_{{\scriptscriptstyle n } } ) } \cdot \left(m_f [ a^f]\right)= j^{{\scriptscriptstyle f } } \cdot w^{(f , f_{{\scriptscriptstyle n } } ) } , \end{aligned } \right .\ ] ] where the dots now refer to the metric @xmath268 and where we have taken into account the homogeneity of both the electromagnetic and the current one - form fields .
they have respectively a degree of @xmath224 and of @xmath225 which correspond to conformal weights of zero and @xmath33 for the components of the fields  : @xmath269 , @xmath270 .
finally , the second equation simplifies and we obtain the announced result on the @xmath94-chart .
then , using the above remark on o@xmath6 invariance , we obtain the result on the whole space .
note that the maxwell equations are obtained independently of the condition @xmath252 , which is used only to obtain the eastwood - singer gauge condition .
note also that , in absence of source ( @xmath271 ) , the eastwood - singer gauge condition is automatically fulfilled as soon as the transversality condition on @xmath31 is fulfilled .
it is important to point out that from a physical perspective the initial conditions and currents are given in the space @xmath14 . in particular ,
the minkowskian currents appearing in ( [ eq - systmink - w ] ) are defined through the equation @xmath272 in order to satisfy the weyl relation .
they have in general no physical meaning in the minkowski space .
besides the results presented in the previous sections , the formalism depicted in the present paper explicits the geometrical nature of the various objects used in our previous works .
this allows us to consider some straightforward generalizations , let us comment briefly about the classical and quantum situations .
the method used for the classical propagation problem considered in @xcite extends naturally to the conformally flat spaces @xmath14 . to summarize , given a set of initial conditions an currents in @xmath14 the problem of finding the @xmath30 solution of the maxwell equations in the eastwood - singer gauge ( or in a gauge contained in it ) in @xmath14 amounts to use in the minkowskian atlas the propagation formula established in @xcite , which uses only the minkowskian scalar green s function .
the quantization scheme used in @xcite for the free electromagnetic field on de  sitter space can be transposed here for the most part .
the definition of an so@xmath5-invariant scalar product on the space of the solutions of ( [ eq - equationsurfibre ] ) and the obtention ( from the known mode solutions of ( [ eq - equationsurfibre ] ) ) of a reproducing kernel for the @xmath273 s with respect to that product can be reproduced almost verbatim .
the result is : @xmath274 , where @xmath275 is the scalar two - point function on @xmath14 . the construction of the fock space , including the determination of the physical subspace is formally identical .
the general form of a covariant two - point function , a bi - tensor on a specific @xmath14 , requires more developments ( note that , such a function has been proposed recently in minkowskian coordinates @xcite ) .
this could be the object of future investigations .
the authors thank c. cachot for valuables discussions and suggestions about the secs.([sec - minkatlas ] ) and ( [ sec - rw ] ) .
for reference , we recall here some properties of an homogeneous @xmath276-tensors @xmath277 of @xmath7 .    the set @xmath7 , together with the scalar product @xmath278 @xmath279 , is naturally endowed with a structure of linear space . a point in @xmath7 is located by a vector @xmath12 whose components in cartesian coordinates are @xmath234 .
the linear structure allows us to identify the tangent spaces at any point through @xmath280 .
this makes sense for expressions as @xmath281 since the tangent space at @xmath12 is identified with the tangent space at @xmath282 .
let us consider the dilatation map in @xmath7 @xmath283 @xmath284 being a positive real number . since @xmath285
this is also the result of a push - forward on the vector @xmath12 . in other words
for a @xmath276-tensor field @xmath277 of the pullback reads @xmath287 .
now , an homogeneous @xmath276-tensor field @xmath277 of degree @xmath288 is defined through the relation @xmath289 thus , homogeneous @xmath276-tensor field @xmath277 of degree @xmath115 satisfy @xmath290    the homogeneity of a tensor field can be related to that of its components in some coordinate basis whose coordinates @xmath291 @xmath292 are each one homogeneous functions of degree @xmath293 .
precisely , let us show that : for @xmath294 a @xmath276-tensor field of @xmath7 , homogeneous of degree @xmath295 one has @xmath296 in particular for @xmath297 @xmath298 this result is obvious in the particular case @xmath297 since the eq .
( [ eq - tenshomoliendefpullback ] ) reads in components @xmath299 now , moving to the more general coordinate basis @xmath291 on has @xmath300 where @xmath301 , is homogeneous of degree @xmath302 .
taking homogeneity into account in the r.h.s .
of the above expression leads to @xmath303 from which the result follows .
the system appearing in sec . [ sec - minkatlas ] has already been used in the literature ( see for instance @xcite ) , it is reminiscent of systems called polyspherical systems ( see for instance @xcite ) .
one may note that they are framed in such a way that the `` extremal''coordinates @xmath304 are functions of the constraints ( defining the cone @xmath15 and the space @xmath148 respectively ) whereas the `` central '' coordinates @xmath205 are those on the background space @xmath14 .
cartesian minkowskian coordinates can be obtained on @xmath14 setting @xmath305 instead of @xmath306 , with this choice conformally invariant equations will appears on their minkowskian form on @xmath14 .    for the de  sitter space
, one may verify that the @xmath94 and @xmath29 minkowski charts correspond to the stereographic projections from the north and south poles from which their names originates .
finally , we note for convenience that the inverse of system ( [ eq - coordnord ] ) reads @xmath307
we summarize here the conventions , main structures and maps used in this paper .
they are discussed in sec .
[ sec - overviewformalism ] . here are the conventions for indices : @xmath9 the coefficients of the metric @xmath308 of @xmath309 are denoted @xmath310 .
the definitions of spaces and maps reads  :        @xmath19 denotes the set of the half lines of the cone @xmath15 , namely the cone modulo the dilations .
one has @xmath314
\in   \mathcal{c'}$ ] iff @xmath314=\{u \in \mathcal{c } : \exists \lambda > 0 , u = \lambda y\}$ ]                      the laplace - beltrami operators on @xmath1 and @xmath14 are denoted respectively by @xmath3 and @xmath38 .
the operator @xmath327 is the scalar laplace - beltrami operator on @xmath14 : it acts on its argument as if it were a scalar .
aaa m. eastwood and m. singer , phys .
107a * , 73 ( 1985 ) .
dirac , ann . math . * 37 * , 429 ( 1936 ) .
g.  mack and a.  salam , ann . phys .
( n.y . ) * 53 * , 174 ( 1969 ) . f.  bayen , m.  flato , c. fronsdal , and a.  haidari , phys . rev .
d * 32 * , 2673 ( 1985 ) .
e. huguet , j. queva , and j. renaud , phys .
d * 73 * , 084025 ( 2006 ) .
e. huguet , j. queva , and j. renaud , phys .
d * 77 * , 044025 ( 2008 ) .
s. faci , e. huguet and j. renaud , phys .
d * 84 * , 124050 ( 2011 ) .
a. higuchi and l. y. cheong , phys .
d * 78 * , 084031 ( 2008 ) .
s. faci , e. huguet , j. queva , and j. renaud , phys .
d * 80 * , 124005 ( 2009 ) .
b.  binegard , c.  fronsdal , and w.  heidenreich , j. math
* 24 * , 2828 ( 1983 ) .
m.  ibison , j. math .
phys . * 48 * , 122501 ( 2007 ) .
s.  faci , arxiv:1205.3804 . h.a .
kastrup , annalen phys .
* 17 * , 631 ( 2008 ) .


the high mass x - ray binary ( hmxb ) system was discovered in 1969 april during a balloon experiment @xcite .
the system consists of an accreting neutron star ( ns ) fed by the surrounding stellar wind of the b type emission line companion wray 977 @xcite .
a recent luminosity estimate derived from atmospheric models puts its distance at @xmath0kpc @xcite , the value utilized in this paper .
the orbital period was established to be @xmath1days @xcite using _ ariel _  5 observations and was refined with the _ burst and transient source experiment ( batse ) _ to @xmath2days with an eccentricity of @xmath3 @xcite .
@xcite discussed a possible orbital evolution and determined an orbital period of @xmath4d , assuming no change in orbital period .
@xcite determined that the mass of the companion was in the range @xmath5 and the radius of wray977 was @xmath662@xmath7 , obtained by fitting atmosphere models .
the x - ray flux is highly variable throughout an individual binary orbit but follows a distinct pattern when averaged over multiple orbits ( see figure  [ fig : orbitlc ] ) . shortly before the periastron passage ,
the x - ray luminosity increases drastically in the energy band above @xmath8kev , as seen in _ rossi x - ray timing explorer ( rxte)_/ all sky monitor ( asm ) data @xcite .
the ns passes closest to the companion at a distance of @xmath9 @xcite . shortly after the periastron passage , @xmath10
, the x - ray luminosity dips for a short period of time .
@xcite demonstrated that neither a simple spherical wind model nor a circumstellar disk model around wray 977 are sufficient to describe the observed variations in the folded _ rxte_/asm data .
an additional stream component is able to account for the sudden increase in x - ray luminosity , as the ns passes trough the stream shortly before periastron and accretes more material .
this model also explains a slightly higher x - ray luminosity around @xmath11 , when the ns passes through the accretion stream a second time .
pulsations with a period of @xmath12s were discovered in the _ ariel-5 _
observations @xcite , making one of the slowest known pulsars .
the pulse period has varied drastically throughout the last @xmath13 years @xcite . prior to 1984 ,
the pulse period stayed relatively constant at 695s@xmath14700s and then spun up between 1985 and 1990 to @xmath15s . from 1993 until the beginning of 2008 , the change in the spin reversed again , showing a decline .
_ fermi_/galactic burst monitor ( gbm ) data have revealed that experienced another spin reversal and briefly spun - up with the pulse period decreasing from @xmath16s to @xmath17s between may 2008 and october 2010 . since october 2010
, the pulse period has shown only small variations around @xmath17s .
most recently , @xcite discovered a peculiar 1ks dip in the luminosity of , where the pulsations disappeared for one spin cycle during the dip .
several such dips have been previously observed in vela  x-1 @xcite , where it is assumed that the accretion on the ns was interrupted for a short period of time .
the pulse phase average spectrum of is described using a power law with a high energy cutoff .
the continuum does not show a strong variation in the intrinsic parameters ( @xmath18 , and @xmath19 ) throughout the orbit @xcite , as seen in two data sets from _ rxte _
, taken in 1996 and 2000 , sampling most phases of the binary orbit .
one of the major characteristics of the x - ray spectrum of is the high and strongly variable column density of its line - of - sight absorber throughout the orbit , indicative of a clumpy stellar wind ( @xmath20 @xmath21 ) .
in addition to the high column density , a very bright fe k@xmath22 emission line can be observed .
this line has shown a strong correlation with the observed luminosity , indicating that the line is produced by local clumpy matter surrounding the neutron star @xcite .
@xcite used the _ rxte _ data set from 2000 to perform phase resolved spectroscopy and showed that an absorbed and partially covered pulsar continuum ( power law with fermi - dirac cutoff ) as well as a reflected and absorbed pulsar continuum were consistent with the data .
a cyclotron resonance scattering feature ( crsf ) at @xmath23kev was first discovered with _
ginga _ @xcite . @xcite found systematic deviations from a power law continuum at @xmath13 and @xmath24kev in _ bepposax _ , where the former could not be confirmed as a crsf due to the proximity of the continuum cutoff .
@xcite excluded the existence of a crsf at @xmath13kev and showed that the crsf centroid energy varies between 3038kev over the pulse rotation of the ns .
furthermore , they showed that the crsf centroid energy and width are correlated .
we report on two observations of performed with the _ suzaku _ satellite in mid 2008 and early 2009 .
this paper is structured as follows : section 2 discusses the observations and data reduction .
section 3 shows the phase averaged results .
section 4 discusses the pulse profiles and phase resolved spectra .
sections 5 and 6 discuss the results and conclusions , respectively .
_ suzaku _ observed on 2008 august 25 with an exposure time of @xmath25ks ( obsid 403044010 ; hereafter obs .  1 ) .
the observation was cut short by a set of target of opportunity observations and was continued on 2009 january 5 , acquiring an additional @xmath26ks exposure time ( obsid 403044020 ; obs .  2 ) .
both main instruments , the x - ray imaging spectrometer ( xis ; * ? ? ?
* ) and the hard x - ray detector ( hxd ; * ? ? ?
* ) were used in these observations .
the two observations correspond to orbital phases of 0.19 ( obs .  1 ) and 0.38 ( obs .  2 ) , where obs .
1 falls into the lowest flux part of the binary orbit ( see figure  [ fig : orbitlc ] ) .
the @xmath27kev absorbed flux was @xmath28erg@xmath21s@xmath29 for obs .  1 and @xmath30erg@xmath21s@xmath29 for obs .  2 ( see table 1 for details ) .
both observations were performed using the hxd nominal pointing to enhance the sensitivity of the hxd detectors .
the xis detectors consist of two front illuminated ( fi ) ccds xis 0 and 3 , and one back illuminated ( bi ) ccd , xis 1 .
all three instruments are sensitive between @xmath31kev and @xmath32kev , although the two fi cameras have a higher sensitivity above @xmath33kev , and the bi chip is more sensitive below @xmath33kev . to minimize possible pile - up ,
the xis instruments were operated with the 1/4 window option with a readout time of 2s .
data were taken in both @xmath34 and @xmath35 editing modes , which were extracted individually with the _ suzaku _ ftools version 16 as part of heasoft 6.9 .
the unfiltered xis data were reprocessed with the most recent calibration files available and then screened with the standard selection criteria as described by the _ suzaku _
abc guide .
the response matrices ( rmfs ) and effective areas ( arfs ) were weighted according to the exposure times of the different editing modes .
the xis data were then grouped with the number of channels per energy bin corresponding to the half width half maximum of the spectral resolution , i.e. grouped by 8 , 12 , 14 , 16 , 18 , 20 , and 22 channels starting at 0.5 , 1 , 2 , 3 , 4 , 5 , 6 , and 7kev , respectively ( m.  a.  nowak , 2010 , private communication ) .
the xis spectral data were used in the energy range of @xmath27kev for reasons described below .
the hxd consists of two non - imaging instruments : the pin silicon diodes ( sensitive between 1260kev ) and the gso / bgo phoswich counters ( sensitive above @xmath36kev ) . to determine the pin background
, the _ suzaku _ hxd team provides the tuned pin non x - ray background ( nxb ) for each individual observation .
in addition , the cosmic x - ray background ( cxb ) was simulated following the example in the abc guide and both backgrounds were added together .
the pin data were grouped by a factor of 5 below and by a factor of 10 above 50kev for obs .  1 . in obs .
2 the channels were grouped by 3 throughout the whole energy range .
gso data were extracted and binned following the _ suzaku _ abc guide .
the pin data energy range of 1560kev and the gso data energy range of @xmath37kev were used for the phase averaged analysis . due to the short exposure times , no gso data were extracted in obs .  1 and in the phase resolved analysis of obs .  2 .
is a source with very large and highly variable photoelectric absorption in the line of sight as well as a very strong fe k@xmath22 emission line . for spectral modeling
, the calculated responses of the xis detectors showed a ` leakage ' emission level @xcite , which stemmed from the instrument characteristics .
that is , a small fraction of an event charge cloud is registered at energies below the peak energy , creating a low energy tail in the spectrum .
this ` tail ' is about three orders of magnitude weaker than the intensity of the peak count rate of the measured photon .
the response to the fe line is therefore described as a gaussian plus a constant at lower energies , approximately three orders of magnitude below the peak value and extending to below 1kev .
due to the strong fe  k@xmath22 line and the very strong absorption in this source , the constant level response of the fe line extends above the measured flux below 2kev .
this results in significant residuals below @xmath38 kev , which seem to be stronger for the bi spectrum then for that of the fi ccd .
figure  [ fig : badspec ] shows the phase averaged spectrum from the second observation of and the best fit models using the heasoft generated xis responses as detailed in section 3 .    together with the _ suzaku _ guest observer facility
( gof ) , which provided for us an experimental response parameter file where the ` low energy tail ' is ignored , we studied this behavior in more detail .
this new response for all three instruments was applied to the individual xis spectra . due to the fact that the new response completely excluded the missing ` tail ' , we now observed an excess of the data with respect to the model below @xmath39kev .
the introduction of additional power law and/or black body components could improve the fits at lower energies .
due to the fact that the true level of the tail is not known , these additional components can not be interpreted physically , and we decided to use the original response matrix and confine the selected xis energy range to above 2kev .
phase averaged broad band spectra in the 260kev ( obs .  1 ) and 290kev ( obs .
2 ) energy ranges were obtained . from the technical description of the xis instruments
, it is known that small discrepancies , e.g. , in fitted power law slope , have been observed between fi and bi xis instruments .
a difference of @xmath40 in the power law index has been observed in calibration data and was also observed here . when modeling the fi and bi instruments with a common power law index , the residuals of the bi instrument show a systematic deviation .
however , only the fi or bi instruments could be modeled together with the hxd instruments . due to this fact and the higher sensitivity of the fi xis instruments above 2kev
, we concentrated our discussion on the results obtained with the two fi instruments .
each data set , fi and bi , was modeled individually with the hxd data to compare the differences in the continuum .
we found that when using the fi and bi instruments individually with the hxd instruments , in both observations the best fit spectral values are consistent with each other within error bars , indicating that the usage of the hxd reduced the observed discrepancies in the power law parameters .    the continuum model for pulsars
can so far only be modeled with empirical models , consisting of a ` powerlaw ` with a cutoff at higher energies , which is typical for this kind of source @xcite .
three empirical models are widely used : the simple high energy cutoff ( ` highecut ` ) and the fermi - dirac cutoff ( ` fdcut ` , * ? ? ?
* ) both are used in combination with a simple power law component . in addition , the negative - positive exponential powerlaw model ` npex ` @xcite is a slightly more complicated model including the power law component . for the data analyzed in this work , the best results have been obtained with the ` fdcut ` model : @xmath41 where @xmath42 is the normalization at 1kev , @xmath43 is the cutoff energy , and @xmath19 is the folding energy of the fermi - dirac cutoff . for the obs .
1 and obs .  2 fi data sets the best fit @xmath44 s were 260 and 542 with 245 and 267 degrees of freedom ( dof ) , respectively ( see figure  [ fig : spec010 ] and [ fig : spec020 ] ) . replacing the smoother fermi - dirac cutoff with the ` highecut ` model in obs .
2,resulted in slightly worse residuals ( @xmath45 ) compared to the best fit results of the fd - cutoff model ( figure  [ fig : spec010 ] d ) and an emission line - like residual at @xmath46kev , which is most likely due to the sharp break of the power law at the cutoff energy in this model . in comparison ,
the npex model shows even worse residuals ( @xmath47 ) and only results in reasonable fits when the exponential curvature at higher energies is independent in the partially and fully covered component .
the best fit cutoff energy is very close to that of the observed crsf feature ( see below ) and these two parameters are rather strongly correlated ( see figure  [ fig : cont ] ) . to avoid a degeneracy of the crsf values due to a changing cutoff energy , @xmath43 was frozen in all the fits to the best fit value from the obs
.  2 fi spectrum ( 29.2kev ) .
previous observations of indicated the existence of clumps in the stellar wind @xcite , which were modeled using partial covering absorption in addition to fully covered photoelectric absorption of the smooth stellar wind . in the present analysis ,
the low energy portion of the spectrum also required a partial covering component ( figure  [ fig : spec010]b and [ fig : spec020]b ) , which was modeled using the ` tbnew ` model ( wilms et al .
, 2011 , in prep . ) , an updated version of the existing ` tbabs ` model @xcite .
in addition , a non - relativistic , optically - thin compton scattering component ` cabs ` was included , as is necessary for column densities @xmath48 @xmath21 , where the plasma becomes compton thick and part of the emission is scattered out of the line of sight .
the @xmath49 values from ` tbabs ` and from the ` cabs ` component were set equal and treated as one model component for the smooth stellar wind ( @xmath50 = ` tbabs1*cabs1 ` ) and a second component for the clumpy partial coverer ( @xmath51 = ` tbabs2*cabs2 ` ) .
the @xmath50 and @xmath51 column densities were left independent of each other . for
spectral fitting the abundances of @xcite and the cross - sections @xcite were used in all data sets .
both observations show residuals to the continuum in the @xmath52kev energy range , which we interpret as the previously observed crsf @xcite . modeling these residuals with an absorption line with a gaussian optical depth
( ` gabs ` ) improved the residuals significantly in both observations ( see figures  [ fig : spec010 ] and [ fig : spec020 ] ) .
the best fit values of @xmath53kev and @xmath54kev for obs .  1 and obs .  2 ,
respectively , are consistent within 90% confidence intervals .
the widths of the crsfs , @xmath55kev ( obs .  1 ) and @xmath56kev ( obs .
2 ) , are also consistent within errors .
an increase of @xmath57 with higher centroid energy , as previously observed in _
data @xcite , could not be detected , although the best fit values hint at such a behavior .
table  [ tab1 ] shows the best fit values for the continuum with a ` gabs ` crsf line where the given errors are 90% confidence values . to calculate the significance of the crsf in the first and weaker observation ,
the null hypothesis approach was applied , where 10,000 spectra were created with monte carlo simulations using the best fit parameters without the crsf .
gaussian uncertainties were used for the individual model parameter .
each model was fitted with and without the crsf line to compare how much an inclusion of the line actually improves the fit . in @xmath58 of all fits ,
no crsf feature were observed with a larger optical depth than the observed lower limits of @xmath59 , making the existence of the feature significant with @xmath60 .
including the crsf in the real data , the best fit improved by a @xmath61 .
a similar improvement of @xmath44 was only observed in @xmath62 ( @xmath63 ) of the simulated spectra , concluding that the observed line in the first observation is indeed real .
for the second observation , the line was even more pronounced , where the exclusion of the crsf increased the @xmath44 by @xmath64 .
lcccc parameter & & + & fi & bi & fi & bi + @xmath65 @xmath21 ] & @xmath67 & 16.6 ( frozen ) & @xmath68 & 20.9 ( frozen ) + & 1.55 ( frozen ) & 1.55 ( frozen ) & @xmath69 & @xmath70 + & 1.17 ( frozen ) & 1.17 ( frozen ) & @xmath71 & @xmath72 + @xmath73 @xmath21 ] & @xmath74 & 76.9 ( frozen ) & @xmath75 & 28.4 ( frozen ) + @xmath76 & @xmath77 & @xmath78 & @xmath79 & @xmath80 + @xmath81 @xmath82 ^ 1 $ ] fully covered & @xmath83 & @xmath84 & @xmath85 & @xmath86 + @xmath87 @xmath82 ^ 1 $ ] part . covered & @xmath88 & @xmath89 & @xmath90 & @xmath91 + @xmath43 [ kev ] & 29.2 ( frozen ) & 29.2 ( frozen ) & @xmath92 & 29.2 ( frozen ) + @xmath19 [ kev ] & @xmath93 & @xmath94 & @xmath95 & @xmath96 + @xmath97 [ kev ] & @xmath98 & @xmath99 & @xmath54 & @xmath100 + @xmath57 [ kev ] & @xmath101 & @xmath102 & @xmath56 & @xmath103 + @xmath104 & @xmath105 & @xmath106 & @xmath107 & @xmath108 + @xmath109 [ kev ] & @xmath110 & @xmath111 & @xmath112 & @xmath113 + @xmath114 ^ 2 $ ] & @xmath115 & @xmath116 & @xmath117 & @xmath118 + eqw@xmath119 [ ev ] & @xmath120 & @xmath121&@xmath122 & @xmath123 + @xmath124 [ kev ] & @xmath125 & @xmath126 & @xmath127 & @xmath128 + @xmath129 ^ 2 $ ] & @xmath130 & @xmath131 & @xmath132 & @xmath133 + eqw@xmath134 [ ev ] & @xmath135 & @xmath136 & @xmath137 & @xmath138 + @xmath139 [ kev ] & @xmath140 & @xmath141 & @xmath142 & @xmath143 + @xmath144 ^ 2 $ ] & @xmath145 & @xmath146 & @xmath147 & @xmath148 + eqw@xmath149 [ ev ] & @xmath150 & @xmath151 & @xmath152 & @xmath153 + @xmath154 [ kev ] & @xmath155 & @xmath156 & @xmath157&@xmath158 + @xmath159 [ ev ] & @xmath160 & @xmath161 & @xmath162 & @xmath163 + @xmath164 ^ 2 $ ] & @xmath165 & @xmath166 & @xmath167 & @xmath168 + eqw@xmath169 [ ev ] & @xmath170 & @xmath171 & @xmath172 & @xmath173 + @xmath174 [ kev ] & @xmath175 & @xmath176 & @xmath177 & @xmath178 + @xmath179 ^ 2 $ ] & @xmath180 & @xmath181 & @xmath182 & @xmath183 + eqw@xmath184 [ ev ] & @xmath185 & @xmath186 & @xmath187 & @xmath188 + @xmath189 [ kev ] & @xmath14 & @xmath14 & @xmath190 & @xmath191 + @xmath192 ^ 2 $ ] & @xmath14 & @xmath14 & @xmath193 & @xmath194 + eqw@xmath195 [ ev ] & @xmath14 & @xmath14 & @xmath196 & @xmath197 + flux@xmath198 absorb .
& @xmath199&@xmath200 & @xmath201 & @xmath202 + flux@xmath198 unabsorb . &
@xmath203 & @xmath204 & @xmath205 & @xmath206 + @xmath207 & @xmath208 /  & / @xmath209 /  & @xmath210 &  / @xmath211 + @xmath44/dof & 260 / 245 & 145 / 130 & 542 / 267 & 328 / 150 +   +    the crsf can alternatively be described with the lorenzian shaped ` cyclabs ` xspec model @xcite . `
cyclabs ` is described by the centroid energy @xmath97 , the width @xmath57 , and the resonance depth @xmath104 , similar to the ` gabs ` parameters .
the best fit continuum parameters are consistent with the best fit values determined with the ` gabs ` component .
the ` cutoff ` energy with a best fit value of @xmath212kev is consistent with the value obtained with the ` gabs ` model .
the observed centroid energies of @xmath213kev and @xmath214kev are of the order of @xmath215 lower than the energies obtained with the ` gabs ` model .
this discrepancy stems from a different calculation of the line centroid energy and is described in detail in @xcite .
the width @xmath216kev and @xmath217kev for obs .  1 and obs .  2 ,
respectively , is bigger than with the ` gabs ` model .
the @xmath44 / dof values of 260 / 245 and 548 / 267 for the two observations using the ` cyclabs ` model and thus were not a better fit when compared to the ` gabs ` model .
for the final discussion we use the values determined by the ` gabs ` model .
several emission features were observed in the residuals of the fits in both observations and were subsequentially modeled with gaussian emission lines ( see figures  [ fig : spec010 ] and [ fig : spec020 ] , inlay ) .
the width of each line was set equal to that of the fe k@xmath22 emission line , while the intensities were left to vary independently .
energies were loosely constrained around the expected values of neutral material to avoid runaway of the line energies .
a constant ( ` const ` ) was applied , taking small instrumental differences in the overall flux normalization into account .
the constant was fixed at 1 for xis 0 and was left free for the other instruments .
@xmath218 , @xmath219 and @xmath220 in table 1 are the cross calibration constants with respect to xis0 for fi fits and with respect to xis1 for the bi fits .
the final model had the form : ` const`*`n_\text{h,1 } ` ( ` pl1`+ ` n_\text{h,2}`*`pl2 ` ) * ` fdcut`*`gabs`+ @xmath221`gaussians ` , where the power law indices of ` pl1 ` and ` pl2 ` are set to be equal to each other .
the individual power law normalizations were independent and were used to calculate the fraction of the partial covering .
table [ tab1 ] summarizes the best fit values for both observations and for the individual fi / bi data sets . for the interpretation we will concentrate on the fi data for both observations , for reasons mentioned above .
the best fit spectral parameters were generally consistent with previous observations , with the exception of the cutoff energy . with a value of @xmath92kev ,
@xmath43 was significantly higher than the @xmath13kev value measured in previous _
rxte _ @xcite and _ bepposax _ @xcite observations , although these observations used slightly different spectral models for the cutoff energy .
@xcite used also the here applied fermi
dirac cutoff , resulting in cutoff energies of @xmath222kev .
the folding energy for obs .
1 , @xmath223kev , is consistent with values obtained for a similar orbital phase with _
@xcite . in obs .  2 ,
the value of @xmath224kev is consistent with the _ rxte _ data for the pre - periastron flare and the periastron passage @xcite .
the observed column densities for the absorption due to the smooth stellar wind ( @xmath50 ) were consistent between both observations .
a larger difference could be observed in the column density associated with the absorption due to the clumped wind ( @xmath51 ) , where the best fit value of obs .
1 is almost a factor three larger than that of obs .  2 .
in addition to the h column density in the ` tbnew ` model , the relative ca and fe abundances were also left independent in obs .  2 .
the best fit values are measured from the ca and fe k edges at @xmath225kev and @xmath226kev , respectively , and were slightly higher than solar abundances : @xmath69 for ca and @xmath71 for fe . for the first observation ,
the ca and fe abundances could not be well constrained and were fixed to 1.55 for ca and 1.17 for fe .
the partial covering fractions can be calculated from the measured normalization values of the two power law components , @xmath227 and @xmath228 : @xmath229 the factor is almost unity in both observations due to the dominance of the @xmath228 value .
similar large covering fractions have been observed in _ xmm - newton _
data @xcite and in the _ rxte _ data @xcite , indicating that the covering fraction does not change significantly for different parts of the orbit .
the spectrum softened slightly from @xmath230 for obs .  1 to @xmath79 for obs .  2 .
the fe k@xmath22 line was detected at @xmath157kev ( obs .
2 ) , which is , given the instrumental gain systematics , consistent with neutral fe .
in addition , the following lines have also been observed in both observations : s k@xmath22 line , ar k@xmath22 line , ca k@xmath22 line , fe k@xmath231 line , and the ni k@xmath22 line , all with energies consistent with neutral material ( e.g. * ? ? ? * ) .
note that the ni k@xmath22 line was not detected in the fainter first observation .
the observed line intensities are summarized in table  [ tab1 ] with 90% confidence errors .
the values are consistent between the fi and bi instruments .
the widths of the lines were set to be equal to the fe k@xmath22 width within each observation , which was found to be @xmath160ev for obs .  1 and had an upper limit of 5ev in obs .  2 ( fi values ) .
the compton shoulder to the fe k@xmath22 line , as observed by the _
observatory @xcite , was not significantly detected in either observation .
the observation of @xcite was performed in the pre - periastron phase , where the luminosity was higher than in the two _ suzaku _ observations .
a compton shoulder has also been observed with _
xmm - newton _ in data taken during another pre - periastron flare @xcite .
barycentric and binary corrected light curves for different energy bands were extracted for both observations using the orbital parameters from @xcite with an updated orbital period and periastron time @xmath232 determined by @xcite . due to
the strong variability of the pulse period on short time scales , individually determined pulse periods were used for each observation .
an accurate pulse period could be established for obs .  2 , but not for obs .  1 due to its short duration of 10ks .
for the second observation the calculated pulse period had a value of @xmath233s .
fermi_/gbm instrument measures the pulse period of on a regular basis , resulting in values of @xmath234s for obs .  1 and @xmath235s for obs .  2 ( finger , 2011 , priv .
comm ) . to create the pulse profiles , the _ fermi_/gbm - provided pulse periods and epoch times for each observation were used .
we studied the energy dependent pulse profiles for obs .
1 and obs .  2 .
both pulse profiles are very similar , showing a double peak shape where the main peak ( p1 ) is broader below 10kev .
the second peak ( p2 ) stays rather constant in width , but increases its relative intensity toward higher energies . as an example , the pulse profile for obs .  2 in different energy bands
is shown in figure  [ fig : prof020 ] . comparing these pulse profiles with previous _
rxte _ and _ bepposax _ data showed that the general shape is consistent through all parts of the orbit .    the data of the longer obs .
2 were divided into 10 equally spaced phase bins and individual spectra were extracted for the xis and pin instruments .
figure  [ fig : prof020 ] shows the individual phase bins in the lowest panel .
this division resulted in an exposure time in each phase bin of @xmath236ks for each individual instrument .
no gso data were used in this analysis due to the reduced exposure time per phase bin . for spectral analysis ,
the phase averaged model was applied to all phase bins , resulting in best fit values summarized in table  [ tab2 ] , as well as figure  [ fig : phasespec1 ] .
again the cutoff energy was frozen to the phase averaged value of 29.2kev to avoid the previously discussed degeneracy with @xmath97 .
lccc cccc ccc parameter & pb1 & pb2&pb3 & pb4&pb5&pb6&pb7 & pb8&pb9&pb10 + @xmath65/@xmath237 & @xmath238 & @xmath239 & @xmath240 & @xmath241 & @xmath242 & @xmath243 & @xmath244 & @xmath245 & @xmath246 & @xmath247 + @xmath73/@xmath237 & @xmath248 & @xmath249 & @xmath250 & @xmath251 & @xmath252 & @xmath253 & @xmath254 & @xmath255 & @xmath256 & @xmath257 + @xmath76 & @xmath258 & @xmath259 & @xmath260 & @xmath261 & @xmath262 & @xmath263 & @xmath264 & @xmath265 & @xmath266 & @xmath267 + @xmath268 ^ 1 $ ] ful .
@xmath269 & @xmath270 & @xmath271&@xmath272&@xmath273&@xmath274 & @xmath275 & @xmath276 & @xmath277 & @xmath278 + @xmath279 ^ 1 $ ] par
@xmath280 & @xmath281 & @xmath282 & @xmath283 & @xmath284 & @xmath285 & @xmath286 & @xmath287 & @xmath288 & @xmath289 + @xmath19 [ kev ] & @xmath290 & @xmath291 & @xmath292 & @xmath293 & @xmath294 & @xmath295 & @xmath296 & @xmath297 & @xmath298 & @xmath299 + @xmath97 [ kev ] & @xmath300 & @xmath301 & @xmath302 & @xmath303 & @xmath304 & @xmath305 & @xmath306 & @xmath307 & @xmath308 & @xmath309 + @xmath57 [ kev ] & @xmath310 & @xmath311 & @xmath312 & @xmath313 & @xmath314 & @xmath315 & @xmath316 & @xmath317 & @xmath318 & @xmath319 + @xmath104 & @xmath320 & @xmath321 & @xmath322 & @xmath323 & @xmath324 & @xmath325 & @xmath326 & @xmath327 & @xmath328 & @xmath329 + @xmath330 ^ 2 $ ] & @xmath331 & @xmath332 & @xmath333 & @xmath334 & @xmath335 & @xmath336 & @xmath337 & @xmath338 & @xmath339 & @xmath340 + flux@xmath341 & @xmath342&@xmath343 & @xmath344 & @xmath345 & @xmath346 & @xmath347&@xmath348&@xmath349 & @xmath350 & @xmath351 + c@xmath352/c@xmath353 & 0.94 / 1.15 & 0.94 / 1.15 & 0.95 / 1.19 & 0.94 / 1.21 & 0.94 / 1.21 & 0.94 / 1.19 & 0.95 / 1.16 & 0.95 / 1.14 & 0.94 / 1.10 & 0.94 / 1.11 + @xmath44/dof & 233 / 243 & 288 / 243 & 353 / 243 & 312 / 243 & 253 /243 & 303 / 243 & 283 / 243 & 292 / 243 & 264 / 243 & 247 / 243 + @xmath44/dof no crsf & 241 / 246 & 307 / 246 & 391 / 246 & 341 / 246 & 273 / 246 & 305 / 246 & 291 / 246 & 344 / 246 & 297 / 246 & 270 / 246 +    the two absorbing components , @xmath50 from the smooth stellar wind and @xmath51 from the partial coverer vary between 20 and 40@xmath354 @xmath21 , but show an anti - correlated trend throughout the first peak .
although the error bars are rather big , one can see an indication that @xmath50 follows the flux in p1 , whereas the @xmath51 value dips at the same time .
in contrast to @xcite , the power law index @xmath76 is varying strongly throughout the pulse .
the values are 0.91.3 throughout the first peak and drop suddenly to 0.60.8 for the second peak ( figure  [ fig : phasespec1 ] ) . measuring the second peak to be significantly harder than the first peak supports the behavior seen in the pulse profile , where the intensity of p2 increases at higher energies , becoming similar to the intensity in p1 .
power law normalizations are relatively small and badly constrained for the first power law component .
the partially covered power law normalization seems to follow the flux , showing a higher value throughout the first peak .
the folding energy does not vary significantly throughout the orbit and shows values between 5 and 6kev ( not shown in figure 7 )    the crsf energy varies between 3040kev , similar to the values observed in the _ rxte _ data @xcite .
table 2 shows that the addition of the crsf in the individual spectra improved the @xmath44 for most of the phase bins , except in phase bin 6 , which falls in the gap between both p1 and p2 .
phase bins 1 and 7 , also minima in the pulse profile , only show a small improvement in the best fit when the crsf was included .
the crsf energy is not very well constrained in phase bins 1 and 6 due to the lack of sufficient statistics .
for the geometrical discussion below , these two phase bins are ignored . in
all other phase bins the addition of the crsf improves the fit .
the energy changes very smoothly throughout the pulse and does not follow explicitly the observed flux in each phase bin .
the crsf energy increases during the main peak and reaches the maximum at the falling flank of p1 . during p2 the crsf energy decreases until it reaches a minimum at the dip between p2 and p1 .
to estimate the significance of the crsf in the phase resolved analysis , again the null hypothesis method was applied in two out of the 8 remaining phase bins , which are good representatives of all data points .
phase bin 2 is a good example for a shallow crsf , whereas phase bin 5 is an example for a phase bin with a higher flux .
similarly as with the phase averaged data , 10,000 spectra were simulated and modeled without and with the crsf .
the interpretation of the observed feature as being due to stochastic fluctuations with a depth similar to the observed values could be ruled out with a probability of over @xmath355 ( @xmath60 ) in both phase bins .
this leads to the conclusion that the crsf can be observed with sufficient significance to allow the interpretation below .
in both observations very strong absorption can be observed .
@xmath50 , which is the absorption due to the smooth stellar wind , is consistent within both observations .
@xmath51 , the absorption from the partial coverer , is significantly stronger in the first observation , establishing the existence of an additional component in the line of sight .
@xcite used _
rxte_/asm and pca data to show that a stellar wind and a stream model component can describe the observed count rate throughout the orbit and that the increase of the column density at orbital phase @xmath356 is mainly due to the stream component .
the findings in both _ suzaku _ observations are supportive of this theoretical picture
.    shows very strong luminosity variability on very short time scales @xcite , where the column density increases to up to 50% in periods of lower activity .
rxte _ and _ xmm - newton _ data showed a much higher column density in the pre - periastron flare than in the two data points observed with _
suzaku_. one possible explanation is that the wind is indeed so variable and clumpy that these variations are just not properly predictable , especially in times of high activity , such as the pre - periastron flare .
@xcite investigated _
rxte_/pca @xmath49 values of the absorbed component from archival observations for different parts of the orbit and found no increased column density in the pre - periastron flare .
observations throughout one full binary orbit could help to understand the variations in @xmath49 on time scales of days .
the observed luminosity dependence of the power law index @xmath76 and the folding energy is similar to other sources , such as v0332 + 53 @xcite and 4u0115 + 63 @xcite .
the power law index @xmath76 shows a hardening for the first observation , whereas the folding energy @xmath19 is slightly higher when compared to obs .  2 .
@xcite observed the @xmath19 variation in her x-1 phase resolved _
heao-1 _ data and concluded that the parameter is dependent on the viewing angle of the accretion column and can be used to describe the plasma temperature of the system . in the case of the smaller @xmath19 value in obs .  2
could indicate a lower plasma temperature which can be interpreted that the x - ray emission region is further up in the accretion column @xcite ( see also section 5.2 ) .
the softening of the power law index @xmath76 with increased luminosity is also in agreement with the basic model of the accretion column that the plasma temperature decreases with increased height . with increased luminosity
, the rate of accreted material would increase and the amount of soft photons created by the lateral walls of the relatively taller column would increase , leading to a softer spectrum , as is observed .    for the second observation
the abundances for fe and ca were left independent and a slight overabundance was observed .
taking into consideration that the abundances used in @xcite are derived from the interstellar medium ( ism ) a small overabundance from an evolved star may be expected .      in both observations the crsf
was clearly observed and the ` gabs ` absorption component improved the overall fit significantly .
strong magnetic fields ( @xmath357 g ) exist close to the ns magnetic poles , where photons at energies close to the landau levels are resonantly scattered from the line of sight and result in an absorption line - like feature in the spectrum .
this feature provides a direct method to measure the magnetic field strength close to the ns surface , where the fundamental centroid energy can be described as @xmath358 where @xmath359 is the gravitational redshift at the scattering site , with @xmath360 for typical ns . using the @xmath97 values obtained in the phase averaged spectra
, the calculated magnetic field has a strength of @xmath361 g ( obs .  1 ) and @xmath362 g ( obs .  2 ) which is consistent within errors for both observations
.    the crsf parameters are consistent with previous results using _ rxte _ data and are lower than the observed _ bepposax _ @xcite and _ integral _ @xcite values of @xmath363kev .
these observations were obtained during pre - periastron outburst , where the luminosity of the source is much higher than in the two _ suzaku _ observations presented in this paper .
a luminosity dependence of the phase - averaged crsf centroid energy is observed in multiple other sources , where an anti - correlation between crsf and luminosity was observed in v0332 + 53 and 4u0115 + 63 @xcite , and a positive correlation was observed in her  x-1 @xcite .
@xcite have confirmed , from pulse to pulse variability studies , such correlations for v0332 + 53 , 4u0115 + 63 and her  x-1 , and also found an anti - correlation for a0535 + 26 . although the crsf values are consistent with a single value in both
_ suzaku _ observations , a hint of an anti - correlation can be seen in the data .
compared with the _ bepposax _ and _ integral _ data , the crsf centroid energy might have a positive correlation .
parallel to the crsf - luminosity correlation , @xcite also observed in pulse to pulse data of four sources that the power law index @xmath76 shows an opposite luminosity dependence than the @xmath97 . in the _ suzaku _ observations discussed here ,
@xmath76 is observed to soften with increased luminosity , another indication that a possible negative crsf - luminosity correlation exists .
the two opposite scenarios seem to depend on the luminosity , i.e. , if the observed luminosity is below or above a critical luminosity ( cl ) which can be derived from the eddington luminosity of a system @xcite .
the cl is of order @xmath364 erg s@xmath29 , but depends also on the accretion geometry ( spherical or disc ) as well as the height and diameter of the accretion column .    above the cl
, the infalling proton density becomes so large that the protons begin to interact and decelerate , creating a radiation pressure dominated shock region above the magnetic pole .
this region of increased density is most likely the region where the crsf is created @xcite .
with increasing luminosity , the shock region moves higher up in the accretion column , where a smaller local magnetic field value results in a lower observed crsf centroid energy , as is observed in v0332 + 53 or 4u0115 + 63 .
the observed crsf - luminosity dependence as well as the @xmath76-luminosity dependence in is very similar to the scenario described here .    below the cl
, the accreting matter slows down via hydrodynamical shock , ` coulomb friction ' or nuclear collision @xcite . with increasing accretion rate , the luminosity increases and the deceleration region
is pushed closer to the ns surface , where the magnetic fields are higher .
this would result in a positive crsf - luminosity correlation , as observed in her  x-1 and a0535 + 26 .    with a distance of 3kpc ,
the intrinsic unabsorbed @xmath27kev luminosity of @xmath365ergs@xmath29 is significantly below the typical cl of @xmath364ergs@xmath29 . although it has been observed that sources with similar luminosities can show the opposite correlation @xcite due to individual critical luminosity values , the large luminosity difference compared to cl puts this source in the same regime as her  x-1 and a0535 + 26 and at odds with the notion of a well - defined cl . the observed variation in @xmath76 and possibly in @xmath97
would then stem from another mechanism .
the existence of multiple fluorescence emission lines , especially at lower energies , where the absorption is very dominant , indicates that the source of the emission originates from a region where the column density is not very large , i.e. the outer layers of the stellar wind in our line of sight @xcite .
if the emission lines would be embedded deeper in the stellar wind , the lines , especially at lower energies would have to be significantly stronger , to be detected at all . on the other hand , if the line emitting region is in the outer layers of the wind , the incident soft x - ray flux would be drastically reduced , making the equivalent width very large , as observed in the low energy emission lines ( table  1 ) .
a possible explanation would be that the emission region is very large for the lines , and maybe spread over the entire surface of the stellar wind .
the most dominant line emission stems from the fe  k@xmath22 transition at @xmath366kev .
the observed energies in both observations are consistent with the emission in neutral material .
the ratio of the intensities of fe  k@xmath231/fe  k@xmath22 for both observations , @xmath367 for obs .  1 and @xmath368 for obs .  2 , is consistent with the fluorescing material being neutral or only slightly ionized @xcite .
the equivalent widths ( eqw ) of the two fe k@xmath22 lines show that the intensity relative to the continuum is reduced by a factor of @xmath38 in the second observation . based on emission line widths
, @xcite estimated in _ asca _ data that the fe emission originates within @xmath369 cm of the continuum emission source .
@xcite use _ xmm - newton _
pulse by pulse observations to observe a correlation with increased flux of the source , showing that the fe emission region is not far from the x - ray source .
in the _ suzaku _ observation , however , the fe k@xmath22 line flux does not directly follow the observed luminosity .
the increase in absorbed @xmath27kev luminosity is a factor of @xmath8 , whereas the increase in the fluorescence line intensities is only @xmath370 .
this difference in intensity change indicates that the distance between continuum and fluorescence emission lines is large , especially when compared to the proposed @xmath371 lt - s from @xcite . in the 6 months time between the two observations , however , the overall emission geometry could have changed , accounting for the difference in the fluorescence intensities .
if the emission originates from a confined region in the stellar wind , observations taken during two different orbital phases would each yield a different distance to this region .
in addition to the strong fe k@xmath22 line , many other emission lines are observed in the xis spectra in both observations .
the observed intensities of these lines show the same behavior , where the flux for the second observation is @xmath372 or less compared to the first observation , although the observed uncertainties for the different lines , especially for s  k@xmath22 and ar  k@xmath22 , are very large .
furthermore , the s  k@xmath22 line falls into an energy regime where calibration uncertainties have been previously observed ( e.g. , * ? ? ?
* ) , so that a detailed analysis is not possible at this time .
historically , shows very strong pulse to pulse variations and an intensive flaring behavior on very short time scales @xcite , e.g. , throughout the pre - periastron flare .
the observations obtained here do not show such flaring behavior , and only show the regular pulsation throughout the observation .
a pulse period could only be established for obs .  2 of the _ suzaku _ data , which was consistent with _ fermi_/gbm data .
with the gbm pulse period for the first observation , pulse profiles for both observations could be produced for different energy bands .
the two peaked pulse profiles show a similar behavior for both observations , where the second peak is getting strong at higher energies .
this behavior is in contrast to many similar sources , such as 4u  0115 + 63 @xcite , v  0332 + 53 @xcite , 4u1909 + 07 @xcite , and 1a1118@xmath1461 @xcite , where the pulse profile turns into a single peak profile at higher energies .
a reduction in the intensity of the second pulse has also been observed in _ bepposax _ data for different parts of the orbit @xcite .
the phase resolved analysis of the second observation showed strong variations in several spectral parameters throughout the pulse profile .
as previously observed by @xcite , the crsf centroid energy varies by @xmath373 where the highest energy is detected at the peak and the falling flank of p1 . such behavior has been observed in multiple other sources , e.g. , vela  x-1 @xcite , 4u  0352 + 309 @xcite , and her  x-1 @xcite . a model for the crsf variations
is based upon the change in the viewing angle throughout the pulse and thus different heights of the accretion column are probed , yielding a different local magnetic field observed for each phase bin .
a simple approach to derive a possible geometry for the neutron star and the magnetic field uses the variation in the observed magnetic field throughout the pulse phase . in this case the very smooth and sinusoidal variation ( figure  8) can be modeled as a simple dipole where the total magnetic moment @xmath374 is calculated from the phase averaged crsf energy value .
the variation of the crsf energy is then fitted by changing the geometrical angles of the system until a best fit converges .
appendix a discusses the model in detail , introducing the three free parameters : @xmath375 , the viewing angle between line of sight and the ns spin axis , the inclination angle @xmath22 of the magnetic moment with respect to the spin axis , and the angle @xmath231 , indicating the lag of the @xmath376-field plane with respect to the ephemeris , i.e. , the observed shift in pulse phase ( see also figure  9 ) the best fit values for @xmath231 were either @xmath377 or @xmath378 , depending on the starting values of the fit ( see appendix a ) , where the latter value is equivalent to the former when the dipole polarity is flipped .
@xmath22 and @xmath375 angles showed a very strong interdependence , and the best fit @xmath379 pairs were : @xmath380 , @xmath381 , @xmath382 , and @xmath383 .
all best fit values have a @xmath384 with 7 dof .
these pairs do show a degeneracy of @xmath385 and @xmath386 , which can be explained with a geometrical symmetry , when calculating the magnetic field for each individual phase bin .
the first and last two pairs can be each treated as the same geometry , rotated by 180 degrees .
figure  [ fig : magnet ] shows the calculated magnetic field values from the crsf centroid energies and the theoretical values for @xmath376 from the best fit angles for a simple dipole model with @xmath387 , @xmath388 and @xmath389 .
the phase resolved energy and width of the crsf show a strong correlation , where the width varies in phase with the crsf energy ( figure .
@xcite observed a similar correlation , where the magnetic pole was observed under different viewing angles @xmath390 , where @xmath390 is the angle between the line of sight and the magnetic axis .
@xcite showed that the anisotropic velocity field of the electrons in the accretion column leads to a fractional line width of : @xmath391    the @xmath392 ratios vary in the 0.16  0.24 range , corresponding to variations of @xmath393 throughout the pulse phase .
the two outliers with large errors in phase bins 1 and 6 were ignored .
the variation of @xmath390 can be estimated as @xmath394 and @xmath395 throughout one pulse phase , where the magnetic pole rotates around the ns spin axis , tilted by @xmath375 with respect to the line of sight .    using the average ratio of 0.2
, we calculated the values of @xmath396 for the two geometries with @xmath397 ( @xmath398 ) and @xmath387 ( @xmath399 ) . for each geometry
, we used the variation of the @xmath392 ratio ( @xmath400 ) to calculate the variation in @xmath390 , resulting in values of @xmath401 and @xmath402 .
these variations are smaller then the values expected from the the geometrical discussion , but show a similar behavior , where the variation of the angle is larger for the smaller value of @xmath403 . from equation 4 , we were also able to estimate a possible plasma temperature for the two geometries from the calculated values of @xmath396 .
the plasma temperature of @xmath404kev for @xmath405 is very similar to the observed folding energy @xmath19 , which is an indication of the plasma temperature ( * ? ? ? * and references therein ) .
the value of @xmath406kev for the @xmath387 geometry is much larger and is not consistent with the observed @xmath19 .
assuming that the ns spin axis is aligned with the inclination of the binary system , a value of @xmath405 would put the inclination at @xmath407 , a value which is only marginally above the upper limit of @xmath408 , as determined by @xcite .
the strong fe line was detected in all 10 phase bins .
the measured line flux did not change significantly throughout the pulse phase , while the @xmath27kev flux varied by more than a factor of 2 . from this
we conclude that the distance to the fe fluorescence region is greater than @xmath12 lt - s ( @xmath409 cm ) , from the ns .
this is in conflict with the conclusion of @xcite that the emission region must be closer than @xmath369 cm from the ns , based upon the width of the emission lines measured by asca .
their assumption is that the emission region is close to the ns and that the material free falls onto the ns surface , encountering much faster velocities , which caused the observed broadening of the lines .
the smaller measured fe k@xmath22 line width in the phase averaged spectra is more consistent with the assumption that the line broadening is due to the terminal velocity of @xmath410 km of the line driven wind @xcite .
we presented results of two _ suzaku _ observations of taken shortly after the periastron passage .
the spectra were modeled with a partially and fully covered power law with a fermi - dirac cutoff and a crsf at @xmath46kev .
the column density of the smooth stellar component only marginally changed and the clumpy wind component significantly changed between both observations .
flux dependencies of the crsf and @xmath76 were not significant , but do hint at similar correlations such as those observed in v0332 + 53 and 4u0115 + 63 , although the observed luminosity is significantly below the calculated critical luminosity of this system .
they are in the realm where one might expect the energy of the crsf to decrease with decreasing flux .
the variations in the pulse profiles and the crsf throughout the pulse phase have a signature of a dipole magnetic field . using a simple dipole model we calculated the expected magnetic field values for different pulse phases and were able to extract a set of geometrical angles , loosely constraining the dipole geometry in the ns .
model constraints derived using the observed ratio of @xmath57 to @xmath97 , together with calculated plasma temperatures , favor a solution wherein the spin axis is tilted @xmath411 to our line of sight .
ss acknowledges support by nasa contract nas5 - 30720 and grant nnx08ax83 g .
ff is supported by dlr grant 50 or 0808 and via a daad fellowship .
ic acknowledges financial support from the french space agency cnes through cnrs . this research has made use of data obtained from the suzaku satellite , a collaborative mission between the space agencies of japan ( jaxa ) and the usa ( nasa )
the variability of the crsf energy throughout the pulse phase can be used to understand the geometry of the neutron star .
we assume a simple dipole magnetic field around a spherical neutron star , however we want to emphasize that these assumptions are very rudimentary , and effects such as the displacement of the dipole field from the neutron star center are not taken into account here .
the total magnetic moment @xmath374 of the neutron star can be calculated with : @xmath412 using the average magnetic field from the phase averaged crsf centroid energy ( @xmath413 g ) and a typical neutron star radius of @xmath414 km . the value of the last part , including the azimuthal angle @xmath390 of the magnetic pole was estimated to the average value of 1.54 .
the total average magnetic moment is divided into three components based on the cartesian coordinate system in figure  [ fig : nsmodel ] : @xmath415 where @xmath22 is the angle between the dipole field and the spin axis and @xmath231 is the rotation angle from the ephemeris , i.e. the x - axis .
next , the vector @xmath416 was calculated for each individual phase bin @xmath417 used in the phase resolved analysis .
@xmath418 the individual components of @xmath416 are the projections of each phase bin on the three axes , as seen from the line of sight of the observer .
the angle @xmath375 corresponds to the angle between the line of sight and the spin axis .
for each of the individual phase bins , the magnetic field was calculated using the position vector @xmath416 indicating the position of the magnetic pole and the vector @xmath419 , which indicates the magnetic moment : @xmath420    the calculated magnetic field values were then fitted to the measured magnetic field , varying the three angles @xmath375 , @xmath22 , and @xmath231 until the best fit converged . to avoid wrong solutions due to local minima ,
the starting parameter of all three angles were varied in steps of @xmath421 before the fit was started .
the best fit solutions resulted in the values that are discussed in detail in section 5.5 .