A&A 458, 523-532 (2006)
DOI: 10.1051/0004-6361:20065241
M. Stejner - J. Madsen
Department of Physics and Astronomy, University of Aarhus, 8000 Århus C, Denmark
Received 21 March 2006 / Accepted 4 August 2006
Abstract
Context. Both superbursters and soft X-ray transients probe the process of deep crustal heating in compact stars. It was recently shown that the transfer of matter from crust to core in a strange star can heat the crust and ignite superbursts provided certain constraints on the strange quark matter equation of state are fulfilled.
Aims. We derive corresponding constraints on the equation of state for soft X-ray transients assuming their quiescent emission is powered in the same way, and further discuss the time dependence of this heating mechanism in transient systems.
Methods. We approach this using a simple parametrized model for deep crustal heating in strange stars assuming slow neutrino cooling in the core and blackbody photon emission from the surface.
Results. The constraints derived for hot frequently accreting soft X-ray transients are always consistent with those for superbursters. The colder sources are consistent for low values of the quark matter binding energy, heat conductivity and neutrino emissivity. The heating mechanism is very time dependent which may help explain cold sources with long recurrence times. Thus deep crustal heating in strange stars can provide a consistent explanation for superbursters and soft X-ray transients.
Key words: stars: neutron - dense matter - equation of state - X-rays: binaries - X-rays: bursts - X-rays: general
Strange stars are a class of compact stars made entirely of absolutely stable strange quark matter except for a possible thin crust of ordinary nuclear matter below the neutron drip density separated from the core by a strong Coulomb barrier (Baym & Chin 1976; Witten 1984; Alcock et al. 1986; Itoh 1970; Haensel et al. 1986). The stability of strange quark matter depends on poorly constrained strong interaction properties and remains to be decided by observation or experiment (see Weber et al. 2006; Madsen 1999; Weber 2005, for reviews). If true, the strange matter hypothesis would imply that all compact stars are likely to be strange stars (because the galactic flux of quark matter lumps from strange star collisions will trigger the hadron-to-quark transition in neutron stars and supernova cores (Caldwell & Friedman 1991; Madsen 1988)), and so strange star models should be able to reproduce observed characteristics of compact stars. Specifically strange stars should be able to reproduce the bursting behavior of low mass X-ray binaries (LMXBs) expected to contain some of the heaviest compact stars due to accretion, and here we test the consistency of the recent model for superburst ignition in strange stars by Page & Cumming (2005) with similar considerations for the quiescent emission from soft X-ray transients; two of the most important LMXB activities.
Superbursts, first discovered by Cornelisse et al. (2000), are a rare
kind of type I X-ray bursts from LMXBs distinguishing themselves by
extreme energies (1042 erg), long durations (4-14 h)
and long recurrence times (
1 year) - see Kuulkers (2004)
for a recent review of their observational properties. The mechanism
behind superbursts is thought to be unstable thermonuclear burning of
carbon in the crust at densities around
(Cumming & Bildsten 2001; Brown 2004; Cooper & Narayan 2005; Cumming & Macbeth 2004; Strohmayer & Brown 2002; Cumming et al. 2005; in't Zand et al. 2004). The
fuel for such an event is provided by the burning of accreted hydrogen
and helium at the surface, which may take place either stably or
unstably. In the latter case it gives rise to ordinary type I X-ray
bursts. The ashes from this process, containing about 10% carbon by
mass, are compressed under the increasing weight of the layers above,
and when enough carbon has been accumulated it may be ignited and burn
unstably. Observed superburst light curves require the ignition to
take place at temperatures around
K and column depths
of approximately
,
and as the heat released by
accretion and burning at the surface is quickly reradiated and does
not penetrate to this depth, these ignition conditions must be set by
the process of deep crustal heating and therefore probe the physics of
the deep crust and core (Cooper & Narayan 2005; Brown 2004).
In neutron stars deep crustal heating refers to the heat released as
the sinking ashes undergo electron captures and most importantly - at
densities above
- pycnonuclear reactions
(Haensel & Zdunik 2003; Brown et al. 1998; Haensel & Zdunik 1990). This amounts to about 1.4 MeV per nucleon and heats the
core to
108 K depending on the accretion rate. It has however
been difficult to achieve superburst ignition at the appropriate
column depth, and in a recent investigation Cumming et al. (2005) showed
that neutrino emission from Cooper pairing of neutrons at densities in
the range
limits the crust
temperature below
K. Thus unstable carbon ignition
seems impossible at the observed depth unless Cooper pair emission is
less efficient than currently thought or some additional heating
source can be found.
In strange stars there can be no pycnonuclear reactions because the crust density is limited by the neutron drip density, but for the same reason there can be no Cooper pair emission from neutrons either. Furthermore the energy released by matter tunneling through the Coulomb barrier at the base of the crust and converting to quark matter could be much greater than that released by pycnonuclear reactions in neutron stars. Thus it seems reasonable to expect that strange stars in LMXBs can be heated to higher temperatures than neutron stars, and Page & Cumming (2005) recently showed that this can lead to superburst ignition at the observed column depths if certain constraints on the strange quark matter core parameters are fulfilled.
Soft X-ray transients are LMXBs which undergo accretion outbursts
lasting days to years separated by long periods of quiescence in which
no or very little accretion takes place and the X-ray luminosity drops to
(Tomsick et al. 2004; Campana et al. 1998; Rutledge et al. 1999). When accretion stops,
the heat deposited at the surface quickly radiates away, and so in the
absence of any residual accretion the surface temperature should
reflect the temperature of the core including the effects of deep
crustal heating. This has been thoroughly investigated for neutron
stars in the literature (Yakovlev et al. 2004; Colpi et al. 2001; Brown et al. 1998; Ushomirsky & Rutledge 2001; Rutledge et al. 2002a; Brown et al. 2002; Yakovlev et al. 2003, and references
therein)
establishing deep crustal heating as an unavoidable heating mechanism
in the core, which should at least provide a "rock bottom''
luminosity for these sources. Residual accretion may still play a role
however - indeed the power law components in the spectra of many of
these sources are often taken as a sign of residual accretion or shock
emission from a rotation powered pulsar - and furthermore the
properties of the core remains uncertain, so the precise nature of the
quiescent emission is not firmly established.
Considering the many possible phases and poorly constrained properties of strange quark matter one would expect - for the same general reasons they were able to fit superbursters - that a wide range of strange star models can reproduce the observed properties of soft X-ray transients. Our approach is therefore to use a simple parametrized model for deep crustal heating in strange stars to discuss the consistency between constraints on strange quark matter derived from superbursters and soft X-ray transients, as such a mutual fit is likely to constrain the strange quark matter equation of state further.
We describe the overall structure of strange stars and their crusts in Sect. 2 and then turn to deep crustal heating in Sect. 3. In Sect. 4 we compare models with different assumptions about the cooling mechanism to observations of soft X-ray transients, and we discuss the importance of time dependence in Sect. 5 before concluding in Sect. 6.
Since the crust has no significant heating sources of its own
(electron captures in the crust amount to less than 40 keV per
nucleon (Haensel & Zdunik 2003)) its thermal profile is set by the
diffusion of heat from the core towards the surface. The surface
temperature, ,
is therefore only a function of the core
surface temperature,
,
and thus probes the process of
deep crustal heating. A detailed account for such
relations can be found in
e.g. Yakovlev et al. (2004), but in a reasonable approximation
(Gudmundsson et al. 1983), which
will be sufficient for our purpose. Here
K and
K is actually the temperature at a
density of
,
but as the crust is largely
isothermal at this depth, one can take this to be the core surface
temperature (Yakovlev et al. 2004; Gudmundsson et al. 1983).
In the discussion above and in what follows we have assumed the conventional description of strange stars with a sharp core surface, an electron filled gap and a nuclear crust. It should be mentioned however that Jaikumar et al. (2006) recently showed, that if the surface tension of strange quark matter is not too high, it may be unstable to phase separation at low pressure leading to a crust of quark nuggets embedded in a degenerate electron gas, which would in many ways resemble the nuclear crust on neutron stars. The consequences of such a crust in an accreting system are not clear at present, and we restrict ourselves to the traditional scenario.
Deep crustal heating in strange stars was recently studied by
Page & Cumming (2005) as a means of achieving superburst ignition conditions
at the appropriate column depth and thus constrain the properties of
strange quark matter (see also the work on continuously accreting
strange stars by Miralda-Escude et al. 1990). It was shown that observed
ignition conditions are achieved for a wide range of crust density,
,
strange quark matter binding energy,
,
neutrino
emissivity,
,
and thermal conductivity, K, provided
that the direct Urca process does not operate in the core implying
that strange quark matter should be a color
superconductor. Furthermore a set of useful analytical approximations
to the core surface temperature were shown to reasonably approximate
numerical calculations. Here we will use those same approximations to
find the corresponding constraints on strange quark matter in soft
X-ray transients and a relationship between the observable parameters
of superbursters and soft X-ray transients. For clarity and to
emphasize the approximations we therefore briefly review the
calculation by Page & Cumming (2005) in the following subsection.
As accretion on the strange star drives the crust towards its maximum
density, the rate at which nuclei (or neutrons if the neutron drip
density can be reached) tunnel through the Coulomb barrier into the
quark matter core will increase until an equilibrium crust density is
reached. This density may vary as a function of the unknown Coulomb
barrier height at the quark matter surface but this would not change
the qualitative features of the model and for our purpose we will
assume that the equilibrium crust density is
.
As the nuclei cross into the strange quark matter core at
the rate
(in units of nucleons cm-2 s-1unless otherwise stated) and are dissolved, they release an energy,
MeV per nucleon (the actual value depends on
uncertain values for strong interaction parameters but is fixed if
these parameters were known). This heats the core until
an equilibrium or steady state temperature is reached as the heating
balances energy loss by neutrino emission in the core. This is assumed
to take place at a (slow) rate of
We will - at first - assume that this equilibrium is simply determined by the average accretion rate, which would be reasonable in systems where the tunneling is more or less continuous. We return to this point in Sect. 5 however as this steady state assumption may not be sufficient for all systems.
Following Page & Cumming (2005) we model the thermal profile
of the core in a plane parallel approximation with heat transport
taking place by conduction and with a heat sink from neutrino
emission. The inward directed heat flux from the conversion of nucleons
to strange quark matter, F, is determined by
![]() |
(4) |
![]() |
(5) |
If however the thermal conductivity of strange quark matter is sufficiently
large the core becomes isothermal with the temperature given by
so
It should be noted that although Eqs. (6) and (7) are certainly oversimplified they were shown by Page & Cumming (2005) to agree reasonably well with numerical calculations of a full general relativistic model. Furthermore they offer a physical transparency, which given the many poorly constrained properties of strange quark matter, is perhaps just as desirable in what follows as full numerical models. Specifically Eqs. (6) and (7) are general enough not to depend on the properties of a particular strange quark matter phase - except for the assumption of slow neutrino emissivity necessary to obtain superburst ignition. They are in fact general enough to also describe neutron stars at the same level of detail, and indeed Eq. (7) corresponds entirely to case (ii.a) in Yakovlev et al. (2003) for low-mass neutron stars with slow neutrino emission. Thus much of what follows would also apply for neutron stars - although in that case it has of course been done in much greater detail.
In deriving Eqs. (6) and (7) we assumed that cooling by
neutrino emission dominated photon emission from the surface. This is
certainly the case in superbursters, where the core is very hot, but
for cold enough cores this may not be the case and one should take the
surface luminosity into account. Approximating the surface with a
blackbody we can estimate when this becomes important by setting,
,
which using
the Gudmundsson et al. (1983) relation,
,
gives a core temperature
![]() |
(8) |
![]() |
(9) |
In the following subsections we shall ignore photon cooling from the
surface, which allows analytical constraints on the core parameters
and simple scaling relations for the surface temperature independent
of assumptions about the core parameters. This approach is justified
for high
,
but we return to discuss the importance of
photon cooling in Sect. 4 as we shall find neutrino cooling unable to
account for the coldest soft X-ray transients.
Assuming that soft X-ray transients are
also powered by deep crustal heating, a
different independent constraint on the quark core parameters
can be obtained by applying Eqs. (6) and (7)
to observations of the accretion rate and surface temperature
of soft X-ray transients using a
relation such as the one in Yakovlev et al. (2004)
or Gudmundsson et al. (1983).
Aql X-1 is an example of a well studied
soft X-ray transient which has
,
(Rutledge et al. 2002a; Yakovlev et al. 2003) and therefore
using a
canonical mass-radius relationship,
,
M/R=0.2, and
.
The relations between the quark matter core parameters
thus obtained for superbursts and for the
soft X-ray transient Aql X-1 are compared in
Fig. 1. If we keep
fixed and increase
then K must also increase to produce the required
temperature until the critical conductivity is reached and the star
must be isothermal. Therefore only the regions below the isothermal curves
contain possible solutions for each class of object
and there can be no common solution for superbursters and
Aql X-1 between the two isothermal curves shown.
It is clear that the range of quark matter parameters permitted for
the soft X-ray transient Aql X-1 is significantly smaller
than that allowed by the superbursters; in particular the quark matter
binding energy per baryon,
,
necessary to fit the soft
X-ray transient must be in the lower part of the range allowed for fitting
superbursters, other parameters being equal.
We further note that the
non-isothermal relations from Eq. (6) agree remarkably well
in the range where common parameters fit superbursters and
Aql X-1, whereas no simultaneous isothermal solutions for the
two classes of objects seem to exist.
![]() |
Figure 1:
Constraints on quark core parameters obtained from
superbursters (solid)
and the soft X-ray transient Aql X-1 (dashed) assuming slow neutrino
cooling for both.
![]() ![]() ![]() ![]() |
Open with DEXTER |
For a fixed (albeit unknown) set of strong interaction parameters
neither the quark matter binding energy
,
nor the
quark core surface density should vary between stars, because both are
given directly by the strong interaction parameters.
The neutrino luminosity is dominated by the core surface contribution
and strange stars are furthermore characterized by an almost constant
density in the core unless we are very close to the maximum mass strange
star, so the slow
variations with density of
and K should be unimportant.
K may be inversely proportional to the
temperature (Haensel 1991) and can therefore vary by a small
factor between the coldest soft X-ray transients with quark core temperature
108 K and superbursters with core temperature
K.
If no density thresholds between phases or emission
mechanisms are crossed when going from superbursters to soft X-ray
transients it
should therefore be reasonable to expect approximate agreement between
the relations in Fig. 1, as seen for the non-isothermal
solutions, but not order of magnitude
differences as found in the isothermal case.
To the extent that
and
do not vary between stars we can utilize
Eqs. (6) and (7) to
compare the other parameters for two arbitrary strange stars with
deep crustal heating. In the following we shall in particular
be interested in cases where one
of the stars is a superburster (SB) and the
other a soft X-ray transient (X).
From Eq. (6) we have in the non-isothermal case
We will use Eqs. (10) and (11) to
compare superbursters and soft X-ray transients so writing
for the left hand side scaled to fit the superburst
data and dropping the subscript on the soft X-ray transient parameters we get
Table 1:
Properties of the observational sample commonly used to test
deep crustal heating. The temperature quoted for Aql X-1 is the
redshifted
.
We use
with the
canonical M/R=0.2. For Cen X-4 we quote the temperatures from both
Rutledge et al. (2001) and Campana et al. (2004).
![]() |
Figure 2:
Left: comparing cooling by a non-isothermal (dashed and
dot-dashed; cf. Eq. (12)) and isothermal (dashed-doubledot
and solid; cf. Eq. (14))
core to the observations in
Table 1. For Cen X-4 the quoted temperature intervals
overlap. Right: the "mixed'' case with non-isothermal superbursts
and isothermal soft X-ray transients (dashed and solid; cf. Eq. (16)) and
photon cooling (dotted) assuming a
pure blackbody with
![]() |
Open with DEXTER |
The observational sample most commonly used to test relations for deep
crustal heating such as Eqs. (14) and (12)
consists of soft X-ray transients in the field with short outburst
durations and recurrence times, for which X-ray observations in
quiescense have determined or limited the surface temperature. The
outburst duration should be short so the crust has not been
significantly heated by the release of nuclear energy
from the accreted matter and the surface temperature is therefore set
by deep crustal heating. The recurrence time should be short enough
that at least a few outbursts have been observed allowing a reasonable
estimate of the time averaged accretion rate - although one should
remember that the last 37 years of X-ray observations are not
necessarily representative of the 104 years of accretion
history needed to reach thermal equilibrium. Similarly since the
accretion history for stars in globular clusters is uncertain one
should only use systems in the field.
Useful lists of such systems and their properties including
references can be found in e.g. Yakovlev et al. (2003) and Tomsick et al. (2004). We
will take this data, summarized in Table 1, as our
observational basis and briefly discuss each system below.
The best examples of such systems are Aql X-1, Cen X-4 and 4U 1608-522
being short outburst field soft X-ray transients with detected thermal
components. Tomsick et al. (2004) discuss the outburst history of these
sources and we cite their accretion rates which are based on the
review of soft X-ray transient activity by Chen et al. (1997) using
For SAX J1808.4-3658, SAX J1810.8-2609 and XTE J2123-058 no thermal components have been detected and so we only cite upper limits to the surface temperature. SAX J1810.8-2609 and XTE J2123-058 have each shown only one outburst so the accretion rate is actually an upper limit. The accretion rate cited for SAX J1808.4-3658 is based on the two outbursts in 1996 and 1998, but further outburst have ocurred in 2000, 2002 and 2005 (see Wijnands 2004, for a review). These systems are therefore not expected to agree well with the theory, but at least the upper limits to their surface temperatures should not lie significantly below the prediction of deep crustal heating if they have accreted for long enough to have reached thermal equilibrium.
XTE J1709-267 lies 9-10 arcmin from NGC 6293 which has a tidal radius of 14.2 arcmin and so may be associated with this globular cluster although distance measurements are inconclusive on this (Jonker et al. 2004). It has shown two outbursts since 1995 and has a varying power law component.
As should be clear from the discussion above the data is rather
uncertain and should be treated with some caution. In particular the
time averaged accretion rates are hardly more than estimates, and the
surface temperatures, when they can be found, are uncertain and may
not only reflect deep crustal heating. Nevertheless deep crustal
heating should provide a minimum temperature for systems that have
accreted long enough and so we should at least test for consistency
with the observations if not perfect agreement. In Fig. 2 we
plot the data from Table 1 and compare with photon cooling
curves as well as Eqs. (12), (14)
and (16), that assume neutrino cooling,
for values of
corresponding to
and
.
The overall agreement with our predictions from neutrino emission in
superbursts is not impressive. Except Aql X-1 and 4U
1608-522 the observations fall significantly below the predicted
temperature for neutrino cooling, and actually pure photon
cooling does better if
MeV although it should not be dominant for
the high accretion rates in Aql X-1 and 4U 1608-522. This
raises two issues.
First it is noteworthy
that the two sources which fall close to the superburst prediction
without photon cooling,
are the hottest and most frequently accreting sources thus resembling
superbursters the most. In turn this raises the question whether the
accretion history may lead to heating with a pronounced time
dependence and whether the time averaged accretion rate is then still
a good predictor for the temperature - we discuss this possibility in
the next section.
![]() |
Figure 3:
Left: comparing combined photon and neutrino cooling of
an isothermal core (solid; cf. Eq. (22)) to the
observations in Table 1 with
![]() ![]() |
Open with DEXTER |
Second we see that if
is sufficiently low, photon
emission from the surface may be very important and so we should
combine the two cooling mechanisms. Thus we now demand - similarly to
the procedure in Yakovlev et al. (2003) -
where
is the
luminosity required to carry away the energy released by deep crustal
heating and
is the usual
blackbody luminosity. Following the calculation that leads to
Eq. (6) the non-isothermal neutrino luminosity is
Brown et al. (1998) similarly found that deep crustal heating of
neutron stars could explain the quiescent emission from Cen X-4 if only 0.1 MeV of the energy released by pycnonuclear
reactions during an outburst was conducted into the
core and deposited there. Yakovlev et al. (2004,2003) instead interpreted Cen X-4 and
SAX J1808.4-3658 as massive neutron stars with enhanced
neutrino emission. In the context of strange stars the energy is
deposited directly in the core, but no enhanced neutrino cooling is
required if the quark matter binding energy
is small.
To fully demonstrate either of these points would require a detailed time dependent model for cooling strange stars, which is beyond our scope here. Thus the calculations below are concerned only with the time dependence of the tunneling rate - and hence the heating - following an outburst, but we note that with respect to its consequences for the surface temperature we may be guided by similar considerations for neutron stars.
In neutron stars heat is released during and after an accretion event at the depth and on the time scale (months for the pycnonuclear reactions) corresponding to the individual electron captures, neutron emissions and pycnonuclear reactions. It then takes about 1 yr for the heat from the pycnonuclear reactions to diffuse to the surface (Ushomirsky & Rutledge 2001; Brown et al. 1998) and since it spreads out along the way, it is seen there as a heat wave lasting approximately the thermal diffusion time to the surface. Ushomirsky & Rutledge (2001) investigated the time dependence of this process for 1-day long outbursts, and found that energy deposited at depths, where the diffusion time to the surface was shorter than the outburst recurrence time, could produce significant variations in the surface luminosity following a pattern in time corresponding to the diffusion time for each nuclear reaction. If for instance the recurrence time was 1 year and assuming a low thermal conductivity in the crust, the luminosity variations were a few percent of the average, and only one heat wave could be seen because the diffusion time from the depth of the pycnonuclear reactions was on the order of the recurrence time. If the recurrence time was 30 years the variations in the luminosity could be up to 30 percent for cores with standard cooling and a factor of 3-4 for cores with accelerated cooling (e.g. pion condensates). In this case two separate heat waves were seen - one corresponding to electron captures in the outer crust visible a few days after the outburst and the other to the pycnonuclear reactions deeper in the crust emerging about a year later. For higher thermal conductivity variations were larger and happened on shorter time scales.
By analogy if the tunneling rate at the bottom of the crust - and hence the heating - is sufficiently time dependent following an accretion event, we would expect the qualitative conclusions in Ushomirsky & Rutledge (2001) to apply, possibly leading to strong variations in the surface temperature following an outburst as the resulting heat waves reach the surface. One difference may be though that in strange stars the heat is released at a lower density in the crust with a shorter diffusion time to the surface so the two heat waves may be closer in time.
The transfer of mass from crust to core by tunneling takes place at a
rate of (Stejner & Madsen 2005; Alcock et al. 1986)
![]() |
Figure 4:
Plot of Eq. (24) assuming
![]() ![]() |
Open with DEXTER |
Equation (2) gives a relation between
and
so we can solve Eq. (23) to find the crust
mass and tunneling rate for given accretion scenarios. Assuming the
same parameters as in Fig. 4 we show a solution in
Fig. 5 with 10 day long accretion outbursts separated by
quiescent intervals of 30 years and an average accretion rate of
- roughly
corresponding to Cen X-4. We start the integration with a
crust mass sufficiently low that the tunneling rate is negligible and
let it build up until the mass transferred to the core during each
cycle matches the average accretion rate. At this point the crust mass
is at its maximum - as determined by the choice of
MeV - during the accretion outbursts and then relaxes to a state with
very little tunneling within about a year after each
outburst. Figure 6 shows the corresponding result for a system
resembling Aql X-1 accreting 10 in every 220 days with an
average accretion rate of
.
Here we more clearly see the crust reach a maximum mass
during the 10 day long outbursts.
![]() |
Figure 5:
A crude model of Cen X-4. Crust mass in a system accreting
with 10 day long accretion outbursts separated by quiescent intervals
of 30 years and an average accretion rate of
![]() |
Open with DEXTER |
![]() |
Figure 6:
A crude model of Aql X-1. Crust mass in a system accreting
with 10 day long accretion outbursts separated by quiescent intervals
of 220 days and an average accretion rate of
![]() |
Open with DEXTER |
Thus practically all the heating takes place during and immediately after the outburst when the transfer rate is essentially the accretion rate, and the heating mechanism is indeed extremely time dependent. In systems such as Aql X-1 and 4U 1608-522 the outburst recurrence time is comparable to the diffusion time to the surface however. The heat wave resulting from an accretion outburst will therefore be smeared out, and the crust temperature will not vary significantly between outbursts unless one assumes a high thermal conductivity in the crust, so a steady state model using the average accretion rate should be able to predict the temperature. Conversely in systems such as Cen X-4 where the recurrence time is significantly longer than the diffusion time to the surface one would expect larger temperature variations following an outburst.
Furthermore in such systems the core - or the surface layers of the
core in the non-isothermal case - may have time to cool significantly
between outbursts. In this case we would not expect
Eqs. (12), (14) and (16) to hold, but
if (speculatively) the tunneling rate late in the cycle is a better
predictor of the temperature at this point, Eq. (12)
assuming neutrino cooling only would actually give a surface
temperature around 76 eV for Cen X-4 as observed. The same
argument for XTE J1709-267 assuming a 6 year recurrence time
would lead to surface temperatures as low as 83 eV. Here the
recurrence time may not be quite long enough however, and furthermore
the temperature we cite was measured in May 2003 not very long after
the March-April 2002 outburst, so the argument may be too simple. The
temperature we cite for Cen X-4 was measured in 1997-2000
long after the 1979 outburst so this should provide a more clear-cut
case. For the temperature of the three remaining sources to reach
their upper bounds in this way with the cited accretion rates, the
outburst recurrence times would have to be very long (50, 100
and 800 years for SAX J1810.8-2609, XTE J2123-058
and SAX J1808.4-3658 respectively).
This shows how the
the time dependence of the heating mechanism may potententially help
explain the coldest soft X-ray transients, even if
is
too high for photon cooling to dominate. More work including
details on the cooling of the core and the transport properties of
both crust and core would be required though to show that it is
possible for the temperature to change this rapidly.
For systems with short recurrence times where the steady state
temperature may be expected to be a good approximation we found
reasonable agreement assuming neutrino cooling only - although the
range of quark matter parameters giving superbursters between the
isothermal curves in Fig. 1 is then excluded by soft X-ray
transients. Assuming a combination of photon and neutrino cooling for
soft X-ray transients, we
saw in Fig. 3 that if
MeV even
the coldest sources may be explained while still fulfilling superburst
ignition conditions at high accretion rates.
We further showed that the tunneling rate near the equilibrium crust
mass depends sensitively on the crust density and consequently that
essentially all the heating takes place during or immediately after an
outburst. This points to the need for a fully time dependent model to
fit systems with recurrence times significantly longer than the
thermal diffusion time to the surface, but if the low level tunneling
rate long after the outburst is a better predictor for such systems,
it would provide a fit for Cen X-4 assuming only neutrino
cooling as well. The upper temperature bounds on SAX J1810.8-2609, XTE J2123-058 and SAX J1808.4-3658
can not be reached for the cited accretion rates in this way unless
the recurrence times are very long. This is consistent with the
outburst histories of SAX J1810.8-2609 and XTE J2123-058 - one known outburst each - but does not compare well
with the 2 yr recurrence time for SAX J1808.4-3658. Similarly to the situation with neutron stars
(Yakovlev et al. 2004; Campana et al. 2002) some form of enhanced cooling would
thus seem to be required if SAX J1808.4-3658 has accreted for
long enough to heat the core to its equilibrium temperature and
is too high for photon cooling to
dominate. XTE J1709-267 has an intermediate recurrence time,
and its temperature was measured shortly after an outburst making
interpretation difficult.
If - as seems reasonable - a combination of the two cooling
mechanisms is required to explain the observations we should further
note that the constraints imposed by Aql X-1 in
Fig. 1 are not relevant since they were derived assuming only
neutrino cooling. The requirement that
MeV
can therefore be met - even if only just - for the quoted
range in
and K assuming the parameters fit superburst
ignition conditions in Fig. 1. Although such fine tuning in
the binding energy can not be ruled out it is rather conspicuous, but
given the simplicity of our model and that we have not dealt with the
uncertainties or the time dependence in detail such quantitative
conclusions should be treated with caution. We have attempted to test
for consistency, not to determine the quark matter parameters. This
would require more detailed modelling, possibly taking the individual
accretion history of each source into account, which may alleviate the
need to keep
small.
Thus even assuming a high
the model seems consistent
for the systems where it should be most credible and with the
exception of SAX J1808.4-3658 no clear inconsistency can be
demonstrated for the other sources studied here. However SAX
J1808.4-3658 is difficult to explain with neutron stars as well, and
its precise nature remains to be decided. At this point strange stars
thus seem to do no worse than neutron stars in explaining the
quiescent emission of soft X-ray transients regardless of the assumed
core parameters, and for a low
they explain even the
coldest sources including SAX J1808.4-3658. Because the two cooling
mechanisms are important in different accretion regimes it finally
seems likely that detailed modelling along with more accurate
temperature measurements at low accretion rates would constrain the
core parameters further - or perhaps demonstrate a clear
inconsistency.
Finally we should not omit a short discussion of the interesting quasi
persistent transients KS 1731-260 and MXB 1659-29 (Wijnands et al. 2004; Cackett et al. 2006; Wijnands et al. 2002; Rutledge et al. 2002b; Wijnands et al. 2001,2003).
These are systems with prolonged accretion outbursts - 12.5 and 2.5 years duration respectively - and a convincing case has been built
for the suggestion that their crusts are significantly heated during
these outbursts and brought out of thermal equilibrium with the
core. The thermal components in their spectrum therefore show the
cooling of the crust to thermal equilibrium and recently
Cackett et al. (2006) using new Chandra and XMM-Newton data
found that both systems are now close to equilibrium allowing the core
temperatures and thermal cooling time scales of the crusts to be
determined. Fitting the cooling curves with an exponential decay
towards a constant temperature presumably reflecting the core
temperature the constant was found to be
eV for KS 1731-260 and
eV for MXB 1659-29 at infinity. Using a
canonical mass-radius relationship this implies surface temperatures
of
eV and
eV respectively, so these
sources are both cold but not unusually so. Only upper limits to
the accretion rates can be established with any certainty however due
to the very limited number of observed outbursts. For KS 1731-260 (one
outburst in 1988-2001) this is
(Yakovlev et al. 2003) using the time-averaged
flux of Rutledge et al. (2002b). For MXB 1659-29 Wijnands et al. (2003)
find bolometric fluxes of
for the 1999-2001 outburst and
for the 1976-1979 outburst
which using the highest of these values in Eq. (18) (with
f=1) means an upper limit of
for a distance of 10 kpc. Referring to Fig. 3 both of these limits are very high
given the surface temperatures, and would require enhanced cooling for
both strange stars and neutron stars - but then the accretion rates
are hardly certain enough to back up such a claim. Perhaps more
interesting the e-folding times for the crust temperature relaxing
back to the core temperature were determined both sources. Both of
these were rather short (
days for KS 1731-260 and
days for MXB 1656-29) and by comparison with the neutron star cooling
curves computed by Rutledge et al. (2002b) for the accretion history of
KS 1731-260 may be taken as evidence for high thermal
conductivity in the crust. Strange stars are missing the deep crust
between 1011 and
and as we have seen
there is very little tunnelling in quiescence, so although it would
go too far to show this in detail here, it should be safe to assume that
they cool faster than neutron stars relaxing the need for high
thermal conductivity in the crust. This new development therefore
further underlines how rewarding a time dependent model for accreting
strange stars might be.
Acknowledgements
This work was supported by The Danish Natural Science Research Council.