A&A 450, 725-733 (2006)
DOI: 10.1051/0004-6361:20053768
K. Werner1 - T. Nagel1 - T. Rauch1 - N. J. Hammer2,1 - S. Dreizler3
1 - Institut für Astronomie und Astrophysik, Universität Tübingen, Sand 1, 72076 Tübingen, Germany
2 -
Max-Planck-Institut für Astrophysik, Karl-Schwarzschild-Straße 1, 85741 Garching, Germany
3 -
Institut für Astrophysik, Universität Göttingen, Friedrich-Hund-Platz 1, 37077 Göttingen, Germany
Received 5 July 2005 / Accepted 2 January 2006
Abstract
Aims. We present new medium-resolution high-S/N optical spectra of the ultracompact low-mass X-ray binaries 4U 0614+091 and 4U 1626-67, taken with the ESO Very Large Telescope. They are pure emission line spectra and the lines are identified as due to C II-IV and O II-III.
Methods. Line identification is corroborated by first results from modeling the disk spectra with detailed non-LTE radiation transfer calculations. Hydrogen and helium lines are lacking in the observed spectra.
Results. Our models confirm the deficiency of H and He in the disks. The lack of neon lines suggests an Ne abundance of less than about 10 percent (by mass), however, this result is uncertain due to possible shortcomings in the model atom. These findings suggest that the donor stars are eroded cores of C/O white dwarfs with no excessive neon overabundance. This would contradict earlier claims of Ne enrichment concluded from X-ray observations of circumbinary material, which was explained by crystallization and fractionation of the white dwarf core.
Key words: accretion, accretion disks - binaries: close - X-rays: binaries - stars: individual: 4U 0614+091 - stars: individual: 4U 1626-67
Low-mass X-ray binaries (LMXBs) consist of a neutron star or black-hole accretor
and a low-mass donor star (M
1
). Of particular interest
are those systems with orbital periods
80 min, which is the
minimum period for LMXBs with hydrogen-rich main sequence donors. In these
ultracompact binaries (UCBs) the mass donor must be a non-degenerate
hydrogen-deficient star or a white dwarf (e.g. Verbunt & van den Heuvel
1995). Currently eight such systems with measured orbital periods (
= 11-50 min) are known (Ritter & Kolb 2003).
In the recent past, the existence of a group of five ultracompact systems with neon-rich white dwarf donors has been claimed based on X-ray spectral properties (Schulz et al. 2001; Juett et al. 2001; Juett & Chakrabarty 2003). The donors are then C/O or even O/Ne/Mg white dwarfs that have transferred a significant fraction of their mass to the neutron star, in contrast to the usual wisdom that the donors are the remains of He white dwarfs. This has caused new explorations of the formation of these systems (e.g. Yungelson et al. 2002). Our motivation for studying these systems is that the stripped donor stars offer the possibility to probe the interior composition of white dwarfs, which depends on the interplay of gravitational settling and crystallization of chemical elements.
Three of these five Ne-rich systems belong to the above-mentioned eight UCBs with measured orbital periods, while two of them are believed to be UCBs because of their similar optical and X-ray properties. The class of Ne-rich UCBs consists of:
(1) 4U 1626-67 (
= 41 min)
(2) 4U 0614+091
(3) 2S 0918-549
(4) 4U 1543-624 (
= 18 min, Wang & Chakrabarty 2004)
(5) 4U 1850-087 (
= 20 min).
4U 1626-67 may be regarded as the prototype of this class. The donor's Ne-rich C/O-WD nature is derived from X-ray and ultraviolet spectra that exhibit double-peaked emission lines that obviously stem from the accretion disk (Schulz et al. 2001; Homer et al. 2002).
The close relation of the other four objects (2)-(5) to 4U 1626-67 was
based on extraordinary high Ne/O abundance ratios (when compared to the ISM value) measured from ASCA spectra of (3)-(5) and a Chandra spectrum of (2). The spectra exhibit photoelectric absorption edges of neutral O and Ne in the interstellar medium (ISM) along the line-of-sight, which is suspected to originate from expelled material close to the binary systems. New X-ray
spectroscopic observations of (3) and (4) obtained with Chandra and XMM-Newton
confirmed the ASCA results, although different values for the Ne/O enrichment
have been obtained from a Chandra and a XMM-Newton spectrum of (4) (Juett &
Chakrabarty 2003). In contrast to the ASCA measurement of (5), recent XMM-Newton
and Chandra observations found no evidence of an unusual Ne/O ratio (Sidoli
et al. 2005; Juett & Chakrabarty 2005). These apparently contradictory results
can be attributed to a variable Ne I/O I ratio due to changes in
the ionisation structure in the measured absorption columns that, however, are
not understood (Juett & Chakrabarty 2005). Hence, this means that the measured
Ne/O ratio does not reflect the donor composition. Although the orbital periods
of these systems are below 80 min (or at least believed to be so small), this
does not necessarily mean that they contain C/O donors. For example, the X-ray
burst properties of the 4U 1820-30 (
= 11 min) suggest an He-WD donor in that ultracompact system (e.g. Strohmayer & Brown 2002). Another
example are the relatively short X-ray bursts observed from (3) that even
suggest H-burning of material accreted onto the neutron star (NS) (Jonker et al. 2001), while optical spectra suggest a C/O WD donor (Nelemans et al. 2004 and this
work). A possible explanation is that the heavy elements (C, O, Ne) undergo
spallation during accretion leaving H and He nuclei for thermonuclear burning on
the NS (Bildsten et al. 1992). In a recent paper In't Zand et al. (2005)
conclude that the Ne/O ratio and the X-ray burst properties are all best
explained with an He-rich donor.
In conclusion, the only way to confirm that the four systems (2)-(5) - besides the prototype (1) - indeed contain C/O-rich donors, perhaps enriched with Ne, is by UV and/or optical spectroscopy. Nelemans et al. (2004) show that the optical spectra of (2)-(4) are devoid of hydrogen and helium emission lines, which are usually seen in H- or He-rich accreting systems. Their spectra exhibit low-ionisation C and O emission lines, most prominent in the brightest of these objects, 4U 0614+091. In this paper we present new optical spectra of this system that cover a larger wavelength interval. They confirm the earlier conclusion by Nelemans et al. (2004) that the emission lines arise from the C/O-dominated material that is probably located in the accretion disk and not in the X-ray heated face of the white-dwarf donor. We also present a first detailed optical spectrum of the prototype 4U 1626-67, which mainly shows weak oxygen emission lines and proves the H- and He-deficiency in this system, too. Nelemans & Jonker (2005) also performed VLT observations of this system (also in spring 2004) with a similar setup. They show a section of their spectrum and emphasize the similarity with 4U 0614+91.
Table 1: Log of the observations.
In addition we present results from first attempts to model the observed spectra with non-LTE accretion-disk models. We derive upper limits for the H and He abundances and investigate the formation of neon lines. Detailed C and O line-formation calculations can already qualitatively explain the observed emission lines. The ultimate goal is to derive detailed abundances and other disk parameters from the observed line profiles.
The paper is organised as follows. We describe our observations in the following section. We then present our line identifications in Sect. 3. Section 4 contains a description of our disk model calculations and we summarise and conclude in Sect. 5.
![]() |
Figure 1: Blue spectra of 4U 0614+091 and 4U 1626-67. H and He lines are lacking. The emission lines are identified as O II-III and C II-III (see also Table 2). All absorption features are of interstellar (i.s.) origin: Ca II 3934/3968 Å, the 4300 Å CH band, and two diffuse interstellar bands (DIB). Plotted over the two observed spectra is a synthetic accretion disk spectrum reddened with E(B-V) = 0.55 and E(B-V) = 0.40, respectively. |
Open with DEXTER |
![]() |
Figure 2: Red spectra; 4U 0614+091 shows emission lines from C II-IV and O II-III, while the 4U 1626-67 spectrum is virtually featureless. All absorption features are either telluric (marked by horizontal bars) or of interstellar origin (Na I 5890/5896 Å and a DIB at 5780 Å). The features at 5575 Å are artifacts ("x'') due to an O I sky emission line. Overplotted is the synthetic accretion disk spectrum like in Fig. 1. |
Open with DEXTER |
We obtained optical medium-resolution long-slit spectra of 4U 1626-67 and
4U 0614+091. The V magnitude of both systems is 18.5 (Ritter & Kolb 2003). We
used the FORS1 spectrograph attached to UT1 of ESO's Very Large Telescope
(VLT) on Paranal in Chile. The slit width was 1
.
We used two grisms (600B
and 600R, the latter in combination with order separation filter GG435) and
obtained spectra covering the regions 3600-6000 Å and 5400-7500 Å with a mean dispersion of 1.20 and 1.07 Å pix-1, respectively. Observations were
performed in service mode between Nov. 2003 and Mar. 2004. Each target was
observed at least twice, see Table 1 for details. For 4U 1626-67
the exposure time covers almost one orbital period so that any respective spectral
variation is smeared out. The data were
processed through the standard ESO reduction pipeline. In addition we performed
a flux calibration using observations of the DA white dwarf EG 274 with the same
instrumental setup. Since this flux standard was observed only once per
wavelength range and since the science targets were mostly observed at different
dates, this provides only a rough absolute flux calibration. For each object the
spectra were co-added to obtain one final spectrum. The blue and red spectra are
displayed in Figs. 1 and 2, respectively.
Archival spectra of 4U 1626-67 taken with the Hubble Space Telescope (HST) and the STIS instrument cover the complete UV/optical wavelength range (1150-10 000 Å). The FUV spectrum is described in detail by Homer et al. (2002) and a first comparison with synthetic accretion-disk spectra was presented by Werner et al. (2004). In Fig. 3 we compare the optical HST spectrum taken from the MAST archive with our VLT spectrum on an absolute flux scale. The blue and red sections of the VLT spectrum were scaled by a factor of 1.65 and 1.45, respectively, to normalise them to the HST flux at 5600 Å. It is obvious that the VLT spectrum is not only weaker but also flatter than the HST spectrum. This can be attributed either to flux calibration problems or to source variability.
![]() |
Figure 3:
Complete VLT spectrum of 4U 1626-67 compared to a spectrum taken with
HST. The VLT spectrum, scaled to the HST spectrum at 5600 Å, is flatter
possibly due to problems in absolute flux calibration. The HST spectrum has
a poorer resolution and S/N-ratio, although some emission features at
![]() |
Open with DEXTER |
Table 2: Emission features observed in the spectra of 4U 1626-67 and 4U 0614+091 and suggested line identifications. A plus sign in parentheses denotes an uncertain detection.
Emission lines from hydrogen and helium are completely lacking in the VLT spectra of both systems. We detect neither Balmer lines nor lines from He I (e.g. 4471 or 5876 Å) or He II (e.g. 4686 Å), which are usually seen in optical spectra of interacting binaries. In the following two sections we discuss the line spectra of both systems in detail. Line identification is performed "by eye'' without an equivalent-width estimation.
![]() |
Figure 4: Detail from three single spectra of 4U 0614+091. Thick line: spectrum taken on Nov. 03, 2003. Two thin lines: spectra taken consecutively on Dec. 15, 2003 (see Table 1). The strength of the two strongest emission features, the C II multiplets at 6580 Å and 7235 Å, show an increase over a time interval of several weeks. |
Open with DEXTER |
The spectrum of 4U 0614+091 shows numerous emission features that can be
assigned to ionized carbon and oxygen, namely C II-IV and O II-III. The detected emission features are listed in
Table 2 together with possible line identifications. Most
features are blends of lines from at least two ions, making interpretation
quite difficult, but some of them are probably unblended. These are the features
at 5590 Å (O III), 5650 Å (C II), 5695 Å (C III),
and 7235 Å (C II). The emission features at 4710 and 4940 Å are
perhaps exclusively due to O II, but they are blends of four and two multiplets from this ion, respectively. Nelemans et al. (2004) present
optical spectra of this star taken with VLT+FORS2 with slightly different
wavelength coverage (4620-8620 Å). We essentially confirm the detection of
their emission features and line identifications, however, our spectra extend to
shorter wavelengths, down to 3600 Å, which allows us to detect some additional
emission lines. On the other hand, we do not see four of the weak emission features
identified in the Nelemans et al. (2004) spectra (at
5140, 5190,
5280, and 6070 Å). This could be the consequence of line variability. In
Fig. 4 we compare our single spectra of 4U 0614+091 in detail, and it
is possible that the emission strength of the strongest C II lines is
variable on a time scale of weeks. In order to assess the significance of this
variability, we simulated several times with different noise the strongest emission features observed on Nov. 3 and Dec. 15, applying the S/N ratio of the VLT spectra. We find a 2
probability that the variability is real. The line widths correspond to a projected
rotational velocity of the order 1200 km s-1.
The spectrum of 4U 1626-67 clearly shows less features than does that of 4U 0614+091. Only five emission lines can be detected in the blue spectrum (see Table 2) and they are seen in both binaries. The red spectrum of 4U 1626-67 is virtually continuous. The strong C II emissions seen in 4U 0614+091 at 6580 Å and 7235 Å are lacking in 4U 1626-67. That could point to a higher ionisation of carbon in 4U 1626-67, although the C IV 5810 Å line is lacking, too. The presence of O II lines in the blue spectrum of 4U 1626-67 (e.g. at 4345 Å) also contradicts a higher ionisation, rather, it appears that the differences in the spectra of both binaries could be assigned to a different C/O ratio, being higher in 4U 0614+091. The entire set of emission lines in 4U 1626-67 could be due to oxygen alone, but in the UV spectrum C is clearly present (Homer et al. 2002). In addition, the line widths appear broader than in 4U 0614+091, however, this is difficult to quantify because of possible line blends and an uncertain identification of the continuum.
All absorption features seen in our optical spectra are either of interstellar origin or telluric. They are marked in Figs. 1 and 2.
We began the construction of accretion disk models to calculate synthetic
spectra and report here on the current state of our work. We use our accretion
disk code AcDc, which is described in detail by Nagel et al. (2004). In essence,
it assumes a radial -disk structure (Shakura & Sunyaev 1973). Then the
disk is divided into concentric annuli. For each annulus we solve the radiation
transfer equation (assuming plane-parallel geometry) together with the non-LTE rate equations for the atomic level populations, plus energy- and hydrostatic
equations, in order to calculate a detailed vertical structure. The integrated
disk spectrum is then obtained by co-adding the specific intensities from the
individual annuli, accounting for inclination and Keplerian rotation.
It is not the aim of this paper to present a detailed fit to the emission line spectra of the two binaries. This requires extensive parameter studies that are extremely time consuming. We rather choose to calculate synthetic spectra from selected disk regions in which we believe the physical properties are representative for the formation of the observed optical spectra. This should at least give a rough idea of the relative strength of lines from different ionisation stages and, thus, is primarily thought to put our line identifications on firm ground. At present we neglect irradiation of the disk by the neutron star, because it would introduce new free parameters. We expect that the radial ionisation structure of the disk will be shifted to larger radii when irradiation is taken into account, but we hope that the relative line strengths are not affected to the extent that our identifications become completely wrong. Although the studied systems are strong X-ray sources, we do not see recombination lines from highly ionised species in the optical spectra as might be expected. It is conceivable that the outer parts of the disks in which the optical spectrum arises is shielded from X-ray irradiation by an inflated inner-disk region.
Table 3: Summary of model atoms used in the disk model calculations. For each ion we list the number of NLTE levels and the number of line transitions. In brackets we give the number of lines in the wavelength range covered by our optical spectra after fine-structure splitting for detailed line-profile calculations.
Since we do not compute the spectrum of the entire disk, we cannot expect to match the overall observed continuum flux, and we have already pointed out that there are problems with the absolute flux calibration. We use interstellar reddening as a free parameter in order to roughly fit the model to the observed flux level. The applied reddening is given in the figure captions.
The synthetic spectra presented here are based on the following choice of
disk-model parameters. The central object is a neutron star with 10 km radius and a mass of 1.4 .
The mass accretion rate is 2
/yr. The disk extends from 1000 to 2000 stellar radii. The corresponding Keplerian velocities at these radii amount to 4300 and 3000 km s-1, and the effective temperatures to 47 000 K and 28 000 K,
respectively. The disk is divided into five annuli such that
decreases
almost linearly with radius steps. The inclination angle is set to 10
.
The
Reynolds number used to parametrise the disk viscosity was set to 10 000, which
corresponds to
.
The chemical composition is C = 50% and O = 50% (by mass).
For the opacity and emissivity calculations, it is essential to solve the non-LTE rate equations with detailed model atoms. Our main emphasis was put on the C II-III and O II-III ions. The number of non-LTE levels and radiative line transitions are summarised in Table 3. Level energies, oscillator strengths, and bound-free cross-sections for
photoionisation are taken from the Opacity Project (Seaton et al. 1994)
TOPbase. Electron
collisional rates for (de-) excitation, ionisation, and recombination are
computed with the usual approximate formulae. The neon model atom was taken from
Dreizler (1993). For the final spectrum synthesis, fine structure splitting of
atomic levels must be considered; level populations were taken from the models
with appropriate statistical weights. Level energies were obtained from the NIST
database. The optical synthetic spectra
finally include a total number of 377 lines from H I, He I-II,
C II-IV, O II-III, and Ne II. They essentially contain most
of the possible C and O line identifications given in Table 2,
plus many others that turn out to be too weak to be seen in the rotationally
broadened spectra.
As an example we present the vertical structure of the accretion disk model at a distance of 2000 stellar radii from the neutron star. The emergent flux at this
location corresponds to
= 28 000 K. Figure 5 shows the
run of several quantities above the disk midplane on a column mass scale m. The Rosseland optical depth reaches unity at
.
Here the
gravity amounts to
= 5.8. Mass and electron densities are 6
10-8 g cm-3 and 5
1015 cm-3, respectively. Hence the
physical conditions in the line-forming regions are comparable to those in the
atmosphere of a subluminous B star. The vertical distribution of ionic fractions
of C and O is shown in Fig. 6. There we also show the H, He, and Ne fractions that were taken from the test models that include these species in an amount of 10%. The dominant ionisation stages of oxygen at
are O II and O III, giving rise to prominent emission lines of these ions in the spectra. In the case of carbon, C III dominates,
closely followed by C II and C IV. This explains why we see lines
from three C ions. The dominant neon ions in the line-forming regions are
Ne II and Ne III. For helium we have He II dominating,
followed by He I, so we would expect strong lines from these ions if neon
and helium were abundant. Hydrogen is ionised by about 99.9%, but
still, prominent emission lines are predicted if H were an abundant species (see below).
![]() |
Figure 5:
Vertical structure of the C/O-disk model at a distance of 20 000 km from the neutron star. The emergent disk flux at this location corresponds to
![]() |
Open with DEXTER |
![]() |
Figure 6:
Vertical ionisation stratification of chemical elements in the disk models at
a distance of 20 000 km from the neutron star. The vertical line drawn at
![]() ![]() |
Open with DEXTER |
![]() |
Figure 7:
Comparison of models with zero and 10% hydrogen and helium
content to the observations. From the lack of H![]() |
Open with DEXTER |
The synthetic spectrum is plotted together with the observed spectra in Figs. 1 and 2. Let us first consider 4U 0614+091. Generally, many of the observed features are also seen in the model, although they do not match the strength. This basically corroborates our line identifications. The C II lines of the model are too weak (e.g. at 6580 Å and 7235 Å), while the C III emissions are in some cases too weak (e.g. 4650 Å) or in other cases too strong (4180 Å). The C IV 5810 Å line matches well. The oxygen lines show a similar behaviour. The O II lines of the model are too weak (e.g. at 4940 Å). Some lines of O III match reasonably well (e.g. at 3720 Å), while others do not (e.g. at 5590 Å, which is much too weak in the model).
The model comparison with 4U 1626-67 confirms our ideas from the first inspection of the spectrum. The lines in the blue region can be explained qualitatively by the mere presence of oxygen lines. This, and the missing C IV 5810 Å line might originate in a relatively low C/O ratio when compared to 4U 0614+091.
It is disappointing that our NLTE disk model obviously gives a poorer fit to the observed line spectrum than the simple isothermal, constant-density, LTE slab model presented by Nelemans et al. (2004). However, that slab model can at best fit a limited spectral range, as it is emitted from a particular emission region with an assumed value for temperature and density and will never be able to simultaneously fit the observed spectra from the X-ray to the UV and optical ranges. In good agreement with the observation, our model exhibits lines from three ionisation stages of carbon, indicating that temperature and density in the line-forming regions are reproduced reasonably well. The fact that the strength of many emission lines of a particular ion is either over- or underestimated may stem from drawbacks in the model atoms. One reason could be that the model ions are still too small and ignore the interlocking effects of neglected energy levels. A further extension of the model ions is hampered by the lack of atomic data, mostly oscillator strengths. Another reason for the poor line fits could be errors in the electron collisional rates. Only few are known from experiments or quantum mechanical calculations.
These problems also might affect our neon line-formation calculations. However,
the employed Ne model atom has been designed for the analysis of sdO stars and
gave successful fits to observed Ne II lines (Dreizler 1993). Therefore
we think that the predicted Ne emission lines in the model are more reliable
than the majority of the C and O lines.
Atomic data for H and He are accurately known and our model atoms have been
employed with success to analyse many classes of stellar spectra, hence,
we regard the computed line strengths of H and He as much more reliable than those
of the metal lines. Another reason for the poor line fits might be that the assumed
-disk does not describe the physics of the emitting region well.
The upper abundance limits we derive for H and He are regarded as realistic because of the above considerations, but the limit for neon is less secure.
![]() |
Figure 8: Comparison of models with zero and 10% neon content to the observations. The model including neon exhibits weak Ne II lines in the region 4350-4450 Å, which cannot be detected in the observed spectra. |
Open with DEXTER |
In addition to the pure C/O disk model we also computed a model including 10%
hydrogen and a model including 10% helium in order to see if this allows the H
and He abundances to be constrained. In Fig. 7 we display the
resulting emission lines of H
and He II 5412 Å. From the lack
of observed emission lines, it is evident that H and He are at most present at
the 10% level. This excludes the possibility that the disks of 4U 1626-67 and
4U 0614+091 are dominated by hydrogen or helium and confirms the suspicion that
they are in fact C/O-dominated.
We also computed a C/O model that includes 10% neon. It exhibits weak Ne II lines, but no such lines are detected in the VLT spectra (Fig. 8). We conclude that the neon abundance cannot be larger than .
We have presented new high-quality optical spectra of the ultracompact low-mass X-ray binaries 4U 0614+091 and 4U 1626-67. They are pure emission line spectra and most probably stem from the accretion disk. The spectral lines are identified as due to C II-IV and O II-III. Line identifications are corroborated by first results from modeling the disk spectra with detailed non-LTE radiation transfer calculations. Hydrogen and helium lines are lacking and our models confirm the deficiency of H and He in the disk. Hence, the donor stars in these systems are in fact the eroded cores of C/O white dwarfs. There are indications that the O/C ratio in 4U 0614+091 is higher than in 4U 1626-67. This could suggest that the stripping process of the WD in 4U 0614+091 is more advanced.
It is hard to estimate the systematic error of the derived upper limit for
the neon abundance so that the following conclusions are at the moment rather
uncertain. From the lack of Ne II lines we find that the Ne abundance
cannot exceed
or, in other words, the Ne/O ratio is at most of
the order 0.2. The much higher Ne/O ratios (
0.7) derived from the ISM
X-ray absorption edges of neutral Ne and O (Schulz et al. 2001; Juett et al. 2001)
would produce detectable Ne II lines in the disk spectra. This confirms
the conclusion of Juett & Chakrabarty (2005), that the determined ISM abundances of Ne and O are affected by ionisation effects and, hence, do not
reflect the abundances of the donor stars.
For an initial solar metallicity the 22Ne abundance in an Ne-enriched
crystallized and fractionated WD core can be estimated from theoretical models
to 0.07, but this value is very uncertain (Yungelson et al. 2002). It
could be even higher by a factor of 3 (Isern et al. 1991). Such an Ne-rich core
would have a mass of
(Yungelson et al. 2002). The mass
of the WD donor in 4U 1626-67 is much smaller (0.01
;
Yungelson et al. 2002). If we accept an upper limit of Ne = 0.1 for the observed neon abundance,
then we may draw the following conclusion. Either the WD core has crystallized
and fractionated, then our observation favors a relatively small Ne-enrichment
as a result of the crystallization process. Or, if one accepts that
fractionation would result in a high Ne abundance of 0.2 (unobserved), then the
WD core in 4U 1626-67 had no time to crystallize which, depending on various
details, lasts several Gyr (e.g. Hernanz et al. 1994).
Future work will concentrate on the disk modeling for the prototype 4U 1626-67, for which the observational database is the best of all such systems. Ultimately, the spectral properties (flux distribution and emission line strengths) over the complete wavelength range comprising spectral observations with Chandra, HST, and VLT, must be explained by a unique disk model.
Acknowledgements
T.R. is supported by the DLR under grant 50 OR 0201. We thank the referee for constructive criticism that helped to improve the paper.