Source author record

J. L. Zdunik

J. L. Zdunik appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

Catalog footprint

What is connected

28works
5topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

28 published item(s)

preprint2021arXiv

Magnetic fields in late-stage proto-neutron stars

We explore the thermal and magnetic-field structure of a late-stage proto-neutron star. We find the dominant contribution to the entropy in different regions of the star, from which we build a simplified equation of state for the hot neutron star. With this, we numerically solve the stellar equilibrium equations to find a range of models, including magnetic fields and rotation up to Keplerian velocity. We approximate the equation of state as a barotrope, and discuss the validity of this assumption. For fixed magnetic-field strength, the induced ellipticity increases with temperature; we give quantitative formulae for this. The Keplerian velocity is considerably lower for hotter stars, which may set a de-facto maximum rotation rate for non-recycled NSs well below 1 kHz. Magnetic fields stronger than around $10^{14}$ G have qualitatively similar equilibrium states in both hot and cold neutron stars, with large-scale simple structure and the poloidal field component dominating over the toroidal one; we argue this result may be universal. We show that truncating magnetic-field solutions at low multipoles leads to serious inaccuracies, especially for models with rapid rotation or a strong toroidal-field component.

preprint2020arXiv

Redshift factor and diffusive equilibrium of unbound neutrons in the single nucleus model of accreting neutron star crust

Using a Wigner-Seitz approximation with spherical cells, we re-analyze a widely used single nucleus model of accreting neutron star crust. We calculate beta disequilibrium within the crust, which is sizable, and implies that neutron and baryon chemical potentials, $μ_n$ and $μ_{\rm b}$, are not equal. We include also non-equilibrium reactions, driven by matter compression, and proceeding in the reaction layers. The constancy of $e^Φμ_n$, where the spacetime metric component $g_{00}=e^{2Φ}$, in the shells between the reaction layers is not applicable, because single electron captures are blocked, so that the neutron fraction is fixed, and therefore neutrons are not an independent component of the crust matter. The absence of neutron diffusion in the shells between the reaction layers, stems from the constancy of the neutron fraction (concentration) in these shells. In the reaction layers, the outward force resulting from neutron fraction gradient is balanced by the inward gravitational force acting on unbound neutrons. Neglecting the thickness of the reaction layers compared to the shell thickness, we obtain condition $e^{Φ(r)}f_Q(r)g(r)=$constant, where $g$ is Gibbs energy per nucleon, undergoing discontinuous drops on the reaction surfaces, and $f_Q(r)g(r)={\widetilde{g}}(r)$ is a continuous function, due to the factor $f_Q(r)$ canceling the discontinuities (drops) in $g(r)$. The function $f_Q(r)$ is calculated using the Tolman-Oppenheimer-Volkov equations from $f_Q(P)$ and $g(P)$ obtained from the equation of state (EOS) with discontinuites. The constancy of of $e^{Φ(r)}{\widetilde{g}}(r)$ is an extension of the standard relation $e^{Φ(r)}{μ}_{\rm b}(r)=$constant, valid in hydrostatic equilibrium for catalyzed crust.

preprint2016arXiv

Analytic approximations of neutrino luminosities and heat capacities of neutron stars with nucleon cores

We derive analytic approximations of neutrino luminosities and heat capacities of neutron stars with nucleon cores valid for a wide class of equations of state of dense nucleon matter. The neutrino luminosities are approximated for the three cases in which they are produced by (i) direct Urca or (ii) modified Urca processes in non-superfluid matter, or (iii) neutrino-pair bremsstrahlung in neutron-neutron collisions (when other neutrino reactions are suppressed by strong proton superfluidity). The heat capacity is approximated for the two cases of (i) non-superfluid cores and (ii) the cores with strong proton superfluidity. The results can greatly simplify numerical simulations of cooling neutron stars with isothermal interiors at the neutrino and photon cooling stages as well as simulations of quasi-stationary internal thermal states of neutron stars in X-ray transients. For illustration, a model-independent analysis of thermal states of the latter sources is outlined.

