XMM-Newton observation of 4U 1820-30:

Xmm-Newton observation of 4U 1820-30:

Broad band spectrum and the contribution of the cold interstellar medium
E. Costantini SRON, Netherlands Institute for Space Research, Sorbonnelaan, 2, 3584 CA, Utrecht, The Netherlands
   C. Pinto SRON, Netherlands Institute for Space Research, Sorbonnelaan, 2, 3584 CA, Utrecht, The Netherlands
   J.S. Kaastra SRON, Netherlands Institute for Space Research, Sorbonnelaan, 2, 3584 CA, Utrecht, The Netherlands
e.costantini@sron.nl Astronomical Institute, Utrecht University, PO Box 80000, 3508 TA Utrecht, The Netherlands
   J.J.M. in’t Zand SRON, Netherlands Institute for Space Research, Sorbonnelaan, 2, 3584 CA, Utrecht, The Netherlands
   M.J. Freyberg Max-Planck-Institut für extraterrestrische Physik, Giessenbachstr. 1, D-85748 Garching bei München, Germany    L. Kuiper SRON, Netherlands Institute for Space Research, Sorbonnelaan, 2, 3584 CA, Utrecht, The Netherlands
   M. Méndez Kapteyn Astronomical Institute, University of Groningen, Postbus 800, 9700 AV, Groningen, The Netherlands    C.P.  de Vries SRON, Netherlands Institute for Space Research, Sorbonnelaan, 2, 3584 CA, Utrecht, The Netherlands
   L.B.F.M. Waters SRON, Netherlands Institute for Space Research, Sorbonnelaan, 2, 3584 CA, Utrecht, The Netherlands
e.costantini@sron.nl Sterrenkundig Instituut Anton Pannekoek, University of Amsterdam, Science Park 904, P.O. Box 94249, 1090 GE Amsterdam, The Netherlands
Key Words.:

Context:We present the analysis of the bright X-ray binary 4U 1820-30, based mainly on XMM-Newton-RGS data, but using complementary data from XMM-Newton-EPIC-pn, INTEGRAL, and Chandra-HETG, to investigate different aspects of the source. The broad band continuum is well fitted by a classical combination of black body and Comptonized emission. The continuum shape and the high flux of the source () are consistent with a ”high state” of the source. We do not find significant evidence of iron emission at energies  keV. The soft X-ray spectrum contain a number of absorption features. Here we focus on the cold-mildly ionized gas. The neutral gas column density is  cm. The detailed study of the oxygen and iron edge reveals that those elements are depleted, defined here as the ratio between dust and the total ISM cold phase, by a factor and , respectively. Using the available dust models, the best fit points to a major contribution of Mg-rich silicates, with metallic iron inclusion. Although we find that a large fraction of Fe is in dust form, the fit shows that Fe-rich silicates are disfavored. The measured Mg:Fe ratio is . Interestingly, this modeling may point to a well studied dust constituent (GEMS), sometimes proposed as a silicate constituent in our Galaxy. Oxygen and iron are found to be slightly over- and under-abundant, respectively (1.23 and 0.85 times the solar value) along this line of sight. We also report the detection of two absorption lines, tentatively identified as part of an outflow of mildly ionized gas () at a velocity of  km s.




Conclusions:Astrochemistry – ISM: dust – X-rays: binaries – X-rays: individuals: 4U 1820-30

1 Introduction

The interstellar medium (ISM) in the plane of our Galaxy is a dynamic and complex environment, composed of mainly neutral matter in both gas and dust form and by a warmer gas phase in the form of diffuse emission in, and above, the Galactic plane. The properties of the cold phase in the diffuse ISM have been extensively studied at long wavelengths, from the far-infrared to the far-UV (e.g. Draine 2003, for a review). A sizeable fraction of the cold phase is locked up in dust grains (e.g. Savage & Sembach 1996; Jenkins 2009, and references therein). Amorphous silicate materials together with graphite and policyclic aromatic carbon should account for the majority of the depleted elements measured in the ISM: C, O, Fe, Mg, Si (e.g. Weingartner & Draine 2001; Wooden 2008). One of the major spectroscopic signatures of the presence of Fe- and Mg- rich amorphous silicates is the 10 m emission feature. Extensive studies of this feature lead to the conclusion that the Fe:Mg ratio should be approximatively 1 (e.g. Li & Draine 2001). Main sources of both Mg and Fe silicates are O-rich asymptotic giant branch stars (AGB) and supernovae. However, the process of amorphization of the dust agglomerates, for instance by rapid cooling of the gas phase (Wooden et al. 2005) or cosmic rays bombardment (Carrez et al. 2002), strongly favors the survival of Mg silicates. In fact, a recent analysis successfully models the 10 m feature in terms of Mg-rich silicates, if non-spherical grain shapes are used (Min et al. 2007). Glassy material, consisting of Mg-rich silicates with metallic iron and sulfide inclusion (called GEMS, Bradley 1994) have also been commonly found during the Stardust mission. Their origin is mostly from the interplanetary environment (e.g. Keller & Messenger 2004), but a fraction have a composition compatible with an ISM origin (Keller & Messenger 2008). Iron is however a highly depleted element (70–99% of Fe is in dust, Wilms et al. 2000; Whittet 2003) whose inclusion into solid grains is not completely understood (e.g. Whittet 2003). This is mainly due to the difficulty of modeling iron emission, which does not display any sharp feature in long-wavelength spectra.

The abundances of the most important metals in the Galactic disk smoothly decrease with the galactocentric distance. The average slope of the distribution is  dex kpc (Chen et al. 2003, and references therein). However, a large scatter in the abundance measurements as a function of the Galactic radius is reported. This is attributed to different factors which contaminate the smooth mixing due to the pure stellar evolution process. Indeed the medium can be locally influenced by e.g. supernovae ejecta, and in-falling metal-poor material into the disk (Lugaro et al. 1999; Nittler 2005).

In recent years it has become clear that the X-ray band could provide an excellent laboratory to study the silicate content of the diffuse ISM, as the absorption K edges of O ( keV), Mg ( keV) and Si ( keV) together with the Fe LII and LIII edges ( and  keV, respectively) fall in the low-energy X-ray band. The method used is to study the absorbed spectrum of bright X-ray binaries, located in different regions in the disk, observed with high-energy resolution instruments. This allows to probe the interstellar dust (ID) content in a variety of environments, with different extinction and with possibly different dust formation history. Previous studies of X-ray spectra taken along different lines of sight led first to the recognition that not only gas but also dust plays a role in shaping the iron and oxygen edges (Takei et al. 2002; Kaastra et al. 2009) and later led to the quantitative modeling of those edges (Lee et al. 2009; de Vries & Costantini 2009; Pinto et al. 2010). There is not yet a clear picture of the chemical composition of the ID as seen in X-rays. Silicates containing andratite (iron-rich silicates) were reported studying the oxygen edge of GS~1826-238 (Pinto et al. 2010), while iron oxides, rather than iron silicates were reported along the line of sight of Cyg~X-1 (Lee et al. 2009). This may point to a chemically inhomogeneous distribution of ID.

4U~1820-30 is an extensively studied source, by virtue of its extraordinary intrinsic properties. It is an ultracompact (orbital period 11.4 minutes, Stella et al. 1987) X-ray binary consisting of a neutron star and a He white-dwarf (Rappaport et al. 1987). The presence of X-ray bursts associated with the neutron star has been early recognized (first by Grindlay et al. 1976). 4U 1820-30 is classified as an atoll source (Hasinger & van der Klis 1989), showing kilohertz quasi-periodic oscillation at different frequencies in its power spectrum (e.g. Smale et al. 1997; Zhang et al. 1998).
The broad band spectrum has been studied with several instruments. The source displays a Comptonized continuum and a soft black body component (e.g. Sidoli et al. 2001).

With the advent of high resolution spectroscopy, the low-energy X-ray spectrum also revealed interesting features. Absorption by ionized ions has been reported in several studies (Futamoto et al. 2004; Yao & Wang 2005; Juett et al. 2006; Cackett et al. 2008a). The absence of blueshift in the lines of this gas and lack of variability generally points to an interstellar origin (but see Cackett et al. 2008a). The hypothesis of a gas intrinsic to the source is intriguing as other X-ray binaries often display absorption by ionized gas either at the rest frame of the source (van Peet et al. 2009) or outflowing (e.g. Miller et al. 2006; Neilsen & Lee 2009). The source has also been used as a backlight to illuminate the interstellar dust along the line of sight. This results in a dust scattering halo, which in this source is moderate (% at 1 keV, Predehl & Schmitt 1995), given the relatively low Galactic column density ( cm, this study).

Absorption by cold interstellar dust has been never studied in detail in this source. In this paper we aim for a comprehensive view of the cold absorbing medium along the line of sight to this source. This study benefits from the combined information provided by the Chandra and XMM-Newton high resolution instruments, which allows us to meaningfully study both the Fe L and O K edges. In addition, we make use of dust models updated with the latest laboratory measurements (e.g. Lee et al. 2008, 2009). We also present the broad band continuum behavior underlying the absorption features.

