# On the nature of the compact star in 4U 1728-34

## 1 Introduction

With the Rossi X-ray Timing Explorer kHz QPOs have been discovered in about 20 neutron star LMXBs (see van der Klis vk98 (1998) for a review). In many cases, two simultaneous kHz peaks ($`2001200`$ Hz) are observed in the power spectra of the X-ray count rate variations, with the separation frequency $`\mathrm{\Delta }\nu `$ roughly constant (e.g., Strohmayer et al. s96 (1996); Ford et al. f97 (1997); Wijnands et al. w97 (1997)). Sometimes a third kHz peak is detected in a few atoll sources during type I X-ray bursts at a frequency $`\nu _\mathrm{b}`$ equal to the separation frequency $`\mathrm{\Delta }\nu `$ of the two peaks or twice that (see Strohmayer, Swank, & Zhang ssz98 (1998) for a review). These two results suggest that a beat-frequency mechanism is at work, with the upper kHz peak at the Keplerian orbital frequency at the inner edge of the accretion disk around the neutron star, the third peak at the neutron star spin frequency (or twice that), and the lower kHz peak at the difference frequency between them (Strohmayer et al. s96 (1996); Miller, Lamb, & Psaltis mlp98 (1998)).

Within the beat-frequency model, the burst oscillations $`\nu _\mathrm{b}`$ and the kHz peak separations $`\mathrm{\Delta }\nu `$ are both believed to be close to the neutron star spin frequency and thus $`\mathrm{\Delta }\nu `$ should remain constant. However, recent observations have challenged this interpretation: in Sco X-1 (van der Klis et al. vk97 (1997)), 4U 1608-52 (Méndez et al. m98 (1998)), and 4U 1735-44 (Ford et al. f98 (1998)), $`\mathrm{\Delta }\nu `$ slowly decreases as the frequencies of both QPOs increase. Moreover, in 4U 1636-53 the frequeny $`\nu _\mathrm{b}`$ of the burst oscillations (or half this value) does not match the frequency difference $`\mathrm{\Delta }\nu `$ between the kHz QPOs (Méndez, van der Klis, & van Paradijs mkp98 (1998)). In the atoll source 4U 1728-34, Méndez & van der Klis (mk99 (1999)) find that $`\mathrm{\Delta }\nu `$ is always significantly smaller than $`\nu _\mathrm{b}`$, even at the lowest inferred mass accretion rate, when $`\mathrm{\Delta }\nu `$ seems to reach its maximum value, and that $`\mathrm{\Delta }\nu `$ decreases significantly as the frequency of the lower kHz QPO increases.

A different approach was recently invoked by Osherovich & Titarchuk (1999a ) and Titarchuk & Osherovich (to99 (1999)), who proposed a unified classification of kHz QPOs and the related low frequency phenomena. In this model, kHz QPOs are modeled as Keplerian oscillations under the influence of the Coriolis force in a rotating frame of reference (magnetosphere). The frequencies $`\nu _\mathrm{u}`$ of the upper kHz QPO branch hold a hybrid frequency relation with the Keplerian frequencies $`\nu _\mathrm{k}`$ (referred to the lower kHz QPO branch): $`\nu _\mathrm{u}^2=\nu _\mathrm{k}^2+(2\nu _\mathrm{m})^2`$, where $`\nu _\mathrm{m}`$ is the rotational frequency of the star’s magnetosphere. For Sco X-1, the QPOs with frequencies $`45`$ and 90 Hz are interpreted as the 1st and 2nd harmonics of the lower branch of the Keplerian oscillations (Osherovich & Titarchuk 1999a ). The same interpretation is applied to the $`35`$ Hz QPOs in the atoll source 4U 1702-42 (Osherovich & Titarchuk 1999b ). The observed low Lorentzian frequency in 4U 1728-34 is suggested to be associated with radial oscillations in the boundary layer of the disk, whereas the observed break frequency is determined by the characteristic diffusion time of the inward motion of the matter in the accretion flow (Titarchuk & Osherovich to99 (1999)). Predictions of this model regarding relations between the QPO frequencies mentioned above compare favorably with recent observations of Sco X-1, 4U 1608-52, 4U 1702-429, and 4U 1728-34.