preprint2016arXiv

Cooling of neutron stars with diffusive envelopes

We study the effects of heat blanketing envelopes of neutron stars on their cooling. To this aim, we perform cooling simulations using newly constructed models of the envelopes composed of binary ion mixtures (H--He, He--C, C--Fe) varying the mass of lighter ions (H, He or C) in the envelope. The results are compared with those calculated using the standard models of the envelopes which contain the layers of lighter (accreted) elements (H, He and C) on top of the Fe layer, varying the mass of accreted elements. The main effect is that the chemical composition of the envelopes influences their thermal conductivity and, hence, thermal insulation of the star. For illustration, we apply these results to estimate the internal temperature of the Vela pulsar and to study the cooling of neutron stars of ages of 0.1 - 1 Myr at the photon cooling stage. The uncertainties of the cooling models associated with our poor knowledge of chemical composition of the heat insulating envelopes strongly complicate theoretical reconstruction of the internal structure of cooling neutron stars from observations of their thermal surface emission.

preprint2016arXiv

Neutron star properties and the equation of state for its core

Few unified equations of state for neutron star matter where core and crust are described using the same nuclear model are available. However the use of non-unified equations of state with a simplified matching between the crust and the core has been shown to introduce uncertainties in the radius determination which can be larger than the expected precision of the next generation of X-ray satellites. We aim at eliminating the dependence of the radius and mass of neutron staron the detailed model for the crust and on the crust-core matching procedure. We solve the approximate equations of the hydrostatic equilibrium for the crust of neutron stars obtaining a precise formula for the radius which depends only on the core mass and radius, and on the baryon chemical potential at the core-crust interface and on the crust surface. For a fully accreted crust one needs additionally the value of the total deep crustal heating per one accreted nucleon. For typical neutron star masses the approximate approach allows to determine the neutron star radius with an error ~0.1% (~ 10 m, equivalent to a 1% inaccuracy in the crust thickness). The formalism applies to neutron stars with a catalyzed or a fully accreted crust. The difference in the neutron star radius between the two models is proportional to the total energy release due to deep crustal heating. For a given model of dense matter describing the neutron star core, the radius of a neutron star can be accurately determined independently of the crust model with a precision much better than the ~5% one expected from the next generation of X-ray satellites. This allows to circumvent the problem of the radius uncertainty which may arise when non-unified equations of state for the crust and the core are used.

preprint2016arXiv

Rotating neutron stars with exotic cores: masses, radii, stability

A set of theoretical mass-radius relations for rigidly rotating neutron stars with exotic cores, obtained in various theories of dense matter, is reviewed. Two basic observational constraints are used: the largest measured rotation frequency (716 Hz) and the maximum measured mass ($2\;M_\odot$). Present status of measuring the radii of neutron stars is described. The theory of rigidly rotating stars in general relativity is reviewed and limitations of the slow rotation approximation are pointed out. Mass-radius relations for rotating neutron stars with hyperon and quark cores are illustrated using several models. Problems related to the non-uniqueness of the crust-core matching are mentioned. Limits on rigid rotation resulting from the mass-shedding instability and the instability with respect to the axisymmetric perturbations are summarized. The problem of instabilities and of the back-bending phenomenon are discussed in detail. Metastability and instability of a neutron star core in the case of a first-order phase transition, both between pure phases, and into a mixed-phase state, are reviewed. The case of two disjoint families (branches) of rotating neutron stars is discussed and generic features of neutron-star families and of core-quakes triggered by the instabilities are considered.

preprint2015arXiv

Neutron drip transition in accreting and nonaccreting neutron star crusts

