A&A 471, 159-164 (2007)
DOI: 10.1051/0004-6361:20077347
M. Revnivtsev1,2 - S. Sazonov1,2
1 - Max-Planck-Institute für Astrophysik,
Karl-Schwarzschild-Str. 1, 85740 Garching bei München,
Germany
2 -
Space Research Institute, Russian Academy of Sciences,
Profsoyuznaya 84/32, 117997 Moscow, Russia
Received 23 February 2007 / Accepted 18 June 2007
Abstract
We analyzed deep Chandra observations of the Galactic
plane region centered at
,
with the
aim to obtain the best possible constraints on the contribution of weak
point sources to the Galactic ridge X-ray emission (GRXE) in this
region. We demonstrate that the vast majority of the detected sources
are Galactic in origin and are probably cataclysmic variables and
coronally active stars. We use the number-flux function of detected
sources to constrain the luminosity function of Galactic X-ray sources
in the range
1030-1032 erg s-1 and find good agreement with
the luminosity function of sources in the Solar vicinity. The fraction
of the total flux at energies 1-7 keV resolved into point sources at
the current sensitivity level is
25%. Excluding the expected
contribution of extragalactic sources,
19% of the GRXE is due
to point Galactic sources with interstellar absorption
corrected fluxes higher than
erg s-1 cm-2 in the energy band
1-7 keV.
Key words: stars: binaries: general - Galaxy: disk - X-rays: general - X-rays: stars
The origin of the Galactic ridge X-ray emission (GRXE) has been a long-standing problem of X-ray astronomy. Since the discovery of this emission in the late 1970s (Cooke et al. 1970; Bleach et al. 1972; Worrall et al. 1982), different explanations have been proposed for its nature. These can be divided into two major branches: i) the GRXE is truly diffuse, i.e. it is generated in the interstellar medium (see e.g. Ebisawa et al. 2005; Koyama et al. 1986,1989; Tanaka 2002), or ii) it is the integrated emission of faint Galactic point sources (Mukai & Shiokawa 1993; Koyama et al. 1986; Ottmann & Schmitt 1992; Worrall & Marshall 1983; Revnivtsev et al. 2006b).
The main problem with the former hypothesis is the difficulty of keeping plasma that is apparently thermal and hot (>5-10 keV) (Koyama et al. 2006,1986,1989; Tanaka 2002) within a thin layer along the Galactic plane (see e.g. Warwick et al. 1985; Revnivtsev et al. 2006b; Worrall et al. 1982; Yamauchi & Koyama 1993).
The latter scenario recently received strong support from work of our team (Revnivtsev et al. 2006a; Sazonov et al. 2006; Revnivtsev et al. 2006b; Krivonos et al. 2006). In particular, we obtained high-quality maps of the GRXE that showed striking similarity with the near-infrared map representing the distribution of stellar mass over the Milky Way. We concluded that throughout the Galaxy the GRXE volume emissivity is proportional to the stellar mass density. Furthemore, we showed that the GRXE emissivity per unit stellar mass is consistent with the cumulative emissivity of known classes of weak X-ray sources in the Solar vicinity - cataclysmic variables and coronally active stars (Sazonov et al. 2006).
We recently used deep Chandra observations of a field in the close vicinity of the Galactic Center (GC) to demonstrate that at least 40-50% of the total X-ray emission from that region is created by point sources (Revnivtsev et al. 2006c). It is known that the properties of the X-ray emission from the GC region are very similar to those of the large-scale GRXE, suggesting that both have the same physical origin (e.g. Tanaka 2002). However, as the space density of stars near the GC is 3-4 orders of magnitude higher than is typical for the Galactic plane (GP), one could expect that truly diffuse X-ray emission, should such exist, might provide a more significant contribution to the GRXE from the GP.
This motivated us to study one of the deepest Chandra
observations of the GP, centered at
,
,
with the aim to estimate the contribution of resolved
Galactic point sources to the total X-ray flux from this
direction. This observation was previously analyzed by
Ebisawa et al. (2005).
We used the Chandra observation (OBSID 2298) of the GP region in
the direction
,
with a
total exposure time of
102 ks. In order to maximize the
sensitivity to point sources, we restricted our analysis to the data
collected within
radius around the aim point of the
telescope, where its angular resolution is better than 1'' (FWHM). At larger distances from the aimpoint the
size of the point spread function of Chadnra significantly
increases (
2.6 times at 4
distance from the aimpoint)
and the sensitivity drops by more then 1.5 times (e.g. Muno et al. 2006), therefore we do not consider
these regions below.
The solid angle of the field of our study is thus
ster. Note that we did not analyze existing
Chandra observations of a nearby GP region (OBSIDs 949 and 1523),
because the central
-area of this field does not overlap
with the region of our study, making it impossible to increase the
point source sensitivity, which is key to us, by adding these extra data.
The data were reduced following a standard procedure fully described in
Vikhlinin et al. (2005). The only difference is that the detector
background was now modeled using the stowed dataset (http://cxc.harvard.edu/contrib/maxim/stowed). Point sources were
detected using the wavelet decomposition package
of
software
(Vikhlinin et al. 1998)
.
Due to the significant drop in the Chandra sensitivity at energies
E>2 keV and especially at E>7 keV, and because of the strong
effect of interstellar photoabsorption in the GP at E<1 keV, we used the
energy band 1-7 keV for our analysis. A total of 6.5 kcnts
in this energy band were accumulated from the studied area of 0.01396 sq. deg over 102 ks. This implies the total observed X-ray
flux
erg s-1 cm-2. For conversion of the count rate into the flux we used
the simplest best fit model to the observed data points - a power law
in energy band 1-7 keV with the photon index
After constructing the number-flux function of deteced sources, we corrected it for the incompleteness near the detection threshold using the procedure developed for our GC study (Revnivtsev et al. 2006c). In essence, Poisson fluctuations in the number of counts detected from a weak source effectively reduce the probability of its detection in a given observation (see e.g. Hasinger et al. 1993; Kenter & Murray 2003). To correct for this "leakage'' of sources under the detection threshold, we derived the allowed range of intrinsic log N-log S functions of sources by simulating mock Chandra images for a large number of such functions and checking which of these images are consistent with the measured number-flux relation. The multifold of allowed intrinsic log N-log Sfunctions was then used to determine the resolved fraction of the GRXE.
The studied Chandra field is characterized by heavy interstellar
absorption. The total column density of atomic and molecular gas through the
Galactic disk in this direction is estimated at
H atoms per cm2 (Ebisawa et al. 2005). The Galactic sources in the field
are expected to be located at widely differing distances, from our
close vicinity all the way to the outer boundary of the Milky Way, whereas the
distribution of the absorbing gas within the sampled cone of the
Galaxy is poorly known. Ebisawa et al. (2005) made an attempt to estimate
typical column densities from stacked Chandra spectra of subsets
of sources detected in this field and found values ranging from
to >
cm-2. In view of the
remaining uncertainty, we assume in our analysis that all
Galactic X-ray sources in the field are characterized by the same
absorbing column density
cm-2, while
all extragalactic ones by
cm-2.
To make a direct comparison of the source counts in the GP region with those in typical extragalactic fields (maximally devoid of
Galactic sources), we also analyzed Chandra observations of the Deep
Fields South (CDFS, OBSID 5021) and North (CDFN, OBSID 3389) in
exactly the same manner as the GP data. To have the same exposure time
as for the GP field, we used only 100 ks of the available
(ultra-deep) observations for each of these extragalactic
fields. However, even in this case the number of detected
extragalactic sources in the GP region is expected to be significatly
reduced compared to the CDFS and CDFN by the strong interstellar absorption.
The majority of weak X-ray sources (
erg s-1) in our Galaxy
are cataclysmic variables and coronally active binary stars
(e.g. Sazonov et al. 2006; Vaiana 1981). These classes of sources are expected to be
distributed over the Milky Way just like ordinary stars (with the
notable exception of globular clusters where the relative fraction of X-ray
sources may be different due to the greatly increased role of
dynamical processes in forming such systems). Such proportionality
between the space density of weak X-ray sources and that of stars has
indeed been observed (see e.g. Revnivtsev et al. 2006c; Muno et al. 2006). This allows
us to predict the distribution of X-ray sources along our
studied light of sight and consequently the expected log N-log Sfunction for a given source luminosity function, e.g. for
that measured in the Solar vicinity (Sazonov et al. 2006).
In estimating the contribution of extragalactic sources to the
number-flux function measured in the GP field, we preferred to carry
out a direct comparison with 100 ks subsets of CDFN and CDFS
observations, rather than relying on published log N-log Sfunctions of extaragalactic sources (e.g. Moretti et al. 2002). This
allows us to use the same energy band (1-7 keV) for detecting sources, thus
avoiding any rescaling and its associated uncertainties. Also,
number-flux relations constructed from these observations should be
affected by the leakage of weak sources (see above) similarly to
the GP observation, faciliating comparison of the measured log N-log S curves.
In order to make a prediction for the number-flux function of Galactic
X-ray sources in the studied GP region, we should take into account
the spatial distribution of stars (which have been shown to be good tracers of
weak X-ray sources) along the line of sight. For this purpose we
should adopt some mass model of the Galaxy.
Since our observation is directed relatively far from the GC (
), we do not expect any significant contribution
from sources residing in the Galactic bulge. Thus, it is sufficient to
consider only the disk component of the Galaxy.
There is general agreement that the space density of stars in the
Galactic disk is exponentially declining (with scale-height
)
with Galactocentric radius R and that there is a similar
decline of density with radius within some minimal radius
due to the presence of the Galactic bulge/bar inside this radius
(e.g. Freudenreich 1998). We therefore adopt below the following
mass model of the Galactic disk:
![]() |
Figure 1:
Profile of the near-infrared brightness of the Galaxy along
its plane as seen by COBE/DIRBE at ![]() |
Open with DEXTER |
Parameter of the adopted model are slightly different from those used in
Sazonov et al. (2006) and Revnivtsev et al. (2006b) due to special requirement to predict
correct near-infrared (NIR) surface brightness value at the region of
our study (
). Good correspondence of the NIR surface
brightness profile to that observed by COBE/DIRBE (corrected for the
interstellar extinction following
Revnivtsev et al. 2006a) is seen in Fig. 1.
The extinction corrected surface brightness of the
Galactic disk at 3.5 m in the direction of our observation
(
,
)
is
MJy/sr (a significant fraction of the quoted uncertainty is associated
with the correction of the COBE/DIRBE data for interstellar
extinction), or
erg s-1 cm-2 deg-2.
Given the area (
deg-2) of the studied Chandra field, the total NIR flux from this region
erg s-1 cm-2.
From our study of the large-scale distribution of the GRXE with RXTE
we found the ratio of the NIR surface brightness to the GRXE intensity
in the 3-20 keV band to be
(Revnivtsev et al. 2006b). Assuming that the GRXE spectrum is a power-law with
photon index
,
we can rescale this factor to the 1-7 keV
band:
.
Multiplying this ratio by the measured NIR flux yields the
expected GRXE intensity in the considered direction:
erg s-1 cm-2 deg-2. Consequently, the GRXE
flux from the entire Chandra field is expected to be
erg s-1 cm-2. This estimate does not
take into account the significant interstellar absorption in the studied
direction. Considering column densities
cm-2 to be typical of Galactic X-ray sources in this
field (see the discussion of this issue in Sect. 2
above), the expected absorbed GRXE flux
erg s-1 cm-2.
The above estimate still does not take into account the
contribution of extragalactic sources. The all-sky average intensity
of the cosmic X-ray background (CXB) is
erg s-1 cm-2 deg-2 (e.g. Revnivtsev et al. 2005,2003; Hickox & Markevitch 2006),
or
erg s-1 cm-2 deg-2 (assuming
a power-law spectrum with
). This translates into a
flux
erg s-1 cm-2 for our Chandra field. Taking into account that the full Galactic column
density in the direction of our study is
cm-2, the CXB contribution to the observed X-ray flux is
expected to be
erg s-1 cm-2. Therefore, the
total (GRXE plus CXB) X-ray flux from the studied region should be
erg s-1 cm-2.
This is compatible with the flux actually measured by Chandra:
erg s-1 cm-2. We should stress here that correcting for the interstellar
absorption is very important for the 1-7 keV energy band but is not
straightforward and subject to significant uncertainties due to
possible variations in the absorption along the line of sight and
across the studied field.
The consistency of cumulative NIR and X-ray characteristics of the
studied Chandra field provides additional support to our adopted Galaxy
model, which is used in the subsequent analysis. We emphasize that it
is crucial that the adopted mass model correctly reproduces the NIR
flux from the considered direction of the GP.
![]() |
Figure 2:
Comparison of the raw Chandra images (1-7 keV) of
three sky regions. From left to right - Chandra Deep Field
North (CDFN), Chandra Deep Field South (CDFS), and Galactic
plane region at
![]() ![]() |
Open with DEXTER |
Despite the long history of studying weak (
erg s-1) Galactic X-ray sources, their luminosity function remained
poorly constrained until lately. We recently used two all-sky X-ray
surveys (the RXTE slew survey in the 3-20 keV band and
the ROSAT all-sky survey below 2 keV) to construct a luminosity
function of X-ray sources located near the Sun in the
broad range
1027-1034 erg s-1 (Sazonov et al. 2006). Below we use
this luminosity function in conjunction with the Galaxy model
described above for predicting the number-flux function of sources in the
studied GP region.
If there were no significant population of weak Galactic X-ray
sources, one could expect the surface density of point sources
to be similar in different regions of the sky observed by Chandra to the same depth. However, if we compare the Chandra
image of the GP region with those of the CDFN and CDFS obtained with
the same exposure (Fig. 2), we clearly see that
the surface density of point sources is much higher in the GP region.
This difference becomes even more pronounced if
we take into account that due to the considerable
interstellar absorption in the direction of
,
,
the number of detected extragalactic sources in that field
is expected to be significantly reduced compared to the CDFN and
CDFS. For extragalactic sources with fluxes
10-15-10-14 erg s-1 cm-2
in the 1-7 keV energy band, whose spectra can typically be described as power
laws with photon indexes
(e.g. Hickox & Markevitch 2006), the Galactic absorption with
cm-2 is expected to suppress the Chandra
count rates in the 1-7 keV band by a large factor of
4.3-5.3 (note
that since absorption mostly removes photons with energies
1-2 keV, the corresponding decrease in the energy flux at 1-7 keV is
much smaller). Therefore, for a given exposure time, Chandra
should detect in the CDFN and CDFS roughly 5 times as many
extragalactic sources as in the GP field.
The much higher surface density of sources in the GP region is further reflected in the log N-log S function of detected sources (see Fig. 5 below). We conclude that in the Galactic plane weak point X-ray sources of Galactic origin greatly outnumber extragalactic ones.
![]() |
Figure 3:
Differential number-flux function of sources measured in the
GP region and corrected for the incompleteness at low count rates
(points with error bars). The hatched region represents the manifold
of number-flux functions allowed by the data (see main text and the
corresponding range of allowed luminosity functions in
Fig. 4). The gray band shows the expected range of
number-flux functions based on the luminosity function of weak X-ray
sources in the Solar vicinity from Sazonov et al. (2006) and including the
predicted contribution of extragalactic sources. The width of this
region is determined by the ![]() ![]() ![]() ![]() ![]() ![]() ![]() ![]() ![]() |
Open with DEXTER |
We now address the number-flux function of detected sources. In Fig. 3 we show the differential log N-log Sdistribution measured by Chandra in the GP region and corrected for the incompleteness near the detection threshold (see Sect. 2). These data are compared with the results of our modelling described below.
To make a prediction for the number-flux function, we convolved the
luminosity function of weak X-ray sources measured in the Solar
vicinity (Sazonov et al. 2006) with the Galactic disk model described in
Sect. 3.1. To convert the 2-10 keV luminosity function
of Sazonov et al. (2006) into the 1-7 keV band, we assumed a power-law
spectrum with
.
In addition, to take into account the interstellar
absorption in the considered direction, we multipled all predicted
1-7 keV count rates by a factor of 0.34, which corresponds to
a power-law spectrum with
absorbed by neutral gas with
column density
cm-2. As we discussed
in Sect. 2, this is only a rough (but reasonable)
approximation of the true situation, as in
reality there is a large scatter in
values from source to
source.
As regards the expected contribution of extragalactic sources, we
estimated it from the log N-log S curves measured in the CDFN and
CDFS by applying an absorption correction as appropriate for our GP
field (
cm-2). Adding this expected
extragalatic contribution (shown by the dashed line in
Fig. 3) to the expected contribution of Galactic
sources yields the expected range of log N-log S functions, shown
in gray in Fig. 3. There is apparently a good
agreement between the expected and observed distributions. The
apparent small difference between them may well be caused by some
inaccuracy in our correction for the interstellar absorption.
On the other hand, we can use the observed number-flux function to
put constraints on the luminosity function of point sources
in the studied GP region. To this end, we assumed that the luminosity
function in a
representation can be described by a set of
constants in the following intervals:
1028-1029,
1029-1030,
1030-1031,
1031-1032, and
1032-1033 erg s-1. By allowing these constants to vary and
convolving the resulting luminosity functions with our
Galactic disk model, we obtained a manifold of trial number-flux
functions of sources. These were then compared in terms of the reduced
value (
per degree of freedom) with the observed log N-log S function (taking into account the estimated contribution of
extragalactic sources). Those trial luminosity functions that resulted
in a reduced
of less than 1.5 were regarded as allowed by the
data. The resulting allowed range of number-flux functions (including
extragalactic sources) is shown as the dashed area in
Fig. 3, and the corresponding range of acceptable luminosity
functions of Galactic sources is shown as the dashed area in
Fig. 4, where it is compared with the luminosity function of local
sources taken from Sazonov et al. (2006).
![]() |
Figure 4: Range of luminosity functions of Galactic X-ray sources in the GP region (hatched area) allowed by the measured number-flux function (see Fig. 3). This is compared with the luminosity function (gray area) of weak X-ray sources in the Solar neighborhood (Sazonov et al. 2006). For the latter we included the uncertainty in the fraction of young single stars in the studied GP region (relative to cataclysmic variables and coronally active binary stars, see Sazonov et al. 2006 for details). Specifically, we allowed this contribution to vary from zero to the value measured near the Sun. Note that we effectively probe only the luminosity range 1030-1032 erg s-1 (as indicated by the wavy line). |
Open with DEXTER |
As follows from Fig. 4, the allowed range of luminosity
functions in the GP region is fully compatible with the luminosity
function of point sources in the Solar vicinity. It is important to
note here that the 100 ks Chandra observation of the
GP region allows us to efficiently constrain the luminosity function
only in the range
1030-1032 erg s-1. The upper bound appears
due to the limited stellar mass within the small solid angle covered
by Chandra (<
). Sources brighter
than
1032 erg s-1 are simply too rare to be found in such a
small area of the Galactic disk. The lower bound is due to finite
Chandra's sensitivity.
In Fig. 5 we compare the cumulative log N-log S distributions constructed for the GP field and for the CDFS and CDFN (all corrected for the incompleteness at low count rates). One can clearly see how Galactic sources become progressively more dominant toward lower fluxes.
![]() |
Figure 5: Cumulative number-flux functions of point sources detected in the GP region (upper histogram) and in extragalactic fields ( Chandra Deep Fields North and South). All data (GP, CDFN and CDFS) are corrected for the incompleteness at low count rates. The dashed line shows an approximation of the extragalactic source counts by a simple analytic model similar to that of Moretti et al. (2002). |
Open with DEXTER |
Fraction of the energy flux resolved into point sources in a relatively broad energy range that we consider (1-7 keV) depends on the spectral shapes both of the Galactic ridge and the cumulative emission of point sources. Due to very strong and variable (for example due to different distances to sources) interstellar absorption of the X-rays in the region of our study it is very hard to predict the cumulative spectrum of point sources independently. Therefore in the following estimations of the fraction of the emission resolved into point sources we can refer either simply to Chandra count rates in the energy band 1-7 keV, or to count rate converted into energy fluxes using some fiducial shape of sources cumulative spectrum.
Integrating the number-flux function measured in the GP region
(Fig. 5) down to our effective detection limit of 5 counts
(over the 102 ks exposure, which corresponds to an absorption corrected
flux of
erg s-1 cm-2 in the 1-7 keV band for Galactic
sources, assuming a power low spectral shape with
and
cm-2) yields
1.6 kcnts
out of the total
6.5 kcnts collected by Chandra. This
implies that
25% of the total flux at energies 1-7 keV from the
GP region is already resolved by Chandra into point
sources. If we similarly integrate the expected log N-log S function of
extragalactic sources (aslo shown in Fig. 5) down to the
same count rate limit (which corresponds to a somewhat higher absorption
corrected flux of
erg s-1 cm-2 for
cm-2), we find that
6% of the total flux
is resolved into extragalactic sources.
This implies that at least 19% of the GRXE is due to Galactic
point sources, presumably cataclysmic variables and coronally active
stars. Moreover, an extrapolation of the measured log N-log S curve to fluxes below the current detection limit based on the
luminosity function of local weak X-ray sources (Sazonov et al. 2006) is
consistent with all of the X-ray flux from the GP being due to point
sources, mostly of Galactic origin. These conclusions are consistent
with our results for the Galactic Center region (Revnivtsev et al. 2006c).
We should say in conclusion a few words about the apparent difference of our results from those of Ebisawa et al. (2005). We can mention two main differences in the approaches taken. First, Ebisawa et al. (2005) assumed some fiducial luminosity functions of weak Galactic X-ray sources and some fiducial mass model of the Galaxy while predicting the number-flux function of possible galactic sources. In our paper we use the luminosity function of weak sources measured in the Solar neighborhood (Sazonov et al. 2006), which is quite different from the fiducial functions of Ebisawa et al. (2005), and demonstrate that this luminosity function is fully consistent with the number-flux function of sources masured in the GP region. Also the mass model of the Galaxy was specially checked to provide correct NIR flux in the direction of the study ensuring the more or less correct overlooked mass of stars. Second, in order to increase the sensitivity to weak point sources, we used an optimized broad energy band, which includes softer energies (1-2 keV), where Chandra has the maximum sensitivity to the galactic point sources, used only the central region of the Chandra field additionally maximizing our sensitivity to point sources and also we took into account the effects of the Possion noise on the observed number-flux function. All together these efforts allowed us to move deeper in sensitivity and to make more accurate comparison with the point sources hypothesys of the formation of the ridge emission.
Acknowledgements
This research made use of data obtained from the High Energy Astrophysics Science Archive Research Center Online Service, provided by the NASA/Goddard Space Flight Center. We acknowledge the use of the Legacy Archive for Microwave Background Data Analysis (LAMBDA). Support for LAMBDA is provided by the NASA Office of Space Science. MR aknowledge partial support of grant RFFI 07-02-00961.