The discovery of kHz QPO features allows to probe not only the interaction between the accretion disk and the stellar magnetosphere, but also the structure of the compact stars involved. In this Letter, based on the work of Osherovich & Titarchuk (1999a ) and Titarchuk & Osherovich (to99 (1999)), we study the mass - radius ($`MR`$) relation of the atoll X-ray source 4U 1728-34, and show that it is possibly a strange star.

## 2 The mass - radius relation for 4U 1728-34

Titarchuk & Osherovich (to99 (1999)) suggest a model for the radial oscillations and diffusion in the viscous boundary layer of the accretion disk in LMXBs. Their dimensional analysis has identified the corresponding frequencies which are consistent with the low Lorentzian and break frequencies for 4U 1728-34 observed by Ford & van der Klis (fk98 (1998)), the predicted values for the break frequency related to the diffusion in the boundary layer are also consistent with the observed ones for the same source. The presence of the break frequency and the correlated Lorentzian frequency suggests the introduction of a new scale in the phenomenon. One attractive feature of the model is the introduction of such a scale in the model through the Reynolds number. The best fit for the observed data was obtained when

$$a_\mathrm{k}=(M/M_{})(R_0/3R_\mathrm{s})^{3/2}(\nu /364\mathrm{Hz})=1.03,$$
(1)   
where $`M`$ is the stellar mass, $`R_0`$ is the inner edge of the accretion disk, $`R_\mathrm{s}=2GM/c^2`$ is the Schwarzschild radius, and $`\nu `$ is the spin frequency of the star. Given the 364 Hz spin frequency of 4U 1728-34 (Strohmayer et al. s96 (1996)), from Eq. (1) we derive the inner disk radius

$$R_09(a_\mathrm{k}/1.03)^{2/3}m^{1/3}\mathrm{km},$$
(2)   
where $`m=M/M_{}`$. Since the innermost radius of the disk must be larger than the radius of the star itself, this leads to a mass-dependent upper bound on the stellar radius,

$$R\stackrel{<}{}9(a_\mathrm{k}/1.03)^{2/3}m^{1/3}\mathrm{km},$$
(3)   
which is plotted in dot-dashed curve in Fig. 1. The dotted curve in Fig. 1 represents the Schwarzschild radius, which presents the lower limit of $`R`$.

Figure 1 compares Eq. (3) with the theoretical $`MR`$ relations (solid curves) for nonrotating neutron stars given by six recent realistic models for the equation of state (EOS) of dense matter. In models UU (Wiringa, Fiks, & Fabrocini wff88 (1988)), BBB1 and BBB2 (Baldo, Bombaci, & Burgio bbb97 (1997)), the neutron star core is assumed to be composed by an uncharged mixture of neutrons, protons, electrons and muons in equilibrium with respect to the weak interaction ($`\beta `$-stable nuclear matter). These models are based on microscopic calculations of asymmetric nuclear matter by use of realistic nuclear forces which fit experimental nucleon-nucleon scattering data, and deuteron properties. In model Hyp (Prakash et al. p97 (1997)), hyperons are considered in addition to nucleons as hadronic constituents of the neutron star core. The curve labeled BPAL12 is a very soft EOS for $`\beta `$-stable nuclear matter (Prakash et al. p97 (1997)), giving a lower limiting mass and a smaller radius with respect to a stiff EOS. Finally, we consider the possibility that neutron stars may possess a core with a Bose - Einstein condensate of negative kaons (e.g. Glendenning & Schaffner-Bielich gs99 (1999)). The main physical effect of the onset of K<sup>-</sup> condensation is a softening of the EOS with a consequent lowering of the neutron star maximum mass and possibly of the radius. This is shown with the curve labeled K<sup>-</sup>. It is clearly seen that only the upper part of EOSs UU, BBB1, BBB2, and BPAL12 are compatible with Eq. (3).

Another constraint on the mass and radius of 4U 1728-34 results from the requirement that the inner radius $`R_0`$ of the disk must be larger than the last stable circular orbit $`R_{\mathrm{ms}}`$ around the star, which is given by (Bardeen, Press, & Teukolsky bpt72 (1972)),

$$R_{\mathrm{ms}}=R_\mathrm{h}(3+Z_2[(3Z_1)(3+Z_1+2Z_2)]^{1/2}),$$
(4)   
where