The neutron-drip transition in the dense matter constituting the interior of neutron stars generally refers to the appearance of unbound neutrons as the matter density reaches some threshold density $ρ_\textrm{drip}$. This transition has been mainly studied under the cold catalyzed matter hypothesis. However, this assumption is unrealistic for accreting neutron stars. After examining the physical processes that are thought to be allowed in both accreting and nonaccreting neutron stars, suitable conditions for the onset of neutron drip are derived and general analytical expressions for the neutron drip density and pressure are obtained. Moreover, we show that the neutron-drip transition occurs at lower density and pressure than those predicted within the mean-nucleus approximation. This transition is studied numerically for various initial composition of the ashes from X-ray bursts and superbursts using microscopic nuclear mass models.

preprint2015arXiv

Neutron stars with hyperon cores: stellar radii and EOS near nuclear density

The existence of 2 Msun pulsars puts very strong constraints on the equation of state (EOS) of neutron stars (NSs) with hyperon cores, which can be satisfied only by special models of hadronic matter. The radius-mass relation for these models is sufficiently specific that it could be subjected to an observational test with future X-ray observatories. We want to study the impact of the presence of hyperon cores on the radius-mass relation for NS. We aim to find out how, and for which particular stellar mass range, a specific relation R(M), where M is the gravitational mass, and R is the circumferential radius, is associated with the presence of a hyperon core. We consider a set of 14 theoretical EOS of dense matter, based on the relativistic mean-field (RMF) approximation, allowing for the presence of hyperons in NSs. We seek correlations between R(M) and the stiffness of the EOS below the hyperon threshold needed to pass the 2 Msun test. For NS masses 1.0<M/Msun<1.6, we get R>13km, because of a very stiff pre-hyperon segment of the EOS. At nuclear density, the pressure is significantly higher than a robust upper bound obtained recently using chiral effective field theory. If massive NSs do have a sizable hyperon core, then according to current models the radii for M=1.0-1.6 Msun are necessarily >13km. If, on the contrary, a NS with a radius R<12 km is observed in this mass domain, then sizable hyperon cores in NSs, as we model them now, are ruled out. Future X-ray missions with <5% precision for a simultaneous M and R measurement will have the potential to solve the problem with observations of NSs. Irrespective of this observational test, present EOS allowing for hyperons that fulfill condition M_max>2 Msun yield a pressure at nuclear density that is too high relative to up-to-date microscopic calculations of this quantity.

preprint2015arXiv

Progenitor neutron stars of the lightest and heaviest millisecond pulsars

The recent mass measurements of two binary millisecond pulsars, PSR J1614-2230 and PSR J0751+1807 with a mass M=1.97+/-0.04 Msun and M= 1.26 +/- 0.14 Msun, respectively, indicate a wide range of masses for such objects and possibly also a broad spectrum of masses of neutron stars born in core-collapse supernovae. Starting from the zero-age main sequence binary stage, we aim at inferring the birth masses of PSR J1614-2230 and PSR J0751+1807 by taking the differences in the evolutionary stages preceding their formation into account. Using simulations for the evolution of binary stars, we reconstruct the evolutionary tracks leading to the formation of PSR J1614-2230 and PSR J0751+1807. We analyze in detail the spin evolution due to the accretion of matter from a disk in the intermediate-mass/low-mass X-ray binary. We consider two equations of state of dense matter, one for purely nucleonic matter and the other one including a high-density softening due to the appearance of hyperons. Stationary and axisymmetric stellar configurations in general relativity are used, together with a recent magnetic torque model and observationally-motivated laws for the decay of magnetic field. The estimated birth mass of the neutron stars PSR J0751+1807 and PSR J1614-2230 could be as low as 1.0 Msun and as high as 1.9 Msun, respectively. These values depend weakly on the equation of state and the assumed model for the magnetic field and its accretion-induced decay. The masses of progenitor neutron stars of recycled pulsars span a broad interval from 1.0 Msun to 1.9 Msun. Including the effect of a slow Roche-lobe detachment phase, which could be relevant for PSR J0751+1807, would make the lower mass limit even lower. A realistic theory for core-collapse supernovae should account for this wide range of mass.