The adopted protosolar abundances follow the prescription given by Lodders & Palme (2009) and discussed in Lodders (2010). The broad band energy spectrum (EPIC-pn and INTEGRAL) is fitted using the minimization method, taking care that at least 20 counts per bin are present in each data set. For the high resolution spectra (RGS and Chandra-HETG) the Cash statistic has been used111http://heasarc.gsfc.nasa.gov/docs/xanadu/xspec/manual/XSappendixCash.html. The errors quoted are for 68% confidence level, corresponding to or . The spectral fitting package used in this paper is SPEX222www.sron.nl/spex (Kaastra et al. 1996). The adopted distance is 7.6 kpc (Kuulkers et al. 2003). The nominal hydrogen column density toward the source is  cm (Kalberla et al. 2005), to be compared with  cm (Dickey & Lockman 1990).

This paper is organized as follows: in Sect. 2, we illustrate the data handling for the different instruments used. Sect. 3 is devoted to the continuum determination, using EPIC-pn and INTEGRAL data. In Sect. 4 we describe in detail the modeling of the different absorption components using RGS and Chandra-MEG. The discussion can be found in Sect. 5 and the conclusions in Sect. 6.

2 The data handling

2.1 Epic-pn

4U 1820-30 was observed by XMM-Newton for 41 ks on April 2, 2009 (revolution 1706). The data reduction was performed using SAS (ver. 9.0). The EPIC-pn (Strüder et al. 2001) was operated in full-frame masked mode. In this mode, the central 13x13 (in RAW coordinates) pixels are flagged as bad on board, masking the central part of the PSF. The full-frame mode was chosen to optimally study the halo of diffuse emission surrounding the source, caused by the scattering of the intervening dust in the line of sight. However, the spectrum of the central source is recovered from the so called Out of Time (OoT) events, which are the photons trailed through the detector during the read out of the frame. The exposure time of the OoT events is 6.3% of the total exposure. We extracted the source spectrum from the RAW coordinate event file. We selected 11 columns from RAWX 32 to RAWX 42, avoiding the central 2 columns, which were affected by pile up. We selected RAWY89 to avoid the innermost region of the PSF, which is still affected by pile up or X-ray loading333http://xmm2.esac.esa.int/docs/documents/CAL-TN-0050-1-0.ps.gz.
The background was selected in a neighboring region. The exposure time of the background was scaled to match the effective exposure of the OoT events. Given the brightness of the source, the nominal background contribution is modest. However the whole area is contaminated by the photons of the scattering halo, which have a different distribution depending on the position in the detector. The OoT events were extracted from a modified event file processed in such a way that the charge transfer inefficiency (CTI) corrections are not applied and all the photons are considered as coming from the central source coordinates.

2.2 Rgs

Due to the brightness of the source, some RGS (den Herder et al. 2001) portions of the grating data were affected by pileup. In addition, in RGS2 the high number of events per CCD frame (twice the exposure time compared to RGS1, due to the single node readout) tended to fill up the buffers too fast. Limits on the data handling in relation to telemetry prevented the buffers to be cleared in time, leading to data losses for some CCD’s. A method to recognize the spectral regions not affected by pile up is to compare the first and the second order for each grating. Calibration shows that in absence of pile up, the ratio of the spectra of the two orders, in the region where they overlap, should be unity444http://xmm.vilspa.esa.es/docs/documents/CAL-TN-0075-1-0.pdf. In Fig. 1, we show the ratio for RGS1. We see that the maximum pile up is recorded between 14 and 19 Å, where the effective area of RGS is the largest. In this region the ratio is about 1.1, which translate in only a 6–7% pile up. In RGS2 the ratio between 2nd and 1st order reaches even a value of five in selected regions. At longer wavelengths ( Å) the combined effect of absorption and decrease of the effective area lowers significantly the count rate. Therefore, this region is not affected by pile up and can be safely used. Absorption features are however in general minimally influenced by this effect. We have chosen to keep the useful data given by RGS2 around the iron L edge, with care of locally fitting a different continuum for RGS1 and RGS2, using “sectors” in SPEX, as defined below (Sect. 4).
We also used archival RGS data (Table 1) in order to improve the signal-to-noise (S/N) ratio. This observation displayed a short period of flaring background, which was filtered out. This resulted in a cut of 3 ks on the total exposure time.

Figure 1: Upper panel: comparison between the 2nd and 1st order for RGS1. Lower panel: the ratio between the 2nd and the 1st order provides an indication of pile up in the spectrum.

2.3 Chandra-Hetg

In order to gain more signal to noise and energy resolution at wavelengths18Å, where the Fe, Mg and Si edges stand, we considered archival Chandra-HETG data, selecting the higher-flux data sets available (Table 1). We obtained the final products, processed in April 2010, from the TGCat archive555http://tgcat.mit.edu/. The three observations displayed similar fluxes and continuum parameters, therefore we combined the data, separately for the HEG and the MEG arms. HEG data are of low signal-to-noise at wavelengths larger than 17Å. In the following we refer to the MEG data only, unless otherwise specified.

2.4 Integral

In order to constrain the continuum spectrum of 4U 1820-30 across an as-large-as-possible energy band, INTEGRAL data can be very useful. Fortunately, the Galactic center region was observed by INTEGRAL for 73 ks from March 30, 2009 to March 31, 2009 during satellite revolution 789, and for 38 ks on April 3, 2009 during revolution 790. The XMM observation ended during the latter observation. Therefore we can consider these INTEGRAL observations as quasi-simultaneous with the XMM-observation. 4U 1820-30 was in the field of view of the INTEGRAL Soft Gamma-Ray Imager ISGRI (15-300 keV; Lebrun et al. 2003) at source angles smaller than for all 21 individual pointings (so called science windows) of the rev. 789 observation and all 10 pointings of rev-790. The source was outside the field of view ( diameter zero response) of the Joint European Monitor for X-rays JEM-X (3-35 keV; Lund et al. 2003) during the rev. 789 observation, but within the field of view during 4 individual pointings of rev. 790.

We generated mosaic maps for 10 logarithmically spaced energy bands across the 20–300 keV range for the ISGRI rev. 789 observation using OSA version 9.0 (distributed by the ISDC; Courvoisier et al. 2003) imaging tools. We checked these maps for significantly detected sources and subsequently derived spectral information for all these sources in 13 pre-defined energy windows across 13–520.9 keV.

For JEM-X (telescope 1) we followed similar imaging and spectral extraction procedures yielding spectra of 4U 1820-30 in 16 pre-defined energy bands across the 3.04–34.88 keV.The observation log is displayed in Table 1.

Inst. orbit/obsid date net exp. rate
(dd/mm/yy) (ks) (c/s)
XMM-pn 1706 02/04/09 39.8 20.9
INTEGRAL-JEMX 790 03/04/09 0.61 28.8
INTEGRAL-ISGRI 789 31/03/09 42.4 15.5
XMM-RGS 1706 02/04/09 41.5 27.6
XMM-RGS 0336 09/10/01 34.7 33.1
Chandra-MEG 6633 13/08/06 25 125.1
Chandra-MEG 6634 20/10/06 25 172.8
Chandra-MEG 7032 05/11/06 46 146.2

Rates refers to the full band for each instrument.
Rate from selected pixels of the OoT events.

Table 1: The multi-instrument observation log for 4U 1820-30.
Figure 2: Light curve of EPIC-pn data as selected in Sect. 2.1. Note that the displayed count rate cannot be used to recover the flux of the source. The bin time is 500 s.
Figure 3: Upper panel: residuals to a simple absorbed power law for EPIC-pn and INTEGRAL data. JEMX data have been highlighted for clarity. Lower panel: residuals to the best fit consisting of a black body plus Comptonized emission, modified by both a neutral and ionized absorbing material. EPIC-pn data have been rebinned for displaying purpose.
Figure 4: Upper panel: The best fit model consists of a black body (dashed line) plus Comptonized emission (dot-dashed line), modified by both a neutral and ionized absorbing material (Table 2). Lower panel: Best fit model superimposed to the EPIC-pn, JEMX and ISGRI spectra.

Units are: cm for column densities, keV for temperatures, Fluxes are measured in erg cm s.
from RGS data
from JEMX data
from JEMX and ISGRI data.

Table 2: Broad band modeling of the source using EPIC-pn and INTEGRAL data.

3 The broad band spectrum

