Source author record

P. Haensel

P. Haensel 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

45works
9topics
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

45 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

Neutron star radii and crusts: uncertainties and unified equations of state

The uncertainties in neutron star (NS) radii and crust properties due to our limited knowledge of the equation of state (EOS) are quantitatively analysed. We first demonstrate the importance of a unified microscopic description for the different baryonic densities of the star. If the pressure functional is obtained matching a crust and a core EOS based on models with different properties at nuclear matter saturation, the uncertainties can be as large as $\sim 30\%$ for the crust thickness and $4\%$ for the radius. Necessary conditions for causal and thermodynamically consistent matchings between the core and the crust are formulated and their consequences examined. A large set of unified EOS for purely nucleonic matter is obtained based on 24 Skyrme interactions and 9 relativistic mean-field nuclear parametrizations. In addition, for relativistic models 17 EOS including a transition to hyperonic matter at high density are presented. All these EOS have in common the property of describing a $2\;M_\odot$ star and of being causal within stable NS. A span of $\sim 3$ km and $\sim 4$ km is obtained for the radius of, respectively, $1.0\;M_\odot$ and $2.0\;M_\odot$ star. Applying a set of nine further constraints from experiment and ab-initio calculations the uncertainty is reduced to $\sim 1$ km and $2$ km, respectively. These residual uncertainties reflect lack of constraints at large densities and insufficient information on the density dependence of the EOS near the nuclear matter saturation point. The most important parameter to be constrained is shown to be the symmetry energy slope $L$ which exhibits a linear correlation with the stellar radius, particularly for masses $\sim 1.0\;M_\odot$. Potential constraints on $L$, the NS radius and the EOS from observations of thermal states of NS are also discussed. [Abriged]

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.

preprint2016arXiv

Theory of neutrino emission from nucleon-hyperon matter in neutron stars: Angular integrals

Investigations of thermal evolution of neutron stars with hyperon cores require neutrino emissivities for many neutrino reactions involving strongly degenerate particles (nucleons, hyperons, electrons, muons). We calculate the angular integrals $I_n$ (over orientations of momenta of $n$ degenerate particles) for major neutrino reactions with $n$ =3, 4, 5 at all possible combinations of particle Fermi momenta. The integrals $I_n$ are necessary ingredients for constructing a uniform database of neutrino emissivities in dense nucleon-hyperon matter. The results can also be used in many problems of physical kinetics of strongly degenerate systems.

preprint2015arXiv

Energy bursts from deconfinement in high-mass twin stars

We estimate the energy reservoir available in the deconfinement phase transition induced collapse of a neutron star to its hybrid star mass twin on the "third family" branch, using a recent equation of state of dense matter. The available energy corresponding to the mass-energy difference between configurations is comparable with energies of the most violent astrophysical burst processes. An observational outcome of such a dynamical transition might be fast radio bursts, specifically a recent example of a FRB with a double-peak structure in its light curve.

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.

preprint2014arXiv

Physics input for modelling superfluid neutron stars with hyperon cores

Observations of massive ($M \approx 2.0~M_\odot$) neutron stars (NSs), PSRs J1614-2230 and J0348+0432, rule out most of the models of nucleon-hyperon matter employed in NS simulations. Here we construct three possible models of nucleon-hyperon matter consistent with the existence of $2~M_\odot$ pulsars as well as with semi-empirical nuclear matter parameters at saturation, and semi-empirical hypernuclear data. Our aim is to calculate for these models all the parameters necessary for modelling dynamics of hyperon stars (such as equation of state, adiabatic indices, thermodynamic derivatives, relativistic entrainment matrix, etc.), making them available for a potential user. To this aim a general non-linear hadronic Lagrangian involving $σωρϕσ^\ast$ meson fields, as well as quartic terms in vector-meson fields, is considered. A universal scheme for calculation of the $\ell=0,1$ Landau Fermi-liquid parameters and relativistic entrainment matrix is formulated in the mean-field approximation. Use of this scheme allow us to obtain numerical tables with the equation of state, Landau quasiparticle effective masses, adiabatic indices, the $\ell=0,1$ Landau Fermi-liquid parameters, and the relativistic entrainment matrix for the selected models of nucleon-hyperon matter. These data are available on-line and suitable for numerical implementation in computer codes modelling various dynamical processes in NSs, in particular, oscillations of superfluid NSs and their cooling.

preprint2014arXiv