preprint2012arXiv

Formation of millisecond pulsars - NS initial mass and EOS constraints

Recent measurement of a high millisecond pulsar mass (PSR J1614-2230, 1.97+-0.04 Msun) compared with the low mass of PSR J0751+1807 (1.26+-0.14 Msun) indicates a large span of masses of recycled pulsars and suggests a broad range of neutron stars masses at birth. We aim at reconstructing the pre-accretion masses for these pulsars while taking into account interaction of the magnetic field with a thin accretion disk, magnetic field decay and relativistic 2D solutions for stellar configurations for a set of equations of state. We briefly discuss the evolutionary scenarios leading to the formation of these neutron stars and study the influence of the equation of state.

preprint2012arXiv

Maximum mass of neutron stars and strange neutron-star cores

Recent measurement of mass of PSR J1614-2230 rules out most of existing models of equation of state (EOS) of dense matter with high-density softening due to hyperonization, based on the recent hyperon-nucleon and hyperon-hyperon interactions, leading to a "hyperon puzzle". We study a specific solution of "hyperon puzzle", consisting in replacing a too soft hyperon core by a sufficiently stiff quark core. We construct an analytic approximation fitting very well modern EOSs of 2SC and CFL color superconducting phases of quark matter. This allows us for simulating continua of sequences of first-order phase transitions from hadronic matter to the 2SC, and then to the CFL state of color superconducting quark matter. We obtain constraints in the parameter space of the EOS of superconducting quark cores, resulting from M_max> 2 M_sol. We also derive constraints that would result from significantly higher measured masses. For 2.4 M_sol required stiffness of the CFL quark core should have been close to the causality limit, the density jump at the phase transition being very small. Condition M_max > 2 M_sol puts strong constraints on the EOSs of the 2SC and CFL phases of quark matter. Density jumps at the phase transitions have to be sufficiently small and sound speeds in quark matter - sufficiently large. A strict condition of thermodynamic stability of quark phase results in the maximum mass of hybrid stars similar to that of purely baryon stars. Therefore, to get M_max>2 M_sol for stable hybrid stars, both sufficiently strong additional hyperon repulsion at high density baryon matter and a sufficiently stiff EOS of quark matter would be needed. However, it is likely that the high density instability of quark matter (reconfinement) indicates actually the inadequacy of the point-particle model of baryons in dense matter at very high densities.

preprint2011arXiv

Compression of matter in the center of accreting neutron stars

To estimate the feasibility of dense-matter phase transition, we studied the evolution of the central density as well as the baryon chemical potential of accreting neutron stars. We compared the thin-disk accretion with and without the magnetic field torque with the spin-down scenario for a selection of recent equations of state. We compared the prevalent (in the recycled-pulsar context) Keplerian thin-disk model, in which the matter is accreted from the marginally-stable circular orbit, with the recent magnetic-torque model that takes into account the influence of stellar magnetic field on the effective inner boundary of the disk. Calculations were performed using a multi-domain spectral methods code in the framework of General Relativity. We considered three equations of state consistent with the recently measured mass of PSR J1614-2230, 1.97 +- 0.04 M_sun (one of them softened by the appearance of hyperons). If there is no magnetic torque and efficient angular momentum transfer from the disk to the star, substantial central compression is limited to the region of initial stellar masses close to the maximum mass. Outside the maximum mass vicinity, accretion-induced central compression is significant only if the angular momentum transfer is inefficient. Accounting for the magnetic field effectively decreases the efficiency of angular momentum transfer and implies a significant central compression. An efficient angular momentum transfer from a thin disk onto a non-magnetized neutron star does not provide a good mechanism for the central compression and possible phase transition. Substantial central compression is possible for a broad range of masses of slowly-rotating initial configurations for magnetized neutron stars. Accretion-induced central compression is particularly strong for stiff equation of state with a high-density softening.