$$Z_1=1+[1(a/R_\mathrm{h})^2]^{1/3}[(1+a/R_\mathrm{h})^{1/3}+(1a/R_\mathrm{h})^{1/3}],$$
(5)   
$$Z_2=[3(a/R_\mathrm{h})^2+Z_1^2]^{1/2},$$
(6)   
with $`R_\mathrm{h}=GM/c^2`$ and $`a=I\mathrm{\Omega }/Mc`$ (where $`I`$ is the moment of inertia and $`\mathrm{\Omega }=2\pi \nu `$ is the angular velocity of the star, respectively). For a given mass $`M`$, from Eq. (2) we can find $`R_0`$. From the inequality $`R_0`$$`\stackrel{>}{}`$$`R_{\mathrm{ms}}`$ and Eqs. (4-6), one can calculate the lower limit of $`a`$, and from there the lower limit of the moment of inertia $`I`$. To account for the inflation of the radius caused by rotation, we adopt the formula

$$I=0.21MR^2(1R_\mathrm{s}/R)^1,$$
(7)   
proposed by Burderi et al. (b99 (1999)), to obtain the lower limit on the stellar radius $`R`$ <sup>1</sup><sup>1</sup>1There are actually two solutions of Eq. (7) with given $`I`$ and $`M`$, one serves as the upper limit of $`R`$ and the other the lower limit of $`R`$. But the former solution gives an unphysically small radius, and is rejected.. The results are shown in the dashed curve in Fig. 1. Combining the constraints $`R\stackrel{<}{}R_0`$ and $`R_0`$$`\stackrel{>}{}`$$`R_{\mathrm{ms}}`$ reveals that the allowed distribution of the mass and radius of 4U 1728-34 must lie between the dot-dashed curve and the dashed curve. None of the neutron star EOSs extends to this region. Thus it seems to rule out the possibility that 4U 1728-34 is a neutron star! Including rotational effects will shift the theoretical neutron star $`MR`$ curves to up-right in Fig. 1, making the contrast between the theoretical neutron star models and the derived mass and radius range of 4U 1728-34 even worse.

The difficulty of the neutron star model for 4U 1728-34 suggests that it may belong to another type of compact objects, that is, strange stars, which are entirely made of deconfined u,d,s quark matter (strange matter). The possible existence of strange stars is a direct consequence of the conjecture (Bodmer b71 (1971); Witten w84 (1984)) that strange matter may be the absolute ground state of strongly interacting matter. Detailed studies have shown that the existence of strange matter is allowable within uncertainties inherent in a strong interaction calculation (Farhi & Jaffe fj84 (1984)); thus strange stars may exist in the universe.

Most of the previous calculations (e.g. Alcock, Farhi, & Olinto afo86 (1986); Haensel, Zdunik, & Schaeffer hzs86 (1996)) of strange star properties used an EOS for strange matter based on the phenomenological nucleonic bag model, in which the basic features of quantum chromodynamics, such as quark confinement and asymptotic freedom are postulated from the beginning, though the deconfinement of quarks at high density is not obvious in the bag model. Recently, Dey et al. (d98 (1998)) derived an EOS for strange matter, which has asymptotic freedom built in, shows confinement at zero baryon density, deconfinement at high density, and gives a stable configuration for chargeless, $`\beta `$-stable strange matter. In this model the quark interaction is described by an interquark vector potential originating from gluon exchange, and by a density dependent scalar potential which restores the chiral symmetry at high density. Using the same model (but different values of the parameters with respect to those employed in Dey et al. d98 (1998)) we calculated the $`MR`$ relations, which are also shown in solid curves labeled ss1 and ss2 in Fig. 1, corresponding to strange stars with maximum masses of $`1.44M_{}`$ and $`1.32M_{}`$ and radii of 7.07 km and 6.53 km, respectively. It is seen that the region confined by the dashed and dot-dashed curves for 4U 1728-34 in Fig. 1 is in remarkable accord with the strange star models with masses $`\stackrel{<}{}1.1M_{}`$. Figure 1 clearly demonstrates that a strange star model is more compatible with 4U 1728-34 than a neutron star one. This seems to be important since the model of Dey et al. (d98 (1998)) employs quark masses which decrease with increasing density, a manifestation of chiral symmetry restoration (CSR), and the parametrization of CSR can be deduced from the mass and radius of the strange star fitted.