We first extracted the light curve from the EPIC-pn OoT events as described in Sect. 2.1, therefore the resulting count rate cannot be used to recover the absolute flux. The 0.5–10 keV light curve is shown in Fig. 2, using a bin size of 500 s. It shows a modest variation of at most 13% during our observation. The rates in the soft (0.5–2 keV) and the hard (2–10 keV) band followed the same pattern.
We fitted the EPIC-pn and INTEGRAL spectra simultaneously. The source was caught at a relatively high flux, with a 2–10 keV flux of  erg cm s (Table 2). In Fig. 4 (upper panel) we show the residuals with respect to a simple powerlaw with photoelectric absorption. From the residuals we see a complex spectrum at energies below 1 keV, probably due to ionized absorbing gas. The spectrum is well detected up to 40 keV.
As a simple powerlaw is not an acceptable model (, where is the number of degrees of freedom), we added first a black body component, to mimic the often observed soft energy emission (). This is parameterized by the black-body temperature () and its normalization. However a black body plus power law model fails to explain the spectral curvature seen at high energies. We then substituted the powerlaw model with a Comptonization model (COMT model in SPEX; Titarchuk 1994), where the parameters are: the temperature of the seed photons (), the optical depth of the electron cloud () where the photons are Compton scattered to the final temperature (). This model provides a satisfactory fit to the broad band continuum (; see Fig. 4, lower panel, Fig. 4, and Table 2).

We repeated the fits using a disk black body model (i.e., emission from a standard Shakura-Sunyaev disk, model DBB in SPEX) instead of a simple black body. We do not obtain an equally good fit using DBB either in addition to a simple powerlaw model () or DBB plus Comptonization ().

In the 6.4 keV region, there is no clear evidence of iron emission lines. A 3  upper limit of 25 eV on the equivalent width is obtained if a delta line is put at a fixed energy of 6.4, 6.70 and 6.97 keV, respectively. We also searched for a relativistically modified line profile (LAOR model in SPEX; Laor 1991). In this model, the line arises from the accretion disk where the intensity follows a profile, where is the distance of the disk gas from the source. Leaving the energy as a free parameter did not lead to a meaningful fit. We then fix the energy of the line at 6.97 keV, following previous studies of 4U 1820-30 (e.g. Cackett et al. 2010). The best fit yelds a very broad, but weak line profile, with parameters and , where is the disk inclination. The flux of the line is  photons cm s.

The complex absorption at soft energies has been modeled following the prescription given by the high-resolution data (Sect 2.2). Therefore we added a minor contribution from ionized absorbers ( cm) and the Galactic neutral absorber, which contributes to most of the low-energy curvature. We note that the best fit neutral column density is about 40% lower than the best fit obtained in the RGS analysis. This may be caused by foreground scattering halo soft-emission (observed for the full exposure) located just in front of the OoT event (observed for 6.3% of the exposure time) which still remains after the subtraction of the local background. The scattering halo appears indeed as diffuse emission which extends on top of both the wings of the point spread function of the source and the OoT events themselves (e.g. Predehl & Schmitt 1995; Costantini et al. 2005). The scattering process is strongly energy dependent (Mathis & Lee 1991). In particular, the halo is brighter at softer energies. The net effect is therefore to add more photons to the soft energy spectrum of the OoT events, reducing the measured absorption toward the source. This effect appears even more enhanced by the reduction of the OoT absolute flux due to the cut of the central two columns, where most of the source photons are.

4 The high resolution spectrum

Since the energy band of the RGS is insufficient to discriminate among different broad-band models, we adopted the model defined in the previous section. The spectrum shows evident O i K and Fe i L edges, due to absorption by neutral material. In addition, absorption by ionized gas is highlighted by the O viii, O vii and O vi absorption lines.
We used a collisionally ionized plasma model (model HOT in SPEX), with a temperature frozen to  keV, in order to mimic a neutral gas. This model fits to first order the edges and resonant lines from the neutral species, leaving however noticeable residuals in the fit. In order to extract more robust results from the absorbtion features, we added archival RGS (Sect. 2.2) and Chandra-MEG data (Sect. 2.3). As these data sets were obtained at different epochs, the continuum shape may differ significantly. To bypass this complication, we used the ”sectors” option is SPEX666See Ch. 5 of the SPEX cookbook: http://www.sron.nl/files/HEA/SPEX/manuals/spex-cookbook.pdf and we left the continuum parameters free to vary for each sector. Note that the RGS1 and RGS2 of the recent observations were treated as different data sets (thus with different continua), as the RGS2 broad band shape was affected by pileup. The parameters of absorber components were coupled together for all the data sets.

4.1 The ionized gas

The RGS spectrum of 4U 1820-30 displays narrow absorption features from ionized gas, mainly from oxygen, iron and neon. In particular Ne ix, O vi–O viii are prominent features of this gas. For this fit we ignored the data below 13 Å for RGS. In this region some wiggles in the local continuum, mainly caused by pile up, could add uncertainties in the determination of the absorption parameters. Here we focus on the neutral or mildly ionized phase of the ISM. Therefore we simply model the higher ionization lines with a phenomenological model (SLAB in SPEX) which calculates the transmission from a thin layer of gas. We defer a more detailed analysis of this gas component to a future publication. The resulting ionic column densities of the main ions, as measured by RGS, are reported in Table 3. Following Yao & Wang (2005), we assume a velocity broadening of  km s. All subsequent fits in the present analysis already take into account this highly ionized component.

ion log
C vi
N vi
N vii
O vi
O vii
O viii

The assumed velocity broadening is  km s.

Table 3: Parameters of the main lines of the more ionized ions as measured by RGS.

4.2 The oxygen edge

The oxygen edge at 0.538 keV has been fitted using the high quality data of the two RGS epochs between 19–36 Å. The neutral material, modeled by a low temperature ( keV) in a collisional ionization equilibrium, fits well the long wavelength curvature due to the ISM absorption as well as the O i 1s-2p absorption line at 23.5 Å (Fig. 7). The edge shape appears modified by several sharp absorption features. Some of them can be easily identified with the O vii and O vi lines, belonging to a highly ionized phase (Sect. 4.1). At the shorter wavelength side of the O i line, we see O ii at  Å and also weaker O iii and O iv lines. Two absorbers with different temperatures are required to fit these lines (Table 4, column (1)). The reported line widths are automatically evaluated in SPEX using a curve of growth analysis applied to multiple lines belonging to the same absorption system (e.g. Spitzer 1978). This allows to accurately evaluate the line width using basic parameters such as the ionic column density and equivalent width. In this case, O i, Fe i, N i transitions are used for the cold gas (comp 1 in Table 4). The strongest transitions are saturated, which in principle may introduce additional uncertainty in determining the line width. However, as shown in Fig. 5, for the lines we study here, the ionic column density as a function of the line equivalent width depends marginally on the velocity width. This behaviour of the curve of growth is due to the large -Voigt parameter for the inner-shell transitions (Mihalas 1978; Kaastra et al. 2008, for a full discussion). We note that the nitrogen region is not modified by dust and the line widths are then easier to evaluate. In addition, many of the the weaker transitions of the same neutral ions are not saturated. For the  eV gas (comp 2), O ii, Fe ii, N ii are among the strongest ions. These are relatively weaker lines where saturation does not play a major role. The relatively limited resolution of the instrument does not allow to distinguish multiple velocity-width components within a same absorption line. Therefore we consider here the total ionic column density along this line of sight.

Figure 5: Line equivalent width as a function of the ionic column density for O i, N i and O ii. The velocity dispersions are 10, 100, 200 km s. Curves of different ions have been shifted for clarity.

Dust should also be present along this line of sight, as 4U 1820-30 is known to display a dust scattering halo (Predehl & Schmitt 1995), which is also evident in the present EPIC-pn data. We therefore fitted the oxygen edge adding the AMOL model in SPEX. We considered 16 oxygen compounds for this analysis (see Appendix A). This list includes silicates and oxides which are common in the ISM plus some lighter (also icy) materials. Absorption features by dust, contrary to gas, are smooth and broadened. All the oxygen compounds considered here have no appreciable effect at  Å (Fig. 10).
Below we describe a simple fit using a single dust component at a time, adding a second component if the fit requires it. This makes the description of the model and of the steps to the final best fit simpler to follow. This approach is however the result of a rigorous fitting procedure for both the oxygen and iron edges, which is fully described below (Sect. 4.4).

We first attempted to fit the whole spectral region only in terms of dust absorption. In the fit we let the hydrogen column densities of the gas producing O i and O ii (comp 1 and 2 in Table 4) as free parameters. We also let the abundance of oxygen in the cold phase (comp 1) as a free parameter, in order to balance the amount of oxygen locked up in dust. In this model, we suppressed the gas component (comp 3) mainly producing O iii ( Å), as this line is the main absorption structure in the region where also dust plays a major role. It is clear that fitting the region around 23 Å only in terms of dust remarkably worsens the fit. The worsening range is indeed , depending on the compound, for , with respect to the pure gas modeling. This indicates that not only dust is absorbing the spectrum in the  Å region. In Fig. 6, upper panel, we show the oxygen region where the pure gas fit is compared with the fit with MgSiO (dashed line) or MgFeSiO (olivine, dashed-dotted line), suppressing comp. 3.