preprint2011arXiv

Formation scenarios and mass-radius relation for neutron stars

Neutron star crust, formed via accretion of matter from a companion in a low-mass X-ray binary (LMXB), has an equation of state (EOS) stiffer than that of catalyzed matter. At a given neutron star mass, M, the radius of a star with an accreted crust is therefore larger, by DR(M), than for usually considered star built of catalyzed matter. Using a compressible liquid drop model of nuclei, we calculate, within the one-component plasma approximation, the EOSs corresponding to different nuclear compositions of ashes of X-ray bursts in LMXB. These EOSs are then applied for studying the effect of different formation scenarios on the neutron-star mass-radius relation. Assuming the SLy EOS for neutron star's liquid core, derived by Douchin & Haensel (2001), we find that at M=1.4 M_sun the star with accreted crust has a radius more than 100 m larger that for the crust of catalyzed matter. Using smallness of the crust mass compared to M, we derive a formula that relates DR(M) to the difference in the crust EOS. This very precise formula gives also analytic dependence of DR on M and R of the reference star built of catalyzed matter. The formula is valid for any EOS of the liquid core. Rotation of neutron star makes DR(M) larger. We derive an approximate but very precise formula that gives difference in equatorial radii, DR_eq(M), as a function of stellar rotation frequency.

preprint2011arXiv

Hyperons in neutron-star cores and two-solar-mass pulsar

Recent measurement of mass of PSR J1614-2230 rules out most of existing models of equation of state (EOS) of dense matter with high-density softening due to hyperonization or a phase transition to quark matter or a boson condensate. We look for a solution of an apparent contradiction between the consequences stemming from up-to-date hypernuclear data, indicating appearance of hyperons at 3 nuclear densities and existence of a two-solar-mass neutron star. We consider a non-linear relativistic mean field (RMF) model involving baryon octet coupled to meson fields. An effective lagrangian includes quartic terms in the meson fields. The values of the parameters of the model are obtained by fitting semi-empirical parameters of nuclear matter at the saturation point, as well as potential wells for hyperons in nuclear matter and the strength of the Lambda-Lambda attraction in double-Lambda hypernuclei. We propose a non-linear RMF model which is consistent with up-to-date semiempirical nuclear and hypernuclear data and allows for neutron stars with hyperon cores and M larger than 2 solar masses. The model involves hidden-strangenes scalar and vector mesons, coupled to hyperons only, and quartic terms involving vector meson fields. Our EOS involving hyperons is stiffer than the corresponding nucleonic EOS (with hyperons artificially suppressed) above five nuclear densities. Required stiffening is generated by the quartic terms involving hidden-strangeness vector meson.

preprint2011arXiv

Implications of the measured parameters of PSR J1903+0327 for its progenitor neutron star

Using the intrinsic PSR J1903+0327 parameters evaluated from radio observations (mass, rotation period and dipole magnetic field deduced from the timing properties) we calculate the mass of its neutron star progenitor, M_i, at the onset of accretion. Simultaneously, we derive constraints on average accretion rate Mdot and the pre-accretion magnetic field B_i. Spin-up is modelled by accretion from a thin disk, using the magnetic-torque disk-pulsar coupling model proposed by Kluzniak and Rappaport (2007), improved for the existence of relativistic marginally-stable circular orbit. Orbital parameters in the disk are obtained using the space-time generated by a rotating neutron star in the framework of General Relativity. We employ an observationally motivated model of the surface magnetic field decay. We also seek for the imprint of the poorly known equation of state of dense matter on the spin-up tracks - three equations of state of dense matter, consistent with the existence of 2 Msun neutron star, are considered. We find that the minimum average accretion rate should be larger than 2-8 10^(-10) Msun/yr, the highest lower bound corresponding to the stiffest equation of state. We conclude that the influence of magnetic field in the "recycling" process is crucial - it leads to a significant decrease of spin-up rate and larger accreted masses, in comparison to the B=0 model. Allowed B_i-dependent values of M_i are within 1.0-1.4 Msun, i.e., much lower than an oversimplified but widely used B=0 result, where one gets M_i>1.55 Msun. Estimated initial neutron-star mass depends on the assumed dense-matter equation of state.