The lower limit on the moment of inertia deduced above does not depend on the EOS. In particular Eq. (7) is an extrapolation from neutron stars and it may be doubted whether it is valid for strange stars. We have checked that the standard calculated moment of inertia of the strange star of mass $`M1.1M_{}`$ agrees with Eq. (7) with an error of about $`20\%`$, and the numerical factor 0.21 can change to 0.28 (ss1) or 0.26 (ss2) depending on the model parameters. The moment of inertia calculation of the strange star with variable density layers may not in fact satisfy the assumptions of the standard calculation. It is satisfying to find that an EOS-independent estimate can already reproduce the inertia to a reasonable accuracy.

## 3 Discussion and conclusion

We have shown that strange stars are more consistent with the properties of 4U 1728-34 compared to neutron stars, if the model of Osherovich & Titarchuk (1999a ) and Titarchuk & Osherovich (to99 (1999)) correctly interprets the high- and low-frequency QPO phenomena in 4U 1728-34. Since our analysis does not involve any unknown physical factors such as the magnetic field strength and topology, radiation processes, and the mass accretion rate, it seems to provide the first definite evidence for the existence of strange stars. It remains to see whether other compact stars in LMXBs would be identified as strange stars from observations of their QPO phenomena. If so, there will be very deep consequences for both physics of strong interactions and astrophysics.

Of course our conclusion is subject to the uncertainty in the parameter $`a_\mathrm{k}`$, which is determined by the best fitting of the observed data. Increase or decrease in $`a_\mathrm{k}`$ will enlarge or reduce the allowed region of the mass and radius of 4U 1728-34 in Fig. 1. For example, if $`a_\mathrm{k}=2`$, nearly all neutron star EOSs are compatible with 4U 1728-34. However, according to Titarchuk & Osherovich (to99 (1999)), the dependence of $`\chi ^2`$ on $`a_\mathrm{k}`$ is quite strong: $`\chi ^2=3802473076a_\mathrm{k}+35732a_\mathrm{k}^235732[(a_\mathrm{k}1.03)^2+0.019]`$. A $`15\%`$ variation of $`a_\mathrm{k}`$ from the preferred value 1.03 will double the value of $`\chi ^2`$, greatly reducing the confidence level of the fit. Thus we expect that the possible change in $`a_\mathrm{k}`$ and in the derived $`MR`$ relation is rather small.

If 4U 1728-34 is a strange star, Fig. 1 reveals that it possesses a mass $`\stackrel{<}{}1.1M_{}`$, which is quite small compared to the $`1.4M_{}`$ masses of standard neutron stars like PSR 1913+16, though still barely consistent with the estimated masses of quite a few radio pulsars within $`68\%`$ confidence limit (see Thorsett & Chakrabarty tc99 (1999)). Since the mass accretion rates in LMXB atoll sources are lower than the Eddington limit accretion rate by $`12`$ orders of magnitude, it is conventionally thought that in such systems almost all of the matter transferred from the companion stars should be accreted by the compact stars. The light mass of 4U 1728-34, however, seems to support the conclusion of Thorsett & Chakrabarty (tc99 (1999)) that there is no evidence for extensive mass accretion occurred in LMXBs, implying that even the simplest evolution of LMXBs is far from being fully understood. However, formation of strange stars may be different from that of neutron stars (i.e. type II/Ib supernova explosions) concerning the glitching phenomena observed in some radio pulsars (e.g. Alpar a87 (1987)), suggesting that some exotic formation channels may be at work. Accreting massive neutron stars in binary systems are not favorable progenitor candidates of strange stars, at least in the case of 4U 1728-34, since it requires too large mass ($`\stackrel{>}{}0.7M_{}`$) ejected during the conversion as the neutron star masses exceed $`\stackrel{>}{}1.8M_{}`$. Accretion-induced collapse of white dwarfs remains an open possibility of strange star formation, and could be responsible for the light mass of 4U 1728-34. The detailed investigation of this subject is beyond the scope of this paper and wil be discussed elsewhere.

We are grateful to Lev Titarchuk for helpful and prompt e-mails. X. L. was supported by National Natural Science Foundation of China. S. R., M. D. and J. D. were supported in part by DST grant no. SP/S2/K18/96, Govt. of India.