Finally, we included both the  eV gas absorption (comp 3) and dust absorption (Table. 4, column (2)). The improvement of the fit is then for . Most of the individual dust components improve the fit. However, we can identify two compounds for which the improvement of the fit is the largest (). These are HO, in the form of crystal or amorphous ice, and MgSiO. In Fig. 6 (lower panel) we show how a fit with pure gas and a gas+MgSiO mixture (comp 5 in Table 4) compare with each other. Note that a fit using HO is non distinguishable in practice from MgSiO as the spectral feature is very similar (Fig. 10, Sect. 4.4 for discussion). Interestingly we note that the compounds with both oxygen and iron provided the least improvement (if not a worsening) of the fit. In particular, the fit rules out the most complex aggregates of oxygen and iron (e.g. magnetite, franklinite, olivine and almandine, see Fig. 10 for the chemical composition).

As noted above, none of the dust compounds considered here have features below 23.7 Å. However, from Fig. 6 we note two clear absorption features (present in both RGS data sets) at and Å (marked with arrows) detected with 6.8 and 3.8 significance, respectively. These features do not belong to the ionized absorber producing the O vii line (Sect. 4.1), but to a lower ionization gas. We tentatively fitted these features with a photoionization model (XABS in XSPEC), obtaining a low column density ( cm) and a ionization parameter777The ionization parameter is defined as , where L is the ionizing luminosity, the gas density and the distance of the gas from the source. log (comp 6 in Table 4). Interestingly, the two major lines, identified as O iv and O v present a systematic blue-shift of about  km s (Table 4, column (3)). From the absorber model we notice that the other lines predicted using these parameters fall either in a lower-resolution and noisier part of the spectrum (e.g. lines of iron) or where the spectrum is heavily absorbed by the neutral gas (e.g. nitrogen lines).
The best fit in the oxygen edge region, including RGS data sets taken in two epochs is displayed in Fig. 7.

Figure 6: Fits in the oxygen region. Here only RGS1 of the most recent observation is displayed for clarity. Upper panel: comparison between a fit only in terms of gas (solid line) and the best fit obtained using dust (dashed line). For comparison also an unacceptable fit using olivine (dashed-dotted line) is displayed. Lower panel: comparison between a fit only in terms of gas (solid line) and the best fit obtained with a mixture of gas and dust (dashed line). See more details in the text.
Figure 7: Best fit transmission spectrum (Table 4, column (3)) in the oxygen edge region. This allows to compare absorption at different epochs, removing the contribution of the continuum. The upper curves display the transmission of the various absorbing components. Note that for clarity only absorption by oxygen is displayed (i.e. no hydrogen absorption). Data have been rebinned for clarity.
(1) (2) (3)
comp 1
0.5 fix. 0.5 fix. 0.5 fix.
20 fix. 20 fix.
comp 2
10 fix. 10 fix.
comp 3
100 fix. 100 fix. 100 fix.
comp 5
comp 6
log 1
5863/4572 5745/4571 5684/4568

Units are  cm for column densities (), eV for temperatures () and km s for line broadening () and the outflow velocity ().
Column(1): model with neutral and mildly ionized gas. Column (2): mildly ionized gas + dust. Column (3): milidly ionized gas + dust + mildly ionized outflowing gas. In all models the highly ionized component is already included.

Table 4: Parameters for the oxygen edge modeling.

4.3 The iron L edge

In Fig. 9, a fit of the Fe LII and LIII edges in terms of absorption by pure gas with solar abundances is shown (dotted line). In this fit two RGS epochs and the MEG data sets were used. There is a clear mismatch not only in the depth, but also in the position of the of the LII and LIII edges. We note that in the fit of the iron edges, it is essential to include the higher resolution MEG data (Fig. 9). A straightforward way to improve the fit is to modify the iron abundance. Letting this parameter free, we obtain a reasonable fit, although the position of the edge in the model still does not match the data. In this case the abundance of iron, assuming that the absorption is only in the gas phase, is about 0.37 times the solar one. However, absorption by dust is known to alter the LII/LIII ratio and to shift the position of the edge with respect to absorption by gas (see e.g. Lee et al. 2009). For the fitting, we considered the compounds measured by Lee et al. (2009): metallic Fe, hematite (FeO), lepidocrocite (FeO(OH)), fayalite (FeSiO) and iron sulfate (FeSO). These transmission models have been recently implemented in the AMOL model in SPEX. As for the modeling of the oxygen edge, we first followed a rigorous approach (Sect. 4.4) which in turns justifies the simpler approach described here, i.e. fitting the dust components one by one, together with the gas model. In the latter, the iron abundance is a free parameter. This procedure provides a measure of the depletion of the gas phase. On the goodness of fit bases, the single dust compound which best models the data is metallic iron (Table 5). Then we tried all possible combinations of gas plus two dust components. We only find a marginally significant presence of FeO, in addition to the metallic iron component (Table 5, Fig. 9, solid line). Hematite (FeO) is a compound for which both iron and oxygen edges’ dust profiles are implemented in the AMOL model. We therefore tested the presence of this compound in the oxygen edge region, obtaining an upper limit on the column density, consistent with the value obtained in the iron region. In Fig. 9 we also show the comparison between the best fit model and a model with an iron rich olivine (fayalite, FeSiO), which does not provide a good fit (Table 5). This suggests that fayalite cannot have a major contribution to the absorption here. We cannot test the same kind of olivine as for the oxygen edge (MgFeSiO) as at present no laboratory measurement of this compound at the iron L edge is available. In Fig. 9, the best fit is shown where data from both RGS and Chandra-MEG are displayed. The spectra have been normalized to their continuum shape for displaying purposes. We also normalized to the total hydrogen column density so that only absorption by iron is visible for the neutral phase. Along with the best fit (solid line), also the different contribution are shown: atomic gas (dotted line) and metallic iron and FeO (light dash-dotted lines). Also the contribution of the ionized gas to the spectrum is shown with a dashed line (e.g. O vii line at 17.76 Å).
Despite the significant improvement of the fit with respect to a gas-only model, the iron edge is clearly more complex that our parameterization. In particular, positive residuals are present on the longer wavelength side of the edge. This effect has been found previously using different instruments (Kaastra et al. 2009; Lee et al. 2009).

Gas ( ) 2377/981
Gas () 1544/980
FeSiO 1505/977
FeSO 1486/977
FeO(OH) 1464/977
FeO 1442/977
Fe met. 1424/977
Fe met.+FeO 1421/976
Table 5: Goodness of fit for the iron edge models.
Figure 8: Detail of the iron L edge region. Here for clarity we display only the MEG and a displaced RGS2 data set. The best fit (solid line, see Table 5) is compared with pure gas fit with solar abundances (dotted line) and with a mixture of gas and iron-rich olivine (thick solid light line).
Figure 9: Best fit transmission spectrum of the region around the iron LII and LIII edges. This allows to compare absorption at different epochs, removing the contribution of the continuum. The upper curves display the transmission of the various absorbing components.

4.4 Further insight in the Fe L and O K edge fitting

The modeling of photoelectric edges modified by dust absorption requires some caution, as in principle many components can contribute to the edge shape, which is in turn smeared, contrary to the atomic transitions features, which are very sharp and relatively easy to identify. Therefore the accuracy of the result depends on the resolution of the instruments, the signal-to-noise ratio, the completeness of the dust models and finally on how many independent parameters there are in the problem.
In this study we combined the higher resolution of Chandra-MEG with the sensitivity of RGS to study the Fe L edges. For the oxygen edge we used the high quality RGS data, as the Chandra-HEG effective area dramatically drops in that spectral region. 4U 1820-30 is also a high flux source, absorbed by a gas and dust column density which is optimal in order to study the O and Fe edges.
The dust edge profiles that we use in our study are a collection of what is available in the literature up to now (e.g. Barrus et al. 1997; Parent et al. 2002; van Aken et al. 1998; Lee et al. 2008, 2009). This data-base is not complete, but the main representatives of oxides and silicates for both oxygen and iron are present (see Appendix A, Lee et al. 2009). This is sufficient for a first order analysis.
The simultaneous analysis of two edges produced by dust on the same line of sight allows to cross-validate our results. In addition, common knowledge on abundances and physical conditions in the ISM allow us to put further constraints to the edges fitting. Here we illustrate this in detail.

In the ISM a mixture of chemical compounds exists. We tried to simulate this by mixing all possible combinations of edges profiles available for a given edge. The AMOL model allows to fit four different components at a time. We therefore had for each edge a number of fits where . Here is the number of available edge’s profiles. Examining the output we notice that the combinations which improve the fits prefer one or two components, putting the remaining column densities to zero. If two components are preferred, one strongly dominates over the other, whose relative value is never more than 20%. This applies to both the oxygen and the iron edge. This justifies the simplified way of illustrating the fitting used in Sect 4.2, 4.3 above.