preprint2010arXiv

Approximate analytic expressions for circular orbits around rapidly rotating compact stars

We calculate stationary configurations of rapidly rotating compact stars in general relativity, to study the properties of circular orbits of test particles in the equatorial plane. We search for simple, but precise, analytical formulae for the orbital frequency, specific angular momentum and binding energy of a test particle, valid for any equation of state and for any rotation frequency of the rigidly rotating compact star, up to the mass-shedding limit. Numerical calculations are performed using precise 2-D codes based on multi-domain spectral methods. Models of rigidly rotating neutron stars and the space-time outside them are calculated for several equations of state of dense matter. Calculations are also performed for quark stars consisting of self-bound quark matter. At the mass-shedding limit, the rotational frequency converges to a Schwarzschildian orbital frequency at the equator. We show that orbital frequency for any orbit outside equator is also approximated by a Schwarzschildian formula. Using a simple approximation for the frame-dragging term, we obtain approximate expressions for the specific angular momentum and specific energy on the corotating circular orbits in the equatorial plane of neutron star, which are valid down to the stellar equator. The formulae recover reference numerical values with typically 1% of accuracy for neutron stars with M > 0.5 M_sun. They are less precise for quark stars consisting of self-bound quark matter.

preprint2009arXiv

Keplerian frequency of uniformly rotating neutron stars and quark stars

We calculate Keplerian (mass shedding) configurations of rigidly rotating neutron stars and quark stars with crusts. We check the validity of empirical formula for Keplerian frequency, f_K, proposed by Lattimer & Prakash, f_K(M)=C (M/M_sun)^1/2 (R/10km)^-3/2, where M is the (gravitational) mass of Keplerian configuration, R is the (circumferential) radius of the non-rotating configuration of the same gravitational mass, and C = 1.04 kHz. Numerical calculations are performed using precise 2-D codes based on the multi-domain spectral methods. We use a representative set of equations of state (EOSs) of neutron stars and quark stars. We show that the empirical formula for f_K(M) holds within a few percent for neutron stars with realistic EOSs, provided 0.5 M_sun < M < 0.9 M_max,stat, where M_max,stat is the maximum allowable mass of non-rotating neutron stars for an EOS, and C=C_NS=1.08 kHz. Similar precision is obtained for quark stars with 0.5 M_sun < M < 0.9 M_max,stat. For maximal crust masses we obtain C_QS = 1.15 kHz, and the value of C_QS is not very sensitive to the crust mass. All our C's are significantly larger than the analytic value from the relativistic Roche model, C_Roche = 1.00 kHz. For 0.5 M_sun < M < 0.9 M_max,stat, the equatorial radius of Keplerian configuration of mass M, R_K(M), is, to a very good approximation, proportional to the radius of the non-rotating star of the same mass, R_K(M) = aR(M), with a_NS \approx a_QS \approx 1.44. The value of a_QS is very weakly dependent on the mass of the crust of the quark star. Both a's are smaller than the analytic value a_Roche = 1.5 from the relativistic Roche model.

preprint2007arXiv

Long Gamma Ray Bursts from Quark Stars

If strange quark matter (SQM) is the true ground state of hadronic matter, then conversion of neutron stars (NS) into quark stars (QS) could release some 10^{53} erg. We describe a scenario of burning of a NS into a hot, differentially rotating QS. Emission of released non-baryonic energy through the QS surface is discussed. The role of magnetobuoyancy of SQM is mentioned. The outflow of gamma e+ e- lasting for up to ~1000 s could be at the origin of long GRBs. Advantages of hot, differentially rotating QS as an inner engine of long GRBs are reviewed.