The Large Observatory For x-ray Timing

The Large Observatory For x-ray Timing (LOFT) was studied within ESA M3 Cosmic Vision framework and participated in the final down-selection for a launch slot in 2022-2024. Thanks to the unprecedented combination of effective area and spectral resolution of its main instrument, LOFT will study the behaviour of matter under extreme conditions, such as the strong gravitational field in the innermost regions of accretion flows close to black holes and neutron stars, and the supra-nuclear densities in the interior of neutron stars. The science payload is based on a Large Area Detector (LAD, 10 m 2 effective area, 2-30 keV, 240 eV spectral resolution, 1 deg collimated field of view) and a WideField Monitor (WFM, 2-50 keV, 4 steradian field of view, 1 arcmin source location accuracy, 300 eV spectral resolution). The WFM is equipped with an on-board system for bright events (e.g. GRB) localization. The trigger time and position of these events are broadcast to the ground within 30 s from discovery. In this paper we present the status of the mission at the end of its Phase A study.

preprint2013arXiv

Transport coefficients of nuclear matter in neutron star cores

We calculate thermal conductivity and shear viscosity of nucleons in dense nuclear matter of neutron star cores in the non-relativistic Brueckner-Hartree-Fock framework. Nucleon-nucleon interaction is described by the Argonne v18 potential with addition of the Urbana IX effective three-body forces. We find that this three body force leads to decrease of the kinetic coefficients with respect to the two-body case. The results of calculations are compared with electron and muon transport coefficients as well as with the results of other authors.

preprint2012arXiv

Disk accretion onto a magnetized star

The problem of interaction of the rotating magnetic field, frozen to a star, with a thin well conducting accretion disk is solved exactly. It is shown that a disk pushes the magnetic field lines towards a star, compressing the stellar dipole magnetic field. At the point of corotation, where the Keplerian rotation frequency coincides with the frequency of the stellar rotation, the loop of the electric current appears. The electric currents flow in the magnetosphere only along two particular magnetic surfaces, which connect the corotation region and the inner edge of a disk with the stellar surface. It is shown that the closed current surface encloses the magnetosphere. Rotation of a disk is stopped at some distance from the stellar surface, which is 0.55 of the corotation radius. Accretion from a disk spins up the stellar rotation. The angular momentum transferred to the star is determined.

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.

preprint2007arXiv

Neutron star cooling after deep crustal heating in the X-ray transient KS 1731-260

We simulate the cooling of the neutron star in the X-ray transient KS 1731-260 after the source returned to quiescence in 2001 from a long (>~ 12.5 yr) outburst state. We show that the cooling can be explained assuming that the crust underwent deep heating during the outburst stage. In our best theoretical scenario the neutron star has no enhanced neutrino emission in the core, and its crust is thin, superfluid, and has the normal thermal conductivity. The thermal afterburst crust-core relaxation in the star may be not over.

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

Accelerated expansion of the Crab Nebula and evaluation of its neutron-star parameters

A model of an accelerated expansion of the Crab Nebula powered by the spinning-down Crab pulsar is proposed, in which time dependence of the acceleration is connected with evolution of pulsar luminosity. Using recent observational data, we derive estimates of the Crab neutron-star moment of inertia. Correlations between the neutron star moment of inertia and its mass and radius allow for rough estimates of the Crab neutron-star radius and mass. In contrast to the previously used constant-acceleration approximation, even for the expanding nebula mass ~7 M_sun results obtained within our model do not stay in conflict with the modern stiff equations of state of dense matter.

preprint2003arXiv

Equation of state of dense matter and maximum mass of neutron stars

Theoretical models of the equation of state (EOS) of neutron-star matter (starting with the crust and ending at the densest region of the stellar core) are reviewed. Apart from a broad set of baryonic EOSs, strange quark matter, and even more exotic (abnormal and Q-matter) EOSs are considered. Results of calculations of M_max for non-rotating neutron stars and exotic compact stars are reviewed, with particular emphasis on the dependence on the dense-matter EOS. Rapid rotation increases M_max, and this effect is studied for both neutron stars and exotic stars. Theoretical results are then confronted with measurements of masses of neutron stars in binaries, and the consequences of such a confrontation and their possible impact on the theory of dense matter are 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.

preprint2003arXiv

Thermal state of transiently accreting neutron stars