We further analyzed all fits that were within 3 sigma the best fit values shown in Sects. 4.2 and 4.3. Interestingly, for both iron and oxygen, no models containing iron-rich silicates are among these selected models.
For iron we could rule out models, within the “3 sigma” set, containing hematite (FeO) with a total molecular column density888Note that this translates in  cm and  cm for O and Fe respectively. larger than  cm. This resulted from the cross-validation we performed with the oxygen edge, as for hematite we have both the O and Fe dust profiles. Indeed, none of the “3 sigma” models for oxygen contained this compound. We therefore added it to the oxygen edge best fit model, obtaining the limit above (and compatible with what found for the iron edge best fit, see Sect. 4.3).
We could also rule out the ”3 sigma” models containing FeS as for this case the fit delivered an amount of sulfur 3.8 times solar. The sulfur edge would have been significantly enhanced, but this is not evident in our data.
For the oxygen edge this multicomponent fitting pointed very clearly to a predominance of MgSiO in all the ”3 sigma” set. As noted in Sect. 4.2, a similarly good fit is obtained if water ice is the predominant component. This is because of the very similar edge profiles of the two compounds (Fig. 10). However, we could readily rule out this option, as along the line of sight of 4U 1820-30 the bulk of cold material is diffuse interstellar medium, where water is virtually absent and silicates are on the contrary abundant (e.g. Wooden 2008).

4.5 The magnesium and silicon edges

The Mg K and Si K edge (at 9.47 and 6.723 Å, respectively) were analyzed using the combined information of Chandra-HEG and MEG, which provide in this band a higher energy resolution. In the case of Si, the energy of the edge is also outside the RGS range. The quality of the data and the relatively low Galactic column density hampered any detailed characterization of the edges. We used the auxiliary information provided by the oxygen edge fitting (Sect. 4.2) where a model including MgSiO was preferred. For that compound we could benefit of the laboratory measurements of both the O and Si edge. For Mg, for which we did not have any measurement, the relative column density was included in the model in the form of an artificial gas-shaped edge. We therefore added to the gas absorption model this specific dust component, fixing the amount of magnesium and silicon locked in dust and letting the amount of Mg and Si in the gas phase as free parameters. This model agrees with the data, within the uncertainties. The column density of Si in dust is  cm, while the Si in gas is  cm. Therefore the inclusion of Si in dust is 85% of the total neutral Si along the line of sight. For the Mg edge, we obtain a  cm, while the value imposed by the oxygen fit to the Mg locked in dust is  cm. The amount of Mg locked up in dust according to this model is about 96%.

4.6 The neon and nitrogen edges

In order to study the neon edge ( Å), we used the combined data of RGS and Chandra, while for the nitrogen edge we used only the RGS data. These elements are only present in gas form in the ISM (e.g. Wilms et al. 2000). Fitting neon with our Galactic absorption model we find a column density of  cm, which converts in a slight formal overabundance with respect to solar .
Next to the nitrogen edge at 30.76 Å, we can also distinguish absorption by the 1s-2p line at 31.28 Å. The measured column density of nitrogen is  cm, corresponding to an abundance of times solar, consistent with the value expected using solar abundances.

5 Discussion

5.1 The broad band continuum and the iron emission line

We have analyzed the simultaneous observation of 4U 1820-30 using XMM-EPIC-pn and INTEGRAL. The fit of the continuum is consistent with black body emission plus a Comptonized component, the latter extending up to 40 keV. The source was observed in a high flux state, accreting at . We have assumed , where the unabsorbed luminosity is  erg s. The Eddington luminosity for ultracompact systems is  erg s, following Kuulkers et al. (2003). The spectral shape is also consistent with a high state scenario. The hard energy portion of the spectrum is well fitted by a Comptonization model with an effective cut-off at  keV, consistent with previous observations of the source in the high-state (Bloser et al. 2000; Sidoli et al. 2001). The black body emission, which fits well the soft energy range of the present observation, has a temperature similar to that of the comptonization seed photons. This may suggest that in ultra-compact X-ray binaries seed photons originate from the accretion disk (Sidoli et al. 2001).

A relatively weak or undetected iron line is a common feature of the confirmed ultra-compact X-ray binaries (as defined e.g. in in’t Zand et al. 2007), as 4U 1820-30 is. This may be simply explained by the nature of the donor star in such systems. Due to the high surface gravity (log ), heavier elements in the white-dwarf companion sink in the internal layer of the star (e.g. Fontaine & Michaud 1979), leaving little high metals for the accretion flow. In 4U 1820-30 however, a faint, variable and relativistically smeared iron line as been reported (e.g. Cackett et al. 2008b, 2010). In our data, no narrow emission features, either from neutral or ionized iron, have been found (Sect. 3). A fit with a relativistically smeared iron line leads to a very elusive 2-sigma detection of the feature. The photon flux of this line is about 16 times fainter than that found in recent data (Cackett et al. 2010).
In exceptional circumstances, like superburst episodes, a large quantity of heavy elements ashes may be ejected during the neutron star photosphere expansion. This could possibly cause detectable features (e.g. in absorption) by e.g. iron in the X-ray spectra (in’t Zand & Weinberg 2010). In 4U 1820-30 the presence of the iron line could then be the tracer of burst activity, not detectable at the epoch of our observation.

5.2 The dust components

The high resolution spectrum of 4U 1820-30 shows evidence of absorption by dust, especially around the Fe L and O K edges (Figs. 7, 9). In the fits we explored all the possible combinations of dust mixtures (Sects. 4.3 and 4.2). The best fit shows a preference for absorption by enstatite (MgSiO), which is the end series of the pyroxene silicate mineral, metallic iron (Fe) and traces of iron oxides in the form of hematite (FeO). In the modeling, 16 and 5 different compounds were considered for oxygen and iron, respectively. The set of models does not cover the whole range of possible compounds present in the ISM. However what are commonly believed to be the main constituents (e.g. olivines, pyroxene, oxides as well as simpler compounds) in the diffuse ISM (e.g. Whittet 2003; Wooden 2008) are represented (see Appendix A, Lee et al. 2009). A further effort in obtaining new measurements has recently been carried out (Lee & Ravel 2005; Lee et al. 2009, and de Vries et al., in prep.). The available X-ray measurements were mostly performed on crystalline materials. The presence of dust in the form of crystals is significantly present only in specific astronomical environments, like comets (Wooden 2008) or the inner regions of protoplanetary disks (van Boekel et al. 2004). In the ISM the amount of crystalline silicates, compared to amorphous grains, is % (Li & Draine 2001), therefore we do not expect to find a sizeable amount of it in our data. However, for the handful of compounds for which X-ray laboratory measurements of both glassy and crystalline forms are available, we tested that the spectral difference is only appreciable if the resolution is of the order of  eV. Therefore, with the present data, our analysis using crystalline grains is a good first order approximation of modeling the chemical composition of ID. In the glassy form, crystalline material loses the ordinate internal structure proper of a crystal. For simplicity, although formally not correct, we may call an amorphous silicate with a certain (for instance enstatite) stoichiometry, glassy-enstatite or glassy-MgSiO.

5.3 Abundances and depletion

As X-ray spectra do not display any sharp H feature, either in emission or absorption, the total hydrogen column density has been evaluated from the low energy curvature of the spectrum, which extends down to 36 Å. This is a reliable method, as starting already at about 25 Å, the transmission is largely dominated by He and H. The best fit value for the hydrogen column density is  cm, slightly larger than the nominal values  cm (Dickey & Lockman 1990) and  cm (Kalberla et al. 2005), which are the average values over a region of 1 radius around the source. The reason could be that the spectral curvature measures the total hydrogen (i.e. H i, H ii and H), which should be larger than H i alone. The amount of H in the diffuse ISM is relatively low (e.g. Takei et al. 2002). We estimated the amount of H ii using the O ii/O i ratio as a proxy (following the expression in Field & Steigman 1971). We derive that 4% of the measured hydrogen column density is in the form of H ii. This would reconcile our value with the H i measurement mentioned above. Besides, it could be that a small fraction (in this case no more than  cm) of the total H resides in the immediate surroundings of the source, like in other sources (e.g. Predehl & Schmitt 1995; van Peet et al. 2009). Such a low column density would be insufficient, however, to produce deep absorption features. We compute the following abundances taking as a reference the value we measure from the X-ray spectrum.

In the iron edge region (Sect. 4.3), a fit in terms of pure gas is clearly unacceptable. This points easily to a further contribution from dust. According to our best fit, the total column density of gaseous iron is then  cm. From the iron edge fit we obtain a total dust column density of  cm. The sum of these iron components provides  cm. When we compare this number to the value predicted from our set of solar abundances,  cm (given the best fit total hydrogen column density of  cm) we obtain that the abundance is times solar. The depletion of iron, meant here as the ratio of dust over the total amount of a given element is instead (Table 6). This value is often reported to be higher (e.g. 0.97, Jenkins 2009). However, other studies report lower values for the Fe depletion (0.7, Wilms et al. 2000). We note that, considering some uncertainties that still remain in the edge fitting (Sect. 4.3) the value we obtain may be considered as a lower limit.
At a temperature of 0.5 eV, Fe is mainly Fe i and only 2% of the total iron in Fe ii. The sum of the mildly ionized phase (Fe ii–Fe iv), excluding therefore the high ionization ions described in Sect. 4.1, is  cm along this line of sight, i.e. Fe/Fe0.3, remembering that neutral iron is % depleted.

