Physics input for modelling superfluid neutron stars with hyperon cores

Observations of massive ($M approx 2.0~M_odot$) neutron stars (NSs), PSRs J1614-2230 and J0348+0432, rule out most of the models of nucleon-hyperon matter employed in NS simulations. Here we construct three possible models of nucleon-hyperon matter consistent with the existence of $2~M_odot$ pulsars as well as with semi-empirical nuclear matter parameters at saturation, and semi-empirical hypernuclear data. Our aim is to calculate for these models all the parameters necessary for modelling dynamics of hyperon stars (such as equation of state, adiabatic indices, thermodynamic derivatives, relativistic entrainment matrix, etc.), making them available for a potential user. To this aim a general non-linear hadronic Lagrangian involving $sigmaomegarhophisigma^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.

Disk accretion onto a magnetized star

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

Formation of millisecond pulsars - NS initial mass and EOS constraints

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

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

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.

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

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.

Keplerian frequency of uniformly rotating neutron stars and quark stars

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

Thermal conductivity of ions in a neutron star envelope

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.