Observations of thermal radiation from neutron stars can potentially provide information about the states of supranuclear matter in the interiors of these stars with the aid of the theory of neutron-star thermal evolution. We review the basics of thi
s theory for isolated neutron stars with strong magnetic fields, including most relevant thermodynamic and kinetic properties in the stellar core, crust, and blanketing envelopes.
We study the plasma correlation effects on nonresonant thermonuclear reactions of carbon and oxygen in the interiors of white dwarfs and liquid envelopes of neutron stars. We examine the effects of electron screening on thermodynamic enhancement of t
hermonuclear reactions in dense plasmas beyond the linear mixing rule. Using these improved enhancement factors, we calculate carbon and oxygen ignition curves in white dwarfs and neutron stars. The energy balance and ignition conditions in neutron star envelopes are evaluated, taking their detailed thermal structure into account. The result is compared to the simplified one-zone model, which is routinely used in the literature. We also consider the effect of strong magnetic fields on the ignition curves in the ocean of magnetars.
The equation of state (EOS) for partially ionized carbon, oxygen, and carbon-oxygen mixtures at temperatures 3times10^5 K <~ T <~ 3times10^6 K is calculated over a wide range of densities, using the method of free energy minimization in the framework
of the chemical picture of plasmas. The free energy model is an improved extension of our model previously developed for pure carbon (Phys. Rev. E, 72, 046402; arXiv:physics/0510006). The internal partition functions of bound species are calculated by a self-consistent treatment of each ionization stage in the plasma environment taking into account pressure ionization. The long-range Coulomb interactions between ions and screening of the ions by free electrons are included using our previously published analytical model, recently improved, in particular for the case of mixtures. We also propose a simple but accurate method of calculation of the EOS of partially ionized binary mixtures based on detailed ionization balance calculations for pure substances.
Contrary to what is claimed by Ferrer et al. [Phys. Rev. C 82, 065802 (2010)], the magnetic field of a neutron star cannot exceed 10^{19} G and the thermodynamic pressure of dense magnetized fermion gas is isotropic.
Some of X-ray dim isolated neutron stars (XDINS) and central compact objects in supernova remnants (CCO) show absorption features in their thermal soft X-ray spectra. It has been hypothesized that these features could be due to the periodic peaks in
free-free absorption opacities, caused by either Landau quantization of electron motion in magnetic fields B<10^{11} G or analogous quantization of ion motion in magnetic fields B>10^{13} G. Here, I review the physics behind cyclotron quantum harmonics in free-free photoabsorption, discuss different approximations for their calculation, and explain why the ion cyclotron harmonics (beyond the fundamental) cannot be observed.
We briefly review analytic approximations of thermodynamic functions of fully ionized nonideal electron-ion plasmas, applicable in a wide range of plasma parameters, including the domains of nondegenerate and degenerate, nonrelativistic and relativis
tic electrons, weakly and strongly coupled Coulomb liquids, classical and quantum Coulomb crystals. We present improvements to previously published approximations. Our code for calculation of thermodynamic functions based on the reviewed approximations is made publicly available.
Recently developed analytic approximation for the equation of state of fully ionized nonideal electron-ion plasma mixtures [Potekhin et al., Phys. Rev. E, 79, 016411 (2009); arXiv:0812.4344], which covers the transition between the weak and strong Co
ulomb coupling regimes and reproduces numerical results obtained in the hypernetted chain (HNC) approximation, is modified in order to fit the small deviations from the linear mixing in the strong coupling regime, revealed by recent Monte Carlo simulations. In addition, a mixing rule is proposed for the regime of weak coupling, which generalizes post-Debye density corrections to the case of mixtures and numerically agrees with the HNC approximation in that regime.
We study the thermal structure and evolution of magnetars as cooling neutron stars with a phenomenological heat source in an internal layer. We focus on the effect of magnetized (B > 10^{14} G) non-accreted and accreted outermost envelopes composed o
f different elements, from iron to hydrogen or helium. We discuss a combined effect of thermal conduction and neutrino emission in the outer neutron star crust and calculate the cooling of magnetars with a dipole magnetic field for various locations of the heat layer, heat rates and magnetic field strengths. Combined effects of strong magnetic fields and light-element composition simplify the interpretation of magnetars in our model: these effects allow one to interpret observations assuming less extreme (therefore, more realistic) heating. Massive magnetars, with fast neutrino cooling in their cores, can have higher thermal surface luminosity.
We develop analytic approximations of thermodynamic functions of fully ionized nonideal electron-ion plasma mixtures. In the regime of strong Coulomb coupling, we use our previously developed analytic approximations for the free energy of one-compone
nt plasmas with rigid and polarizable electron background and apply the linear mixing rule (LMR). Other thermodynamic functions are obtained through analytic derivation of this free energy. In order to obtain an analytic approximation for the intermediate coupling and transition to the Debye-Hueckel limit, we perform hypernetted-chain calculations of the free energy, internal energy, and pressure for mixtures of different ion species and introduce a correction to the LMR, which allows a smooth transition from strong to weak Coulomb coupling in agreement with the numerical results.