preprint2005arXiv

Mixed-phase induced core-quakes and the changes in neutron star parameters

We present approximate formulae describing the changes in neutron star parameters caused by the first-order phase transition to an ``exotic'' state (pion or kaon condensate, quark matter) resulting in formation of a mixed-phase core. The analytical formulae for the changes in radius, moment of inertia and the amount of energy released during the core-quake are derived using the theory of linear response of stellar structure to the core transformation. Numerical coefficients in these formulae are obtained for two realistic EOSs of dense matter. The problem of nucleation of the exotic phase as well as possible astrophysical scenarios leading to a core-quake phenomenon is also discussed.

preprint2003arXiv

Nuclear composition and heating in accreting neutron-star crusts

Nuclear reactions in accreting neutron-star crusts and the heat release accompanying them are studied, under different assumptions concerning the composition of the outermost layer formed of the ashes of X-ray bursts. Particular examples of ashes containing nuclides with A ~ 90-110 are considered and compared with a standard A=56 case. In all cases, evolution of a crust shell is followed from 10^8 g/cm^3 to a few times 10^{13} g/cm^3. The total crustal heating produced in the non-equilibrium processes in the accreting crust is 1.1-1.5 MeV per one accreted nucleon. The composition of the accreted crust at densities exceeding the threshold for the pycnonuclear fusion is essentially independent of the assumed initial composition of the X-ray burst ashes.

preprint2001arXiv

Recycling strange stars to millisecond periods

Recycling strange stars to millisecond periods is studied within the framework of general relativity. We employ equations of state of strange quark matter based on the MIT Bag Model, with massive strange quarks and lowest order QCD interactions. The presence of the crust of normal matter is taken into account, with a bottom density assumed to be equal to the neutron-drip one. The calculations are performed by solving the exact 2-D equations for rigidly rotating stationary configurations in general relativity. Evolutionary tracks of accreting strange stars are computed, and their dependence on the initial strange star mass and on the fraction of the angular momentum transferred to the star from the infalling matter, is studied. The differences between recycling strange stars and neutron stars are pointed out.

preprint2001arXiv

Small strange stars and marginally stable orbit in Newtonian theory

It is shown that for very rapidly rotating low mass strange stars the marginally stable orbit is located above the stellar surface. This effect is explained by the very important role of the oblateness of the rotating strange star. The comparison with some ``academic'' examples is presented. This feature is purely Newtonian in its nature and has nothing to do with relativistic marginally stable orbit. The effect is very large and cannot be treated in a perturbative way. It seems that strange stars as a very dense self-bound objects are the only possibility in Nature to represent these toy models.

preprint2001arXiv

The crust of rotating strange quark stars

Calculations of the properties of rotating strange stars with crusts are performed within the framework of general relativity. We employ an equation of state (EOS) of strange quark matter based on the MIT Bag Model with massive strange quarks and lowest order QCD interactions. The crust is described by the BPS equation of state. A significant increase of the stellar radius is found close to the Keplerian (mass-shedding limit) configuration. This leads to the disappearance of the gap between the stellar surface and the innermost stable circular orbit (ISCO) at very high rotation rates, for a rather broad range of stellar masses. The Keplerian configuration for the strange star with crust corresponds to values of J, T/W which are about 20% smaller than in the case of bare strange stars. Because the Keplerian configuration is achieved due to the increase of the stellar oblateness, the Keplerian frequency (of the rotation) remains almost unaltered. The lack of the gap close to the Keplerian rotation could imply a more stringent limit on the frequency of the ISCO, if the existence of such a gap is supported by observations, as in the source 4U 1820-30 with the upper QPO frequency 1.07 kHz. During the spin-down of a magnetized strange quark star with crust, the crust matter is absorbed in the equatorial region by the strange matter core. The deconfinement of absorbed crust matter is a strongly exothermic process, which would influence the cooling curve of this compact object.