In the oxygen region the evidence for dust is not as striking as in the iron region (Sect. 4.2). From the best fit including cold ( eV) gas, the amount of oxygen in the gas form is consistent to be solar ( cm). However, the amount of oxygen in dust compounds is  cm. The total amount of oxygen is then  cm, showing about 23% overabundance with respect to solar (Table 6). With this reference value, the amount of depletion is . Oxygen depletion and abundances, based on Chandra data have been previously presented by Yao et al. (2006). As expected, the column density of O i () that we measure is similar to theirs. However, taking as a reference the Anders & Grevesse (1989) abundances (), Yao et al. (2006) find an underabundance of oxygen of about 30%. In the reference list that we adopted (Lodders & Palme 2009), the absolute oxygen abundance is significantly lower ().
Finally, the amount of mildly ionized oxygen (O ii–O v, with total column density  cm) over the total oxygen in gas form is O/O0.05.

As shown in Sect. 4.5 the inclusion of Si in dust is 85%, while for Mg is 97% of the total ISM. These values are comparable to what has been reported in the literature (e.g. 80–92 and 90–97% for Mg and Si respectively, Wilms et al. 2000; Whittet 2003). However, Mg and Si could not be studied in detail in this case because of the relatively low column density toward 4U 1820-30 which imprints shallow absorption edges. Moreover, our fitting of the Mg and Si edges relies on the oxygen modeling. Therefore, if other Mg or Si compounds are present, here they are difficult to detect. Keeping in mind these limitations we report the depletion values of Mg and Si in Table 6, for completeness. The abundance estimate of both Mg and Si are again driven by the oxygen edge fitting and are formally slightly above the solar values.
In Table 6 we list for each element, the ratio with Proto-solar abundances and the amount of depletion defined as the ratio between dust and total ISM abundance.
A further test on the abundances derived from the X-ray data would of course come from high-resolution () UV data, where the elemental gas phase can be accurately studied (e.g. Savage & Sembach 1996). In the case of 4U 1820-30, high-resolution UV data (either from HST or FUSE) are not available. 4U 1820-30 has been observed only with the HST-STIS-G140L spectrograph (). We find that the signal to noise ratio of the data was insufficient for a quantitative study on the absorption lines.

5.4 The location of the cold matter

It has been established that a gradient in the abundances in our Galaxy exists for the most abundant metals. An average slope of 0.06 dex kpc should roughly apply to the gradient for both O and Fe (see Sect.1).
In our analysis we estimate a slight overabundance of oxygen (by a factor , see Table 6), while iron is a factor of the solar value. This appears to be in contradiction with the expected abundances at the distance of the source, where according to the gradient above, some overabundance is expected. However, 4U 1820-30 is located at latitude , i.e. about 1 kpc below the Galactic disk, therefore our line of sight intercepts a relatively small fraction of the cold ISM near the source. The cold phase that we detect towards this source is rather due to absorption in the environment close to the Sun.
For iron, a gradient as a function of both the radial distance from the Galactic center and the height above the disk has been estimated (e.g. Chen et al. 2003, based on open clusters measurements). Our value of [Fe/H]999defined as log(Fe/H)–log(Fe/H) is –0.07. This is consistent, within the errors, with absorption in the disk at the distance of the Sun rather than absorption local to the source, at the height of 4U 1820-30 below the disk (Chen et al. 2003).
The location of the gas based on the oxygen abundance is difficult to determine, as there is a large scatter in the measurements (e.g. Rudolph et al. 2006, for a compilation of results). Our oxygen abundance fits well with both absorption far or near the observer. It is therefore not straightforward to understand where the bulk of the absorption takes place on the basis of gradient measurements. However, a systematic analysis of dust scattering halos points out that most of the scattering (and therefore the absorption by dust) should happen within from the Galactic plane (Predehl & Schmitt 1995). This further support the intuitive idea that most of the cold absorbing material is located close to the observer.

elem. dust/ISM
N 0.0 0.0
Ne 0.0 0.0

Notes: Abundances are referred to Lodders & Palme (2009)
Column densities are in units of  cm

Table 6: Relative abundances and depletion values from the present analysis.

5.5 The silicates

In the present study, we find that the fit of both the iron and oxygen edges in general rejects silicate models containing iron (e.g. Fig. 69). In the iron region, the shift of the L III edge is consistent with absorption by metallic iron. Also a modest quantity of iron in the form of FeO is allowed by the fit. As discussed above, in the oxygen region the evidence for dust is not striking. Our results may be still partially contaminated by instrumental effects, such as bad pixels in the RGS, which are however included in the response files using the most updated calibration. We keep limitations in mind when discussing the physical implications of this result.

In summary, we find that absorption by dust is mainly caused by metallic iron and glassy-enstatite. This kind of composition is reminiscent of the composition of GEMS (Glass with Embedded Metal and Sulfides), that are small grains abundant among the interplanetary dust particles (Bradley 1994). Most of GEMS particles should not have an ISM origin, but they rather reside in interplanetary environments (Keller & Messenger 2008). However, GEMS with anomalous composition may have been processed in the ISM (Matzel et al. 2008, and references therein). In particular, some of those particles show a low amount of sulfur relative to silicon (Keller & Messenger 2008, S/Si), which is more similar to what is found in the diffuse ISM (Sofia 2004) rather than in the solar neighborhood (Anders & Ebihara 1982). However, we cannot yet test the contribution of sulfides in a typical GEMS in our data, as FeS laboratory measurements around both the S K- and the Fe L-edge are not available. From the model, the glassy-enstatite MgSiO provides the best fit. Therefore the Mg/Si ratio is 1 by definition. This is also the value found for anomalous composition GEMS which are possibly of ISM origin, (Keller & Messenger 2008). Unfortunately we have not yet the means of testing compounds with a varying amount of Mg (or Fe). Therefore, we cannot test wether a Mg/Si ratio of 0.6 (typical of the average GEMS in interplanetary dust Keller & Messenger 2004; Ishii et al. 2008) would be still allowed by the fit.
We find that the Fe/Si ratio ranges from 0.42–0.55, which matches this ratio in any type of GEMS where Fe/Si=0.43–0.54 (Keller & Messenger 2004, 2008). This is different from what is expected in the diffuse ISM (Fe/Si=, Sofia 2004). However, the Fe/Si ratio is currently under debate. Recent studies have revised this value, lowering it in fact significantly to roughly Fe/Si (see Min et al. 2007, for a discussion).
Finally, the amount of Fe along this line of sight is about two times less than Mg (Mg/Fe). This estimate may be biassed because we cannot tune the amount of Mg, which is bound to be equal to Si in MgSiO. Another caveat that should be kept in mind is that large grains would be grey to X-ray radiation (Whittet 2003). Therefore, at least in absorption, the large grain population (which may contain more iron) can be under-represented (Lee et al. 2009).

Even taking into account the limitations imposed by our models, we find that, from both the O and Fe edges fits, Mg-rich rather than Fe-rich silicates are present along the line of sight to 4U 1820-30. Such a result strengthens previous studies in both IR (Min et al. 2007) and X-rays (Costantini et al. 2005) which point to the same conclusion. Interestingly, a similar spectral analysis, but restricted to the Fe L edge only, shows that oxides, rather than olivines or pyroxene are responsible for the absorption along the line of sight to Cyg X-1 (Lee et al. 2009). The column density towards 4U 1820-30 is about 4.5 times smaller than for Cyg X-1. This result may therefore point out a chemical homogeneity on different path lengths within the diffuse ISM.

5.6 A fast outflowing gas?