We study thermal states of transiently accreting neutron stars (with mean accretion rates $\dot{M} \sim 10^{-14}-10^{-9}$ M$_\odot$ yr$^{-1}$) determined by the deep crustal heating of accreted matter sinking into stellar interiors. We formalize a direct correspondence of this problem to the problem of cooling neutron stars. Using a simple toy model we analyze the most important factors which affect the thermal states of accreting stars: a strong superfluidity in the cores of low-mass stars and a fast neutrino emission (in nucleon, pion-condensed, kaon-condensed, or quark phases of dense matter) in the cores of high-mass stars. We briefly compare the results with the observations of soft X-ray transients in quiescence. If the upper limit on the quiescent thermal luminosity of the neutron star in SAX J1808.4-3658 (Campana et al. 2002) is associated with the deep crustal heating, it favors the model of nucleon neutron-star cores with switched-on direct Urca process.

preprint2001arXiv

Bulk viscosity in superfluid neutron star cores. III. Effects of $Σ^-$ hyperons

Bulk viscosity of neutron star cores containing hyperons is studied taking into account non-equilibrium weak process $n+n \rightleftharpoons p+Σ^-$. Rapid growth of the bulk viscosity within the neutron star core associated with switching on new reactions (modified Urca process, direct Urca process, hyperon reactions) is analyzed. The suppression of the bulk viscosity by superfluidity of baryons is considered and found out to be very important.

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

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

Neutrino Emission from Neutron Stars

We review the main neutrino emission mechanisms in neutron star crusts and cores. Among them are the well-known reactions such as the electron-positron annihilation, plasmon decay, neutrino bremsstrahlung of electrons colliding with atomic nuclei in the crust, as well as the Urca processes and neutrino bremsstrahlung in nucleon-nucleon collisions in the core. We emphasize recent theoretical achievements, for instance, band structure effects in neutrino emission due to scattering of electrons in Coulomb crystals of atomic nuclei. We consider the standard composition of matter (neutrons, protons, electrons, muons, hyperons) in the core, and also the case of exotic constituents such as the pion or kaon condensates and quark matter. We discuss the reduction of the neutrino emissivities by nucleon superfluidity, as well as the specific neutrino emission produced by Cooper pairing of the superfluid particles. We also analyze the effects of strong magnetic fields on some reactions, such as the direct Urca process and the neutrino synchrotron emission of electrons. The results are presented in the form convenient for practical use. We illustrate the effects of various neutrino reactions on the cooling of neutron stars. In particular, the neutrino emission in the crust is critical in setting the initial thermal relaxation between the core and the crust. Finally, we discuss the prospects of exploring the properties of supernuclear matter by confronting cooling simulations with observations of the thermal radiation from isolated neutron stars.

preprint1999arXiv

Fast rotation of strange stars

Exact models of uniformly rotating strange stars, built of self bound quark matter, are calculated within the framework of general relativity. This is made possible thanks to a new numerical technique capable to handle the strong density discontinuity at the surface of these stars. Numerical calculations are done for a simple MIT bag model equation of state of strange quark matter. Evolutionary sequences of models of rotating strange stars at constant baryon mass are calculated. Maximally rotating configurations of strange stars are determined, assuming that the rotation frequency is limited by the mass shedding and the secular instability with respect to axisymmetric perturbations. Exact formulae which give the dependence of the maximum rotation frequency, and of the maximum mass and corresponding radius of rotating configurations, on the value of the bag constant, are obtained. The values of T/W for rapidly rotating massive strange stars are significantly higher than those for ordinary neutron stars. This might indicate particular susceptibility of rapidly rotating strange stars to triaxial instabilities.

preprint1999arXiv

Neutrino emission due to electron bremsstrahlung in superfluid neutron-star cores

We study neutrino energy emission rates (emissivities) due to electron bremsstrahlung produced by $ee$ and $ep$ collisions in the superfluid neutron star cores. The neutrino emission due to $ee$ collisions is shown to be the dominant neutrino reaction at not too high temperatures ($T \la 10^8$ K) in dense matter if all other neutrino reactions involving nucleons are strongly suppressed by neutron and proton superfluidity. Simple practical expressions for the $ee$ and $ep$ neutrino emissivities are obtained. The efficiency of various neutrino reactions in the superfluid neutron-star cores is discussed for the cases of standard neutrino energy losses and the losses enhanced by the direct Urca process.

preprint1999arXiv

Neutrino emission due to proton pairing in neutron stars