preprint2000arXiv

Innermost stable circular orbits around strange stars and kHz QPOs in low-mass X-ray binaries

Exact calculations of innermost stable circular orbit (ISCO) around rotating strange stars are performed within the framework of general relativity. Equations of state (EOS) of strange quark matter based on the MIT Bag Model with massive strange quarks and lowest order QCD interactions, are used. The presence of a solid crust of normal matter on rotating, mass accreting strange stars in LMXBs is taken into account. It is found that, contrary to neutron stars, above some minimum mass a gap always separates the ISCO and stellar surface, independently of the strange star rotation rate. For a given baryon mass of strange star, we calculate the ISCO frequency as function of stellar rotation frequency, from static to Keplerian configuration. For masses close to the maximum mass of static configurations the ISCO frequencies for static and Keplerian configurations are similar. However, for masses significantly lower than the maximum mass of static configurations, the minimum value of the ISCO frequency is reached in the Keplerian limit. Presence of a solid crust increases the ISCO frequency for the Keplerian configuration by about ten percent compared to that for a bare strange star of the same mass.

preprint2000arXiv

Rotating compact strange stars

We compute numerical models of uniformly rotating strange stars (SS) in general relativity for the recently proposed QCD-based equation of state (EOS) of strange quark matter (Dey et al. 1998). Static models based on this EOS are characterised by a larger surface redshift than strange stars within the MIT bag model. The frequencies of the fastest rotating configurations described by Dey model are much higher than these for neutron stars (NS) and for the simplest SS MIT bag model. We determine a number of physical parameters for such stars and compare them with those obtained for NS. We construct constant baryon mass equilibrium sequences both normal and supramassive. Similarly to the NS a supramassive SS, prior to collapse to a black hole, spins up as it loses angular momentum. We find the upper limits on maximal masses and maximal frequencies of the rotating configurations. We show that the maximal rotating frequency for each of considered evolutionary sequences is never the Keplerian one. A normal and low mass supramassive strange stars gaining angular momentum always slows down just before reaching the Keplerian limit. For a high mass supramassive SS sequence the Keplerian configuration is the one with the lowest rotational frequency in the sequence. The value of $T/W$ for rapidly rotating SS of any mass is significantly higher than those for ordinary NS. For Keplerian configurations it increases as mass decreases. The results are robust for all linear self-bound equations of state.

preprint2000arXiv

Strange stars - linear approximation of the EOS and maximum QPO frequency

The equation of state of strange matter build of the u,d,s quarks in the framework of the MIT bag model is studied. The scaling relations with bag constant are discussed and applied to the determination of maximum frequency of a particle in stable circular orbit around strange star. The highest QPO frequency of 1.33 kHz observed so far is consistent with the strange stars models for which the maximum QPO frequency is 1.7-2.4 kHz depending on the strange quark mass and the QCD coupling constant. The linear approximation of the equation of state is found and the parameters of this EOS are determined as a functions of strange quark mass, QCD coupling constant and bag constant. This approximation reproduces exact results within an error of the order of 1% and can be used for the complete study of the properties of strange stars including microscopic stability of strange matter and determination of the total baryon number of the star.

preprint1995arXiv

Damping of GRR instability by direct URCA reactions

The role of direct URCA reactions in damping of the gravitational radiation driven instability is discussed. The temperature at which bulk viscosity suppresses completely this instability is calculated. The results are obtained analytically using recent calculations performed in the case of bulk viscosity due to the modified URCA processes (Lindblom 1995; Yoshida & Eriguchi 1995). The bulk viscosity caused by direct URCA reactions is found to reduce significantly the region of temperatures and rotation frequencies where a neutron star is subject to GRR instability.