The oxygen region of the RGS spectrum of 4U 1820-30 displays evidence, clearly detected in two set of RGS observations, of two absorption lines, which are consistent with O iv and O v outflowing at  km s(Fig. 7). In the 4U 1820-30 accretion flow heavy elements should not be abundant, as the companion is classified as a He-white dwarf, on the basis of X-ray binary evolutionary models (Rappaport et al. 1987). From optical spectra analysis Nelemans et al. (2010) found that elements heavier than Ne are absent in ultracompact systems. Oxygen is a relatively light element which could still be present in the outer envelope of the companion and be transferred in the accretion disk. Moreover, a way of producing oxygen could be through triple- burning of He in the white dwarf. Qualitatively, this could justify the detection of an outflow containing oxygen. However this result needs to be tested with additional observations.
Here we fitted the absorption in terms of a photoionized gas. Other lines are predicted by the model, but they are too weak or in a noisy part of the spectrum to be significantly detected. This system is well detected in both RGS data sets, taken 8 years apart. The outflow did not change its physical parameters during that time. We note that this models could not be tested against the Chandra-HETG data, because of the noise affecting the oxygen region. The model predicts also a substantial amount of C iv. We qualitatively checked STIS-G140L low-resolution data (taken in 1998, about 3.5 years before the first RGS measurement) for the presence of a blueshifted C iv doublet. Fixing the outflow velocity to the X-ray value we obtain an upper limit for the C iv column density of  cm, which is almost two orders of magnitude lower than that required by the X-ray model. Therefore, in the hypothesis of an outflow, this must have been absent three and a half years before the first RGS pointing, when the STIS observation was carried out. Fast outflows are not unusual in X-ray binaries and are interpreted as accretion disk winds (e.g. Miller et al. 2006). However, the ions involved in the outflow are generally more highly ionized (e.g. O vii, Ne ix-Ne x).
In terms of goodness of fit, this outflowing system can be equally well fitted by a collisionally ionized plasma with temperature  eV and  cm. The amount of C iv predicted by this model would be  cm cm, less than what is predicted by the photoionized gas, but still in disagreement with the upper limit derived from the STIS data. However, higher-quality STIS data are necessary for a reliable comparison between the UV and X-ray band outflowing absorber.
Fast moving collisionally ionized clouds in the line of sight with such high velocities are not reported by the UV surveys (Savage et al. 2004). Therefore the phenomenon should arise in the proximity of the source, possibly in a less ionized impact region, that would cause ionization by collisions. A candidate is the spot where the accretion flow impacts the accretion disk, which should be less ionized than the disk itself (Boirin et al. 2005; van Peet et al. 2009). The absence of emission lines from this gas implies that, regardless of the absorbing mechanism, the flow pointing toward the observer must be very collimated. Although plausible, the scenarios depicted above need support from further (multiwavelength) observations.
Finally, absorption by dust seems unlikely, as dust features are generally smoother (Appendix A).

6 Conclusions

In this paper we present the X-ray analysis of the continuum of 4U 1820-30 (using a quasi-simultaneous observation of XMM-Newton and INTEGRAL) and of the absorption features due to the cold matter in the line of sight (using XMM-Newton-RGS and Chandra-MEG data).

The continuum shape and the Eddington ratio show that the source was caught in a high-state. The continuum is well fitted by black body emission plus a Comptonization component which extends up to 40 keV. We do not find evidence of iron emission, either from neutral or ionized matter. This may be naturally explained by the metal-poor accretion stream expected from the white-dwarf companion.

The absorption spectrum shows the presence of many components with different ionization. We focused on the cold and mildly ionized phase only. Oxygen has been found slightly overabundant by a factor 1.23 times the solar value. Iron is on the contrary slightly underabundant (0.85 times solar). The abundance values are not dramatically deviating from the solar ones and do not allow us to assign a precise location of the absorbing gas. However, a location close to the observer seems likely.

Thanks to the simultaneous study of absorption by dust and gas we measured also the element depletion. Oxygen is mildly depleted by a factor about 0.20. The depletion of iron is more evident, as the depletion factor is 0.87. The depletion of Mg and Si are more difficult to determine. We find that they are depleted of a factor and , respectively

We modeled the dust contribution with the currently availables absorption profiles of dust compounds. Our conclusions bear the uncertainty due to the still limited dust data-base and a lower sensitivity in selected spectral regions. However we clearly find that both the oxygen and iron edges cannot be fitted by iron-rich silicates. On the contrary, the oxygen edge is consistent to be mostly absorbed by enstatite (MgSiO, possibly in a glass-form). Metallic iron should be the main absorber in the iron edge. This leads to the interesting possibility that a GEMS-like form (Mg-rich silicates with metallic iron inclusion) of grain may be absorbing along this line of sight. A fraction of the studied GEMS, in particular the sulfur-poor grains, are believed to be of ISM origin and have also been proposed as constituents of ISM. For the first time an X-ray absorption analysis provides a tentative confirmation of this scenario.

Finally, we report the tentative detection of a mildly ionized outflow ( km s), highlighted by the O iv and O v absorption lines. Both a photo- or collisional- ionizing process could fit the lines, leaving open the interpretation on the nature of this gas.