We calculate the neutrino energy emission rate due to singlet-state pairing of protons in the neutron star cores taking into account the relativistic correction to the non-relativistic rate. The non-relativistic rate is numerically small, and the relativistic correction appears to be about 10 -- 50 times larger. It plays thus the leading role, reducing great difference between the neutrino emissions due to pairing of protons and neutrons. The results are important for simulations of neutron star cooling.

preprint1999arXiv

Transport properties and neutrino emissivity of dense neutron-star matter with localized protons

As pointed out by Kutschera and W{ó}jcik, very low concentration of protons combined with a specific density dependence of effective neutron-proton interaction could lead to a localization of ``proton impurities'' in neutron medium at densities exceeding four times normal nuclear matter density. We study consequences of the localization of protons for transport processes in dense neutron star cores, assuming random distribution of proton impurities. Kinetic equations, relevant for the transport of charge, heat and momentum, are solved using variational method. Localization of protons removes a T^{-2} factor from the transport coefficients, which leads, at lower temperatures, to a strong decrease of thermal conductivity, electrical conductivity and shear viscosity of neutron star matter, as compared to the standard case, where protons form a Fermi liquid. Due to the localization of protons a number of conventional neutrino emission processes (including modified URCA process) become inoperative in neutron star cores. On the other hand, the energy loss rate from neutrino-antineutrino pair bremsstrahlung due to electron and neutron scattering off (localized) protons, will have a specific T^6 dependence, which could modify the cooling of the neutron star core, as compared to the standard case. Possible astrophysical implications of the localization of protons for neutron star evolution and dynamics are discussed.

preprint1999arXiv

Transport properties of degenerate electrons in neutron star envelopes and white dwarf cores

New calculations of the thermal and electrical electron conductivities are performed for a broad range of physical parameters typical for envelopes of neutron stars and cores of white dwarfs. We consider stellar matter composed of astrophysically important chemical elements from H to Fe in the density range from 10^2-10^4 g/ccm up to 10^7-10^{10} g/ccm, where atoms are fully ionized and electrons are strongly degenerate. We have used modified ion structure factors suggested in physics/9811052. In the ion liquid, these modifications take into account, in an approximate way, instantaneous electron-band structures that reduce the electron-ion scattering rate. In crystallized matter, the new structure factors include multi-phonon processes important at temperatures not very much lower than the melting temperature T_m. The transport coefficients obtained differ significantly from those derived earlier in the important temperature range T_m/5 < T < 5 T_m. The results of our numerical calculations are fitted by analytical expressions convenient for astrophysical applications.

preprint1996arXiv

Electron $ν\barν$ Bremsstrahlung in a Liquid Phase of Neutron Star Crusts

Neutrino emissivity from the electron $ν\barν$ bremsstrahlung in the liquid layers of the neutron star crusts is studied. Nuclear composition of matter in neutron star crusts is considered for various scenarios of neutron star evolution. For the deep layers of the crust, the compositions of cold catalyzed matter, accreted matter and hot matter ($T \ga 5 \times 10^9$~K) are shown to be very different, and this implies differences in the neutrino emissivity at given density and temperature. Neutrino-pair bremsstrahlung, due to collisions of relativistic degenerate electrons in a Coulomb liquid of atomic nuclei, is considered. The neutrino energy loss rate is expressed in a simple form through a Coulomb logarithm $L$ -- a slowly varying function of density, temperature, and nuclear composition. Non-Born corrections, thermal width of electron Fermi sphere, and finite sizes of nuclei are taken into account. Implications for cooling of neutron stars are discussed.

preprint1996arXiv

Ferromagnetism of dense matter and magnetic properties of neutron stars

Possible consequences of ferromagnetic transition in dense matter suggested recently by Kutschera and W{ó}jcik, for the magnetic properties of neutron stars, are studied. Specific model of dense matter, in which a small admixture of protons is completely polarized due to their interaction with neutrons, is considered. Magnetic field of neutron stars with a ferromagnetic core is calculated within the framework of general relativity. Two types of boundary conditions at the ferromagnetic core edge are considered, corresponding to normal and superconducting liquid envelope, respectively. Numerical results for the neutron star magnetic dipole moment are confronted with pulsar timing. To be consistent with observations, ferromagnetic cores surrounded by a non-superconducting envelope, should consist of weakly ordered ferromagnetic domains. If domains are highly ordered, ferromagnetic core should be screened by a superconducting envelope.