Normalized to: Haensel, P.

**117** article(s) in total. **413** co-authors, from **1** to **53** common article(s). Median position in authors list is **2,0**.

[1] oai:arXiv.org:2002.12827 [pdf] - 2056411

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

Comments: To be submitted for publication. Remarks and comments on the present
paper will be welcome by the authors

Submitted: **2020-02-28**

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, $\mu_n$ and $\mu_{\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^{\Phi}\mu_n$, where
the spacetime metric component $g_{00}=e^{2\Phi}$, 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^{\Phi(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^{\Phi(r)}{\widetilde{g}}(r)$ is an extension of the standard
relation $e^{\Phi(r)}{{\mu}}_{\rm b}(r)=$constant, valid in hydrostatic
equilibrium for catalyzed crust.

[2] oai:arXiv.org:1911.08407 [pdf] - 2001943

Bulk viscosity in neutron stars with hyperon cores

Comments: 21 pages, 13 figures, 4 tables, accepted in PRD (selected as PRD
editors' suggestion)

Submitted: **2019-11-19**

It is well-known that r-mode oscillations of rotating neutron stars may be
unstable with respect to the gravitational wave emission. It is highly unlikely
to observe a neutron star with the parameters within the instability window, a
domain where this instability is not suppressed. But if one adopts the
`minimal' (nucleonic) composition of the stellar interior, a lot of observed
stars appear to be within the r-mode instability window. One of the possible
solutions to this problem is to account for hyperons in the neutron star core.
The presence of hyperons allows for a set of powerful (lepton-free)
non-equilibrium weak processes, which increase the bulk viscosity, and thus
suppress the r-mode instability. Existing calculations of the instability
windows for hyperon NSs generally use reaction rates calculated for the
$\Sigma^-\Lambda$ hyperonic composition via the contact $W$ boson exchange
interaction. In contrast, here we employ hyperonic equations of state where the
$\Lambda$ and $\Xi^-$ are the first hyperons to appear (the $\Sigma^-$'s, if
they are present, appear at much larger densities), and consider the meson
exchange channel, which is more effective for the lepton-free weak processes.
We calculate the bulk viscosity for the non-paired $npe\mu\Lambda\Xi^-$ matter
using the meson exchange weak interaction. A number of viscosity-generating
non-equilibrium processes is considered (some of them for the first time in the
neutron-star context). The calculated reaction rates and bulk viscosity are
approximated by simple analytic formulas, easy-to-use in applications. Applying
our results to calculation of the instability window, we argue that accounting
for hyperons may be a viable solution to the r-mode problem.

[3] oai:arXiv.org:1809.08609 [pdf] - 1767678

Bulk viscosity in a neutron star mantle

Comments: 7 pages, 1 figure, MNRAS, accepted

Submitted: **2018-09-23**

We study the bulk viscosity in two (anti-spaghetti and Swiss cheese) phases
of non- spherical nuclei in the mantle of an oscillating neutron star near the
boundary with the stellar core. The bulk viscosity is produced by
non-equilibrium Urca neutrino emis- sion processes. In the mantle, the direct
Urca process may be open (Gusakov et al., 2004) if neutrons and protons move in
a periodic potential created by a lattice of non- spherical nuclei (which
allows the nucleons to have large quasi-momenta and satisfy direct Urca
momentum-conservation). This bulk viscosity can dominate over that due to the
modified Urca process in the outer stellar core and over the shear viscosity.
The bulk viscosity depends strongly on temperature, oscillation frequency and
nucleon superfluidity. The enhanced bulk viscosity in the mantle can control
propagation and damping of neutron star oscillations.

[4] oai:arXiv.org:1807.06569 [pdf] - 1755060

In-medium enhancement of the modified Urca neutrino reaction rates

Comments: 8 pages, 3 figures; accepted for publication in PLB. In v.2 misprint
in eq.(9) corrected and discussion of cooling curves expanded

Submitted: **2018-07-17**, last modified: **2018-09-18**

We calculate modified Urca neutrino emission rates in the dense nuclear
matter in neutron star cores. We find that these rates are strongly enhanced in
the beta-stable matter in regions of the core close to the direct Urca process
threshold. This enhancement can be tracked to the use of the in-medium nucleon
spectrum in the virtual nucleon propagator. We describe the in-medium nucleon
scattering in the non-relativistic Bruckner-Hartree-Fock framework taking into
account two-body as well as the effective three-body forces, although the
proposed enhancement does not rely on a particular way of the nucleon
interaction treatment. Finally we suggest a simple approximate expression for
the emissivity of the n-branch of the modified Urca process that can be used in
the neutron stars cooling simulations with any nucleon equation of state of
dense matter.

[5] oai:arXiv.org:1806.03861 [pdf] - 1795724

Crustal heating in accreting neutron stars from the nuclear
energy-density functional theory. I. Proton shell effects and neutron-matter
constraint

Comments: 20 pages, 13 figures, submitted to A&A

Submitted: **2018-06-11**

Observations of soft X-ray transients in quiescence suggest the existence of
heat sources in the crust of accreting neutron stars. Heat is thought to be
released by electroweak and nuclear processes triggered by the burying of ashes
of X-ray bursts. The heating is studied using a fully quantum approach taking
consistently into account nuclear shell effects. We have followed the evolution
of ashes made of $^{56}$Fe employing the nuclear energy-density functional
theory. Both the outer and inner crusts are described using the same
functional, thus ensuring a unified and thermodynamically consistent treatment.
To assess the role of the neutron-matter constraint, we have employed the set
of accurately calibrated Brussels-Montreal functionals BSk19, BSk20, and BSk21
and for comparison the SLy4 functional. Due to nuclear shell effects, the fully
accreted crust is found to be much less stratified than in previous studies. In
particular, large regions of the inner crust contain clusters with the magic
number $Z=14$. The heat deposited in the outer crust is tightly constrained by
experimental atomic mass data. The shallow heating we obtain does not exceed
$0.2$~MeV and is therefore not enough to explain the cooling of some soft X-ray
transients. The total heat released in the crust is very sensitive to details
of the nuclear structure and is predicted to lie in the range from $1.5$~MeV to
$1.7$~MeV. The evolution of an accreted matter element and therefore the
location of heat sources are governed to a large extent by the existence of
nuclear shell closures. Ignoring these effects in the inner crust, the total
heat falls to $\sim 0.6$~MeV. The neutron-matter constraint is also found to
play a key role. The large amount of heat obtained by Steiner et al. (2012)
could thus be traced back to unrealistic neutron-matter equations of state.

[6] oai:arXiv.org:1711.02414 [pdf] - 1656172

Accurate ray tracing of realistic neutron star atmospheres for
constraining their parameters

Comments: 17 pages, 13 figures; accepted in ApJ ; important updates in Section
2.2 and Fig. 2

Submitted: **2017-11-07**, last modified: **2018-02-21**

Thermal dominated X-ray spectra of neutron stars in quiescent transient X-ray
binaries and neutron stars that undergo thermonuclear bursts are sensitive to
mass and radius. The mass-radius relation of neutron stars depends on the
equation of state that governs their interior. Constraining this relation
accurately is thus of fundamental importance to understand the nature of dense
matter. In this context we introduce a pipeline to calculate realistic model
spectra of rotating neutron stars with hydrogen and helium atmospheres. An
arbitrarily fast rotating neutron star with a given equation of state generates
the spacetime in which the atmosphere emits radiation. We use the
Lorene/nrotstar code to compute the spacetime numerically and the ATM24 code to
solve the radiative transfer equations self-consistently. Emerging specific
intensity spectra are then ray-traced through the neutron star's spacetime from
the atmosphere to a distant observer with the Gyoto code. Here, we present and
test our fully relativistic numerical pipeline. To discuss and illustrate the
importance of realistic atmosphere models we compare our model spectra to
simpler models like the commonly used isotropic color-corrected blackbody
emission. We highlight the importance of considering realistic model-atmosphere
spectra together with relativistic ray tracing to obtain accurate predictions.
We also insist on the crucial impact of the star's rotation on the observables.
Finally, we close a controversy that has been appearing in the literature in
the recent years regarding the validity of the ATM24 code.

[7] oai:arXiv.org:1709.04855 [pdf] - 1626310

Thermal states of neutron stars with a consistent model of interior

Comments:

Submitted: **2017-09-14**

We model the thermal states of both isolated neutron stars and accreting
neutron stars in X-ray transients in quiescence and confront them with
observations. We use an equation of state and superfluid baryon gaps, which are
consistently calculated. We conclude that the direct Urca process is required
to be consistent with low-luminous accreting neutron stars. In addition, proton
superfluidity and sufficiently weak neutron superfluidity are necessary to
explain the cooling of middle-aged neutron stars and to obtain a realistic
distribution of neutron star masses.

[8] oai:arXiv.org:1708.08272 [pdf] - 1587557

Neutrino luminosities and heat capacities of neutron stars in analytic
form

Comments: 15 pages, 9 figures, 5 tables

Submitted: **2017-08-28**

We derive analytic approximations for the neutrino luminosities and the heat
capacities of neutron stars with isothernal nucleon cores as functions of the
mass and radius of stars. The neutrino luminosities are approximated for the
three basic neutrino emission mechanisms, and the heat capacities for the five
basic combinations of the partial heat capacities. The approximations are valid
for for a wide class of equations of state of dense nucleonmatter. The results
significantly simplify the theoretical interpretation of observations of
cooling neutron stars as well as of quasistationary thermal states of neutron
stars in X-ray transients. For illustration, we present an analysis of the
neutrino cooling functions of nine isolated neutron stars taking into account
the effects of their magnetic fields and of the presence of light elements in
their heat blanketing envelopes. These results allow one to investigate the
superfluid properties of neutron star cores.

[9] oai:arXiv.org:1608.07049 [pdf] - 1553870

Consequences of a strong phase transition in the dense matter equation
of state for the rotational evolution of neutron stars

Comments: 9 pages, 7 figures, Astronomy and Astrophysics accepted

Submitted: **2016-08-25**, last modified: **2017-01-08**

We explore the implications of a strong first-order phase transition region
in the dense matter equation of state in the interiors of rotating neutron
stars, and the resulting creation of two disjoint families of neutron-star
configurations (the so-called high-mass twins). We numerically obtained
rotating, axisymmetric, and stationary stellar configurations in the framework
of general relativity, and studied their global parameters and stability. The
instability induced by the equation of state divides stable neutron star
configurations into two disjoint families: neutron stars (second family) and
hybrid stars (third family), with an overlapping region in mass, the high-mass
twin-star region. These two regions are divided by an instability strip. Its
existence has interesting astrophysical consequences for rotating neutron
stars. We note that it provides a natural explanation for the rotational
frequency cutoff in the observed distribution of neutron star spins, and for
the apparent lack of back-bending in pulsar timing. It also straightforwardly
enables a substantial energy release in a mini-collapse to another neutron-star
configuration (core quake), or to a black hole.

[10] oai:arXiv.org:1612.04672 [pdf] - 1533395

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

Comments: 10 pages, 5 figures

Submitted: **2016-12-13**

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.

[11] oai:arXiv.org:1611.01357 [pdf] - 1508846

Neutron star properties and the equation of state for its core

Comments: 9 pages, 8 figures, submitted to Astronomy and Astrophysics

Submitted: **2016-11-04**

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.

[12] oai:arXiv.org:1608.08091 [pdf] - 1475519

Cooling of neutron stars with diffusive envelopes

Comments: 8 pages, 4 figures, MNRAS, accepted

Submitted: **2016-08-26**

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.

[13] oai:arXiv.org:1607.05265 [pdf] - 1447876

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

Comments: 11 pages, 4 figures, Astrophysics and Space Science, accepted

Submitted: **2016-07-16**

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.

[14] oai:arXiv.org:1604.01944 [pdf] - 1494668

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

Comments: Submitted to Phys. Rev. C. Supplemental material not included

Submitted: **2016-04-07**

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]

[15] oai:arXiv.org:1601.05368 [pdf] - 1392805

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

Comments: Matches published version. Minor modifications and reference added

Submitted: **2016-01-20**, last modified: **2016-03-29**

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.

[16] oai:arXiv.org:1409.1120 [pdf] - 1351406

Progenitor neutron stars of the lightest and heaviest millisecond
pulsars

Comments: 13 pages, 10 figures, accepted in A&A

Submitted: **2014-09-03**, last modified: **2015-12-31**

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.

[17] oai:arXiv.org:1506.08645 [pdf] - 1239015

Energy bursts from deconfinement in high-mass twin stars

Comments: 4 pages, 2 figures

Submitted: **2015-06-29**

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.

[18] oai:arXiv.org:1408.3052 [pdf] - 1047933

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

Comments: 10 pages, 10 figures, published in A&A

Submitted: **2014-08-13**, last modified: **2015-05-23**

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.

[19] oai:arXiv.org:1504.04537 [pdf] - 1055963

Neutron drip transition in accreting and nonaccreting neutron star
crusts

Comments: 24 pages, accepted for publication in Physical Review C

Submitted: **2015-04-17**

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 $\rho_\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.

[20] oai:arXiv.org:1408.6526 [pdf] - 864022

The Large Observatory For x-ray Timing

Comments: Proc. SPIE 9144, Space Telescopes and Instrumentation 2014:
Ultraviolet to Gamma Ray, 91442T

Submitted: **2014-08-27**, last modified: **2014-08-29**

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.

[21] oai:arXiv.org:1402.6911 [pdf] - 1203614

EMMI Rapid Reaction Task Force Meeting on 'Quark Matter in Compact Star'

Comments: 31 pages, summary of the EMMI Rapid Reaction Task Force on "Quark
Matter in Compact Stars", October 7-10, 2013, FIAS, Goethe University,
Frankfurt, Germany

Submitted: **2014-02-27**

The recent measurement of two solar mass pulsars has initiated an intense
discussion on its impact on our understanding of the high-density matter in the
cores of neutron stars. A task force meeting was held from October 7-10, 2013
at the Frankfurt Institute for Advanced Studies to address the presence of
quark matter in these massive stars. During this meeting, the recent
oservational astrophysical data and heavy-ion data was reviewed. The
possibility of pure quark stars, hybrid stars and the nature of the QCD phase
transition were discussed and their observational signals delineated.

[22] oai:arXiv.org:1401.2827 [pdf] - 791880

Physics input for modelling superfluid neutron stars with hyperon cores

Comments: 21 pages, 8 figures, 10 tables, accepted for publication in MNRAS

Submitted: **2014-01-13**

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
$\sigma\omega\rho\phi\sigma^\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.

[23] oai:arXiv.org:1311.4278 [pdf] - 1180774

Transport coefficients of nuclear matter in neutron star cores

Comments: 38 pages, 17 figures, accepted in PRC

Submitted: **2013-11-18**

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.

[24] oai:arXiv.org:1307.3995 [pdf] - 747588

On the Maximum Mass of Neutron Stars

Comments: 25 pages, 4 figures; revised version with additions, corrections and
updated references

Submitted: **2013-07-15**, last modified: **2013-11-18**

One of the most intringuing questions about neutron stars concerns their
maximum mass. The answer is intimately related to the properties of matter at
densities far beyond that found in heavy atomic nuclei. The current view on the
internal constitution of neutron stars and on their maximum mass, both from
theoretical and observational studies, are briefly reviewed.

[25] oai:arXiv.org:1211.6236 [pdf] - 1158074

Disk accretion onto a magnetized star

Comments: 7 pages, 2 figures

Submitted: **2012-11-27**

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.

[26] oai:arXiv.org:1211.1231 [pdf] - 778041

Maximum mass of neutron stars and strange neutron-star cores

Comments: 8 pages, 10 figures, submitted to A&A

Submitted: **2012-11-06**

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.

[27] oai:arXiv.org:1210.8331 [pdf] - 1152611

Formation of millisecond pulsars - NS initial mass and EOS constraints

Comments: 4 pages, 2 figures; proceedings of IAUS 290 "Feeding Compact Objects:
Accretion on All Scales", C. M. Zhang, T. Belloni, M. Mendez & S. N. Zhang
(eds.)

Submitted: **2012-10-31**

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.

[28] oai:arXiv.org:1210.0682 [pdf] - 1151796

Lev Landau and the conception of neutron stars

Comments: 12 pages, 4 figures

Submitted: **2012-10-02**

We review the history of neutron star physics in the 1930s that is related to
L. Landau. According to recollections of Rosenfeld (1974, Proc. 16th Solvay
Conference on Physics, p. 174), Landau improvised the concept of neutron stars
in a discussion with Bohr and Rosenfeld just after the news of the discovery of
the neutron reached Copenhagen in February 1932. We present arguments that the
discussion took place in March 1931, before the discovery of the neutron, and
that they in fact discussed the paper written by Landau in Zurich in February
1931 but not published until February 1932 (Phys. Z. Sowjetunion, 1, 285). In
his paper Landau mentioned the possible existence of dense stars which look
like one giant nucleus; this can be regarded as an early theoretical prediction
or anticipation of neutron stars, prior to the discovery of the neutron. The
coincidence of the dates of the neutron's discovery and the paper's publication
has led to an erroneous association of the paper with the discovery of the
neutron. In passing, we outline the contribution of Landau to the theory of
white dwarfs and to the hypothesis of stars with neutron cores.

[29] oai:arXiv.org:1106.2432 [pdf] - 450591

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

Comments: 11 pages, 10 figures; A&A accepted

Submitted: **2011-06-13**, last modified: **2011-12-13**

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.

[30] oai:arXiv.org:1109.1179 [pdf] - 450136

Compression of matter in the center of accreting neutron stars

Comments: 6 pages, 7 figures, Astronomy and Astrophysics in print

Submitted: **2011-09-06**, last modified: **2011-12-11**

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.

[31] oai:arXiv.org:1111.6942 [pdf] - 805843

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

Comments: 7 pages, 5 figures. Main results of this paper were already presented
at the MODE-SNR-PWN Workshop in Bordeaux, France, November 15-17, 2010, and
in a poster at the CompStar 2011 Workshop in Catania, Italy, May 9-12, 2011.
The paper is being submitted to Astronomy & Astrophysics

Submitted: **2011-11-29**

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.

[32] oai:arXiv.org:1104.0385 [pdf] - 1053117

Formation scenarios and mass-radius relation for neutron stars

Comments: 6 pages, 4 figures. Accepted for publication in Astronomy and
Astrophysics

Submitted: **2011-04-03**

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.

[33] oai:arXiv.org:1008.0384 [pdf] - 371968

Approximate analytic expressions for circular orbits around rapidly
rotating compact stars

Comments: 6 pages, 6 figures, A&A in press

Submitted: **2010-08-02**, last modified: **2010-09-27**

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.

[34] oai:arXiv.org:0907.0010 [pdf] - 25830

The relativistic entrainment matrix of a superfluid nucleon-hyperon
mixture. II. Effect of finite temperatures

Comments: 10 pages, 1 figure, accepted for publication in Phys. Rev. C

Submitted: **2009-06-30**

We calculate the important quantity of superfluid hydrodynamics, the
relativistic entrainment matrix for a nucleon-hyperon mixture at arbitrary
temperature. In the nonrelativistic limit this matrix is also termed the
Andreev-Bashkin or mass-density matrix. Our results can be useful for modeling
the pulsations of massive neutron stars with superfluid nucleon-hyperon cores
and for studies of the kinetic properties of superfluid baryon matter.

[35] oai:arXiv.org:0904.3467 [pdf] - 23577

The relativistic entrainment matrix of a superfluid nucleon-hyperon
mixture at zero temperature

Comments: 14 pages, 4 figures, accepted for publication in Phys. Rev. C

Submitted: **2009-04-22**

We calculate the relativistic entrainment matrix Y_ik at zero temperature for
nucleon-hyperon mixture composed of neutrons, protons, Lambda and Sigma^-
hyperons, as well as of electrons and muons. This matrix is analogous to the
entrainment matrix (also termed mass-density matrix or Andreev-Bashkin matrix)
of non-relativistic theory. It is an important ingredient for modelling the
pulsations of massive neutron stars with superfluid nucleon-hyperon cores. The
calculation is done in the frame of the relativistic Landau Fermi-liquid theory
generalized to the case of superfluid mixtures; the matrix Y_ik is expressed
through the Landau parameters of nucleon-hyperon matter. The results are
illustrated with a particular example of the sigma-omega-rho mean-field model
with scalar self-interactions. Using this model we calculate the matrix Y_ik
and the Landau parameters. We also analyze stability of the ground state of
nucleon-hyperon matter with respect to small perturbations.

[36] oai:arXiv.org:0901.3819 [pdf] - 1001355

Dynamic migration of rotating neutron stars due to a phase transition
instability

Comments: 23 pages, 18 figures, minor modifications

Submitted: **2009-01-24**, last modified: **2009-04-09**

Using numerical simulations based on solving the general relativistic
hydrodynamic equations, we study the dynamics of a phase transition in the
dense core of isolated rotating neutron stars, triggered by the back bending
instability reached via angular momentum loss. In particular, we investigate
the dynamics of a migration from an unstable configuration into a stable one,
which leads to a mini-collapse of the neutron star and excites sizeable
pulsations in its bulk until it acquires a new stable equilibrium state. We
consider equations of state with softening at high densities, a simple analytic
one with a mixed hadron-quark phase in an intermediate pressure interval and
pure quark matter at very high densities, and a microphysical one that has a
first-order phase transition, originating from kaon condensation. Although the
marginally stable initial models are rigidly rotating, we observe that during
the collapse (albeit little) differential rotation is created. We analyze the
emission of gravitational radiation, which in some models is amplified by mode
resonance effects, and assess its prospective detectability by interferometric
detectors. We expect that the most favorable conditions for dynamic migration
exist in very young magnetars. We find that the damping of the post-migration
pulsations strongly depends on the character of the equation of state
softening. The damping of pulsations in the models with the microphysical
equation of state is caused by dissipation associated with matter flowing
through the density jump at the edge of the dense core. If at work, this
mechanism dominates over all other types of dissipation, like bulk viscosity in
the exotic-phase core, gravitational radiation damping, or numerical viscosity.

[37] oai:arXiv.org:0811.1090 [pdf] - 18288

Luminosity of a quark star undergoing torsional oscillations and the
problem of gamma ray bursts

Comments: 16 pages, accepted by Astronomy and Astrophysics

Submitted: **2008-11-07**, last modified: **2009-01-22**

We discuss whether the winding-up of the magnetic field by differential
rotation in a new-born quark star can produce a sufficiently-high, energy,
emission rate of sufficiently long duration to explain long gamma-ray bursts.
In the context of magnetohydrodynamics, we study the torsional oscillations
and energy extraction from a new-born, hot, differentially rotating quark star.
The new-born compact star is a rapid rotator that produces a relativistic,
leptonic wind. The star's torsional oscillation modulates this wind emission
considerably when it is odd and of sufficient amplitude, which is relatively
easy to reach. Odd oscillations may occur just after the formation of a quark
star. Other asymmetries can cause similar effects. The buoyancy of wound-up
magnetic fields is inhibited, or its effects are limited, by a variety of
different mechanisms. Direct electromagnetic emission by the torsional
oscillation in either an outside vacuum or the leptonic wind surrounding the
compact object is found to be insignificant. In contrast, the twist given to
the outer magnetic field by an odd torsional oscillation is generally
sufficient to open the star's magnetosphere. The Poynting emission of the star
in its leptonic environment is then radiated from all of its surface and is
enhanced considerably during these open episodes, tapping at the bulk
rotational energy of the star. This results in intense energy shedding in the
first tens of minutes after the collapse of magnetized quark stars with an
initial poloidal field of order of 10**14 Gauss, sufficient to explain long
gamma-ray bursts.

[38] oai:arXiv.org:0901.1268 [pdf] - 1001290

Keplerian frequency of uniformly rotating neutron stars and quark stars

Comments: 6 pages, 6 color figures, submitted to A&A

Submitted: **2009-01-09**

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.

[39] oai:arXiv.org:0812.3955 [pdf] - 1001229

Physics of Neutron Star Crusts

Comments: 182 pages, published version available at
<http://www.livingreviews.org/lrr-2008-10>

Submitted: **2008-12-20**

The physics of neutron star crusts is vast, involving many different research
fields, from nuclear and condensed matter physics to general relativity. This
review summarizes the progress, which has been achieved over the last few
years, in modeling neutron star crusts, both at the microscopic and macroscopic
levels. The confrontation of these theoretical models with observations is also
briefly discussed.

[40] oai:arXiv.org:0805.1814 [pdf] - 12565

Deformation and crustal rigidity of rotating neutron stars

Comments: Accepted by A&A. 11 pages, 11 figures, 1 table. Title changed.
Section on precession of RX J0720-3125 removed due to the withdrawal of
observational evidence by Haberl et al. Section 6 substantially expanded.
Minor addditions to Sections 2 and 10, top margin corrected

Submitted: **2008-05-13**, last modified: **2008-10-02**

We calculate parameters A and B of the Baym-Pines model of the hydro-elastic
equilibrium of rotating neutron stars. Parameter A determines the energy
increase of a non-rotating star due to a quadrupolar deformation of its shape.
Parameter B determines residual quadrupolar deformation due to the crustal
shear strain, in a neutron star that spun-down to a non-rotating state. The
calculations of A are based on precise numerical 2-D calculations for rotating
neutron stars with realistic equations of state (EOSs) of dense matter. An
approximate, but quite precise, formula for B is used, which allows us to
separate the contribution of the crust from the dependence on the stellar mass
M and radius R. The elastic shear strain distribution within the crust is
modeled following Cutler et al. (2003). A(M) and B(M) are calculated for
0.2Msun < M < 0.9M_max for seven EOSs of neutron star core, combined with
several crust models. A standard formula based on the incompressible fluid
model is shown to severely underestimate the value of A. For M<0.7Msun A(M) is
nearly EOS independent and are given (within a few percent) by a universal
formula A=3.87(M/Msun)^7/3 [10^53 erg]. We derive the scaling of B with respect
to R and M, valid also for a thick crust. We show that B for accreted crust
strongly depends on pycnonuclear fusions at rho>10^12 g/cm^3.

[41] oai:arXiv.org:0805.1820 [pdf] - 12566

Fast rotation of neutron stars and equation of state of dense matter

Comments: 7 pages, 7 figures, proceedings of the "Jean-Pierre Lasota, X-ray
binaries, accretion disks and compact stars" conference, edited by M.
Abramowicz and O. Straub, New Astronomy Reviews (Elsevier 2008), in press

Submitted: **2008-05-13**

Fast rotation of compact stars (at submillisecond period) and, in particular,
their stability, are sensitive to the equation of state (EOS) of dense matter.
Recent observations of XTE J1739-285 suggest that it contains a neutron star
rotating at 1122 Hz (Kaaret et al. 2007). At such rotational frequency the
effects of rotation on star's structure are significant. We study the interplay
of fast rotation, EOS and gravitational mass of a submillisecond pulsar. We
discuss the EOS dependence of spin-up to a submillisecond period, via mass
accretion from a disk in a low-mass X-ray binary.

[42] oai:arXiv.org:0707.4614 [pdf] - 3565

Thermal conductivity of ions in a neutron star envelope

Comments: 12 pages, 5 figures; to appear in MNRAS

Submitted: **2007-07-31**, last modified: **2008-03-24**

We analyze the thermal conductivity of ions (equivalent to the conductivity
of phonons in crystalline matter) in a neutron star envelope.
We calculate the ion/phonon thermal conductivity in a crystal of atomic
nuclei using variational formalism and performing momentum-space integration by
Monte Carlo method. We take into account phonon-phonon and phonon-electron
scattering mechanisms and show that phonon-electron scattering dominates at not
too low densities. We extract the ion thermal conductivity in ion liquid or gas
from literature.
Numerical values of the ion/phonon conductivity are approximated by
analytical expressions, valid for T>10^5 K and 10^5 g cm^-3 < \rho < 10^14 g
cm^-3. Typical magnetic fields B~10^12 G in neutron star envelopes do not
affect this conductivity although they strongly reduce the electron thermal
conductivity across the magnetic field. The ion thermal conductivity remains
much smaller than the electron conductivity along the magnetic field. However,
in the outer neutron star envelope it can be larger than the electron
conductivity across the field, that is important for heat transport across
magnetic field lines in cooling neutron stars. The ion conductivity can greatly
reduce the anisotropy of heat conduction in outer envelopes of magnetized
neutron stars.

[43] oai:arXiv.org:0710.5010 [pdf] - 1579263

EOS of dense matter and fast rotation of neutron stars

Comments: 8 pages, 5 eps figures, proceedings of the International Symposium on
Exotic States of Nuclear Matter, EXOCT07, to be published by World Scientific
Publishing Co

Submitted: **2007-10-26**

Recent observations of XTE J1739-285 suggest that it contains a neutron star
rotating at 1122 Hz. Such rotational frequency would be the first for which the
effects of rotation are significant. We study the consequences of very fast
rotating neutron stars for the potentially observable quantities as stellar
mass and pulsar period.

[44] oai:arXiv.org:0710.1528 [pdf] - 1579259

The energy release--stellar angular momentum independence in rotating
compact stars undergoing first-order phase transitions

Comments: 4 pages, 2 figures, proceedings of the International Symposium on
Exotic States of Nuclear Matter, EXOCT07, to be published by World Scientific
Publishing Co

Submitted: **2007-10-08**

We present the general relativistic calculation of the energy release
associated with a first order phase transition (PT) at the center of a rotating
neutron star (NS). The energy release, E_rel, is equal to the difference in
mass-energies between the initial (normal) phase configuration and the final
configuration containing a superdense matter core, assuming constant total
baryon number and the angular momentum. The calculations are performed with the
use of precise pseudo-spectral 2-D numerical code; the polytropic equations of
state (EOS) as well as realistic EOSs (Skyrme interactions, Mean Field Theory
kaon condensate) are used. The results are obtained for a broad range of
metastability of initial configuration and size of the new superdense phase
core in the final configuration. For a fixed ``overpressure'', dP, defined as
the relative excess of central pressure of a collapsing metastable star over
the pressure of the equilibrium first-order PT, the energy release up to
numerical accuracy does not depend on the stellar angular momentum and
coincides with that for nonrotating stars with the same dP. When the equatorial
radius of the superdense phase core is much smaller than the equatorial radius
of the star, analytical expressions for the E_rel can be obtained: E_rel is
proportional to dP^2.5 for small dP. At higher dP, the results of 1-D
calculations of E_rel(dP) for non-rotating stars reproduce with very high
precision exact 2-D results for fast-rotating stars. The energy release-angular
momentum independence for a given overpressure holds also for the so-called
``strong'' PTs (that destabilise the star against the axi-symmetric
perturbations), as well as for PTs with ``jumping'' over the energy barrier.

[45] oai:arXiv.org:0708.3996 [pdf] - 4368

Models of crustal heating in accreting neutron stars

Comments: 8 pages, 5 figures, Submitted to A&A

Submitted: **2007-08-29**

Heating associated with non-equilibrium nuclear reactions in accreting
neutron-star crusts is reconsidered, taking into account suppression of
neutrino losses demonstrated recently by Gupta et al. Two initial compositions
of the nuclear burning ashes, A=56 and A=106, are considered. Dependence of the
integrated crustal heating on uncertainties plaguing pycnonuclear reaction
models is studied.
One-component plasma approximation is used, with compressible liquid-drop
model of Mackie and Baym to describe nuclei. Evolution of a crust shell is
followed from 10^8 g/cm^3 to 10^(13.6) g/cm^3
The integrated heating in the outer crust agrees nicely with results of
self-considtent multicomponent plasma simulations of Gupta et al.; their
results fall between our curves obtained for A=56 and A=106. Total crustal heat
per one accreted nucleon ranges between 1.5 MeV to 1.9 MeV for A=106 and A=56,
respectively. The value of total crustal heat per nucleon depends weakly on the
presence of pycnonuclear reactions at densities 10^(12)-10^(13) g/cm^3.
Remarkable insensitivity of the total crustal heat on the details of the
distribution of nuclear processes in accreted crust is explained.

[46] oai:arXiv.org:0708.0086 [pdf] - 3609

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

Comments: 5 pages, 2 figures, accepted by MNRAS. In v.2, two references added
and typos corrected

Submitted: **2007-08-01**, last modified: **2007-08-24**

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.

[47] oai:arXiv.org:0707.3691 [pdf] - 3404

Strong first-order phase transition in a rotating neutron star core and
the associated energy release

Comments: 9 pages, 12 figures, submitted to A&A

Submitted: **2007-07-25**

We calculate the energy release associated with a strong first-order phase
transition, from normal phase N to an "exotic" superdense phase S, in a
rotating neutron star. Such a phase transition, accompanied by a density jump
rho_N --> rho_S, is characterized by rho_S/rho_N > 3/2(1+P_0/rho_N c^2), where
P_0 is the pressure, at which phase transition occurs. Configurations with
small S-phase cores are then unstable and collapse into stars with large
S-phase cores. The energy release is equal to the difference in mass-energies
between the initial (normal) configuration and the final configuration
containing an S-phase core, total stellar baryon mass and angular momentum
being kept constant. The calculations of the energy release are based on
precise numerical 2-D calculations. Polytropic equations of state (EOSs) as
well as realistic EOS with strong first-order phase transition due to kaon
condensation are used. For polytropic EOSs, a large parameter space is studied.
For a fixed "overpressure", dP, defined as the relative excess of central
pressure of collapsing metastable star over the pressure of equilibrium
first-order phase transition, the energy release E_rel does not depend on the
stellar angular momentum. It coincides with that for nonrotating stars with the
same dP. Therefore, results of 1-D calculations of E_rel(dP) for non-rotating
stars can be used to predict, with very high precision, the outcome of much
harder to perform 2-D calculations for rotating stars with the same dP. This
result holds also for dP_min < dP < 0, corresponding to phase transitions with
climbing over the energy barrier separating metastable N-phase configurations
from those with an S-phase core. Such phase transitions could be realized in
the cores of newly born, hot, pulsating neutron stars.

[48] oai:arXiv.org:0705.4594 [pdf] - 1763

Two branches of neutron stars - reconciling a 2M_sun pulsar and SN1987A

Comments: 4 pages, 5 figures, submitted to A&A Letters

Submitted: **2007-05-31**

The analysis of SN1987A led Brown and Bethe (1995) to conclusion, that the
maximum mass of cold neutron stars is low, M_max ~ 1.5M_sun. Such a low M_max,
due to a kaon condensation in the stellar core, implies collapse of a too
massive deleptonized protoneutron star into a black hole. This would naturally
explain the lack of a neutron star in the SN1987A remnant. On the other hand,
recent evaluation of mass of PSR J0751+1807 gives M_max > 2M_sun. This
contradicts the original Bethe-Brown model, but can be reconciled within
scenarios proposed in the present Letter. We consider two types of dense matter
models with high-density softening, due to a transition from a non-strange
N-phase of matter to a strangeness carrying phase S: kaon condensation and
deconfinement of quarks. Two scenarios of neutron star formation in stellar
core collapse are considered. In the first scenario, realized in sufficiently
hot and dense supernova cores, nucleation of an S-phase is sufficiently rapid
so as to form an S-phase core, and implying M_max = M^S_max =~ 1.5M_sun. In the
second scenario, nucleation of the S-phase at neutron star birth is too slow to
materialize, and the star becomes cold without forming an S-phase core. Then,
stellar mass can increase via accretion, until central density \rho_crit is
reached, and the S phase forms. This N branch of neutron stars ends at
M=M_crit. We select several models of N-phase satifying the necessary condition
M^N_max > 2M_sun and combine them with models of kaon condensation and quark
deconfinement. For kaon condensation, we get M_crit =~ M^S_max =~ 1.5M_sun,
which is ruled out by PSR J0751+1807. On the contrary, for the EOSs with quark
deconfinement we get M_crit =~ M^N_max > 2M_sun, which reconciles SN1987A and
PSR J0751+1807.

[49] oai:arXiv.org:astro-ph/0612216 [pdf] - 87554

Rotation at 1122 Hz and the neutron star structure

Comments: 4 pages, 3 figures, new paragraph on the constraints from the spin-up
by accretion added. Accepted for publication in A&A Letters

Submitted: **2006-12-08**, last modified: **2007-01-30**

Recent observations of XTE J1739-285 suggest that it contains a neutron star
rotating at 1122 Hz. Such rotation imposes bounds on the structure of neutron
star in XTE J1739-285. These bounds may be used to constrain poorly known
equation of state of dense matter. One-parameter families of stationary
configurations rotating rigidly at 1122 Hz are constructed, using a precise 2-D
code solving Einstein equations. Hydrostatic equilibrium solutions are tested
for stability with respect to axi-symmetric perturbations. A set of ten diverse
EOSs of neutron stars is considered. Hypothetical strange stars are also
studied. For each EOS, the family of possible neutron star models is limited by
the mass shedding limit, corresponding to maximum allowed equatorial radius,
R_max, and by the instability with respect to the axi-symmetric perturbations,
reached at the minimum allowed equatorial radius, R_min. We get R_min \simeq
10-13km, and R_max \simeq 16-18km, with allowed mass 1.4-2.3 M_\odot. Allowed
stars with hyperonic or exotic-phase core are supramassive and have a very
narrow mass range. Quark star with accreted crust might be allowed, provided
such a model is able to reproduce X-ray bursts from XTE J1739-285.

[50] oai:arXiv.org:astro-ph/0701258 [pdf] - 260846

Long Gamma Ray Bursts from Quark Stars

Comments: To be published in the proceedings of the conference "SWIFT and GRBs:
Unveiling the Relativistic Universe", Venice, June 5-9, 2006. To appear in
"Il Nuovo Cimento"

Submitted: **2007-01-09**

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.

[51] oai:arXiv.org:astro-ph/0611115 [pdf] - 86501

Nucleon superfluidity versus thermal states of isolated and transiently
accreting neutron stars

Comments: 8 pages, 7 Figs, to be published in the proc. of "Isolated Neutron
Stars: From the Interior to the Surface", eds. D.Page, S.Zane, R.Turolla

Submitted: **2006-11-03**

The properties of superdense matter in neutron star (NS) cores control NS
thermal states by affecting the efficiency of neutrino emission from NS
interiors. To probe these properties we confront the theory of thermal
evolution of NSs with observations of their thermal radiation. Our
observational basis includes cooling isolated NSs (INSs) and NSs in quiescent
states of soft X-ray transients (SXTs). We find that the data on SXTs support
the conclusions obtained from the analysis of INSs: strong proton superfluidity
with T_{cp,max} >= 10^9 K should be present, while mild neutron superfluidity
with T_{cn,max} =(2*10^8 -- 2*10^9) K is ruled out in the outer NS core. Here
T_{cn,max} and T_{cp,max} are the maximum values of the density dependent
critical temperatures of neutrons and protons. The data on SXTs suggest also
that: (i) cooling of massive NSs is enhanced by neutrino emission more powerful
than the emission due to Cooper pairing of neutrons; (ii) mild neutron
superfluidity, if available, might be present only in inner cores of massive
NSs. In the latter case SXTs would exhibit dichotomy, i.e. very similar SXTs
may evolve to very different thermal states.

[52] oai:arXiv.org:astro-ph/0610549 [pdf] - 85973

Equation of state of neutron star cores and spin down of isolated
pulsars

Comments: 7 pages, 11 figures, to appear in Astrophysics and Space Science, in
the proceedings of "Isolated Neutron Stars: from the Interior to the
Surface", edited by D. Page, R. Turolla and S. Zane

Submitted: **2006-10-18**

We study possible impact of a softening of the equation of state by a phase
transition, or appearance of hyperons, on the spin evolution of of isolated
pulsars. Numerical simulations are performed using exact 2-D simulations in
general relativity. The equation of state of dense matter at supranuclear
densities is poorly known. Therefore, the accent is put on the general
correlations between evolution and equation of state, and mathematical
strictness. General conjectures referring to the structure of the one-parameter
families of stationary configurations are formulated. The interplay of the back
bending phenomenon and stability with respect to axisymmetric perturbations is
described. Changes of pulsar parameters in a corequake following instability
are discussed, for a broad choice of phase transitions predicted by different
theories of dense matter. The energy release in a corequake, at a given initial
pressure, is shown to be independent of the angular momentum of collapsing
configuration. This result holds for various types of phases transition, with
and without metastability. We critically review observations of pulsars that
could be relevant for the detection of the signatures of the phase transition
in neutron star cores.

[53] oai:arXiv.org:astro-ph/0610188 [pdf] - 85612

Energy release associated with a first-order phase transition in a
rotating neutron star core

Comments: 8 pages, 8 figures, submitted to A&A

Submitted: **2006-10-06**

We calculate energy release associated with a first order phase transition at
the center of a rotating neutron star. The results are based on precise
numerical 2-D calculations, in which both the polytropic equations of state
(EOS) as well as realistic EOS of the normal phase are used. Presented results
are obtained for a broad range of metastability of initial configuration and
size of the new superdense phase core in the final configuration. For small
radii of the superdense phase core analytical expressions for the energy
release are obtained. For a fixed "overpressure" dP (the relative excess of
central pressure of collapsing metastable star over the pressure of equilibrium
first-order phase transition) the energy release remarkably does not depend on
the stellar angular momentum and coincides with that for nonrotating stars with
the same dP. The energy release is proportional to dP^2.5 for small dPs, when
sufficiently precise brute force 2-D numerical calculations are out of
question. For higher dPs, results of 1-D calculations of energy release for
non-rotating stars are shown to reproduce, with very high precision, the exact
2-D results for rotating stars.

[54] oai:arXiv.org:nucl-th/0603018 [pdf] - 118778

Entrainment parameters in cold superfluid neutron star core

Comments: 10 pages, 7 figures. Revised version

Submitted: **2006-03-07**, last modified: **2006-09-13**

Hydrodynamical simulations of neutron star cores, based on a two fluid
description in terms of a neutron-proton superfluid mixture, require the
knowledge of the Andreev-Bashkin entrainment matrix which relates the momentum
of one constituent to the currents of both constituents. This matrix is derived
for arbitrary nuclear asymmetry at zero temperature and in the limits of small
relative currents in the framework of the energy density functional theory. The
Skyrme energy density functional is considered as a particular case. General
analytic formulae for the entrainment parameters and various corresponding
effective masses are obtained. These formulae are applied to the liquid core of
a neutron star, composed of an homogeneous plasma of nucleons, electrons and
possibly muons in beta equilibrium.

[55] oai:arXiv.org:astro-ph/0408083 [pdf] - 260559

Entrainment coefficient and effective mass for conduction neutrons in
neutron star crust: II Macroscopic treatment

Comments: 21 pages Latex. Part II of article whose Part I (Simple microscopic
models) is given by nucl-th/0402057. New version extended to include figures

Submitted: **2004-08-04**, last modified: **2006-01-30**

Phenomena such as pulsar frequency glitches are believed to be attributable
to differential rotation of a current of ``free'' superfluid neutrons at
densities above the ``drip'' threshold in the ionic crust of a neutron star.
Such relative flow is shown to be locally describable by adaption of a
canonical two fluid treatment that emphasizes the role of the momentum
covectors constructed by differentiation of action with respect to the
currents, with allowance for stratification whereby the ionic number current
may be conserved even when the ionic charge number Z is altered by beta
processes. It is demonstrated that the gauge freedom to make different choices
of the chemical basis determining which neutrons are counted as ``free'' does
not affect their ``superfluid'' momentum covector, which must locally have the
form of a gradient (though it does affect the ``normal'' momentum covector
characterising the protons and those neutrons that are considered to be
``confined'' in the nuclei). It is shown how the effect of ``entrainment''
(whereby the momentum directions deviate from those of the currents) is
controlled by the (gauge independent) mobility coefficient K, estimated in
recent microscopical quantum mechanical investigations, which suggest that the
corresponding (gauge dependent) ``effective mass'' m* of the free neutrons can
become very large in some layers. The relation between this treatment of the
crust layers and related work (using different definitions of ``effective
mass'') intended for the deeper core layers is discussed.

[56] oai:arXiv.org:astro-ph/0509806 [pdf] - 76387

Phase transitions in rotating neutron stars cores: back bending,
stability, corequakes and pulsar timing

Comments: 13 pages, 15 figures, accepted by A&A

Submitted: **2005-09-27**, last modified: **2006-01-25**

The back-bending phenomenon for compact stars is studied by means of
analytical equations of state, for both constant-pressure phase transitions and
the transitions through the mixed-phase region. We restrict ourselves to the
case of normal rotating configurations, with baryon mass below the maximum
allowable baryon mass for non-rotating stars. We use high-precision 2-D
multi-domain spectral code LORENE to search the parameter space for possible
instability regions, and possible changes in the stability character of
rotating stars with phase transitions in their cores. Conditions on the density
jump in constant-pressure phase transitions, leading to the existence of the
unstable segments in the evolutionary sequences of spinning down isolated
normal neutron stars, are derived. Conjectures concerning the existence of two
disjoint families of non-rotating and rotating stationary configurations of
neutron stars are formulated. Particular case of EOSs leading to marginal
instability of static and rotating configurations is also studied: marginal
instability point in non-rotating configurations continues to exist in all
evolutionary spin-down tracks. The fate of rotating stars entering the region
of instability is discussed. The change in radius, energy release, and spin-up
associated with the corequake in rotating neutron star, triggered by the
instability, are calculated. The energy release is found to be very weakly
dependent on the angular momentum of collapsing star.

[57] oai:arXiv.org:astro-ph/0504572 [pdf] - 72669

Non-equilibrium beta processes in superfluid neutron star cores

Comments: 11 pages, 6 Figs., published, minor changes, subroutines can be found
on line at http://luth2.obspm.fr/~etu/villain/Micro/Resolution.html

Submitted: **2005-04-26**, last modified: **2005-12-05**

The influence of nucleons superfluidity on the beta relaxation time of
degenerate neutron star cores, composed of neutrons, protons and electrons, is
investigated. We numerically calculate the implied reduction factors for both
direct and modified Urca reactions, with isotropic pairing of protons or
anisotropic pairing of neutrons. We find that due to the non-zero value of the
temperature and/or to the vanishing of anisotropic gaps in some directions of
the phase-space, superfluidity does not always completely inhibit beta
relaxation, allowing for some reactions if the superfluid gap amplitude is not
too large in respect to both the typical thermal energy and the chemical
potential mismatch. We even observe that if the ratio between the critical
temperature and the actual temperature is very small, a suprathermal regime is
reached for which superfluidity is almost irrelevant. On the contrary, if the
gap is large enough, the composition of the nuclear matter can stay frozen for
very long durations, unless the departure from beta equilibrium is at least as
important as the gap amplitude. These results are crucial for precise
estimation of the superfluidity effect on the cooling/slowing-down of pulsars
and we provide online subroutines to be implemented in codes for simulating
such evolutions.

[58] oai:arXiv.org:astro-ph/0508105 [pdf] - 74994

Constraints on the dense matter EOS from the measurements of PSR
J0737-3039A moment of inertia and PSR J0751+1807 mass

Comments: 6 pages, 3 figures, accepted for publication in MNRAS Main Journal

Submitted: **2005-08-03**, last modified: **2005-09-13**

The moment of inertia of the pulsar A in the neutron star binary J0737-3039
will soon be measurable through detailed measurements of the periastron
advance. We present the calculation of the moment of inertia of neutron stars
with the masses of the components of the binary J0737-3039 for a broad range of
equations of state of dense matter and discuss the implications of such
measurement for constraining the equation of state. An observational
determination of the moment of inertia of the pulsar A in J0737-3039 with the
accuracy of 10% shall narrow down considerably the range of viable equations of
state. We also show that limits on maximal mass of a neutron star provide a
complementary set of constraints on the properties of dense nuclear matter.

[59] oai:arXiv.org:astro-ph/0502297 [pdf] - 71088

Gamma-Ray Bursts from quark stars

Comments: 5 pages, (v3) MNRAS Letters, in press

Submitted: **2005-02-15**, last modified: **2005-08-27**

Long gamma-ray bursts (GRBs) are believed to be related to the explosion of
type Ic supernovae, which have been stripped of their hydrogen and helium
envelopes. There appear to be two types of these explosions: those which are
approximately spherical (GRB980425/1998bw), and which are associated with weak
bursts, and the classical GRBs which generate ultrarelativistic jets
(GRB030329/SN2003dh). If this bimodality is real Swift will provide a clear
evidence for it.
We propose that classical powerful GRBs, which generate ultrarelativistic
outflows, are a result of a formation of quark stars. Quark stars may provide
an additional energy for the explosion of SN Ic, but far more important is a
creation of a surface which acts as a membrane which cannot be penetrated by
baryons. A surface of a quark star allows only ultrarelativistic matter to
escape: photons, neutrinos, electron -- positron pairs and magnetic fields. The
formation of a quark star follows the initial core collapse in several minutes.
Possible evidence for this time delay is provided by BATSE precursors to GRBs,
as analyzed by Lazzati (2005).

[60] oai:arXiv.org:astro-ph/0406228 [pdf] - 65385

Effect of BCS pairing on entrainment in neutron superfluid current in
neutron star crust

Comments: 30 pages, revised published version

Submitted: **2004-06-09**, last modified: **2005-08-23**

The relative current density $n^i$ of conduction neutrons in a neutron star
crust beyond the neutron drip threshold can be expected to be related to the
corresponding particle momentum covector $p_i$ by a linear relation of the form
$n^i=K^{ij}p_j$ in terms of a physically well defined mobility tensor $K^{ij}$.
This result is describable as an ``entrainment'' whose effect - wherever the
crust lattice is isotropic - will simply be to change the ordinary neutron mass
m to an effective mass $m_\star$ such that in terms of the relevant number
density n of unconfined neutrons we shall have $K^{ij}=n/m_\star\gamma^{ij}$.
In a preceding analysis based on an independent particle treatment using Bloch
boundary conditions to obtain the distribution of energy $E_k$ and associated
group velocity $v_k^i=\partial E_k/\partial\hbar k_i$ as a function of
wavenumber $k_i$, it was shown that the mobility tensor would be given by
$K^{ij}\propto\int d^3 k v_k^i v_k^j\delta\{E_k-\mu\}$, where $\mu$ is the
Fermi energy. Using the approach due to Bogoliubov, it is shown here that the
effect of BCS pairing with an energy gap \Delta_F and corresponding
quasiparticle energy function \EE_k=\sqrt{(E_k-\mu)^2+\Delta_F^2}$ will just be
to replace the Dirac distributional integrand by a smoother distribution :
$K^{ij}\propto \int d^3 k v_k^i v_k^j \Delta_F^2 /\EE^3_k$. It is also shown
how the pairing condensation gives rise to superfluidity in the technical sense
of providing (meta) stability against resistive perturbations for a current
that is not too strong (its momentum $p_i$ must be small enough to give $2|p_i
v_k^i|<\EE_k^2/|E_k-\mu|$ for all modes). It is concluded that the prediction
of a very large effective mass enhancement in the middle layers of the crust
will not be significantly affected by the pairing mechanism.

[61] oai:arXiv.org:astro-ph/0508104 [pdf] - 74993

The entrainment matrix of a superfluid neutron-proton mixture at a
finite temperature

Comments: 19 pages, 1 figure

Submitted: **2005-08-03**

The entrainment matrix (also termed the Andreev-Bashkin matrix or the
mass-density matrix) for a neutron-proton mixture is derived at a finite
temperature in a neutron star core. The calculation is performed in the frame
of the Landau Fermi-liquid theory generalized to account for superfluidity of
nucleons. It is shown, that the temperature dependence of the entrainment
matrix is described by a universal function independent on an actual model of
nucleon-nucleon interaction employed. The results are presented in the form
convenient for their practical use. The entrainment matrix is important, e.g.,
in kinetics of superfluid nucleon mixtures or in studies of the dynamical
evolution of neutron stars (in particular, in the studies of star pulsations
and pulsar glitches).

[62] oai:arXiv.org:astro-ph/0502140 [pdf] - 70931

Spin-up of the hyperon-softened accreting neutron stars

Comments: 5 pages, 5 figures, accepted by Astronomy & Astrophysics

Submitted: **2005-02-07**, last modified: **2005-06-20**

We study the spin-up of the accreting neutron stars with a realistic
hyperon-softened equation of state. Using precise 2-D calculations we study the
evolutionary tracks of accreting neutron stars in the angular-momentum -
frequency plane. In contrast to the case of spinning-down solitary
radio-pulsars, where a strong back-bending behavior has been observed, we do
not see back-bending phenomenon in the accretion-powered spinning-up case. We
conclude that in the case of accretion-driven spin-up the back-bending is
strongly suppressed by the mass-increase effect accompanying the
angular-momentum increase.

[63] oai:arXiv.org:gr-qc/0407081 [pdf] - 111588

Inertial modes in stratified rotating neutron stars : An evolutionary
description

Comments: 34 pages, 24 figures, published version

Submitted: **2004-07-21**, last modified: **2005-04-20**

With (non-barotropic) equations of state valid even when the neutron, proton
and electron content of neutron star cores is not in beta equilibrium, we study
inertial and composition gravity modes of relativistic rotating neutron stars.
We solve the relativistic Euler equations in the time domain with a three
dimensional numerical code based on spectral methods, in the slow rotation,
relativistic Cowling and anelastic approximations. Principally, after a short
description of the gravity modes due to smooth composition gradients, we focus
our analysis on the question of how the inertial modes are affected by
non-barotropicity of the nuclear matter. In our study, the deviation with
respect to barotropicity results from the frozen composition of non-superfluid
matter composed of neutrons, protons and electrons, when beta equilibrium is
broken by millisecond oscillations. We show that already for moderatly fast
rotating stars the increasing coupling between polar and axial modes makes
those two cases less different than for very slowly rotating stars. In
addition, as we directly solve the Euler equations, without coupling only a few
number of spherical harmonics, we always found, for the models that we use, a
discrete spectrum for the $l = m = 2$ inertial mode. Finally, we find that, for
non-barotropic stars, the frequency of this mode, which is our main focus,
decreases in a non-negligible way, whereas the time dependence of the energy
transfer between polar and axial modes is substantially different due to the
existence of low-frequencies gravity modes.

[64] oai:arXiv.org:astro-ph/0502348 [pdf] - 880823

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

Comments: 9 pages, 6 figures, accepted for publication in MNRAS

Submitted: **2005-02-17**

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.

[65] oai:arXiv.org:nucl-th/0402057 [pdf] - 118684

Entrainment coefficient and effective mass for conduction neutrons in
neutron star crust: simple microscopic models

Comments: 27 pages, latex. Revised published version

Submitted: **2004-02-17**, last modified: **2005-01-26**

In the inner crust of a neutron star, at densities above the ``drip''
threshold, unbound ``conduction'' neutrons can move freely past through the
ionic lattice formed by the nuclei. The relative current density $n^i= n \bar
v^i$ of such conduction neutrons will be related to the corresponding mean
particle momentum $p_i$ by a proportionality relation of the form $n^i= {\cal
K}p^i$ in terms of a physically well defined mobility coefficient $\cal K$
whose value in this context has not been calculated before. Using methods from
ordinary solid state and nuclear physics, a simple quantum mechanical treatment
based on the independent particle approximation, is used here to formulate
$\cal K$ as the phase space integral of the relevant group velocity over the
neutron Fermi surface. The result can be described as an ``entrainment'' that
changes the ordinary neutron mass m to a macroscopic effective mass per neutron
that will be given -- subject to adoption of a convention specifying the
precise number density n of the neutrons that are considered to be ``free'' --
by $m_\star=n/{\cal K}$. The numerical evaluation of the mobility coefficient
is carried out for nuclear configurations of the ``lasagna'' and ``spaghetti''
type that may be relevant at the base of the crust. Extrapolation to the middle
layers of the inner crust leads to the unexpected prediction that $m_\star$
will become very large compared with m.

[66] oai:arXiv.org:astro-ph/0406234 [pdf] - 65391

Impact of the nuclear equation of state on the last orbits of binary
neutron stars

Comments: 11 pages, 10 EPS figures, minor revisions, accepted for publication
in A&A

Submitted: **2004-06-09**, last modified: **2004-10-18**

We present calculations of quasiequilibrium sequences of irrotational binary
neutron stars based on realistic equations of state (EOS) for the whole neutron
star interior. Three realistic nuclear EOSs of various softness and based on
different microscopic models have been joined with a recent realistic EOS of
the crust, giving in this way three different EOSs of neutron-star interior.
Computations of quasiequilibrium sequences are performed within the
Isenberg-Wilson-Mathews approximation to general relativity. For all
evolutionary sequences, the innermost stable circular orbit (ISCO) is found to
be given by mass-shedding limit (Roche lobe overflow). The EOS dependence on
the last orbits is found to be quite important: for two 1.35 M_sol neutron
stars, the gravitational wave frequency at the ISCO (which marks the end of the
inspiral phase) ranges from 800 Hz to 1200 Hz, depending upon the EOS. Detailed
comparisons with 3rd order post-Newtonian results for point-mass binaries
reveals a very good agreement until hydrodynamical effects (dominated by
high-order functions of frequency) become important, which occurs at a
frequency ranging from 500 Hz to 1050 Hz, depending upon the EOS.

[67] oai:arXiv.org:astro-ph/0408324 [pdf] - 66815

Analytical representations of unified equations of state of neutron-star
matter

Comments: 7 pages, 5 figures, 3 tables, accepted by A&A. In v.2, auxiliary fits
(15) and (16) are corrected

Submitted: **2004-08-18**, last modified: **2004-09-07**

Analytical representations are derived for two equations of state (EOSs) of
neutron-star matter: FPS and SLy. Each of these EOSs is unified, that is, it
describes the crust and the core of a neutron star using the same physical
model. Two versions of the EOS parametrization are considered. In the first
one, pressure and mass density are given as functions of the baryon density. In
the second version, pressure, mass density, and baryon density are given as
functions of the pseudo-enthalpy, which makes this representation particularly
useful for 2-D calculations of stationary rotating configurations of neutron
stars.

[68] oai:arXiv.org:astro-ph/0404165 [pdf] - 64073

Direct Urca Process in a Neutron Star Mantle

Comments: 7 pages, 3 figures, submitted to A&A

Submitted: **2004-04-07**

We show that the direct Urca process of neutrino emission is allowed in two
possible phases of nonspherical nuclei (inverse cylinders and inverse spheres)
in the mantle of a neutron star near the crust-core interface. The process is
open because neutrons and protons move in a periodic potential created by
inhomogeneous nuclear structures. In this way the nucleons acquire large
quasimomenta needed to satisfy momentum-conservation in the neutrino reaction.
The appropriate neutrino emissivity in a nonsuperfluid matter is about 2--3
orders of magnitude higher than the emissivity of the modified Urca process in
the stellar core. The process may noticeably accelerate the cooling of low-mass
neutron stars.

[69] oai:arXiv.org:astro-ph/0403550 [pdf] - 63741

Surface gravity of neutron stars and strange stars

Comments: Accepted by A&A

Submitted: **2004-03-23**

The upper bound on the value of the surface gravity, g_s, for neutron stars
with equations of state respecting v_sound <= c, is derived. This bound is
inversely proportional to the maximum allowable mass M_max, and it reads g_s <=
1.411 x 10^15 (M_sun/M_max) cm/s^2. It implies an absolute upper bound 7.4 x
10^14 cm/s^2 if one uses the 2sigma lower bound on the neutron mass measured
recently in 4U1700-37, 1.9 M_sun. A correlation between g_s and the compactness
parameter 2GM/Rc^2 for baryonic stars is analyzed. The properties of g_s of
strange quark stars and its upper bounds are discussed using the scaling
properties of the strange-star models.

[70] oai:arXiv.org:nucl-th/0403066 [pdf] - 118690

Static response of Fermi liquids with tensor interactions

Comments:

Submitted: **2004-03-22**

We use Landau's theory of a normal Fermi liquid to derive expressions for the
static response of a system with a general tensor interaction that conserves
the total spin and the total angular momentum of the quasiparticle-quasihole
pair. The magnetic susceptibility is calculated in detail, with the inclusion
of the center of mass tensor and cross vector terms in addition to the exchange
tensor one. We also introduce a new parametrization of the tensor Landau
parameters which significantly reduces the importance of high angular harmonic
contributions. For nuclear matter and neutron matter we find that the two most
important effects of the tensor interaction are to give a contribution from
multipair states and to renormalize the magnetic moments. Response to a weak
probe may be calculated using similar methods, replacing the magnetic moments
with the matrix elements of the weak charges.

[71] oai:arXiv.org:astro-ph/0311470 [pdf] - 61050

Hyperon softening of the EOS of dense matter and the spin evolution of
isolated neutron stars

Comments: 11 pages, 13 figures. Accepted for publication in A&A

Submitted: **2003-11-20**

The effect of the hyperon softening of the equation of state (EOS) of dense
matter on the spin evolution of isolated neutron stars is studied for a broad
set of hyperonic EOSs. We use a multidomain 2-D code based on a spectral
method, and show how important the precision of solving the equations of
stationary motion is for the stability analysis. For some EOSs, the hyperon
softening leads to spin-up by the angular momentum loss, in the form of the
back bending phenomenon, for a rather broad range of stellar baryon mass. We
show that large segments of the evolutionary tracks exhibiting the back bending
behaviour in the moment-of-inertia - rotation-frequency plane, are actually
unstable and therefore not astrophysically relevant. We show also that during
the spin-up - angular-momentum-loss epoch, an isolated neutron star (e.g., a
radio pulsar) can lose a sizable part of its initial angular momentum without
significantly changing its rotation period. We propose also simple arguments
and criteria allowing to connect the presence of a back bending epoch with the
mass-radius relations and the stiffness and/or softness of the nucleon and
hyperon EOSs of the neutron star core.

[72] oai:arXiv.org:astro-ph/0311128 [pdf] - 60708

Are rotating strange quark stars good sources of gravitational waves?

Comments: 15 pages, 10 figures, to appear in Astronomy and Astrophysics

Submitted: **2003-11-05**

We study the viscosity driven (Jacobi-like) bar mode instability of rapidly
rotating strange stars in general relativity. A triaxial, "bar shaped" compact
star could be an efficient source of continuous wave gravitational radiation in
the frequency range of the forthcoming interferometric detectors. We locate the
secular instability point along several constant baryon mass sequences of
uniformly rotating strange stars described by the MIT bag model. Contrary to
neutron stars, strange stars with T/|W| (the ratio of the rotational kinetic
energy to the absolute value of the gravitational potential energy) much lower
than the corresponding value for the mass-shed limit can be secularly unstable
to bar mode formation if shear viscosity is high enough to damp out any
deviation from uniform rotation. The instability develops for a broad range of
gravitational masses and rotational frequencies of strange quark stars. It
imposes strong constraints on the lower limit of the frequency at the innermost
stable circular orbit around rapidly rotating strange stars. The above results
are robust for all linear self-bound equations of state assuming the growth
time of the instability is faster than the damping timescale. We discuss
astrophysical scenarios where triaxial instabilities (r-mode and viscosity
driven instability) could be relevant in strange stars described by the
standard MIT bag model of normal quark matter. Taking into account actual
values of viscosities in strange quark matter and neglecting the magnetic field
we show that Jacobi-like instability cannot develop in any astrophysicaly
interesting temperature windows. The main result is that strange quark stars
described by the MIT bag model can be accelerated to very high frequency in Low
Mass X-ray binaries if the strange quark mass is ~ 200 MeV or higher.

[73] oai:arXiv.org:astro-ph/0209027 [pdf] - 880607

Thermal state of transiently accreting neutron stars

Comments: 7 pages, 2 figures, revised after the referee remarks, to appear in
A&A

Submitted: **2002-09-02**, last modified: **2003-07-10**

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.

[74] oai:arXiv.org:astro-ph/0209026 [pdf] - 51432

What we can learn from observations of cooling neutron stars

Comments: 5 pages, 2 figures, revised after the referee remarks

Submitted: **2002-09-02**, last modified: **2003-07-10**

A generic toy model of a cooling neutron star (NS) is used to analyze cooling
of NSs with nucleon and exotic composition of the cores. The model contains the
parameters which specify the levels of slow and enhanced neutrino emission as
well as the lower and upper densities of the layer where the slow emission
transforms into the enhanced one. The prospects to constrain these parameters
from the present and future observations of isolated middle-aged NSs are
discussed.

[75] oai:arXiv.org:astro-ph/0305220 [pdf] - 363031

Nuclear composition and heating in accreting neutron-star crusts

Comments: 5 pages, 2 figures, accepted for publication in A&A Letters

Submitted: **2003-05-13**

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.

[76] oai:arXiv.org:astro-ph/0301071 [pdf] - 362931

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

Comments: to be submitted to Astronomy & Astrophysics

Submitted: **2003-01-06**, last modified: **2003-03-19**

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.

[77] oai:arXiv.org:astro-ph/0301073 [pdf] - 1233039

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

Comments: 35 pages, 11 figures, uses eas.cls, to be published in ``Final Stages
of Stellar Evolution'', edts. J.-M. Hameury and C. Motch, EAS Publications
Series (EDP Sciences, 2003)

Submitted: **2003-01-06**

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.

[78] oai:arXiv.org:astro-ph/0209151 [pdf] - 51557

Moments of inertia for neutron and strange stars: limits derived for the
Crab pulsar

Comments: 5 pages, 2 figures included - accepted for publication in A&A

Submitted: **2002-09-09**, last modified: **2002-10-09**

Recent estimates of the properties of the Crab nebula are used to derive
constraints on the moment of inertia, mass and radius of the pulsar. To this
purpose, we employ an approximate formula combining these three parameters. Our
"empirical formula" I = a(x)MR^2, where x = (M/M_sun)(km/R), is based on
numerical results obtained for thirty theoretical equations of state of dense
matter. The functions a(x) for neutron stars and strange stars are
qualitatively different. For neutron stars a_NS(x) = x/(0.1+2x) for x <= 0.1
(valid for M > 0.2 M_sun) and a_NS(x) = (2/9)(1+5x) for x > 0.1. For strange
stars a_SS(x) = (2/5)(1+x) (not valid for strange stars with crust and M < 0.1
M_sun). We obtain also an approximate expression for the maximum moment of
inertia I_max,45 = (-0.37 + 7.12x_max) (M_max/M_sun)(R_(M_max)/(10 km))^2,
valid for both neutron stars and strange stars. Applying our formulae to the
evaluated values of I_Crab, we derive constraints on the mass and radius of the
pulsar. A very conservative evaluation of the expanding nebula mass, M_neb = 2
M_sun, yields M_Crab > 1.2 M_sun and R_Crab = 10-14 km. Setting the most recent
evaluation ("central value") M_neb = 4.6 M_sun rules out most of the existing
equations of state, leaving only the stiffest ones: M_Crab > 1.9 M_sun, R_Crab
= 14-15 km.

[79] oai:arXiv.org:astro-ph/0208078 [pdf] - 50894

Adiabatic Index of Dense Matter and Damping of Neutron Star Pulsations

Comments: 5 pages, 2 figures, accepted for publication in A&A

Submitted: **2002-08-03**

The adiabatic index Gamma_1 for perturbations of dense matter is studied
under various physical conditions which can prevail in neutron star cores. The
dependence of Gamma_1 on the composition of matter (in particular, on the
presence of hyperons), on the stellar pulsation amplitude, and on the baryon
superfluidity is analyzed. Timescales of damping of stellar pulsations are
estimated at different compositions, temperatures, and pulsation amplitudes.
Damping of pulsations by bulk viscosity in the neutron-star cores can prevent
the stars to pulsate with relative amplitudes > (1-15)% (depending on the
composition of matter).

[80] oai:arXiv.org:astro-ph/0204233 [pdf] - 48778

The Cooling Neutron Star in 3C 58

Comments: 3 pages, 2 figures, AA Letters (submitted)

Submitted: **2002-04-15**, last modified: **2002-05-13**

The upper limit of the effective surface temperature of the neutron star (NS)
PSR J0205+6449 in the supernova remnant 3C 58 obtained recently by Slane et al.
(2002) is analyzed using a modern theory of NS cooling (Kaminker et al. 2002).
The observational limit can be explained by cooling of a superfluid NS with the
core composed of neutrons, protons and electrons, where direct Urca process is
forbidden. However, combined with the data on the surface temperatures of other
isolated NSs, it gives evidence (emphasized by Slane et al.) that direct Urca
process is open in the inner cores of massive NSs. This evidence turns out to
be less stringent than the evidence provided by the well known observations of
Vela and Geminga.

[81] oai:arXiv.org:astro-ph/0201434 [pdf] - 47378

Equation of state of dense matter and the minimum mass of cold neutron
stars

Comments: 7 pages, 5 figures, uses aa.cls, accepted in Astronomy and
Astrophysics

Submitted: **2002-01-26**

Equilibrium configurations of cold neutron stars near the minimum mass are
studied, using the recent equation of state SLy, which describes in a unified,
physically consistent manner, both the solid crust and the liquid core of
neutron stars. Results are compared with those obtained using an older FPS
equation of state of cold catalyzed matter. The value of M_min\simeq 0.09M_sun
depends very weakly on the equation of state of cold catalyzed matter: it is
0.094 M_sun for the SLy model, and 0.088 M_sun for the FPS one. Central density
at M_min is significantly lower than the normal nuclear density: for the SLy
equation of state we get central density 1.7 10^{14} g/cm^3, to be compared
with 2.3 10^{14} g/cm^3 obtained for the FPS one. Even at M_min, neutron stars
have a small liquid core of radius of about 4 km, containing some 2-3% of the
stellar mass. Neutron stars with 0.09 M_sun <M<0.17 M_sun are bound with
respect to dispersed configuration of the hydrogen gas, but are unbound with
respect to dispersed Fe^56. The effect of uniform rotation on the minimum-mass
configuration of cold neutron stars is studied. Rotation increases the value of
M_min; at rotation period of 10 ms the minimum mass of neutron stars increases
to 0.13 M_sun, and corresponds to the mass-shedding (Keplerian) configuration.
In the case of the shortest observed rotation period of radio pulsars 1.56 ms,
minimum mass of uniformly rotating cold neutron stars corresponds to the
mass-shedding limit, and is found at 0.61 M_sun for the SLy EOS and 0.54 M_sun
for the FPS EOS.

[82] oai:arXiv.org:astro-ph/0110575 [pdf] - 45632

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

Comments: LaTeX, 9 pages, added reference, version accepted by Astron.
Astrophys

Submitted: **2001-10-26**, last modified: **2001-11-15**

Bulk viscosity of neutron star cores containing hyperons is studied taking
into account non-equilibrium weak process $n+n \rightleftharpoons p+\Sigma^-$.
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.

[83] oai:arXiv.org:astro-ph/0105485 [pdf] - 42714

Apparent radii of neutron stars and equation of state of dense matter

Comments: 4 pages, 2 figures, latex, uses aa.cls, to be published in Astronomy
and Astrophysics

Submitted: **2001-05-28**, last modified: **2001-11-12**

Apparent (radiation) radius of neutron star,R_\infty, depends on the star
gravitational mass in quite a different way than the standard coordinate radius
in the Schwarzschild metric, R. We show that, for a broad set of equations of
state of dense matter,
R_\infty(M_max) for the configurations with maximum allowable masses is very
close to the absolute lower bound on R_\infty at fixed M, resulting from the
very definition of R_\infty.
Also, the value of R_infty at given M, corresponding to the maximum
compactness (minimum R) of neutron star consistent with general relativity and
condition v_sound<c, is only 0.6% higher than this absolute lower bound.
Theoretical predictions for R_\infty are compared with existing observational
estimates of the apparent radii of neutron stars.

[84] oai:arXiv.org:astro-ph/0111092 [pdf] - 45858

A unified equation of state of dense matter and neutron star structure

Comments: 17 pages, 10 figures, latex, uses aa.cls, to be published in
Astronomy and Astrophysics

Submitted: **2001-11-05**, last modified: **2001-11-12**

An equation of state (EOS) of neutron star matter, describing both the
neutron star crust and the liquid core, is calculated. It is based on the
effective nuclear interaction SLy of the Skyrme type, which is particularly
suitable for the application to the calculation of the properties of very
neutron rich matter (Chabanat et al. 1997, 1998). The structure of the crust,
and its EOS, is calculated in the T=0 approximation, and under the assumption
of the ground state composition. The crust-core transition is a very weakly
first-order phase transition, with relative density jump of about one percent.
The EOS of the liquid core is calculated assuming (minimal) n-p-e-mu
composition. Parameters of static neutron stars are calculated and compared
with existing observational data on neutron stars. The minimum and maximum
masses of static neutron stars are 0.094 M_sun and 2.05 M_sun, respectively.
Effects of rotation on the minimum and the maximum mass of neutron stars are
briefly discussed.

[85] oai:arXiv.org:astro-ph/0111162 [pdf] - 397389

Recycling strange stars to millisecond periods

Comments: 8 pages, 10 figures, accepted for publication in A&A

Submitted: **2001-11-08**

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.

[86] oai:arXiv.org:astro-ph/0105047 [pdf] - 42276

Nucleon Superfluidity vs Observations of Cooling Neutron Stars

Comments: 4 pages, 2 figures, AA Letters, accepted

Submitted: **2001-05-03**, last modified: **2001-05-23**

Cooling simulations of neutron stars (NSs) are performed assuming that
stellar cores consist of neutrons, protons and electrons and using realistic
density profiles of superfluid critical temperatures $T_{cn}(\rho)$ and
$T_{cp}(\rho)$ of neutrons and protons. Taking a suitable profile of
$T_{cp}(\rho)$ with maximum $\sim 5 \times 10^9$ K one can obtain smooth
transition from slow to rapid cooling with increasing stellar mass. Adopting
the same profile one can explain the majority of observations of thermal
emission from isolated middle--aged NSs by cooling of NSs with different masses
either with no neutron superfluidity in the cores or with a weak superfluidity,
$T_{cn} < 10^8$ K. The required masses range from $\sim 1.2 M_\odot$ for (young
and hot) RX J0822-43 and (old and warm) PSR 1055-52 and RX J1856-3754 to
$\approx 1.45 M_\odot$ for the (colder) Geminga and Vela pulsars. Observations
constrain the $T_{cn}(\rho)$ and $T_{cp}(\rho)$ profiles with respect to the
threshold density of direct Urca process and maximum central density of NSs.

[87] oai:arXiv.org:astro-ph/0105105 [pdf] - 42334

Thermal Conductivity of Neutrons in Neutron Star Cores

Comments: 13 pages, 4 figures, accepted to AA

Submitted: **2001-05-07**

The diffusive thermal conductivity of neutrons in dense matter [$\rho \sim (1
- 8) \times 10^{14}$ g cm$^{-3}$] of neutron star cores is calculated. The
contribution from neutron--neutron and neutron--proton collisions is taken into
account. We use the transition probabilities calculated for symmetric dense
nucleon matter on the basis of the Dirac--Brueckner approach to the in-medium
effects and the Bonn model of bare nucleon--nucleon interaction. The diffusive
thermal conductivity of neutrons in the presence of neutron and proton
superfluidities is analyzed in a microscopic manner; the effects of
superfluidity are shown to be significant. The low temperature behavior of the
thermal conductivity appears to be extremely sensitive to the relation between
critical temperatures of neutrons and protons. The results are fitted by simple
analytic expressions. In combination with the formulae for the electron and
muon thermal conductivities, obtained earlier, the present expressions provide
a realistic description of the full diffusive thermal conductivity in the
neutron star cores for normal and various superfluid phases.

[88] oai:arXiv.org:astro-ph/0104116 [pdf] - 880498

The crust of rotating strange quark stars

Comments: 9 pages, 10 figures, accepted for publication in A&A

Submitted: **2001-04-06**

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.

[89] oai:arXiv.org:astro-ph/0103290 [pdf] - 41507

Bulk viscosity in superfluid neutron star cores. II. Modified Urca
processes in $npe\mu$ matter

Comments: 8 pages, 4 figures, accepted by A&A

Submitted: **2001-03-19**

We study the bulk viscosity in neutron star cores due to modified Urca
processes involving nucleons, electrons and muons and analize its reduction by
singlet-state or triplet-state superfluidity of nucleons. In combination with
the results of our previous paper on the bulk viscosity due to direct Urca
processes, a realistic description of the bulk viscosity in superfluid neutron
star cores is obtained. Switching off direct Urca processes with decreasing
density in a nonsuperfluid matter lowers the bulk viscosity by 3-5 orders of
magnitude. The presense of muons opens additional source of bulk viscosity due
to muon Urca processes and lowers the threshold density of the electron direct
Urca process. The superfluidity may strongly reduce the bulk viscosity and
affect thus damping of neutron star vibrations.

[90] oai:arXiv.org:astro-ph/0012543 [pdf] - 40132

Protoneutron stars and neutron stars

Comments: ASP Conference Series "Pulsar Astronomy - 2000 and Beyond", eds
Kramer, Wex, Wielebinski, Vol 202, 663

Submitted: **2000-12-29**

We find constraints on minimum and maximum mass of ordinary neutron stars
imposed by their early evolution (protoneutron star stage). We calculate models
of protoneutron stars using a realistic standard equation of state of hot,
dense matter valid for both supranuclear and subnuclear densities. Results for
different values of the nuclear incompressibility are presented.

[91] oai:arXiv.org:astro-ph/0012541 [pdf] - 40130

Protoneutron stars and constraints on maximum and minimum mass of
neutron stars

Comments: 4 pages, 2 figures, ASP Conference Series, Vol 138, 131, eds K. L.
Chan, K. S. Cheng, H. P. Singh

Submitted: **2000-12-29**

Constraints on minimum and maximum mass of ordinary neutron stars are imposed
by the consideration of their early evolution (protoneutron star stage).
Calculations are performed for a realistic standard model of hot, dense matter
(Lattimer & Swesty 1991) valid for both supranuclear and subnuclear densities.
Various assumptions concerning protoneutron star interior (large trapped lepton
number, no trapped lepton number, isentropic, isothermal) are taken into
account.

[92] oai:arXiv.org:astro-ph/0012122 [pdf] - 387036

Neutrino Emission from Neutron Stars

Comments: review, 165 pages, Physics Reports, 2001 in press

Submitted: **2000-12-05**

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.

[93] oai:arXiv.org:astro-ph/0009282 [pdf] - 38124

Rapidly rotating strange stars

Comments: ASP Conference Series "Pulsar Astronomy - 2000 and Beyond", eds
Kramer, Wex, Wielebinski, Vol 202, 661

Submitted: **2000-09-18**

We study effects of the strange quark mass and of the QCD coupling constant
in the MIT bag model of quark matter with lowest order quark-gluon interactions
on the rapid rotation of strange stars. The influence of rotation on global
parameters of strange stars is greater than in the case of the neutron star
models. We show that independently of the strange quark mass and the coupling
constans the ratio of the rotational kinetic energy to the absolute value of
the gravitational potential energy for rotating strange stars is significantly
higher than for ordinary neutron stars. This might indicate that rapidly
rotating strange stars could be important sources of gravitational waves.

[94] oai:arXiv.org:astro-ph/0006135 [pdf] - 36515

Inner edge of neutron-star crust with SLy effective nucleon-nucleon
interactions

Comments: 14 pages, 5 figures, latex, to be published in Physics Letters B

Submitted: **2000-06-09**

The boundary between the solid crust, assumed to be in complete thermodynamic
equilibrium (cold catalyzed matter), and the liquid core of a neutron star is
studied using Skyrme SLy effective N-N interactions. An approximate value of
the density at the inner edge of the crust is obtained from the threshold for
the instability of homogeneous $npe$ matter with respect to small periodic
density perturbations. Calculations of the bottom layer of the crust are
performed within the Compressible Liquid Drop Model. Spherical nuclei are
energetically preferred over exotic ones (cylinders, plates, tubes, bubbles),
down to the inner edge of the crust, found at the density 0.08 fm^{-3}.

[95] oai:arXiv.org:astro-ph/0004278 [pdf] - 35637

On the mass of moderately rotating strange stars in the MIT bag model
and LMXBs

Comments: 5 pages, 3 figures, A&A in press

Submitted: **2000-04-19**

We compute the maximum mass of moderately rotating strange stars as a
function of the strange quark mass, of the QCD coupling constant, $\alpha_c$,
and of the bag constant (vacuum energy density), $B$, in the MIT bag model of
quark matter with lowest order quark-gluon interactions. For a fixed value of
$B$, the maximum stellar mass depends only weakly on $\alpha_c$, and is
independent of this coupling in the limit of massless quarks. However, if it is
the value of the chemical potential of quark matter at zero pressure which is
held constant, for example at the value corresponding to the stability limit of
nucleons against conversion to quark matter, the maximum mass of the strange
star is higher by up to 25% for $\alpha_c=0.6$, than for non-interacting
quarks, and this may be relevant in the discussion of kHz QPO sources. The
maximum mass of a non-rotating strange star could be sufficiently high to allow
an orbital frequency as low as 1.0 kHz in the marginally stable orbit. However,
for all $\alpha_c<0.6$, the stellar mass cannot exceed $2.6M_\odot$ at any
rotational period $\ge1.6 $ms.

[96] oai:arXiv.org:astro-ph/0004185 [pdf] - 35544

Cooling neutron stars with localized protons

Comments: 5 pages, 3 figures, latex, uses aa.cls, to be published in Astronomy
and Astrophysics

Submitted: **2000-04-13**

We analyze cooling of neutron stars, assuming the presence of localized
protons in the densest region of their cores. Choosing a single threshold
density for proton localization and adjusting neutron star mass, we reproduce
the observational data on effective surface temperatures of Vela and PSR
0656+14, with or without an accreted hydrogen envelope. However, the presence
of a tiny hydrogen envelope is mandatory, in this model, for reproducing the
Geminga data.

[97] oai:arXiv.org:astro-ph/0004183 [pdf] - 35542

Bulk viscosity in superfluid neutron star cores. I. Direct Urca
processes in npe\mu matter

Comments: 14 pages, 7 figures, latex, uses aa.cls, to be published in Astronomy
and Astrophysics

Submitted: **2000-04-13**

The bulk viscosity of the neutron star matter due to the direct Urca
processes involving nucleons, electrons and muons is studied taking into
account possible superfluidity of nucleons in the neutron star cores. The cases
of singlet-state pairing or triplet-state pairing (without and with nodes of
the superfluid gap at the Fermi surface) of nucleons are considered. It is
shown that the superfluidity may strongly reduce the bulk viscosity. The
practical expressions for the superfluid reduction factors are obtained. For
illustration, the bulk viscosity is calculated for two models of dense matter
composed of neutrons, protons,electrons and muons. The presence of muons
affects the bulk viscosity due to the direct Urca reactions involving electrons
and produces additional comparable contribution due to the direct Urca
reactions involving muons. The results can be useful for studying damping of
vibrations of neutron stars with superfluid cores.

[98] oai:arXiv.org:astro-ph/0002394 [pdf] - 265299

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

Comments: 8 pages, 4 figures, accepted for publication in A&A

Submitted: **2000-02-21**

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.

[99] oai:arXiv.org:astro-ph/9908249 [pdf] - 332247

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

Comments: 24 pages, 4 figures, latex, uses appb.sty

Submitted: **1999-08-23**

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.

[100] oai:arXiv.org:astro-ph/9907225 [pdf] - 265353

Fast rotation of strange stars

Comments: 13 pages, 7 figures, latex, uses aa.cls, epsf.tex, to be published in
Astronomy and Astrophysics

Submitted: **1999-07-16**

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.

[101] oai:arXiv.org:astro-ph/9906312 [pdf] - 332246

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

Comments: 27 pages, no figures, uses appb.sty (included)

Submitted: **1999-06-18**

As pointed out by Kutschera and W{\'o}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.

[102] oai:arXiv.org:astro-ph/9903127 [pdf] - 363263

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

Comments: 10 pages, 7 figures, LaTeX using aa.cls and epsf.sty. Corrected typos
in Eq.(16) (v2) and in Eq.(A3) (v3)

Submitted: **1999-03-09**, last modified: **1999-05-25**

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.

[103] oai:arXiv.org:astro-ph/9905036 [pdf] - 106362

Maximum redshift and minimum rotation period of neutron stars

Comments: 5 pages, latex, uses texas.sty. Presented at the 19th Texas
Symposium, 14-19 December 1998, Paris, France. To appear in the 19th Texas
Symposium Proceedings (CD-ROM)

Submitted: **1999-05-05**

The lower bound on the period of uniform rotation of neutron stars (NS) with
causal equation of state (EOS), P_ min, is shown to be determined, with a high
precision, by the upper bound on the surface redshift of static neutron star
models, z_ max. We obtain z_max=0.851, which turns out to be about 5% lower
than upper bound obtained by Lindblom (1984). Using then the empirical formula
for the maximum rotation frequency of uniform rotation of neutron star models
with realistic causal EOS, derived by Lasota et al. (1996), we establish an
approximate relation between z_max and P_min, in which P_min is proportional to
the maximum mass of static NS, M_s^ max. In this way we reproduce, within 2%,
the exact formula of Koranda et al. (1997) based on extensive exact numerical
calculations. Replacing M_s^max by the maximum measured NS mass, 1.442 M_sun,
we get a lower bound P_min = 0.288 ms.

[104] oai:arXiv.org:astro-ph/9904166 [pdf] - 363275

Neutrino emission due to proton pairing in neutron stars

Comments: LaTeX, uses aa.cls, 3 pages, 2 figures, published in Astron.
Astrophys

Submitted: **1999-04-13**

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.

[105] oai:arXiv.org:astro-ph/9901118 [pdf] - 104712

On the minimum period of uniformly rotating meutron stars

Comments: 3 pages, latex, uses l-aa.sty, accepted for publication in Astronomy
and Astrophysics

Submitted: **1999-01-11**

We show that the neutron star minimum period of rotation is determined, for
causal equations of state, by the maximum value of the relativistic
(compactness) parameter 2GM/Rc^2, allowed by causality for static neutron
stars, 0.7081, and by the largest measured mass of a neutron star. The relation
between these three quantities, resulting from the extrapolation of the
empirical formulae of Lasota et al. (1996), yields the minimum period of 0.288
ms, only 2% higher than an absolute lower bound obtained in extensive exact
numerical calculations of Koranda et al. (1997).

[106] oai:arXiv.org:astro-ph/9808030 [pdf] - 102374

Bounds on the existence of neutron rich nuclei in neutron star interiors

Comments: Invited talk presented by P. Haensel at the international conference
``Nuclear Physics Close to the Barrier'' Warsaw, Poland, 30.06-4.07.1998. 12
pages, latex, including 2 figures, to appear in Acta Physica Polonica

Submitted: **1998-08-04**

We address the question concerning the maximum density, at which nuclei (and
more generally -- nuclear structures) can exist in neutron star interiors. An
absolute upper bound to such a density is obtained using the bulk
approximation, in which surface and Coulomb effects are neglected. A very good
approximation to the limiting density is given by the threshold for the
instability of a uniform n-p-e plasma with respect to density modulations; this
threshold is calculated using the Extended Thomas-Fermi approximation for the
Skyrme energy functionals. For recent SLy Skyrme forces, which are particularly
suitable for the description of very neutron rich nucleon systems, the maximum
density at which nuclear structures can exist is 0.08 nucleons/fm^3 (i.e., half
of the normal nuclear density); at this density protons constitute only 4\% of
nucleons.

[107] oai:arXiv.org:astro-ph/9711347 [pdf] - 99480

Rapid differential rotation of protoneutron stars and constraints on
radio pulsars periods

Comments: 13 pages, 8 figures, needs l-aa.sty, epsf.sty, accepted for
publication in Astronomy and Astrophysics

Submitted: **1997-11-28**

Models of differentially rotating protoneutron stars are calculated, using
realistic equations of state of dense hot matter. Various conditions within the
stellar interior, corresponding to different stages of protoneutron star
evolution, are considered.
Families of differentially rotating models of a given baryon mass are
calculated, using a two-parameter formula describing the angular velocity
profile within a rotating protoneutron star. Maximum angular momentum, which
can be accommodated by a protoneutron star is calculated, for various thermal
conditions in stellar interior, for a baryon mass of $1.5 M_\odot$. In the case
of a thermally homogeneous (isentropic or isothermal) neutrino-opaque interior
this maximum angular momentum turns out to be somewhat higher than that of a
cold neutron star of the same baryon mass, rotating uniformly at the mass
shedding angular velocity. However, if the protoneutron star has a thermal
structure characteristic of initial state, with a low entropy (unshocked) core,
and a high entropy (shocked) outer half of baryon mass, the maximum angular
momentum is significantly lower. This leads to a minimum period of uniform
rotation of cold neutron stars of baryon mass $\sim 1.5~M_\odot$, formed
directly (i.e. without a subsequent significant accretion of mass) from
protoneutron stars with shocked envelope, of about 1.7 ms and strengthens the
hypothesis that millisecond pulsars are accretion accelerated neutron stars.

[108] oai:arXiv.org:astro-ph/9708181 [pdf] - 98359

Neutrino synchrotron emission from dense magnetized electron gas of
neutron stars

Comments: 11 pages, latex, including 8 figures, accepted by Astronomy and
Astrophysics, Main Journal, uses l-aa.sty

Submitted: **1997-08-20**

We study the synchrotron emission of neutrino pairs by relativistic,
degenerate electrons in strong magnetic fields. Particular attention is paid to
the case in which the dominant contribution comes from one or several lowest
cyclotron harmonics. Calculations are performed using the exact quantum
formalism and the quasiclassical approach. Simple analytic fits to the neutrino
synchrotron emissivity are obtained in the domain of magnetized, degenerate and
relativistic electron gas provided the electrons populate either many Landau
levels or the ground level alone. The significance of the neutrino synchrotron
energy losses in the interiors of cooling neutron stars is discussed.

[109] oai:arXiv.org:astro-ph/9705157 [pdf] - 97419

Radial pulsations and stability of protoneutron stars

Comments: 11 pages, 11 figures, LaTeX, uses l-aa.sty and epsfig.sty, to appear
in Astronomy and Astrophysics CAMK preprint no: 315

Submitted: **1997-05-20**

Radial pulsations of newborn neutron stars (protoneutron stars) are studied
for a range of internal temperatures and entropies per baryon predicted by the
existing numerical simulations. Protoneutron star models are constructed using
a realistic equation of state of hot dense matter, and under various
assumptions concerning stellar interior (large trapped lepton number, zero
trapped lepton number, isentropic, isothermal). Under prevailing conditions,
linear oscillations of a protoneutron star can be treated as adiabatic, and
evolutionary effects can be neglected on dynamic timescale. The dynamic
behavior is governed by the adiabatic index, which in turn depends on the
physical state of the stellar interior. The eigenfrequencies of the lowest
radial modes of linear, adiabatic pulsations are calculated. Stability of
protoneutron stars with respect to small radial perturbations is studied, and
the validity of the static stability criteria is discussed.

[110] oai:arXiv.org:astro-ph/9702155 [pdf] - 96674

Neutrino pair emission due to scattering of electrons off fluxoids in
superfluid neutron star cores

Comments: 11 pages, 2 figures, needs l-aa.sty, epsf.sty, accepted for
publication in Astronomy and Astrophysics

Submitted: **1997-02-18**

We study the emission of neutrinos, resulting from the scattering of
electrons off magnetic flux tubes (fluxoids) in the neutron star cores with
superfluid (superconducting) protons. In the absence of proton superfluidity
(T> T_{cp}), this process transforms into the well known electron synchrotron
emission of neutrino pairs in a locally uniform magnetic field B, with the
neutrino energy loss rate Q proportional to B^2 T^5. For temperatures T not
much below T_{cp}, the synchrotron regime (Q \propto T^5) persists and the
emissivity Q can be amplified by several orders of magnitude due to the
appearance of the fluxoids and associated enhancement of the field within them.
For lower T, the synchrotron regime transforms into the bremsstrahlung regime
(Q \propto T^6) similar to the ordinary neutrino-pair bremsstrahlung of
electrons which scatter off atomic nuclei. We calculate Q numerically and
represent our results through a suitable analytic fit. In addition, we estimate
the emissivities of two other neutrino-production mechanisms which are usually
neglected -- neutrino-pair bremsstrahlung processes due to electron-proton and
electron-electron collisions. We show that the electron-fluxoid and
electron-electron scattering can provide the main neutrino production
mechanisms in the neutron star cores with highly superfluid protons and
neutrons at T < 5 10^8 K. The electron-fluxoid scattering is significant if the
initial, locally uniform magnetic field B > 10^{13} G.

[111] oai:arXiv.org:astro-ph/9610265 [pdf] - 95789

Rapid uniform rotation of protoneutron stars

Comments: 14 pages, LaTeX, 12 eps figures included, uses l-aa.sty and
epsfig.sty, accepted for publication in A&A (10/1996)

Submitted: **1996-10-31**

Rapid uniform rotation of newborn neutron stars (protoneutron stars) is
studied for a range of internal temperatures and entropies per baryon predicted
by the existing numerical simulations. Calculations are performed using general
relativistic equations of hydrostatic equilibrium of rotating, axially
symmetric stars. Stability of rotating configurations with respect to mass
shedding and the axially symmetric perturbations is studied. Numerical
calculations are performed for a realistic dense matter equation of state,
under various assumptions concerning neutron star interior (large trapped
lepton number, no trapped lepton number, isentropic, isothermal). For
configurations with baryon mass well below the maximum one for the non-rotating
models, the mass shedding limit depends quite sensitively on the position of
the ``neutrinosphere''. The absolute upper limit on rotation frequency is, to a
good approximation, obtained for the maximum baryon mass of rotating
configurations. Empirical formula for the maximum rotation frequency of
uniformly rotating protoneutron stars is shown to be quite precise; it actually
coincides with that used for cold neutron stars. Evolutionary sequences at
fixed baryon mass and angular momentum, which correspond to evolution of
protoneutron stars into cold neutron stars are studied, and resulting
constraints on the maximum rotation frequency of solitary pulsars are
discussed.

[112] oai:arXiv.org:astro-ph/9608090 [pdf] - 95208

Neutrino-electron scattering in dense magnetized plasma

Comments: Physical Review D, to be published, 18 pages, using REVTEX, without
figures. Figures (hardcopy) available upon request from one of the authors
(haensel@camk.edu.pl)

Submitted: **1996-08-14**

We derive exact expressions for the cross section of neutrino scattering on
electrons in dense, hot stellar matter, in the presence of strong magnetic
fields. Numerical calculations of the scattering cross sections at various
densities, temperatures and magnetic fields, are performed. Strong, quantizing
magnetic fields modify significantly the angular and energy dependence of the
scattering cross section.

[113] oai:arXiv.org:astro-ph/9605164 [pdf] - 94720

Solid interiors of neutron stars and gravitational radiation

Comments: 26 pages, Latex, uses l-school.sty, epsf.tex. Four figures. Lectures
delivered at Les Houches Summer School 1995 ``Astrophysical Sources of
Gravitational Radiation''. To be published in the Proceedings (eds.
J.-A.Marck and J.-P.Lasota)

Submitted: **1996-05-27**

Contents:
1.Introduction
2. Formation and composition of the solid crust
3. The bottom layers of the crust: exotic nuclear shapes and topologies
4. Crust contribution to stellar mass and moment of inertia
5. Solid cores
6. Deformation, elastic strain and mountains
7. Shape and energy of rotating, partially solid neutron star
8. Oblateness and precession
9. Gravitational radiation from precessing pulsars
10. Gravitational radiation from precessing pulsars and pulsar timing
11. Detectability

[114] oai:arXiv.org:astro-ph/9605149 [pdf] - 363153

Ferromagnetism of dense matter and magnetic properties of neutron stars

Comments: 8 pages, Latex, uses l-aa.sty, epsf.tex. Two figures, fig1.eps,
fig2.eps. To be published in Astronomy and Astrophysics

Submitted: **1996-05-24**

Possible consequences of ferromagnetic transition in dense matter suggested
recently by Kutschera and W{\'o}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.

[115] oai:arXiv.org:astro-ph/9604073 [pdf] - 363146

Electron $\nu \bar{\nu}$ Bremsstrahlung in a Liquid Phase of Neutron
Star Crusts

Comments: 14 pages, 8 eps-files of figures, uses l-aa.sty, accepted for
publication in Astronomy and Astrophysics

Submitted: **1996-04-13**

Neutrino emissivity from the electron $\nu\bar\nu$ 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.

[116] oai:arXiv.org:astro-ph/9508118 [pdf] - 93174

Fast Rotation of Neutron Stars

Comments: 12 pages, LateX, uses aaspp.sty file, figures (Postscript files)
available from haensel@camk.edu.pl; to be published in the Astrophysical
Journal

Submitted: **1995-08-25**

We show that for realistic equations of state of dense matter, the universal
proportionality factor relating the maximum rotation rate of neutron stars due
to mass-shedding limit to the mass and radius of maximum allowable mass
configuration of non-rotating models results from a universal proportionality
between masses and radii of static maximum-mass neutron stars and those of
maximally rotating configurations. These empirical relations cannot be obtained
in the slow rotation approximation.

[117] oai:arXiv.org:nucl-th/9310003 [pdf] - 118846

Experimental Nuclear Masses and the Ground State of Cold Dense Matter

Comments: 15 pages in standard TeX in this file, 5 tables in standard LaTeX in
a separate file and 4 figures in an attached uuencoded tar.Z file. Preprint
number not available, Accepted for publication in Astron. Astrophys

Submitted: **1993-10-04**

We study the consequences of recent progress in the experimental
determination of masses of neutron rich nuclei for our knowledge of the ground
state of cold dense matter. The most recent experimental data determine the
ground state of cold dense matter up to $ \rho \simeq 10^{11}~{\rm g~cm^{-3}}
$. The composition and the equation of state of the ground state of matter, in
this density interval, are calculated.