This research made use of the Chandra Transmission Grating Catalog and archive (http://tgcat.mit.edu). We also made use of the Multimission Archive at the Space Telescope Science Institute (MAST). STScI is operated by the Association of Universities for Research in Astronomy, Inc., under NASA contract NAS5-26555. The Space Research Organization of the Netherlands is supported financially by NWO, the Netherlands Organization for Scientific Research. XMM-Newton and INTEGRAL are ESA science missions with instruments and contributions directly funded by ESA Members States and the USA (NASA). Thanks to O. Madej, M. Min, P. Predehl and E. Ratti for useful discussion. Thanks also to V. Beckmann for pointing out the presence of INTEGRAL data taken quasi-simultaneously to our XMM-Newton data.


  • Anders & Ebihara (1982) Anders, E., & Ebihara, M. 1982, Geochim. Cosmochim. Acta., 46, 2363
  • Anders & Grevesse (1989) Anders, E., & Grevesse, N. 1989, Geochim. Cosmochim. Acta., 53, 197
  • Barret et al. (2000) Barret, D., Olive, J. F., Boirin, L., Done, C., Skinner, G. K., & Grindlay, J. E. 2000, ApJ, 533, 329
  • Barrus et al. (1997) Barrus, D. M., Blake, R. L., Burek, A. J., Chambers, K. C., & Pregenzer, A. L. 1979, Phys. Rev. A, 20, 1045
  • Bloser et al. (2000) Bloser, P. F., Grindlay, J. E., Kaaret, P., Zhang, W., Smale, A. P., & Barret, D. 2000, ApJ, 542, 1000
  • Boirin et al. (2005) Boirin, L., Méndez, M., Díaz Trigo, M., Parmar, A. N., & Kaastra, J. S. 2005, A&A, 436, 195
  • Bowen et al. (2008) Bowen, D. V., et al. 2008, ApJS, 176, 59
  • Bradley (1994) Bradley, J. P. 1994, Science, 265, 925
  • Cackett et al. (2008a) Cackett, E. M., Miller, J. M., Raymond, J., Homan, J., van der Klis, M., Méndez, M., Steeghs, D., & Wijnands, R. 2008a, ApJ, 677, 1233
  • Cackett et al. (2008b) Cackett, E. M., et al. 2008b, ApJ, 674, 415
  • Cackett et al. (2010) Cackett, E. M., et al. 2010, ApJ, 720, 205
  • Carrez et al. (2002) Carrez, P., Demyk, K., Cordier, P., Gengembre, L., Grimblot, J., D’Hendecourt, L., Jones, A. P., & Leroux, H. 2002, Meteoritics and Planetary Science, 37, 1599
  • Chen et al. (2003) Chen, L., Hou, J. L., & Wang, J. J. 2003, AJ, 125, 1397
  • Costantini et al. (2005) Costantini, E., Freyberg, M. J., & Predehl, P. 2005, A&A, 444, 187
  • Courvoisier et al. (2003) Courvoisier, T.J.-L., Walter, R., Beckmann, V., et al. 2003, A&A, 411, L53
  • den Herder et al. (2001) den Herder, J. W., et al. 2001, A&A, 365, L7
  • de Vries & Costantini (2009) de Vries, C. P., & Costantini, E. 2009, A&A, 497, 393
  • Dickey & Lockman (1990) Dickey, J. M., & Lockman, F. J. 1990, ARA&A, 28, 215
  • Draine (2003) Draine, B. T. 2003, ARA&A, 41, 241
  • Field & Steigman (1971) Field, G. B., & Steigman, G. 1971, ApJ, 166, 59
  • Fontaine & Michaud (1979) Fontaine, G., & Michaud, G. 1979, ApJ, 231, 826
  • Futamoto et al. (2004) Futamoto, K., Mitsuda, K., Takei, Y., Fujimoto, R., & Yamasaki, N. Y. 2004, ApJ, 605, 793
  • Grindlay et al. (1976) Grindlay, J., Gursky, H., Schnopper, H., Parsignault, D. R., Heise, J., Brinkman, A. C., & Schrijver, J. 1976, ApJ, 205, L127
  • Hasinger & van der Klis (1989) Hasinger, G., & van der Klis, M. 1989, A&A, 225, 79
  • in’t Zand et al. (2007) in’t Zand, J. J. M., Jonker, P. G., & Markwardt, C. B. 2007, A&A, 465, 953
  • in’t Zand & Weinberg (2010) in’t Zand, J. J. M., & Weinberg, N. N. 2010, A&A, 520, A81
  • Ishii et al. (2008) Ishii, H. A., Bradley, J. P., Dai, Z. R., Chi, M., Kearsley, A. T., Burchell, M. J., Browning, N. D., & Molster, F. 2008, Science, 319, 447
  • Jenkins (2009) Jenkins, E. B. 2009, ApJ, 700, 1299
  • Juett et al. (2006) Juett, A. M., Schulz, N. S., Chakrabarty, D., & Gorczyca, T. W. 2006, ApJ, 648, 1066
  • Kaastra et al. (1996) Kaastra, J. S., Mewe, R., & Nieuwenhuijzen, H. 1996, UV and X-ray Spectroscopy of Astrophysical and Laboratory Plasmas, 411
  • Kaastra et al. (2008) Kaastra, J. S., Paerels, F. B. S., Durret, F., Schindler, S., & Richter, P. 2008, Space Sci. Rev., 134, 155
  • Kaastra et al. (2009) Kaastra, J. S., de Vries, C. P., Costantini, E., & den Herder, J. W. A. 2009, A&A, 497, 291
  • Kalberla et al. (2005) Kalberla, P. M. W., Burton, W. B., Hartmann, D., Arnal, E. M., Bajaja, E., Morras,R., Poumppel, W. G. L. 2005, A&A, 440, 775
  • Keller & Messenger (2004) Keller, L. P., & Messenger, S. 2004, Lunar and Planetary Institute Science Conference Abstracts, 35, 1985
  • Keller & Messenger (2008) Keller, L. P., & Messenger, S. 2008, Lunar and Planetary Institute Science Conference Abstracts, 39, 2347
  • Kuulkers et al. (2003) Kuulkers, E., den Hartog, P. R., in’t Zand, J. J. M., Verbunt, F. W. M., Harris, W. E., & Cocchi, M. 2003, A&A, 399, 663
  • Laor (1991) Laor, A. 1991, ApJ, 376, 90
  • Lebrun et al. (2003) Lebrun, F., Leray, J. P., Lavocat, P., et al. 2003, A&A, 411, L141
  • Lee & Ravel (2005) Lee, J. C., & Ravel, B. 2005, ApJ, 622, 970
  • Lee et al. (2009) Lee, J. C., Xiang, J., Ravel, B., Kortright, J., & Flanagan, K. 2009, ApJ, 702, 970
  • Lee et al. (2008) Lee, S. K. et al.  2008, PNAS vol. 105 no. 23 7925-7929
  • Li & Draine (2001) Li, A., & Draine, B. T. 2001, ApJ, 550, L213
  • Lodders & Palme (2009) Lodders, K., & Palme, H. 2009, Meteoritics and Planetary Science Supplement, 72, 5154
  • Lodders (2010) Lodders, K. 2010, Principles and Perspectives in Cosmochemistry, 379
  • Lugaro et al. (1999) Lugaro, M., Zinner, E., Gallino, R., & Amari, S. 1999, ApJ, 527, 369
  • Lund et al. (2003) Lund, N., et al. 2003, A&A, 411, L231
  • Mathis & Lee (1991) Mathis, J. S., & Lee, C.-W. 1991, ApJ, 376, 490
  • Matzel et al. (2008) Matzel, J., Dai, Z. R., Teslich, N., Hutcheon, I. D., Weber, P., & Bradley, J. P. 2008, Lunar and Planetary Institute Science Conference Abstracts, 39, 2525
  • Mihalas (1978) Mihalas, D. 1978, San Francisco, W. H. Freeman and Co., 1978. 650 p.,
  • Miller et al. (2006) Miller, J. M., et al. 2006, ApJ, 646, 394
  • Min et al. (2007) Min, M., Waters, L. B. F. M., de Koter, A., Hovenier, J. W., Keller, L. P., & Markwick-Kemper, F. 2007, A&A, 462, 667
  • Neilsen & Lee (2009) Neilsen, J., & Lee, J. C. 2009, Nature, 458, 481
  • Nelemans et al. (2010) Nelemans, G., Yungelson, L. R., van der Sluys, M. V.,& Tout, C. A. 2010, MNRAS, 401, 1347
  • Nittler (2005) Nittler, L. R. 2005, ApJ, 618, 281
  • Parent et al. (2002) Parent, P., Laffon, C., Mangeney, C., Bournel, F., & Tronc, M. 2002, J. Chem. Phys., 117, 10842
  • Pinto et al. (2010) Pinto, C., Kaastra, J. S., Costantini, E., & Verbunt, F. 2010, A&A, 521, A79
  • Predehl & Schmitt (1995) Predehl, P., & Schmitt, J. H. M. M. 1995, A&A, 293, 889
  • Rappaport et al. (1987) Rappaport, S., Ma, C. P., Joss, P. C., & Nelson, L. A. 1987, ApJ, 322, 842
  • Rudolph et al. (2006) Rudolph, A. L., Fich, M., Bell, G. R., Norsen, T., Simpson, J. P., Haas, M. R., & Erickson, E. F. 2006, ApJS, 162, 346
  • Savage & Sembach (1996) Savage, B. D., & Sembach, K. R. 1996, ARA&A, 34, 279
  • Savage et al. (2004) Savage, B. D., Wakker, B. P., Sembach, K. R., Richter, P., & Meade, M. 2004, Recycling Intergalactic and Interstellar Matter, 217, 147
  • Sidoli et al. (2001) Sidoli, L., Parmar, A. N., Oosterbroek, T., Stella, L., Verbunt, F., Masetti, N., & Dal Fiume, D. 2001, A&A, 368, 451
  • Smale et al. (1997) Smale, A. P., Zhang, W., & White, N. E. 1997, ApJ, 483, L119
  • Sofia (2004) Sofia, U. J. 2004, Astrophysics of Dust, 309, 393
  • Spitzer (1978) Spitzer, L. 1978, New York Wiley-Interscience, 1978. 333 p.,
  • Stella et al. (1987) Stella, L., Priedhorsky, W., & White, N. E. 1987, ApJ, 312, L17
  • Strüder et al. (2001) Strüder, L., et al. 2001, A&A, 365, L18
  • Takei et al. (2002) Takei, Y., Fujimoto, R., Mitsuda, K., & Onaka, T. 2002, ApJ, 581, 307
  • Tarana et al. (2007) Tarana, A., Bazzano, A., Ubertini, P., & Zdziarski, A. A. 2007, ApJ, 654, 494
  • Titarchuk (1994) Titarchuk, L. 1994, ApJ, 434, 570
  • van Aken et al. (1998) van Aken, P. A., Liebscher, B., & Styrsa, V. J. 1998, Physics and Chemistry of Minerals, 25, 494
  • van Boekel et al. (2004) van Boekel, R., et al. 2004, Nature, 432, 479
  • van Peet et al. (2009) van Peet, J. C. A., Costantini, E., Méndez, M., Paerels, F. B. S., & Cottam, J. 2009, A&A, 497, 805
  • Weingartner & Draine (2001) Weingartner, J. C., & Draine, B. T. 2001, ApJ, 548, 296
  • Whittet (2003) Whittet, D. C. B. 2003, Dust in the galactic environment, 2nd ed.  by D.C.B. Whittet. Bristol: Institute of Physics (IOP) Publishing, 2003 Series in Astronomy and Astrophysics, ISBN 0750306246.
  • Wilms et al. (2000) Wilms, J., Allen, A., & McCray, R. 2000, ApJ, 542, 914
  • Wooden et al. (2005) Wooden, D. H., Harker, D. E., & Brearley, A. J. 2005, Chondrites and the Protoplanetary Disk, 341, 774
  • Wooden (2008) Wooden, D. H. 2008, Space Sci. Rev., 138, 75
  • Yao & Wang (2005) Yao, Y., & Wang, Q. D. 2005, ApJ, 624, 751
  • Yao et al. (2006) Yao, Y., Schulz, N., Wang, Q. D., & Nowak, M. 2006, ApJ, 653, L121
  • Zhang et al. (1998) Zhang, W., Smale, A. P., Strohmayer, T. E., & Swank, J. H. 1998, ApJ, 500, L171

Appendix A Oxygen compounds

Here we show the profiles of the absorption around the oxygen edge for the compounds used in this analysis. These and other profiles are included in the AMOL model, implemented in SPEX. We refer to Pinto et al. (2010) for a complete list and description of the oxygen compounds (with the exeption of MgSiO).

Figure 10: Transmission of the dust absorption models included in the present analysis for the oxygen region. The oxygen column density has been set here to  cm for all compounds. a) Barrus et al. (1997), b) Parent et al. (2002), c) van Aken et al. (1998), d) Lee et al. (2008). See also Pinto et al. (2010) for details.
Comments 0
Request Comment
You are adding the first comment!
How to quickly get a good reply:
  • Give credit where it’s due by listing out the positive aspects of a paper before getting into which changes should be made.
  • Be specific in your critique, and provide supporting evidence with appropriate references to substantiate general statements.
  • Your comment should inspire ideas to flow and help the author improves the paper.

The better we are at sharing our knowledge with each other, the faster we move forward.
The feedback must be of minimum 40 characters and the title a minimum of 5 characters
Add comment
Loading ...
This is a comment super asjknd jkasnjk adsnkj
The feedback must be of minumum 40 characters
The feedback must be of minumum 40 characters

You are asking your first question!
How to quickly get a good answer:
  • Keep your question short and to the point
  • Check for grammar or spelling errors.
  • Phrase it like a question
Test description