based on their broad - band spectral and variability properties , six of the persistently bright neutron star low - mass x - ray binaries ( lmxbs ) were classified as z sources @xcite , after the z - like tracks they trace out in x - ray color - color ( cds ) and hardness - intensity diagrams ( hids ) .
these sources are gx 17 + 2 , cyg x-2 , gx 5 - 1 , gx 340 + 0 , sco x-1 , and gx 349 + 2 .
the z - like tracks consist of three branches , which , from top to bottom , are referred to as the horizontal branch ( hb ) , the normal branch ( nb ) , and the flaring branch ( fb ) .
it is generally believed that the parameter that determines the position along the z track is the mass accretion rate , increasing from the hb to the fb .
note that the definition of the mass accretion rate in this respect is rather vague
it is often defined as the mass accretion rate through the inner disk and onto the neutron star surface .
in addition to spectral changes along the z - track , some of the z sources show long term changes in the shape and position of the z - track in the cd and hid .
these secular changes , as they are referred to , are clearly observed in cyg x-2 @xcite , gx 5 - 1 @xcite , and gx 340 + 0 @xcite ( the cyg - like sources ; the other z sources are referred to as the sco - like sources ) , and more recently also in gx 17 + 2 @xcite .
it has been suggested that they are related to the relatively high inclination at which these sources are seen @xcite , or to different magnetic field strengths of the neutron stars @xcite .
the power spectra of the z sources show several types of quasi - periodic oscillations ( qpos ) and noise components ( see * ? ? ?
* for a review ) .
it was found that their presence and properties are very well correlated with the position of the source along the z track @xcite , even when the z - tracks show secular changes ( e.g. * ? ? ?
three types of low frequency ( @xmath1100 hz ) qpos are seen in the z sources : the horizontal branch ( hbos ) , normal branch ( nbos ) and flaring branch qpos ( fbos ) .
their names derive from the branches on which they were originally found .
the hbo is found on the hb and nb with a frequency ( 1560 hz ) that gradually increases along the hb towards the nb .
when the sources move from the hb onto the nb the frequency increase flattens off . in gx 17 + 2 and
cyg x-2 it was found that when the source passes a certain point on the nb , the hbo frequency starts to decrease @xcite .
the nbo and fbo are most likely the same phenomenon .
they are found on the nb and fb ( near the nb / fb vertex ) but not on the hb . on the nb
the qpo has a frequency of @xmath057 hz , which rapidly increases to @xmath020 hz when the source moves across the nb / fb vertex @xcite . in recent years
two high frequency ( or khz ) qpo were found in the z sources ( @xcite , see @xcite , for a review ) .
they are often observed simultaneously , with a frequency difference of @xmath0300 hz , and have frequencies between 215 hz and 1130 hz , which increase from the hb to the nb . in sco x-1 the frequency difference was found to decrease with increasing qpo frequency @xcite .
the frequency difference in the other z sources is both consistent with the behavior seen in sco x-1 and with being constant @xcite .
three types of noise are seen in the z sources .
they are the very low frequency noise ( vlfn ) , the low frequency noise ( lfn ) , and the high frequency noise ( hfn ) .
the vlfn and hfn are found on all branches , whereas the lfn is only observed on the hb and nb .
the vlfn , which is found at frequencies below 1 hz , can be described by a power law .
the hfn and lfn are both band limited components , with cutoff frequencies of , respectively , 10100 hz , and 210 hz .
many competing models have been proposed for the origin of the qpos and noise components .
it is beyond the scope of this introduction to mention these models in detail - most of them will be discussed in section [ gx17 + 2_discuss_sec ] .    in this paper
we present a study , based on data acquired with the _ rossi x - ray timing explorer _ ( _ rxte _ ) , of the correlated spectral and variability properties of the z source gx 17 + 2 .
it is a continuation of the work by @xcite .
the current paper constitutes the first report on the very large observing campaign of 1999 , which more than doubled the total coverage of the source .
this campaign was undertaken with the express purpose of investigating if a non - monotonic relation exists between the frequency of the khz qpos and the hbo in gx 17 + 2 .
section [ 17 + 2_obs_sec ] deals with the observations and analysis .
the spectral results are presented in section [ spectral_sec ] and the results for each power spectral component in sections [ timing_sec ] and [ sec : res - khz ] .
a number of qualitatively new results is found in our greatly expanded data set . in particular
, we find that when the upper khz qpo frequency is below 1030 hz , it correlates with the hbo frequency , whereas above 1030 hz they anti - correlate .
we also find that khz qpo frequency difference is not constant and that the q values of the hbo and khz qpos are probably not explained by life time broadening . these and other results are compared to observations of other z sources and discussed in terms of current models in section [ gx17 + 2_discuss_sec ] .
finally , we argue that changes in the spectral and variability properties of the z sources may not be not driven by changes in the mass accretion rate .
the data used for the analysis in this paper were all obtained with the proportional counter array ( pca , @xcite ) on board _ rxte _  @xcite .
the pca consists of five xenon - filled proportional counter units ( pcus ) , each with an effective area of @xmath01250 @xmath2 ( at 10 kev ) .
although the five pcus are in principle identical they all have a slightly different energy response .
these responses change continuously due to slow processes such as gas leakage and the aging of the electrodes .
in addition , the high voltage settings of the instruments are occasionally altered ( gain changes ) , resulting in rather more drastic changes in the detector response .
these changes have , at the time of writing , been applied three times during the life time of _ rxte _ thereby defining four gain epochs .
occasionally one or more pcus are not operational .
they can be switched off by an internal safety mechanism , or by the ground control crew , for reasons of detector preservation .
therefore , the number of active detectors varies between the observations .
all our _ rxte_/pca observations of gx17 + 2 were done between 1997 february 2 and 2000 march 31 .
a log of the observations is given in table [ obs_tab ] .
we do not include the observations done in february 1996 , which were used by @xcite .
the reasons for this are the limited time resolution in the energy range of interest and difficulties with scaling the @xmath3  ( see below ) due to an incomplete z track .
data taken during satellite slews and earth occultations were removed .
the nine type i x - ray bursts that were observed are the subject of a separate article by @xcite ; two of those were also studied in this paper .
the total amount of good data that remained was @xmath0600 ks .
data were collected in several modes with different time and spectral resolutions .
two of these modes , ` standard 1 ' and ` standard 2 ' , were always operational .
in addition to these two modes , other modes were active that varied between the observations .
the properties of all modes are given in table [ gx17 + 2_modes_tab ] .
the standard 2 data were used to perform a broad - band spectral analysis .
the data were background subtracted , but no dead time corrections were applied ; these were in the order of 25% .    for each 16 s data segment ( i.e. the intrinsic resolution of the standard 2 mode ) we defined two colors , which are the ratios of count rates in two different energy bands , and an intensity , which is simply the count rate in one energy band .
the energy bands used for the colors ( soft and hard color ) and intensity , are given in table [ colors_tab ] .
the lower energy boundaries for the soft color and intensity were chosen relatively high in order to avoid the lowest and least reliable energy channels . by plotting the two colors against each other a color - color diagram ( cd ) was produced . a hardness - intensity diagram ( hid )
was produced by plotting the hard color against the intensity .
to produce the cds and hids we only used data obtained with those pcus that ( for each gain epoch ) were always on . for gain epoch 3
these were pcus 0 , 1 and 2 , and for gain epoch 4 pcus 0 and 2 .
due to the different numbers of detectors and the differences in the detector settings we decided to produce the cds and hids separately for the two gain epochs .
however , in choosing the energy channels we tried to take channels whose energy boundaries were as close as possible .
due to the aging processes mentioned earlier , observations with the same source spectrum that are made more than a few weeks apart end up at a different location in the cd . to correct for this effect we analyzed a number of _
rxte_/pca observations of the crab pulsar / nebula ( which is assumed to have a constant spectrum , see also @xcite ) that were taken around the time of our gx 17 + 2 observations . for all crab observations we produced the colors in the same energy bands as we used for gx 17 + 2 .
we found that the observed colors of the crab indeed changed .
for each observation we calculated multiplicative scaling factors , for both colors , with respect to those of the first crab observation .
we used these factors to scale the colors of the gx 17 + 2 observations back to those of the first gx 17 + 2 observation .
this procedure could only be applied to the epoch 3 observations , where , as expected , we found it to result in a narrower track in the cd .
no corrections were applied for the intensity . since no crab observations were available for the march 2000 observations , no corrections could be applied to epoch 4 data .
the shifts in that data set appeared to be small anyway .
our power spectral analysis was based on selecting observations as a function of the position along the z track in the cd .
we used a method that is based on the rank number and @xmath3 parameterization methods introduced by @xcite and @xcite , that has been gradually refined in similar studies @xcite . in this method in its current form , all points in the cd are projected onto a bicubic spline ( see * ? ? ?
* ) whose normal points are placed by hand in the middle of the z track ( see figure [ cd - hid_fig ] ) .
the hb / nb and nb / fb vertices are given the values @xmath3=1 and @xmath3=2 , respectively .
the rest of the z track is scaled according to the length of the nb .
problems arise when sharp curves are present in the track , in our study most notably around the nb / fb vertex . due to scatter , points which have fb properties end up on the nb and vice versa .
this is not a problem that is only intrinsic to the @xmath3 parameterization ; it is a limitation that applies to all selection methods based on colors .
no observable parameter , apart from the power spectra , has been identified that could be used to distinguish nb and fb observations around the vertex better than the x - ray colors .    in order to improve on the results of @xcite we wished to combine the epoch 3 and epoch 4 data sets . unfortunately , the tracks traced out in the cds of epoch 3 and 4 are not the same and two different splines had to be used for the @xmath3 parameterization .
since the normal points for these splines were drawn by hand , the @xmath3 scales for the epoch 3 and epoch 4 cds were unlikely to be exactly the same .
therefore , we first transformed the @xmath3 scale of epoch 4 to that of epoch 3 . to accomplish this , we measured the frequency of either the hbo or nbo at several places along the z track of epoch 3 , and determined the @xmath3 interval corresponding to the same frequency in the epoch 4 data .
we found that @xmath4 , showing that , although the scales are the same ( within the errors ) , a small shift is present .
we subsequently scaled the @xmath3 values of epoch 4 using the above expression .
the above scaling method assumes that the hbo and nbo frequencies are strongly related to @xmath3 .
previous studies of z sources ( see e.g. * ? ? ?
* ; * ? ? ?
* ; * ? ? ?
* ; * ? ? ?
* ) , as well as the fact that only little scatter was present around the linear relation seem to confirm this .
power spectra were created from the data in modes with high time resolution ( @xmath5 ; see table [
gx17 + 2_modes_tab ] ) using standard fast fourier transform techniques ( see * ? ? ? * and references therein ) .
the data were not background subtracted and no dead time corrections were applied prior to the fourier transformations .
we made power spectra in several energy bands , with several frequency resolutions and nyquist frequencies .
we finally settled on 0.06254096 hz power spectra in the 5.160 kev ( epoch 3 ) and 5.860 kev ( epoch 4 ) bands , since the qpos were most significantly detected in those bands ( note that the source does actually not contribute significantly to the power spectra above @xmath040 kev , however , our data modes did not allow us to discard data above this energy ) .
this choice of frequency range means that the properties of the vlfn , which dominates the power spectrum below 0.1 hz , could not always be measured satisfactorily , but it does allow us to follow the power spectral evolution on time scales down to 16 s. note that the power spectra of epoch 4 were produced in a slightly higher energy band than those of epoch 3 .
due to the gain changes and the limited choice of energy channels this was the best possible match .
the power spectra were selected on the basis of @xmath3 , as determined from the cd .
the @xmath3 selections usually had a width of 0.1 and did not overlap - no power spectrum was represented in more than one @xmath3-selection .
different widths were used in cases where the power spectrum changed rapidly as a function of @xmath3 ( narrower selections ) or when the powers were weak ( wider selections ) .
all the power spectra in a selection were averaged , and the resulting power spectrum was rms normalized according to a procedure described in @xcite .    the properties of the power spectra were quantified by fitting functional forms to them .
the low frequency ( 0.0625256 hz ) and high frequency ( 1004096 hz ) parts of the power spectrum were fitted separately .
the high frequency part was fitted with one or two lorentzians ( for the khz qpo(s ) ) and with the function @xmath6 for the dead time modified poisson level @xcite ; no separate term was used for the contribution of the very large events ( vle ) count rate since it was absorbed by this fit function . at low frequencies the noise that was present
could not be fitted consistently , and varying combinations of functionals had to be used :    * @xmath3@xmath10.0 : a lorentzian ( lfn+vlfn ) * @xmath3=0.00.1 : a cut - off power law ( lfn+vlfn ) * @xmath3=0.11.4 : a power law ( vlfn ) and a cut - off power law ( lfn ) * @xmath3=1.41.6 : a power law ( vlfn ) and a lorentzian ( lfn+nbo ) * @xmath3=1.65.0 : a power law ( vlfn ) and a lorentzian ( lfn ) * @xmath3@xmath75.0 : a power law ( vlfn )    the expression for a power law is @xmath8 , that for a cut - off power law is @xmath9 ( where @xmath10 is the cut - off frequency ) , and for a lorentzian @xmath11 $ ] ( where @xmath12 is the centroid frequency and @xmath13 is the full - width - at - half - maximum ) .
the hbo , its second harmonic and the nbo / fbo were each fitted with a lorentzian .
below the hbo , at about half its frequency , a broad bump was present that was also fitted with a lorentzian , whose frequency was sometimes fixed to zero . at low frequencies the dead time modified poisson level was fitted with a constant .
an example of a fit to a power spectrum ( @xmath3=0.40.5 ) , displaying the contributing of the several components , is shown in figure [ pds_comp_fig ] .
note that changes in the fit function , such as using a cut - off power law instead of a lorentzian or adding an extra component , may lead to minor changes in the values of other parameters .
errors on the fit parameters were determined using @xmath14 .
upper limits were determined by fixing some or all of the parameters of a component , except the rms amplitude , to values similar to those obtained in the closest @xmath3  selection where it was found to be significant , leaving all other fit parameters free , and using @xmath15 ( 95% confidence ) .
a study of the energy dependence and time lag properties of the qpos and noise components will be presented elsewhere .
the cds and hids of both epochs are shown in figure [ cd - hid_fig ] . in both cds the hb / nb vertex is not well defined , since the hb is almost a continuation of the nb .
this is mainly due to the relatively high energies that we chose for our soft color ( see section [ 17 + 2_obs_sec ] ) . at lower energies
the turn - over is much clearer . in figure [ sz - count_fig ]
we show the count rate in several energy bands as a function of @xmath3 , for epoch 4 .
it shows that the hb / nb vertex in the hid is entirely due to the count rates at low energies . at high energies
the hb is a perfect continuation of the nb , and no vertex is present .    in the hid of epoch 3
the z track appears to be segmented , which was already noted by @xcite .
although no corrections for the slow detector aging processes were applied to the intensity in this hid , we note that a quick look showed that they would only have made the segmentation more apparent .
the shifts in the hid might be due to the secular motion that has been observed in other z sources .
the shifts , of up to @xmath05% in intensity , do not show up in the cd .
this is because colors are ratios of intensities and are therefore not very sensitive to overall intensity changes ; exactly for this reason cds are preferable over hids for our purpose .
moreover , the width of the tracks in the cd is about 5% , so changes smaller than this are hard to observe .
the bottom panel of figure [ hist - vz_fig ] shows the distribution of the time spent by the source in each part of the z track .
the source spent 28% of the time on the hb ( @xmath3@xmath161 ) , 44.2% on the nb ( 1@xmath1@xmath3@xmath162 ) , and 27.8% on the fb ( @xmath3@xmath72 ) .
the average speed along the z track as a function of @xmath3(@xmath17 ) is shown in the top panel of figure [ hist - vz_fig ] .
the speed at a given @xmath18 is defined as @xmath19 ( see also * ? ? ?
* ) , where @xmath20 is used to number the points in order of time .
as expected the @xmath17 increases considerably when the source enters the fb , but also at the top of the hb ( @xmath21 ) .
combined with the small amount of time spent in the upper hb we can conclude that the source reaches to @xmath3  values this low only in brief , quick dashes .    for a more detailed analysis of the spectral changes as a function of the position along the z track in gx 17 + 2 we refer to et al .
( in prep . ) .
figure [ pds_fig ] shows the low frequency part ( 0.0625256 hz ) of the power spectrum for nine different @xmath3 selections . the power spectra shown in figure [ pds_fig ]
are selected from epoch 3 and epoch 4 and are therefore a combination of 5.160 kev and 5.860 kev data .
the percentage of epoch 3 and epoch 4 data varies between the @xmath3-selections .
the contribution of epoch 3 data is highest at the extremes of the @xmath3  range , reaching 100% , and gradually decreases from both ends to @xmath010% around @xmath3=1 .
several components can be seen , which change in strength and shape as a function of @xmath3 .
the different components are identified in figure [ pds_fig ] ( see also figure [ pds_comp_fig ] ) and the @xmath3ranges in which they were detected are given in table [ sz - comp_tab ] .
the component identified as the sub - harmonic of the hbo is rather broad and does not have the appearance of a qpo .
however , based on the frequency ratios ( see below ) it is referred to as the hbo sub - harmonic .    some difficulties were experienced while fitting the low frequency part of the power spectrum between @xmath3=1.4 and @xmath3=1.6 , where two or more components with similar frequencies were simultaneously present . above @xmath3=1.4
the nbo appeared , as a broad feature on top of the lfn .
we were not able to distinguish the two components and decided to fit them together with a single lorentzian .
the fit values are not used in figures and tables , since they do not represent any of the individual components .
above @xmath3=1.6 , the nbo and lfn could be distinguished more easily .
the fit function used for the lfn , which was underlying the nbo , was changed to a lorentzian to be more consistent with fits at higher @xmath3  ( fits with a cut - off power law gave equally good @xmath22 at this @xmath3 ) .    in the following sections the results for each of the components will be presented .
the hbo was detected between @xmath3=@xmath230.6 and @xmath3=2.1 and its second harmonic between @xmath3=@xmath230.6 and @xmath3=1.0 .
the second and third columns of figure [ hbo_fig ] show their rms amplitudes , fwhm , and frequencies as a function of @xmath3(see also table [ low_tab ] ) .
the frequency of the hbo increased from 21.3 hz at @xmath3=@xmath230.43 to 60.3 hz at @xmath3=1.45 and then decreased to 48.2 hz at @xmath3=2.04 .
the second harmonic of the hbo had a frequency that was on average 1.941@xmath240.007 times that of the first harmonic , which is significantly different from the expected value of 2 .
an explanation for this discrepancy is proposed in section [ disc - hbo_sec ] . between @xmath3=@xmath230.6 and @xmath3=1.0 both the hbo and its second
harmonic decreased in strength , with the decrease of the second harmonic occurring somewhat faster than that of the hbo . in that same @xmath3interval
the fwhm of the hbo and its second harmonic were fairly constant ( showing a slight increase ) , although considerable scatter was present around the average values ( which were , respectively , 7.9@xmath240.3 hz and 14.4@xmath240.6 hz ) .
the q values of the hbo and its harmonic are shown in figure [ q - values_fig ] . in the @xmath3 range where both were detected their q values were consistent with each other .
when the second harmonic was not significantly detected anymore ( @xmath3@xmath71.0 ) the rms amplitude and the fwhm of the hbo both decreased from , respectively , 2.45@xmath240.08% and 9.7@xmath247 hz ( @xmath3=0.51.0 ) to 1.48@xmath240.06% and 5.1@xmath240.4 hz ( @xmath3=1.01.5 ) . between @xmath3=1.0 and @xmath3=1.5
the relation between the @xmath3 and the hbo frequency started to flatten . above
@xmath3=1.5 the frequency of the hbo clearly dropped , initially quite abruptly and later on more smoothly .
this frequency drop coincided with an increase in the fwhm to 15.7@xmath241.5 ( @xmath3=1.52.1 ) ; the rms amplitude showed a small increase to 2.1@xmath240.1% ( @xmath3=1.76 ) , followed by a decrease to 1.2@xmath240.2% ( @xmath3=2.04 ) .    underlying the hbo and its second harmonic , but with a central frequency lower than that of the hbo , we found a broad feature that was also fitted with a lorentzian .
it was significantly detected between @xmath3=@xmath230.6 and @xmath3=1.3 .
the properties of the broad feature are shown in the first column of figure [ hbo_fig ] ( see also table [ low_tab ] ) .
below @xmath3=0 the frequency of the lorentzian was fixed to zero . between @xmath3=0.5 and @xmath3=1.5
the frequency of the broad feature is on average 0.539@xmath240.015 times that of the hbo .
this could mean that the broad feature is the sub - harmonic of the hbo , certainly when one takes into account that the frequency of a broad feature is rather sensitive to the shape of the continuum . whereas the rms amplitude and the frequency of the sub - harmonic both change strongly with @xmath3 the fwhm remains more or less constant , with an average value of 29.5@xmath241.5 hz .      the nbo and fbo were detected between @xmath3=1.6 and @xmath3=2.7 .
the fit results for the nbo / fbo are shown in figure [ nbo_fig ] ( see also table [ low_tab ] ) .
below @xmath3=2.1 the nbo ( represented by the filled circles in figure [ nbo_fig ] ) had a fairly constant frequency , with values between 6.3 hz and 7.0 hz .
its rms amplitude increased from 1.7% to 3.2% , and the fwhm was approximately @xmath02.1 or @xmath03.8 . in the @xmath3=2.02.1 selection both the nbo and fbo were present .
this is likely an artifact of the @xmath3selection method , since a careful inspection of dynamical power spectra showed no evidence for simultaneous presence of both qpos .
it is interesting to note though , that in the @xmath3=2.02.1 selection the frequency of the fbo was 2.0@xmath240.1 times that of the nbo , which could mean that the two nbo and fbo are harmonically related ( however , see below ) .
the fbo increased in frequency from 13.9 hz ( @xmath3=2.04 ) to 23.1 hz ( @xmath3=2.65 ) , while its fwhm increased from @xmath05 hz to @xmath013 hz .
the rms amplitude of the fbo initially continued the trend of the nbo rms amplitude ; it increased from 2.8% ( @xmath3=2.04 ) to 5.6% ( @xmath3=2.23 ) , but then decreased to 1.8% ( @xmath3=2.65 ) .    to study the transitions between the nbo and fbo more carefully , we inspected all dynamical power spectra of observations with @xmath3values around 2 .
although no clear transitions were found , mainly due to the limited quality of the dynamical power spectra , we did in some cases see qpos with intermediate frequencies ( @xmath0 10 hz ) , suggesting that the frequency does not jump directly from @xmath07 hz to @xmath014 hz .
fitting the power spectra in consecutive time intervals ( longer than 16 s ) , rather than inspecting dynamical power spectra with the eye , will probably lead to more conclusive results , but is beyond the scope of the current paper .
the time scales on which the nbo / fbo frequency changed from @xmath07 hz to @xmath014 hz and back were as short as a few tens of seconds .
we also studied the behavior of the nbo / fbo during two long type i x - ray bursts .
the first one started on 1998 november 18 at 08:51:26 utc , the second one on 1999 oct 10 at 09:10:47 utc .
their exponential decay times were , respectively , 189@xmath242 s and 144@xmath242 s @xcite .
both bursts occurred near the nb / fb vertex and in the power spectra of their respective observations the nbo / fbo is clearly detected .
no other qpos were detected in the power spectra of these bursts .
figure [ burst - nbo_fig ] shows the dynamical power spectra of both bursts , together with their 260 kev light curves . during the brightest part of the bursts the nbo / fbo seemed to disappear
apparently , the burst flux was not modulated at the nbo / fbo frequencies with the same amplitude as the persistent flux . to quantify this , we determined upper limits on the nbo / fbo strength and compared those with the values outside the bursts .
the results are shown in table [ burst - nbo_tab ] .
clearly , during the brightest part of the bursts the rms amplitude of the nbo / fbo was significantly weaker , not only as a fraction of the total flux , but also in absolute terms ( and hence as a fraction of the persistent flux , if the persistent emission is assumed to continue during the bursts ) .
these measurements constitute the first determination of the effect of x - ray bursts on the nbo / fbo .
the fact that the bursts suppress the qpo can have important consequences for our understanding of its formation mechanism ( see section [ disc - nbo_sec ] ) .
lfn was detected between @xmath3=@xmath230.6 and @xmath3=5.0 . as mentioned before ,
the appearance of the nbo around @xmath3=1.4 , kept us from putting firm constraints on the lfn parameters between @xmath3=1.4 and @xmath3=1.6 .
the fit results for the lfn are shown in figure [ lfn_fig ] ( see also table [ tab : noise ] ) , in two separate columns : one for the fits with a cut - off power law fit ( @xmath3=0.01.4 ) and one for the fits with a lorentzian ( @xmath3@xmath10.0 and @xmath3=1.65.0 ) .
the strength of the lfn was in both cases defined as the integrated power spectral density between 1 hz and 100 hz .
we note that below @xmath3=0.1 the vlfn was not fitted separately from the lfn .
an inspection of 1/2564096 hz power spectra below @xmath3=0.1 showed that a weak power law component was present at frequencies below 0.1 hz .
this component was probably vlfn ; its power in the 1100 hz range was much smaller than that of the lfn , so , although some vlfn power was absorbed by the lfn , this did not affect the lfn rms amplitudes significantly .
the strength of the lfn changed considerably as a function of @xmath3 .
it showed a narrow peak between @xmath3=@xmath230.6 and @xmath3=0.2 of @xmath05% rms . between @xmath3=0.2 and @xmath3=1.7 it gradually decreased from @xmath04.5% rms to @xmath02.5% rms .
another decrease was observed between @xmath3=2.7 and @xmath3=5.0 , from @xmath03.8% rms to @xmath00.9% rms .
the behavior between @xmath3=1.7 and @xmath3=2.7 was quite erratic , probably due to interactions with the fit functions of nbo and fbo .
the centroid frequency of a lorentzian and the cut - off frequency of a cut - off power law can not be directly compared .
since we wanted to see how the typical frequency of the lfn evolved with @xmath3 , we chose to plot @xmath25 , which is the maximum in a @xmath26 plot and the frequency at which most of the power is concentrated per logarithmic frequency interval ( belloni , psaltis , & van der klis 2001 , in prep . ) .
for a lorentzian @xmath25 is @xmath27 and for a cut - off power law @xmath25 is @xmath28 ( see section [ 17 + 2_obs_sec ] for analytical expressions for a lorentzian and cut - off power law ) . as can be seen from figure [ lfn_fig ] , between @xmath3=@xmath230.6 and @xmath3=1.7 @xmath25 smoothly increased from @xmath01.2 hz to @xmath014.7 hz .
above @xmath3=1.7 the errors on @xmath25 were larger and the behavior was less clear
. between @xmath3=1.7 and @xmath3=2.3 @xmath25 decreased to @xmath06.5 hz , and above @xmath3=2.3 it increased again to @xmath015 hz .
we tested whether the change of the fit functions at @xmath3=0.0 affected the values for @xmath25 and the other power spectral parameters , by swapping the fit functions used below and above @xmath3=0.0 ; no significant changes were found .
vlfn was detected over almost the whole @xmath3  range .
although we only started fitting the vlfn separately from the lfn above @xmath3=0.1 ( in the 16 s power spectra ) , vlfn was present at frequencies below 0.1 hz in the 256 s power spectra below @xmath3=0.1 . in the 16 s power spectra
it was only significantly detected above @xmath3=0.3 .
the fit results for the vlfn are shown in figure [ vlfn_fig ] ( see also table [ tab : noise ] )
. between @xmath3=0.3 and @xmath3=1.4 the vlfn strength decreased from @xmath00.9% rms to @xmath00.4% rms . above @xmath3=1.4 its strength increased , to a peak value of @xmath01.3% rms at the nb / fb vertex .
on the lower fb the strength decreased again , to a value of @xmath00.5% rms , and above @xmath3=3.0 it increased to @xmath01.6% rms at the top of the fb .
the index of the vlfn slowly increased from @xmath00.6 to @xmath01.0 between @xmath3=0.1 and @xmath3=2.4 , with a small peak between @xmath3=1.6 and @xmath3=1.9 , where the index had a value of @xmath01.5 . between @xmath3=2.4 and @xmath3=3.2
the vlfn was much steeper .
the indices were not well constrained and had values between 1.7 and 4.0 . above @xmath3=3.2 , where the error bars are much smaller
, the index slowly decreased from 2.1 to 1.8 .
our fits did not require an additional component to fit possible hfn ; no such component was found .
possible explanations for this are discussed in section [ disc : noise ] .
both the lower and upper khz qpo were clearly detected ; the lower khz qpo between @xmath3=0.5 and @xmath3=1.5 , the upper khz qpo between @xmath3=0.3 and @xmath3=1.7 .
the results can be found in table [ khz_tab ] and are shown in figure [ khz_fig ] .
the qpo frequencies showed a clear increase with @xmath3 , although both relations flattened at their low frequency ends .
the fwhm of the lower khz qpo was consistent with being constant at @xmath0100 hz , whereas the rms amplitude showed a peak near @xmath3=1.05 with a value of 3.6@xmath240.3 .
the rms amplitude and fwhm of the upper khz qpo both decreased with @xmath3 .
figure [ khz - diff_fig]a shows the frequency difference of the two khz qpos as a function of the frequency of the upper khz qpo .
some fits were made to the data ; they are also shown in figure [ khz - diff_fig]a .
the best fit to the frequency difference with a constant gave a value of 282@xmath244 hz . the @xmath29 for this fit was 18.6/9 , which means that at a 97% confidence level the frequency difference was not constant .
this is the first time that this is observed for gx 17 + 2 .
fits with first and second order polynomials ( taking into account errors in both coordinates ) resulted in , respectively , @xmath29=15.0/8 and @xmath29=3.7/7 .
although the latter two fits show that with 99.76 % confidence ( 3.2@xmath30 ) a decrease towards higher frequency is not monotonic , it is not clear either whether the decrease towards lower frequencies is significant or not .
this would be the first time that such a decrease towards lower frequencies is observed in any khz qpo source with a non - constant frequency difference . to test this , two fits with a broken line ( not shown )
were performed , where in one case the slope of the low frequency part was fixed to zero .
they resulted in @xmath29=6.28/7 ( slope fixed ) and @xmath29=2.52/6 ( slope free ) , which suggests that the frequency decrease towards lower frequency is only significant at a 97.6% confidence level ( 2.4@xmath30 ) .    in figure
[ khz - diff_fig]b we show additional fits to the data with theoretical curves for the radial epicyclic frequency vs.azimuthal frequency , which in the relativistic precession models are associated with the frequency difference and the upper khz qpo frequency .
the solid line is the best fit for circular orbtis @xcite , which yields a @xmath22 of 7.3 ( d.o.f = 9 ) for a mass of 2.07 m@xmath31 ; the dashed line is the best fit for eccentric orbits @xcite , yielding a @xmath22 of 1.14 ( d.o.f = 8) for a mass of 1.86 m@xmath31 .
both curves were calculated for non - rotating stars ; rotation of the neutron star gives only insignificant improvements . for comparison
we also plotted the data for sco x-1 from @xcite .
clearly , both fits ( especially the solid line ) do not fit the decrease of the frequency difference towards lower frequencies , which occurs in a frequency range for the upper khz qpo that was not observed in sco x-1 , very well .
the q values ( frequency / fwhm ) of the two khz qpos were consistent with each other and with that of the hbo and its harmonic ( see figure [ q - values_fig ] ) ; between @xmath3=1.0 and @xmath3=1.5 the q values of the hbo and khz qpos increased from @xmath05 to @xmath010 .
this unexpected finding may provide a key to understanding the formation mechanism of the three qpos  this will be further discussed in section [ q - values_sec ] .    when comparing figures [ khz_fig ] and [ hbo_fig ] one can see that above @xmath3=1.5 , where the hbo frequency starts to decrease , the upper khz qpo frequency still increases . in figure [ khz - hbo_fig ] we plot the frequencies of both qpos against each other . for values of the upper khz frequency
lower than 1030 hz the frequencies are well correlated , but above there is a clear anti - correlation .
the solid line is the best power law fit to the points below 1000 hz .
the power law index is 2.08@xmath240.07 .
this is the first time in any atoll or z source showing khz qpos that an anti - correlation is observed between the frequencies of the low and high frequency qpos ( see section [ disc - hbo_sec ] for a discussion ) .
we performed a detailed study of the low and high frequency power spectral features of the neutron star low - mass x - ray binary and z source gx 17 + 2 .
as was found in previous studies of gx 17 + 2 and other z sources , the properties of most power spectral features correlated well with the position of the source along the z track in the x - ray color - color diagram , suggesting that the spectral and variability changes are caused by one and the same parameter .
some new results were found , the most interesting being the fact that the frequency separation of the khz qpos was not constant and probably not monotonic ( fig .
[ khz - diff_fig ] ) , that their q values were consistent with those of the hbo and its harmonic ( fig.[q - values_fig ] ) , and that the frequency of the upper khz qpo was anti - correlated with that of the hbo when the latter started to decline on the nb ( fig .
[ khz - hbo_fig ] ) .
we also found a sub - harmonic of the hbo and showed that the absolute amplitude of the nbo / fbo is suppressed during type i x - ray bursts ( fig .
[ burst - nbo_fig ] ) . in the remainder of this section
these and other results will be discussed in more detail . for a more detailed discussion , including a comparison of our results with those obtained from with other satellites ( like _ exosat _ and _ ginga _ ) and of other z sources
, we refer to @xcite .        in our discussion of the hbo and khz qpo results
, we will only consider models that try to explain the behavior of both types of qpos .
currently these models are the relativistic precession model , the sonic point beat frequency model and the two - oscillator model . in the latter model
( see * ? ? ?
* ) an angle between the magnetospheric equator and the disc plane is defined , which depends on the frequencies of the hbo and khz qpos and is supposed to be constant .
however , @xcite show , by using our gx 17 + 2 results and theirs of gx 5 - 1 , that this is definitely not the case . for this reason
, we will only consider the first two models . in the relativistic precession model ( rpm , @xcite )
three frequencies are predicted ; they correspond to the fundamental frequencies of motion of a test particle in the vicinity of a compact object and are identified with the hbo and the two khz qpos . the lowest predicted frequency , which should explain the hbo , is a factor 2 lower than the low frequency qpos in the atoll sources and the sub - harmonic of the hbo in the z sources .
hence , a way to produce second harmonics is required .
note that no actual mechanism is proposed for the production of qpos themselves .
however , in a recent model by @xcite the accretion disk acts as low band - pass filter with resonances close to the three frequencies predicted in the rpm . it is unclear to what extent the details of the rpm ( i.e. the turnover in the relation between hbo and upper khz frequency ) carry over to the @xcite description .
we also note that recent investigations of the rpm by @xcite have cast some doubt on the frequencies predicted by the rpm and relations between them . in the sonic point beat frequency model ( spbfm , * ? ? ?
* ; * ? ? ?
* ) the hbo is explained by the original beat frequency model of @xcite and @xcite .
the observed khz frequencies are close to the orbital frequency at the sonic point of this disk and the beat of this frequency with neutron star spin frequency .
the model requires the presence of a magnetosphere and a solid surface , and can therefore not explain the low and high frequency qpo in the black hole candidates .
the decrease of the hbo frequency on the nb in gx 17 + 2 was first observed in _ rxte _
data that were not used in our analysis @xcite .
there were interesting features in the hbo properties which were reproduced in that data and ours . in both cases
the @xmath3 dependence of the hbo frequency was not symmetric around its peak value ; below , the frequency slowly flattened off to its maximum value , but above it showed a fast decrease by a few hz that was followed by a more or less linear decrease with @xmath3 .
also , simultaneously with the frequency drop , the fwhm increased by a factor of more than two .
possible explanations for the hbo frequency decrease are discussed in section [ hbo - khz_disc ] .    in all z sources where
an harmonic of the hbo has been observed , the frequency ratio of the harmonic and the fundamental is generally less than two .
this can be explained , at least in gx 17 + 2 , by looking at the @xmath3 dependence of the rms amplitudes of both components . in a certain @xmath3 selection
both frequencies tend to be dominated by those with the higher rms amplitudes , and hence by lower @xmath3 and lower frequencies ; this effect will be stronger for the harmonic than for the fundamental , since the @xmath3  dependence of its rms amplitude is steeper ( see figure [ hbo_fig ] ) . as a result
the frequency ratio will end up to be less than two . to test this hypothesis we calculated the weighted averages of two harmonically related linear functions , representing the hbo frequency-@xmath3relation between @xmath3=0.0 and @xmath3=1.0 .
as weight we used linear functions that were fit to the hbo rms-@xmath3  relation between @xmath3=0.0 and @xmath3=1.0 .
the ratio of the weighted frequency averages was 1.93 , which is very close to the observed value of 1.94
.    peaked features with frequencies intermediate to those of the lfn and the hbo have been reported for gx 340 + 0 @xcite , sco x-1 @xcite , gx 5 - 1 @xcite , and now also for gx 17 + 2 .
based on the average frequency ratio of this sub - hbo component and the hbo in gx 17 + 2 , 0.539@xmath240.015 , we suggest that , similar to what was suggested for gx 340 + 0 , it may be the sub - harmonic of the hbo .
a comparative study by @xcite of low frequency ( hbo - like ) qpos and band limited noise components in the neutron star z and atoll sources and black holes revealed strong correlations between the typical frequencies of those components . when plotting the frequency of the low frequency qpo versus the break frequency of the band limited noise they found that the branch traced out by the z sources , although lying parallel to that of the atoll source and black holes , was slightly offset .
it is interesting to see that this discrepancy between the z sources on the one hand and the atoll sources and black holes on the other disappears @xcite , at least for gx 17 + 2 , sco x-1 and gx 340 + 0 , when one uses the frequencies of the sub - hbo component , instead of the hbo frequency ( see also * ? ? ?
if the low frequency qpo in the atoll sources and the sub - hbo component of the three z sources would indeed be the same qpo , this would mean that the(se ) z sources do not anymore follow the relation between the low frequency ( hbo - like ) qpo and the lower khz qpo as defined by atoll and black hole systems @xcite .
it is interesting to see that the three z sources would then end up near or slightly above a second relation that is traced out by atoll sources like 4u 1728 - 34 and 4u 160852 . at the moment
, however , it seems that the relation found by @xcite and the main relation of @xcite are mutually exclusive for the two types of sources .
note that @xcite found that if in their plot for sco x-1 the frequencies of the sub - hbo component and the hbo were used , instead of the lfn and hbo frequencies , the discrepancy between sco x-1 and the non - z sources also disappeared .
this reinterpretation has the advantage of preserving the consistency with the @xcite relation as well .
the dominance of the even harmonics over the odd harmonics in the low frequency ( hbo - like ) qpos , as suggested by the presence of a sub - harmonic of the hbo , is not only observed in the z sources , but also a common feature in bhcs such as gs 112468 @xcite and xte j1550564 @xcite .
it suggests that a two - fold symmetry is present in the production mechanism of these low frequency qpos . for
the spbfm one could think of an asymmetry in the magnetic field , which results in different areas for the polar caps ( see e.g. * ? ? ?
* ) . in the rpm
a two - fold symmetry is introduced by the two points at which the slightly inclined orbits of test particles cross the plane of the accretion disk .
the khz qpos in gx 17 + 2 were discovered with _
rxte _ @xcite . in our data
the lower khz qpo was found between 517 hz and 794 hz ( @xmath3=0.51.5 ) and the upper khz qpo between 618 hz and 1087 hz ( @xmath3=-0.31.7 ) , similar to what was reported by @xcite .
our findings of the khz qpos in gx 17 + 2 are compatible with the general properties of khz in other z sources @xcite .
an important new result is that the frequency difference is no longer consistent ( at a 97% confidence level ) with being constant ( see figure [ khz - diff_fig ] ) .
the frequency separation reached a maximum value of 308 hz , when the upper peak frequency was about 915 hz , and fell off in both directions to values of 263 hz and 239 hz , when the frequencies of the upper peak were , respectively , 780 hz and 1033 hz .
note that the decrease in both directions , when fitted with a parabola , resulted in significantly ( 3.1 @xmath30 ) better fits than that with a constant .
the decrease in the frequency separation at the low frequency side of the upper khz qpo , although itself only significant at a 2.3@xmath30 level , has not been seen before in any khz qpo source . in the other lmxbs that show an non - constant frequency separation , sco x-1 @xcite , 4u 160852 @xcite , 4u 173544 @xcite , and 4u 172834 @xcite , the frequency separation was not measured at upper peak frequencies as low as in gx 17 + 2 .
although both the spbfm and the rpm are capable of explaining the observed decrease in frequency difference towards higher upper peak frequencies , they have severe problems with a decrease towards lower frequencies .
such a decrease is not possible in the spbfm , and , although it is in principle predicted by the rpm , that model only yields reasonale fits to our data when one allows for highly eccentric orbits ( see figure [ khz - diff_fig]b ) .
the principal aim of our observing campaign in 1999 was to check on the relation between hbo frequency ( whose increase was known to turn into a decrease towards higher @xmath3 ) and upper khz qpo frequency .
clearly , a decrease , or even flattening , of the upper khz qpo frequency , is not observed to accompany the decrease in hbo frequency ( see figures [ hbo_fig ] , [ khz_fig ] , and [ khz - hbo_fig ] ) . in a pure rpm , this can not be explained .
however , if classical precession due to oblateness of the neutron star also plays a role , such an explanation is possible .
for certain combinations of neutron star parameters , the prediction is that above a certain value of the upper khz qpo frequency , the frequency of the hbo in fact starts to decrease ( see * ? ? ? * ) , due to the increasing importance of the classical precession . in gx 17 + 2
this turn - over of the hbo frequency occurs at a value of @xmath01 khz for the upper khz qpo with a peak frequency of @xmath060 hz for the hbo .
judging from figure 2 in @xcite , the relation between the hbo and upper khz qpo frequencies in gx 17 + 2 would require a high neutron star spin frequency ( @xmath0500700 hz ) , a neutron star mass of more than 2 @xmath32 , and a hard equation of state ( assuming the observed hbo frequency is four times the predicted frequency ) . in the rpm ,
well below the turn - over , the hbo frequency is expected to scale quadratically with the upper khz qpo frequency .
the best power law fit to the data below 1000 hz yielded an index of 2.08@xmath240.07 ( see figure [ khz - hbo_fig ] ) , which is consistent with a quadratic relation .    in the spbfm
a turnover in the hbo - upper khz qpo relation is not predicted , but as the two frequencies do not directly depend on each other , the model is flexible enough to account for it ; in case the turnover in hbo frequency can be explained by the two flow model put forward by @xcite , the steady increase of the upper khz qpo frequency could in principle be explained by assuming that most of the radial flow falls back to the disk and rejoins the disk flow before the sonic point is reached , so that @xmath33  at that location , and hence the upper khz qpo frequency keep increasing .
however , since the two radii at which the different qpos are generated are very close , this seems a rather unlikely option .
the q values of the khz qpos are remarkably similar to each other , and to those of the hbo and its harmonic ( see figure [ q - values_fig ] ) .
they also show an almost simultaneous increase above @xmath3=1 , when the hbo harmonic is no longer detected . if the qpos were artificially broadened by our selection method , averaging together peaks with different centroid frequencies
, one would expect correlations of the fwhm ( and q value ) with the width of the @xmath3 selection and with @xmath34 .
such correlations , if present at all , are far too weak to explain the observed similarities
.    the fwhm of a qpo can be caused by several mechanisms ; intrinsic frequency variations , phase jumps , the simultaneous presence of several frequencies , and a finite lifetime of the signal will all cause a signal to appear as a broadened peak in a power spectrum .
the data quality did not allow us to test the phase jump option .
lifetime broadening can probably not explain our results .
in the spbfm there is no reason why the lifetimes of the low and high frequency oscillations , which are determined by different processes , would lead to similar q values . in the rpm , where all frequencies are generated by a single blob of matter
, the lifetime for all variations should to first order be the same
. however , in that case the q values of the khz qpos should be much higher ( by a factor 1020 ) than those of the hbo , since more khz cycles fit in a blob s lifetime .
of course , it is possible that the lifetime of the oscillations is not determined by the lifetime of the blob , but rather by a damping factor . assuming that these factors are not the same for the different qpos , the similarity of the q values could in principle be explained .
however , there is no obvious reason why the ( independent ) damping factors would be fine tuned in a such a way that they lead to similar q values .    explaining the q values by broadening due to the presence of several frequencies
does not work either .
below the turnover in the hbo - upper khz qpo frequency relation , the frequency of the hbo scales quadratically with the frequency of the upper khz qpo
. a given range of upper khz frequencies would therefore lead to a q value for the hbo that is half the q value for the upper khz qpo .
this is not consistent with what we find .
frequency modulation could in principle explain the similar q values , but only if the modulation is caused by a different parameter than the one that determines the frequency changes along the z track .
frequency modulation only produces similar q values if the qpo frequencies vary with the modulating parameter by the same factor , i.e.  proportionally to each other . the observed quadratic relation between the hbo frequency and that of the upper khz qpo as well as the subsequent turnover
, shows that this is not the case for the parameter that causes the motion along the z track .
however , assuming that the modulating parameter is not responsible for motion along the z ( see section [ disc - spec ] ) , we discuss the possibility of frequency modulation in a bit more detail . for broadening by frequency modulation to be observable in our power spectra , the time scale of these variations should be longer than that of the hbo . additionally , the strength of those variations should anti - correlate with the q values of the qpo . assuming that the parameter that causes the frequency modulations also causes some changes in the count rate we could try to identify it in the power spectrum .
the only reasonably strong components that we see in our power spectra and whose typical frequencies fulfill the above frequency requirement ( @xmath35 hz ) are the vlfn , the lfn and the nbo / fbo . the vlfn might be associated with the motion along the z ( see section [ disc : noise ] ) ; disentangling that effect from a possible additional one on qpo frequencies is beyond the scope of this work .
the strength of the lfn , which is of interest below @xmath3=1.5 , anti - correlates with the q values below @xmath3=1.5 . above @xmath3=1.5
the q values decrease again to a value between 2 and 5 . in that case
the broadening might be due to the appearance of the nbo / fbo and the increase in the lfn strength .
it is interesting to note that in sco x-1 it has already been found that the khz qpo properties are clearly modulated by the nbo @xcite .
for the first time we were able to study the properties of the nbo / fbo during type i x - ray bursts .
we found that the absolute rms amplitude significantly decreased when the radiation from the neutron star surface increased .
this is clearly different from the behavior of the @xmath01 hz qpos in the dipping lmxbs 4u 132362 @xcite , exo 0748676 @xcite , and 4u 174637 @xcite , whose absolute rms amplitudes increased during type i x - ray bursts .
it is thought that these qpos are not related to the nbo / fbo in the z sources .
the behavior of the nbo / fbo in gx 17 + 2 shows that nbo / fbo mechanism is very sensitive to the radiation field , thereby lending strong support for radiation feedback mechanisms ( see below ) . in the framework of such mechanisms
it is easy to understand that the nbo / fbo disappears during bursts : the nbo / fbo is only observed in a narrow @xmath3 range , suggesting that it requires a delicate balance of certain parameters , among which the radiation field .
the observations during the burst tails indicate that the nbo / fbo mechanism does not switch on instantly but rather becomes stronger in a more gradual way .
spectral fits during the bursts @xcite suggest that the spectral properties of the accretion flow itself are not affected much by the increase in radiation .
the frequency of the nbo / fbo changed from being more or less constant below @xmath3=2 to being strongly correlated with @xmath3 , above @xmath3=2 . in our data
set we found the nbo between @xmath3=1.6 and @xmath3=2.1 , although it may already have been present as a very broad feature as early as @xmath3=1.4 . between @xmath3=2.0 and @xmath3=2.1 ,
the nbo ( 6.86 hz ) was detected simultaneously with the fbo ( 13.9 hz ) . above @xmath3=2.1
the fbo frequency increased , to @xmath023 hz at @xmath3=2.65 .
note that the simultaneous presence of the nbo and fbo between @xmath3=2.0 and @xmath3=2.1 is likely an artifact of the @xmath3  selection method .
inspection of dynamical power spectra showed that they were never present at the same time .
it is interesting though , that in that @xmath3  selection their frequencies differed by a factor of 2 , which is suggestive of mode switching of the nbo / fbo .
however , there is evidence for qpos with intermediate frequencies , which undermines that idea . the only way to see whether the frequencies really gradually transform into each other is to study the power spectral changes as a function of time .
however , the quality of the gx 17 + 2 data is not sufficient to perform such a study .
the presence of the nbo / fbo in the z sources has often been related to the fact that these sources accrete at near eddington mass accretion rates .
it is thought that at these high mass accretion rates a significant fraction of the accretion flow is in the form of a thick and perhaps spherical flow .
also , the effects of radiation pressure are thought to play an important role in this regime .
the model for the nbo by @xcite is based on a radiation feedback mechanism in a spherical accretion flow .
another model for the nbo was proposed by @xcite , in which the nbo frequency is basically that of sound waves in a thickened accretion disk .
this model does not explain how the nbo changes into the fbo
. a major problem for both models could be the recently discovered nbo / fbo - like qpos in the two atoll sources 4u 182030 @xcite and xte j1806246 @xcite .
certainly in the first source the luminosity is believed to be only @xmath02040 % of the eddington luminosity , while near - eddington luminosities are required in both models .
the lack of hfn in our data , compared to , e.g. , in the _
exosat_@xcite data can be explained by instrumental effects and different approaches to fitting the data .
a study of _ exosat _ data of several sources by @xcite showed that a component with an rms amplitude of @xmath03% is always present in the frequency range where the hfn was often found in the _ exosat _  data .
when hfn was found at lower frequencies , it was , based on strength and frequency , most likely the component that we identified as the sub - harmonic of the hbo .
past measurements of the energy spectrum of the vlfn suggest that it is caused by motion of the source along the z @xcite .
the strength of the vlfn ( fig . [ vlfn_fig ] ) and the speed along the z track ( @xmath36 , see figure [ hist - vz_fig ] ) should therefore have a similar dependence on @xmath3 .
our results are inconclusive ; although some correlations are observed ( the increase above @xmath3=3 , and the steep increase around @xmath3=2 ) there are also some differences .
most notable is the dip in the vlfn strength around @xmath3=3 , which is not observed as a decrease in @xmath17 .
however , it is possible that due to the steepness of the vlfn around @xmath3=2.5 most of its power was outside the range in which we determined the strength .
the lfn followed the behavior of the hbo , both in strength and frequency , until the latter disappeared ( @xmath37 ) . between @xmath3=2.1 and @xmath3=2.8
the behavior is rather complex and probably considerably influenced by the presence of the nbo / fbo .
above @xmath3=3.0 the behavior was rather clean , with a decrease in strength and an increase in frequency .
it is not clear , however , whether the lfn on the fb is the same as that on the hb and nb .
unfortunately no hbo was observed above @xmath3=2.1 to compare its behavior with .
van der klis et al .
( 1997 ) suggested that in sco x-1 the nbo / fbo emerged from the lfn .
this seems not to be the case in gx 17 + 2 , since the rms amplitude of the lfn already started to decrease well before the nbo appeared .
moreover , @xmath25 shows a gradual increase between @xmath3=-0.6 and @xmath3=1.4 to values well above the nbo frequency .
a comparative study by @xcite strongly suggests that lfn has a similar origin as the band - limited noise component in atoll sources , bhcs , and the millisecond x - ray pulsar sax 1808.4 - 3658 .
they suggest that this noise component is produced in the accretion disk , and does not require the presence of magnetosphere or solid surface .      in our observations
gx 17 + 2 traced out the well known z - like tracks in the cd and hid .
in addition to motion along this track , we also observed secular motion of the z track itself in the hid of epoch 3 ( see figure [ cd - hid_fig ] ) , as was already reported by @xcite ( see also * ? ? ? * for a possible shift in the exosat hid of gx 17 + 2 ) .
below @xmath3=1 the fluxes in the hard and soft energy bands behaved quite differently ( fig.[sz - count_fig ] ) ; unlike the soft flux , the hard flux continued the trend that is observed above @xmath3=1 .
hence , the hb / nb vertex in the cd is entirely due to the turn over of the soft spectral flux .
it is not clear what causes the soft spectral flux to decrease below @xmath3=1 .
although it is partly an instrumental effect
the peak of the soft component slowly shifts out of the pca effective energy range as it becomes softer most likely intrinsic source changes occur at @xmath3=1 .
most variability components increase in strength below @xmath3=1 , suggesting that they are related to the hard spectral component .
it is commonly believed that the motion of the source along and its position on the z track are determined by the mass accretion rate ( @xmath33 ) through the inner accretion disk and onto the neutron star , increasing from the hb to the fb .
even though the count rates actually decrease on the nb and in some sources also on the fb ( implying an anti - correlation between @xmath33 and flux ) there are several observational results that support this view : ( 1 ) the frequencies of the hbo and khz qpos , that in many models are directly related to @xmath33 , gradually increase from hb to nb ( and in sco x-1 even to the fb , see * ? ? ?
( 2 ) the optical and uv flux increase from hb to fb @xcite ; in x - ray binaries they are thought to be due to reprocessing of x - rays from the central source in the outer accretion disk .
there are several reasons to believe that the optical and uv flux are betters tracers of the mass accretion rate than the x - ray flux @xcite .
this would then solve the apparent anti - correlation between @xmath33 and the ( x - ray ) flux @xcite .
( 3 ) motion of a source is always along the z track ; jumps between branches are never observed , in accordance with the assumption that @xmath33 varies continuously .
it has been found that secular motion of the z tracks in the cds of the cyg - like sources does not affect the relation between the rapid variability properties and source position along the track ( see , e.g. @xcite , however see @xcite ) . in view of the @xmath33-driven picture
mentioned above , these changes can therefore not be due to changes in @xmath33 , but should have a different origin .
it has been suggested that the cyg - like sources are viewed at a higher inclination than the sco - like sources @xcite .
the differences in shapes of the z tracks and some of the variability properties are also explained in this model @xcite .
@xcite suggest that the difference among the z sources might be explained by a difference in the magnetic field strength , with the sco - like sources having a lower magnetic field .    in the last few years results
have been obtained that may challenge the above described @xmath33-driven picture .
these results are mainly due to the arrival of _ rxte _  and _ bepposax _ , whose high quality data revealed new phenomena and allowed for a better comparison between different classes of x - ray binaries . in the following we discuss how some of our results of gx 17 + 2 conflict with the @xmath33-driven picture .
the frequency of the hbo in gx17 + 2 starts to decrease when the source is half way through the nb .
although there are several explanations for this ( see section [ disc - hbo_sec ] ) , it means that the first assumption given above ( i.e.  qpo frequencies increase with @xmath33 ) is not as solid as it seemed a few years ago .
the theory of type i x - ray bursts predicts that their properties change and correlate with @xmath33   @xcite . a study of the type i x - ray bursts in gx 17 + 2 by @xcite showed that their properties did not correlate well with the position along the z track at which they occurred , and therefore not well with the inferred @xmath33 .
although other explanations are possible ( see * ? ? ?
* ) , this could mean that @xmath33  does not determine the position of the source along the z track .
however , note that there does seem to be a correlation between burst occurrence and position along the z track , which is usually explained by a changes in @xmath33  onto the neutron star ( along the z track ) .
recently , @xcite reported the discovery with bepposax of a hard tail in the energy spectrum of gx 17 + 2 .
it increased in strength from the hb / nb vertex ( where only upper limits on its strength could be determined ) onto the hb .
they also showed that the total 0.1200 kev flux increased monotonically from the nb / fb vertex ( @xmath38 ) to the top of the hb ( @xmath39 ) .
this is different from the behavior of the 2.920.1 kev count rate ( see figure [ sz - count_fig]a ) , which increased from nb / fb vertex to the hb / nb vertex , but decreases afterwards .
as can be seen from figure [ sz - count_fig]b the behavior of the hard count rate confirms the findings of @xcite .
also , there is no evidence for a change in its behavior near the hb / nb vertex .
unfortunately no bepposax data were taken on the fb to compare to our rxte data .
hard power law tails have also been found in the z sources gx 5 - 1 @xcite , cyg x-2 @xcite , sco x-1 @xcite and gx 349 + 2 @xcite . except for sco x-1 , where no clear correlation with spectral state was found @xcite , in all these sources
the strength of the hard power law tail seems to decrease in the same direction as was observed in gx 17 + 2 , i.e. from hb to fb .
the behavior of the 0.1200 kev flux on the nb and hb suggest that @xmath33  actually increases in the opposite direction from what is commonly assumed .
however , since the flux changes are rather moderate ( @xmath020% ) , other options are that @xmath33  does not change at all along the z track , or at least that it shows no good correlation with the position along the z. in that case motion along the z track ( i.e. spectral change ) is caused by another parameter , for example the inner disk radius @xcite .
a similar possibility was suggested by @xcite to explain the observed behavior of the black hole transient xte j1550564 .
they found that spectral transitions occurred at different levels of the inferred @xmath33 .
the spectral changes were ascribed to an unknown parameter , which they suggested to be the inner disk radius ( based on variability ) , the size of the comptonizing region ( based on spectral hardness ) , and/or the accretion flow geometry ( based on radio brightness ) .
note that options like the latter one may lead to changes in the @xmath33  onto the neutron star , while the @xmath33  through the inner disk remains constant ; in that way the occurrence of type x - rays bursts might become dependent on the position along the z track . if the above also applies to the z sources , the secular motion that is observed could perhaps be explained by changes in the @xmath33 .
the apparent lack of a correlation between the type i burst properties and the position along the z track could then also be explained , since @xmath33  would determine the properties of the bursts but no longer the spectral properties .
it is interesting to see that gx 17 + 2 and the z sources , show in fact several similarities with xte j1550564 and other black hole sources .
the first was already briefly mentioned above , and concerns the independent behavior of the hard and soft flux on the hb .
this independence of the hard and soft spectral components is commonly observed in black hole candidates ( bhcs , see * ? ? ?
moreover , there are several properties that correlate well with the spectral hardness in bhcs and also with the spectral hardness in z sources , which increases from the nb to the hb ( for the moment we will exclude the fb from the comparison ; it will be discussed later ) .
these properties are the radio flux , the optical flux , and to a certain degree also qpo frequencies .
( 1 ) the radio flux in bhcs is correlated with the spectral hardness ( * ? ? ?
* ; * ? ? ?
* and references therein ) .
this is also true for gx 17 + 2 , where the radio flux increases from the nb to the hb @xcite .
sco x-1 @xcite and cyg x-2 @xcite show similar behavior , whereas gx 5 - 1 only showed one radio flare on the nb @xcite .
the behavior of the radio flux suggests that as the spectral hardness increases , part of the flow onto the neutron star evolves into an outflow .
( 2 ) @xcite found that in xte j1550564 the optical flux decreased during several small hard x - ray flares .
although the optical counterpart of gx 17 + 2 @xcite has never been studied in much detail , this behavior is similar to the anti - correlation between the hard color and optical / uv flux on the nb and hb in the z sources cyg x-2 @xcite and sco x-1 @xcite .
( 3 ) in xte j1550564 the frequencies of the low and high frequency qpo anti - correlate with the hard flux @xcite .
this is also found in gx 17 + 2 for the khz qpos and hbo ( compare figure [ sz - count_fig]b with figures [ hbo_fig ] and [ khz_fig ] ) and other z sources .
the only exception to this , together with the odd behavior of the hard tail in sco x-1 ( see above ) is the correlation between the hbo frequency and spectral hardness that we observed on the lower nb . in that respect
it is interesting to mention that @xcite report a switch from a correlation to an anti - correlation between the qpo frequency and spectral hardness in several bhcs , when the latter increases beyond a certain value .
in addition to the qpo frequencies also the band - limited noise behaves similarly .
it becomes stronger with spectral hardness , both in bhcs and z sources @xcite    it is tempting to compare the different black hole states @xcite with those of the z sources .
based on the above similarities , we suggest that the hb corresponds to the low state , the nb to the intermediate and/or very high state state , and the fb to the high state . however , there are still considerable differences : the low state spectra are much harder than those of the hb , and also the noise in the low state is much stronger . moreover , flaring is not observed in the black hole high state . however , in both the high state and fb the variability at high frequencies ( @xmath750 hz ) is extremely weak .
it is not clear if and how the presence of a solid surface can account for these differences .
our most important findings are the variable frequency difference of the khz qpo , the non - monotonic relation between the hbo and the upper khz qpo , the similar q values of the hbo , its second harmonic , and the khz qpos and the suppression of the nbo amplitude during type i x - ray bursts .    in the following
we summarize what the implications of each of our new findings are for some of the proposed models .
the models we considered for the hbo and khz qpos were the relativistic precession model ( rpm ) , the disk - filter model , the sonic point beat frequency mode ( spbfm ) , and the two - oscillator model .
the latter model can be discarded on the basis of our results ( see * ? ? ?
* ) . for the disk - filter model we assumed that it produces the same frequencies as the rpm and it will therefore not be discussed it separately .
\(1 ) hbo - khz qpo frequencies : the rpm can only explain a turnover of the hbo - khz qpo frequency relation when one assumes that the observed hbo is four times the predicted frequency and then only for rather extreme neutron star parameters : a high neutron star spin frequency ( @xmath0500700 hz ) , a neutron star mass of more than 2 @xmath32 , and a hard equation of state .
the spbfm needs some extension ; for example , decoupling of part of the accretion flow from the disk flow outside the magnetosphere ( to explain the hbo frequency decrease ) and a subsequent recoupling of this flow with the disk flow before the sonic point radius is reached ( to explain the increase in the khz qpo frequencies ) .
how this should exactly work is not clear .
\(2 ) variable khz qpo separation : although not significant at a 3 sigma level , the khz qpo separation seems to decrease both towards low and high values for the upper khz qpo frequency .
the rpm basically predicts such behavior , but in its current form it is inconsistent with our results , as is apparent from figure [ khz - diff_fig ] .
the spbfm can not explain the possible decrease of the frequency separation at the low frequency side .
\(3 ) similar q values : in the spbfm the low and high frequency qpos are produced independently from each other and there is no a priori reason why their q values should be similar . in the rpm ,
all the oscillations are produced by a single test particle .
lifetime broadening and the presence of multiple test particles both lead to different q values for the low and high frequency qpo .
frequency modulation could work for both models , but only if the modulating parameter is not the same as the one that causes the motion along the z track .
the observed behavior of the nbo during two type i x - ray bursts showed that this qpo is rather sensitive to the strength of the radiation field from the neutron star .
this seems to confirm the basic idea behind the feedback mechanism model that was proposed by @xcite , i.e. it requires a delicate balance between the inwards directed forces and the radiation field .
finally , we compared the behavior of gx 17 + 2 ( and the other z sources ) with that of ( some ) black hole lmxbs .
many similarities exist , particularly with regard to the relation between the hard spectral component on the one hand , and variability properties ( qpo frequency , broad band noise strength ) and radio flux on the other .
based on recent findings in black hole lmxbs we conclude that the mass accretion might not be the parameter that determines the position along the z track after all .
this work was supported in part by the netherlands organization for scientific research ( nwo ) grant 614 - 51 - 002 and by nwo spinoza grant 08 - 0 to e.p.j .
van den heuvel .
rw was supported by nasa through the chandra postdoctoral fellowship grant number pf9 - 10010 awarded by the chandra x - ray center , which is operated by the smithsonian astrophysical observatory for nasa under contract nas8 - 39073 .
mm is a fellow of the consejo nacional de investigaciones cientficas y tcnicas de la repblica argentina .
whgl is grateful for support from nasa .
the authors thank the referee for his / her thorough reading and usefull remarks .
finally , we would like to thank draza markovi for helpfull comments and for providing theoretical fits to our data .
, m. & barret , d. 2000 , in proceedings of the fifth compton symposium , held in portsmouth , nh , usa , september 1999 .
melville , ny : american institute of physics ( aip ) , 2000 . edited by mark
l. mcconnell and james m. ryan aip conference proceedings , vol . 510 .
, p. 157
, y. & lewin , w. h.  g. 1995 , in x - ray binaries ( cambridge astrophysics series , cambridge , ma : cambridge university press , |c1995 , edited by lewin , walter h.g . ; van paradijs , jan ; van den heuvel , edward p.j . ) , p. 126
, m. 1989 , in proceedings of the nato advanced study institute on timing neutron stars , held in eme , .
izmir , turkey , april 415 , 1988 .
editors , h. gelman and e.p.j .
van den heuvel ; publisher , kluwer academic , dordrecht , the netherlands , boston , massachusetts , p. 27
, m. 1995 , in proceedings of the nato advanced study institute on the lives of the neutron stars , held in kemer , turkey , august 29-september 12 , 1993 .
editor(s ) , m.a .
alpar , u. kiziloglu , j. van paradijs ; publisher , kluwer academic , dordrecht , the netherlands , boston , massachusetts , 1995 . , p. 301
( epoch 4 data only ) .
the hb / nb vertex in the 2.920.1 kev band ( a ) is caused by the contribution from the low energies ( b ) .
this vertex does not show up at high energies ( c ) .
[ sz - count_fig],width=377 ]      = 0.40.5 power spectrum , showing the relative contribution of the individual components .
the combined fit is represented by the thick solid line , the three lorentzians / qpos by the thin solid lines , the cut - off power law / lfn by the dashed line , and the power law / vlfn by the dotted line .
the poisson level was subtracted .
the reduced @xmath40 of this fit was 1.03 for 204 degrees of freedom .
[ pds_comp_fig],width=377 ]            .
the left column shows the results for the fits with a lorentzian , the right column those for fits with a cut - off power law . in both cases
the rms amplitude is the integrated power in the 1100 hz range .
@xmath25 is the frequency at which most of the power is concentrated ( see text for the expressions for @xmath25 . ) for reasons of comparison we also plotted the values of the other fit function ( gray ) in the panels for the rms amplitude and @xmath25.[lfn_fig],width=377 ]        23 . for comparison
we also plotted the data for sco x-1 ( gray ) from @xcite . for additional details
see text[khz - diff_fig],title="fig:",width=245 ] 23 . for comparison we also plotted the data for sco x-1 ( gray ) from @xcite . for additional details
see text[khz - diff_fig],title="fig:",width=245 ]      .
a log of all _ rxte_/pca observations used in this paper .
mode 7 was not always active during the november 1998 observations .
note that none of the observations represents an uninterrupted interval .
each is a collection of observations that were done around the same time .
these observations were separated in time from each other for various reasons , such as earth occultations , passages of the south atlantic anomaly , or observations of other sources.[obs_tab ] [ cols="<,<,^,^,^",options="header " , ]