Ages and thermal luminosities of neutron stars, inferred from observations, can be interpreted with the aid of the neutron star cooling theory to gain information on the properties of superdense matter in neutron-star interiors. We present a survey of estimated ages, surface temperatures and thermal luminosities of middle-aged neutron stars with relatively weak or moderately strong magnetic fields, which can be useful for these purposes. The catalogue includes results selected from the literature, supplemented with new results of spectral analysis of a few cooling neutron stars. The data are compared with the theory. We show that overall agreement of theoretical cooling curves with observations improves substantially for models where neutron superfluidity in stellar core is weak.

We study long-term thermal evolution of neutron stars in soft X-ray transients (SXTs), taking the deep crustal heating into account consistently with the changes of the composition of the crust. We collect observational estimates of average accretion rates and thermal luminosities of such neutron stars and compare the theory with observations. We perform simulations of thermal evolution of accreting neutron stars, considering the gradual replacement of the original nonaccreted crust by the reprocessed accreted matter, the neutrino and photon energy losses, and the deep crustal heating due to nuclear reactions in the accreted crust. We test and compare results for different modern theoretical models. We update a compilation of the observational estimates of the thermal luminosities in quiescence and average accretion rates in the SXTs and compare the observational estimates with the theoretical results. Long-term thermal evolution of transiently accreting neutron stars is nonmonotonic. The quasi-equilibrium temperature in quiescence reaches a minimum and then increases toward the final steady state. The quasi-equilibrium thermal luminosity of a neutron star in an SXT can be substantially lower at the minimum than in the final state. This enlarges the range of possibilities for theoretical interpretation of observations of such neutron stars. The updates of the theory and observations leave unchanged the previous conclusions that the direct Urca process operates in relatively cold neutron stars and that an accreted heat-blanketing envelope is likely present in relatively hot neutron stars in the SXTs in quiescence. The results of the comparison of theory with observations favor suppression of the triplet pairing type of nucleon superfluidity in the neutron-star matter.

Magnetars are believed to host the strongest magnetic fields in the present universe ($B\gtrsim10^{14}$ G) and the study of their persistent emission in the X-ray band offers an unprecendented opportunity to gain insight into physical processes in the presence of ultra-strong magnetic fields. Up to now, most of our knowledge about magnetar sources came from spectral analysis, which allowed to test the resonant Compton scattering scenario and to probe the structure of the star magnetosphere. On the other hand, radiation emitted from magnetar surface is expected to be strongly polarized and its observed polarization pattern bears the imprint of both scatterings onto magnetospheric charges and QED effects as it propagates in the magnetized vacuum around the star. X-ray polarimeters scheduled to fly in the next years will finally allow to exploit the wealth of information stored in the polarization observables. Here we revisit the problem of assessing the spectro-polarimetric properties of magnetar persistent emission. At variance with previous investigations, proper account for more physical surface emission models is made by considering either a condensed surface or a magnetized atmosphere. Results are used to simulate polarimetric observations with the forthcoming Imaging X-ray Polarimetry Explorer (IXPE). We find that X-ray polarimetry will allow to detect QED vacuum effects for all the emission models we considered and to discriminate among them.

We generalize our earlier work on neutron stars (arXiv:1903.04981), which assumed spherical Wigner-Seitz cells in the inner crust, to admit the possibility of pasta phases, i.e., non-spherical cell shapes. Full fourth-order extended Thomas-Fermi calculations using the density functional BSk24 are performed for cylindrical and plate-like cells. Unlike in our spherical-cell calculations, we do not include shell and pairing corrections, but there are grounds for expecting these corrections for pasta to be significantly smaller. It is therefore meaningful to compare the ETF pasta results with the full spherical-cell results, i.e., with shell and pairing corrections included. However, in view of the many previous studies in which shell and pairing corrections were omitted entirely, it is of interest to compare our pasta results with the ETF part of the corresponding spherical calculations. Making this latter comparison we find that as the density increases the cell shapes pass through the usual sequence sphere $\to$ cylinder $\to$ plate before the transition to the homogeneous core. The filling fractions found at the phase transitions are in close agreement with expectations based on the liquid-drop model. On the other hand, when we compare with the full spherical-cell results, we find the sequence to be sphere $\to$ cylinder $\to$ sphere $\to$ cylinder $\to$ plate. In neither case do any "inverted", i.e., bubble-like, configurations appear. We provide accurate fitting formulas to all our essential numerical results for each of the three phases, designed especially for the density range where the nonspherical shapes are expected, which enable one to capture not only the general behavior of the fitted functions, but also the differences between them in different phases.

We report a 72 ks XMM-Newton observation of the Be/X-ray pulsar (BeXRP) RX J0812.4-3114 in quiescence ($L_X \approx 1.6 \times 10^{33}~\mathrm{erg~s^{-1}}$). Intriguingly, we find a two component spectrum, with a hard power-law ($\Gamma \approx 1.5$) and a soft blackbody-like excess below $\approx 1~\mathrm{keV}$. The blackbody component is consistent in $kT$ with a prior quiescent Chandra observation reported by Tsygankov et al. and has an inferred blackbody radius of $\approx 10~\mathrm{km}$, consistent with emission from the entire neutron star (NS) surface. There is also mild evidence for an absorption line at $\approx 1~\mathrm{keV}$ and/or $\approx 1.4~\mathrm{keV}$. The hard component shows pulsations at $P \approx 31.908~\mathrm{s}$ (pulsed fraction $0.84 \pm 0.10$), agreeing with the pulse period seen previously in outbursts, but no pulsations were found in the soft excess (pulsed fraction $\lesssim 31\%$). We conclude that the pulsed hard component suggests low-level accretion onto the neutron star poles, while the soft excess seems to originate from the entire NS surface. We speculate that, in quiescence, the source switches between a soft thermal-dominated state (when the propeller effect is at work) and a relatively hard state with low-level accretion, and use the propeller cutoff to estimate the magnetic field of the system to be $\lesssim 8.4 \times 10^{11}~\mathrm{G}$. We compare the quiescent thermal $L_X$ predicted by the standard deep crustal heating model to our observations and find that RX J0812.4-3114 has a high thermal $L_X$, at or above the prediction for minimum cooling mechanisms. This suggests that RX J0812.4-3114 either contains a relatively low-mass NS with minimum cooling, or that the system may be young enough that the NS has not fully cooled from the supernova explosion.

We present the results of an XMM-Newton observation of the slowly rotating ($P = 3.4$ s), highly magnetized ($B \approx 3 \times 10^{13}$ G) radio pulsar PSR J0726-2612. A previous X-ray observation with the Chandra satellite showed that some of the properties of PSR J0726-2612 are similar to those of the X-ray Dim Isolated Neutron Stars (XDINSs), a small class of nearby slow pulsars characterized by purely thermal X-ray spectra and undetected in the radio band. We confirm the thermal nature of the X-ray emission of PSR J0726-2612, which can be fit by the sum of two blackbodies with temperatures $kT_1 = 0.074_{-0.011}^{+0.006}$ keV and $kT_2 = 0.14_{-0.02}^{+0.04}$ keV and emitting radii $R_1 = 10.4_{-2.8}^{+10.8}$ km and $R_2 = 0.5_{-0.3}^{+0.9}$ km, respectively (assuming a distance of 1 kpc). A broad absorption line modeled with a Gaussian profile centred at $0.39_{-0.03}^{+0.02}$ keV is required in the fit. The pulse profile of PSR J0726-2612 is characterized by two peaks with similar intensity separated by two unequal minima, a shape and pulsed fraction that cannot be reproduced without invoking magnetic beaming of the X-ray emission. The presence of a single radio pulse suggests that in PSR J0726-2612 the angles that the dipole axis and the line of sight make with the rotation axis, $\xi$ and $\chi$ respectively, are similar. This geometry differs from that of the two radio-silent XDINSs with a double peaked pulse profile similar to that of PSR J0726-2612, for which $\xi \sim 90^\circ$ and $\chi \sim 45^\circ$ have been recently estimated. These results strengthen the similarity between PSR J0726-2612 and the XDINSs and support the possibility that the lack of radio emission from the latter might simply be due to an unfavourable viewing geometry.

The role of the symmetry energy on the internal constitution and the global structure of a cold nonaccreted neutron star is studied using a set of unified equations of state. Based on the nuclear energy-density functional theory, these equations of state provide a thermodynamically consistent treatment of all regions of the star and were calculated using the four different Brussels-Montreal functionals BSk22, BSk24, BSk25 and BSk26. Our predictions are compared to various constraints inferred from astrophysical observations including the recent detection of the gravitational wave signal GW170817 from a binary neutron-star merger.

The theory of the nuclear energy-density functional is used to provide a unified and thermodynamically consistent treatment of all regions of cold non-accreting neutron stars. In order to assess the impact of our lack of complete knowledge of the density dependence of the symmetry energy on the constitution and the global structure of neutron stars, we employ four different functionals. All of them were precision fitted to essentially all the nuclear-mass data with the Hartree-Fock-Bogoliubov method and two different neutron-matter equations of state based on realistic nuclear forces. For each functional, we calculate the composition, the pressure-density relation, and the chemical potentials throughout the star. We show that uncertainties in the symmetry energy can significantly affect the theoretical results for the composition and global structure of neutron stars. To facilitate astrophysical applications, we construct analytic fits to our numerical results.

PSR B0943+10 is a mode-switching radio pulsar characterized by two emission modes with different radio and X-ray properties. Previous studies, based on simple combinations of blackbody and power law models, showed that its X-ray flux can be decomposed in a pulsed thermal plus an unpulsed non-thermal components. However, if PSR B0943+10 is a nearly aligned rotator seen pole-on, as suggested by the radio data, it is difficult to reproduce the high observed pulsed fraction unless magnetic beaming is included. In this work we reanalyze all the available X-ray observations of PSR B0943+10 with simultaneous radio coverage, modeling its thermal emission with polar caps covered by a magnetized hydrogen atmosphere or with a condensed iron surface. The condensed surface model provides good fits to the spectra of both pulsar modes, but, similarly to the blackbody, it can not reproduce the observed pulse profiles, unless an additional power law with an ad hoc modulation is added. Instead, the pulse profiles and phase-resolved spectra are well described using the hydrogen atmosphere model to describe the polar cap emission, plus an unpulsed power law. For the X-ray brighter state (Q-mode) we obtain a best fit with a temperature kT~0.09 keV, an emitting radius R~260 m, a magnetic field consistent with the value of the dipole field of 4x10^12 G inferred from the timing parameters, and a small angle between the magnetic and spin axis, $\xi$=5. The corresponding parameters for the X-ray fainter state (B-mode) are kT~0.08 keV and R~170 m.

As a first step toward a multi-phase equation of state for dense water, we develop a temperature-dependent equation of state for dense water covering the liquid and plasma regimes and extending to the super-ionic and gas regimes. This equation of state covers the complete range of conditions encountered in planetary modeling. We use first principles quantum molecular dynamics simulations and its Thomas-Fermi extension to reach the highest pressures encountered in giant planets several times the size of Jupiter. Using these results, as well as the data available at lower pressures, we obtain a parametrization of the Helmholtz free energy adjusted over this extended temperature and pressure domain. The parametrization ignores the entropy and density jumps at phase boundaries but we show that it is sufficiently accurate to model interior properties of most planets and exoplanets. We produce an equation of state given in analytical form that is readily usable in planetary modeling codes and dynamical simulations {\bf (a fortran implementation can be found at http://www.ioffe.ru/astro/H2O/)}. The EOS produced is valid for the entire density range relevant to planetary modeling, {\bf for densities where quantum effects for the ions can be neglected, and for temperatures below 50,000K. We use this equation of state to calculate the mass-radius relationship of exoplanets up to 5,000M_Earth, explore temperature effects in ocean and wet Earth-like planets, and quantify the influence of the water EOS for the core on the gravitational moments of Jupiter.

We study the relative importance of several recent updates of microphysics input to the neutron star cooling theory and the effects brought about by superstrong magnetic fields of magnetars, including the effects of the Landau quantization in their crusts. We use a finite-difference code for simulation of neutron-star thermal evolution on timescales from hours to megayears with an updated microphysics input. The consideration of short timescales ($\lesssim1$ yr) is made possible by a treatment of the heat-blanketing envelope without the quasistationary approximation inherent to its treatment in traditional neutron-star cooling codes. For the strongly magnetized neutron stars, we take into account the effects of Landau quantization on thermodynamic functions and thermal conductivities. We simulate cooling of ordinary neutron stars and magnetars with non-accreted and accreted crusts and compare the results with observations. Suppression of radiative and conductive opacities in strongly quantizing magnetic fields and formation of a condensed radiating surface substantially enhance the photon luminosity at early ages, making the life of magnetars brighter but shorter. These effects together with the effect of strong proton superfluidity, which slows down the cooling of kiloyear-aged neutron stars, can explain thermal luminosities of about a half of magnetars without invoking heating mechanisms. Observed thermal luminosities of other magnetars are still higher than theoretical predictions, which implies heating, but the effects of quantizing magnetic fields and baryon superfluidity help to reduce the discrepancy.

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 develop an analytical equation of state (EOS) for magnetized fully ionized plasmas, which cover a wide range of temperatures and densities, from low-density classical plasmas to relativistic, quantum plasma conditions. This EOS directly applies to calculations of structure and evolution of strongly magnetized white dwarfs and neutron stars. We review available analytical and numerical results for thermodynamic functions of the nonmagnetized and magnetized Coulomb gases, liquids, and solids. We propose a new analytical expression for the free energy of solid Coulomb mixtures. Based on recent numerical results, we construct analytical approximations for the thermodynamic functions of harmonic Coulomb crystals in quantizing magnetic fields. The analytical description ensures a consistent evaluation of all astrophysically important thermodynamic functions based on the first, second, and mixed derivatives of the free energy. Our numerical code for calculation of thermodynamic functions based on these approximations is made publicly available. Using this code, we calculate and discuss the effects of electron screening and magnetic quantization on the position of the melting point in a range of densities and magnetic fields relevant to white dwarfs and outer envelopes of neutron stars. We consider also the thermal and mechanical structure of a magnetar envelope and argue that it can have a frozen surface which covers the liquid ocean above the solid crust.

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 this 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 review the current status of the theory of thermal emission from the surface layers of neutron stars with strong magnetic fields $B\sim 10^{10}-10^{15}$ G, including formation of the spectrum in a partially ionized atmosphere and at a condensed surface. In particular, we describe recent progress in modeling partially ionized atmospheres of central compact objects in supernova remnants, which may have moderately strong fields $B\sim 10^{10}-10^{11}$ G. Special attention is given to polarization of thermal radiation emitted by a neutron star surface. Finally, we briefly describe applications of the theory to observations of thermally emitting isolated neutron stars.

We construct new models of outer heat blanketing envelopes of neutron stars composed of binary ion mixtures (H - He, He - C, C - Fe) in and out of diffusive equilibrium. To this aim, we generalize our previous work on diffusion of ions in isothermal gaseous or Coulomb liquid plasmas to handle non-isothermal systems. We calculate the relations between the effective surface temperature Ts and the temperature Tb at the bottom of heat blanketing envelopes (at a density rhob= 1e8 -- 1e10 g/cc) for diffusively equilibrated and non-equilibrated distributions of ion species at different masses DeltaM of lighter ions in the envelope. Our principal result is that the Ts - Tb relations are fairly insensitive to detailed distribution of ion fractions over the envelope (diffusively equilibrated or not) and depend almost solely on DeltaM. The obtained relations are approximated by analytic expressions which are convenient for modeling the evolution of neutron stars.

The goal of this work is to develop a simple analytic description of the emission properties (spectrum and polarization) of the condensed, strongly magnetized surface of neutron stars. We have improved the method of van Adelsberg et al. (2005) (arXiv:astro-ph/0406001) for calculating the spectral properties of condensed magnetized surfaces. Using the improved method, we calculate the reflectivity of an iron surface at magnetic field strengths B \sim (10^{12} - 10^{14}) G, with various inclinations of the magnetic field lines and radiation beam with respect to the surface and each other. We construct analytic expressions for the emissivity of this surface as functions of the photon energy, magnetic field strength, and the three angles that determine the geometry of the local problem. Using these expressions, we calculate X-ray spectra for neutron stars with condensed iron surfaces covered by thin partially ionized hydrogen atmospheres. We develop simple analytic descriptions of the intensity and polarization of radiation emitted or reflected by condensed iron surfaces of neutron stars with strong magnetic fields typical for isolated neutron stars. This description provides boundary conditions at the bottom of a thin atmosphere, which are more accurate than previously used approximations. The spectra calculated with this improvement show absorption features different from those in simplified models. The approach developed in this paper yields results that can facilitate modeling and interpretation of the X-ray spectra of isolated, strongly magnetized, thermally emitting neutron stars.

The early 21st century witnesses a dramatic rise in the study of thermal radiation of neutron stars. Modern space telescopes have provided a wealth of valuable information which, when properly interpreted, can elucidate the physics of superdense matter in the interior of these stars. This interpretation is necessarily based on the theory of formation of neutron star thermal spectra, which, in turn, is based on plasma physics and on the understanding of radiative processes in stellar photospheres. In this paper, the current status of the theory is reviewed with particular emphasis on neutron stars with strong magnetic fields. In addition to the conventional deep (semi-infinite) atmospheres, radiative condensed surfaces of neutron stars and "thin" (finite) atmospheres are considered.

Confronting theoretical models with observations of thermal radiation emitted by neutron stars is one of the most important ways to understand the properties of both, superdense matter in the interiors of the neutron stars and dense magnetized plasmas in their outer layers. Here we review the theory of thermal emission from the surface layers of strongly magnetized neutron stars, and the main properties of the observational data. In particular, we focus on the nearby sources for which a clear thermal component has been detected, without being contaminated by other emission processes (magnetosphere, accretion, nebulae). We also discuss the applications of the modern theoretical models of the formation of spectra of strongly magnetized neutron stars to the observed thermally emitting objects.

Context. A unified equation of state (EoS) should describe the crust and the core of a neutron star using the same physical model. The Brussels-Montreal group has recently derived a family of such EoSs based on the nuclear energy-density functional theory with generalized Skyrme effective forces, fitted to the available mass data. At the same time, these forces were constrained to reproduce microscopic calculations of homogeneous neutron matter based on realistic two- and three-nucleon forces. Aims. We represent basic physical characteristics of the latest Brussels-Montreal EoS models by analytical expressions to facilitate their inclusion in astrophysical simulations. Methods. We consider three EoS models, which significantly differ by stiffness: BSk19, BSk20, and BSk21. For each of them we constructed two versions of the EoS parametrization. In the first version, pressure P and gravitational mass density \rho are given as functions of the baryon number density n_b. In the second version, P, \rho, and n_b are given as functions of pseudo-enthalpy, which is useful for two-dimensional calculations of stationary rotating configurations of neutron stars. In addition to the EoS, we derived analytical expressions for several related quantities that are required in neutron-star simulations: number fractions of electrons and muons in the stellar core, nucleon numbers per nucleus in the inner crust, and equivalent radii and shape parameters of the nuclei in the inner crust. Results. We obtain analytical representations for the basic characteristics of the models of cold dense matter, which are most important for studies of neutron stars. We demonstrate the usability of our results by applying them to calculations of neutron-star mass-radius relations, maximum and minimum masses, thresholds of direct Urca processes, and the electron conductivity in the neutron-star crust.

The observed rapid cooling of the Cassiopeia A neutron star can be interpreted as being caused by neutron and proton transitions from normal to superfluid and superconducting states in the stellar core. Here we present two new Chandra ACIS-S Graded observations of this neutron star and measurements of the neutron star mass M and radius R found from consistent fitting of both the X-ray spectra and cooling behavior. This comparison is only possible for individual nuclear equations of state. We test phenomenological superfluid and superconducting gap models which mimic many of the known theoretical models against the cooling behavior. Our best-fit solution to the Cassiopeia A data is one in which the (M,R) = (1.44 Msun,12.6 km) neutron star is built with the BSk21 equation of state, strong proton superconductor and moderate neutron triplet superfluid gap models, and a pure iron envelope or a thin carbon layer on top of an iron envelope, although there are still large observational and theoretical uncertainties.

There is observational evidence that central compact objects (CCOs) in supernova remnants have moderately strong magnetic fields $B\sim10^{11}$ G. Meanwhile, available models of partially ionized hydrogen atmospheres of neutron stars with strong magnetic fields are restricted to $B\gtrsim10^{12}$ G. We extend the equation of state and radiative opacities, presented in previous papers for $10^{12}\mbox{ G}\lesssim B \lesssim 10^{15}$ G, to weaker fields. An equation of state and radiative opacities for a partially ionized hydrogen plasma are obtained at magnetic fields $B$, temperatures $T$, and densities $\rho$ typical for atmospheres of CCOs and other isolated neutron stars with moderately strong magnetic fields. The first- and second-order thermodynamic functions, monochromatic radiative opacities, and Rosseland mean opacities are calculated and tabulated, taking account of partial ionization, for $3\times10^{10}\mbox{ G}\lesssim B\lesssim 10^{12}$ G, $10^5$ K $\lesssim T\lesssim 10^7$ K, and a wide range of densities. Atmosphere models and spectra are calculated to verify the applicability of the results and to determine the range of magnetic fields and effective temperatures where the incomplete ionization of the hydrogen plasma is important.

Using 1D and 2D cooling codes we study thermal emission from neutron stars with steady state internal heaters of various intensities and geometries (blobs or spherical layers) located at different depths in the crust. The generated heat tends to propagate radially, from the heater down to the stellar core and up to the surface; it is also emitted by neutrinos. In local regions near the heater the results are well described with the 1D code. The heater's region projects onto the stellar surface forming a hot spot. There are two heat propagation regimes. In the first, conduction outflow regime (realized at heat rates $H_0 \lesssim 10^{20}$ erg cm$^{-3}$ s$^{-1}$ or temperatures $T_\mathrm{h} \lesssim 10^9$ K in the heater) the thermal surface emission of the star depends on the heater's power and neutrino emission in the stellar core. In the second, neutrino outflow regime ($H_0 \gtrsim 10^{20}$ erg cm$^{-3}$ s$^{-1}$ or $T_\mathrm{h} \gtrsim 10^9$ K) the surface thermal emission becomes independent of heater's power and the physics of the core. The largest (a few per cent) fraction of heat power is carried to the surface if the heater is in the outer crust and the heat regime is intermediate. The results can be used for modeling young cooling neutron stars (prior to the end of internal thermal relaxation), neutron stars in X-ray transients, magnetars and high-$B$ pulsars, as well as merging neutron stars.

We study the impact of plasma correlation effects on nonresonant thermonuclear reactions for various stellar objects, namely in the liquid envelopes of neutron stars, and the interiors of white dwarfs, low-mass stars, and substellar objects. We examine in particular the effect of electron screening on the enhancement of thermonuclear reactions in dense plasmas within and beyond the linear mixing rule approximation as well as the corrections due to quantum effects at high density. In addition, we examine some recent unconventional (Yukawa-potential and "quantum-tail") theoretical results on stellar thermonuclear fusions and show that these scenarios do not apply to stellar conditions.

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 thermonuclear 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 3\times10^5 K <~ T <~ 3\times10^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.

The data from all observations of RBS 1223 (1RXS J130848.6+212708) conducted by XMM-Newton EPIC pn with the same instrumental setup in 2003-2007 were combined to form spin-phase resolved spectra. A number of complex models of neutron stars with strongly magnetized ($B_{pole} 10^{13}-10^{14}$ G) surface, with temperature and magnetic field distributions around magnetic poles, and partially ionized hydrogen thin atmosphere above it have been implemented into the X-ray spectral fitting package XSPEC for simultaneous fitting of phase-resolved spectra. A Markov-Chain-Monte-Carlo (MCMC) approach is also applied to verify results of fitting and estimating in multi parameter models. The spectra in different rotational phase intervals and light curves in different energy bands with high S/N ratio show a high complexity. The spectra can be parameterized with a Gaussian absorption line superimposed on a blackbody spectrum, while the light curves with double-humped shape show strong dependence of pulsed fraction upon the energy band (13%-42%), which indicates that radiation emerges from at least two emitting areas. A model with condensed iron surface and partially ionized hydrogen thin atmosphere above it allows us to fit simultaneously the observed general spectral shape and the broad absorption feature observed at 0.3 keV in different spin phases of RBS 1223. It allowed to constrain some physical properties of X-ray emitting areas, i.e. the temperatures ($T_{p1} ~ 105$ eV, $T_{p2} ~ 99$ eV), magnetic field strengths ($B_{p1}= B_{p2} ~ 8.6x10^{13}$G) at the poles, and their distributions parameters ($a_{1} 0.61, a_{2} 0.29$, indicating an absence of strong toroidal magnetic field component). In addition, it puts some constraints on the geometry of the emerging X-ray emission and gravitational redshift ($z=0.16^{+0.03}_{-0.01}$) of RBS 1223.

Topical problems in the physics of and basic facts about neutron stars are reviewed. The observational manifestations of neutron stars, their core and envelope structure, magnetic fields, thermal evolution, and masses and radii are briefly discussed, along with the underlying microphysics.

Zaghloul [Phys. Plasmas 17, 062701 (2010); arXiv:1010.1161v1] reconsiders the occupation probability formalism in plasma thermodynamics and claims inconsistencies in previous models. I show that the origin of this incorrect claim is an omission of the configurational factor from the partition function. This arXiv version is supplemented with two appendices, where I add remarks and comments on two more recent publications of the same author on the same subject: on his response to this Comment [Phys. Plasmas 17, 124705 (2010)] and on his criticism towards the Hummer and Mihalas's (1988) formalism [Phys. Plasmas 17, 122903 (2010); arXiv:1010.1102v1].

We explore the thermal state of the neutron star in the Cassiopeia A supernova remnant using the recent result of Ho & Heinke (Nature, 462, 71 (2009)) that the thermal radiation of this star is well-described by a carbon atmosphere model and the emission comes from the entire stellar surface. Starting from neutron star cooling theory, we formulate a robust method to extract neutrino cooling rates of thermally relaxed stars at the neutrino cooling stage from observations of thermal surface radiation. We show how to compare these rates with the rates of standard candles -- stars with non-superfluid nucleon cores cooling slowly via the modified Urca process. We find that the internal temperature of standard candles is a well-defined function of the stellar compactness parameter $x=r_g/R$, irrespective of the equation of state of neutron star matter ($R$ and $r_g$ are circumferential and gravitational radii, respectively). We demonstrate that the data on the Cassiopeia A neutron star can be explained in terms of three parameters: $f_\ell$, the neutrino cooling efficiency with respect to the standard candle; the compactness $x$; and the amount of light elements in the heat blanketing envelope. For an ordinary (iron) heat blanketing envelope or a low-mass ($\lesssim 10^{-13}\,M_\odot$) carbon envelope, we find the efficiency $f_\ell \sim 1$ (standard cooling) for $x \lesssim 0.5$ and $f_\ell \sim 0.02$ (slower cooling) for a maximum compactness $x\approx 0.7$. A heat blanket containing the maximum mass ($\sim 10^{-8}\,M_\odot$) of light elements increases $f_\ell$ by a factor of 50. We also examine the (unlikely) possibility that the star is still thermally non-relaxed.

The X-ray spectra of some magnetized isolated neutron stars (NSs) show absorption features with equivalent widths (EWs) of 50 - 200 eV, whose nature is not yet well known. To explain the prominent absorption features in the soft X-ray spectra of the highly magnetized (B ~ 10^{14} G) X-ray dim isolated NSs (XDINSs), we theoretically investigate different NS local surface models, including naked condensed iron surfaces and partially ionized hydrogen model atmospheres, with semi-infinite and thin atmospheres above the condensed surface. We also developed a code for computing light curves and integral emergent spectra of magnetized neutron stars with various temperature and magnetic field distributions over the NS surface. We compare the general properties of the computed and observed light curves and integral spectra for XDINS RBS\,1223 and conclude that the observations can be explained by a thin hydrogen atmosphere above the condensed iron surface, while the presence of a strong toroidal magnetic field component on the XDINS surface is unlikely. We suggest that the harmonically spaced absorption features in the soft X-ray spectrum of the central compact object (CCO) 1E 1207.4-5209 (hereafter 1E 1207) correspond to peaks in the energy dependence of the free-free opacity in a quantizing magnetic field, known as quantum oscillations. To explore observable properties of these quantum oscillations, we calculate models of hydrogen NS atmospheres with B ~ 10^{10} - 10^{11} G (i.e., electron cyclotron energy E_{c,e} ~ 0.1 - 1 keV) and T_eff = 1 - 3 MK. Such conditions are thought to be typical for 1E 1207. We show that observable features at the electron cyclotron harmonics with EWs \approx 100 - 200 eV can arise due to these quantum oscillations.

In the X-ray spectra of most X-ray dim isolated neutron stars (XDINSs) absorption features with equivalent widths (EWs) of 50 -- 200 eV are observed. We theoretically investigate different models to explain absorption features and compare their properties with the observations. We consider various theoretical models for the magnetized neutron star surface: naked condensed iron surfaces and partially ionized hydrogen model atmospheres, including semi-infinite and thin atmospheres above a condensed surface. The properties of the absorption features (especially equivalent widths) and the angular distributions of the emergent radiation are described for all models. A code for computing light curves and integral emergent spectra of magnetized neutron stars is developed. We assume a dipole surface magnetic field distribution with a possible toroidal component and corresponding temperature distribution. A model with two uniform hot spots at the magnetic poles can also be employed. Light curves and spectra of highly magnetized neutron stars with parameters typical for XDINSs are computed using different surface temperature distributions and various local surface models. Spectra of magnetized model atmospheres are approximated by diluted blackbody spectra with one or two Gaussian lines having parameters, which allow us to describe the model absorption features. To explain the prominent absorption features in the soft X-ray spectra of XDINSs a thin atmosphere above the condensed surface can be invoked, whereas a strong toroidal magnetic field component on the XDINS surfaces can be excluded.

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 relativistic 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 Coulomb 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 present a new computer code for modeling magnetized neutron star atmospheres in a wide range of magnetic fields (10^{12} - 10^{15} G) and effective temperatures (3 \times 10^5 - 10^7 K). The atmosphere is assumed to consist either of fully ionized electron-ion plasmas or of partially ionized hydrogen. Vacuum resonance and partial mode conversion are taken into account. Any inclination of the magnetic field relative to the stellar surface is allowed. We use modern opacities of fully or partially ionized plasmas in strong magnetic fields and solve the coupled radiative transfer equations for the normal electromagnetic modes in the plasma. Using this code, we study the possibilities to explain the soft X-ray spectra of isolated neutron stars by different atmosphere models. In particular, the outgoing spectrum using the "sandwich" model (thin atmosphere with a hydrogen layer above a helium layer) is constructed. Thin partially ionized hydrogen atmospheres with vacuum polarization are shown to be able to improve our understanding of the observed spectrum of the nearby isolated neutron star RBS 1223 (RX J1308.8+2127).

We present an equation of state and radiative opacities for a strongly magnetized hydrogen plasma at magnetic fields B, temperatures T, and densities \rho typical for atmospheres of isolated neutron stars. The first- and second-order thermodynamic functions, monochromatic radiative opacities, and Rosseland mean opacities are calculated and tabulated, taking account of partial ionization, for 8x10^{11} G < B < 3x10^{13} G, 2x10^5 K < T < 10^7 K, and a wide range of \rho. We show that bound-bound and bound-free transitions give an important contribution to the opacities at T < (1-5)x10^6 K in the considered range of B in the outer neutron-star atmosphere layers, which may substantially modify the X-ray spectrum of a typical magnetized neutron star. In addition, we re-evaluate opacities due to free-free transitions, taking into account the motion of both interacting particles, electron and proton, in a strong magnetic field. Compared to the previous neutron-star atmosphere models, the free-free absorption is strongly suppressed at photon frequencies below the proton cyclotron frequency. The latter result holds for any field strength, which prompts a revision of existing theoretical models of X-ray spectra of magnetar atmospheres. Note added in v4: for free-free photoabsorption (Appendix B), equivalent but much simpler expressions are available in Appendix B of Potekhin & Lai, arXiv:astro-ph/0701285.

Observed X-ray spectra of some isolated magnetized neutron stars display absorption features, sometimes interpreted as ion cyclotron lines. Modeling the observed spectra is necessary to check this hypothesis and to evaluate neutron star parameters.We develop a computer code for modeling magnetized neutron star atmospheres in a wide range of magnetic fields (10^{12} - 10^{15} G) and effective temperatures (3 \times 10^5 - 10^7 K). Using this code, we study the possibilities to explain the soft X-ray spectra of isolated neutron stars by different atmosphere models. The atmosphere is assumed to consist either of fully ionized electron-ion plasmas or of partially ionized hydrogen. Vacuum resonance and partial mode conversion are taken into account. Any inclination of the magnetic field relative to the stellar surface is allowed. We use modern opacities of fully or partially ionized plasmas in strong magnetic fields and solve the coupled radiative transfer equations for the normal electromagnetic modes in the plasma. Spectra of outgoing radiation are calculated for various atmosphere models: fully ionized semi-infinite atmosphere, thin atmosphere, partially ionized hydrogen atmosphere, or novel "sandwich" atmosphere (thin atmosphere with a hydrogen layer above a helium layer. Possibilities of applications of these results are discussed. In particular, the outgoing spectrum using the "sandwich" model is constructed. Thin partially ionized hydrogen atmospheres with vacuum polarization are shown to be able to improve the fit to the observed spectrum of the nearby isolated neutron star RBS 1223 (RX J1308.8+2127).

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 of 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-component 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.

We construct partially ionized hydrogen atmosphere models for magnetized neutron stars in radiative equilibrium with fixed surface fields between B=10^12 and 2x10^13 G and effective temperatures logT_eff=5.5-6.8, as well as with surface B and T_eff distributions around these values. The models are based on the latest equation of state and opacity results for magnetized, partially ionized hydrogen plasmas. The atmospheres directly determine the characteristics of thermal emission from the surface of neutron stars. We also incorporate these model spectra into XSPEC, under the model name NSMAX, thus allowing them to be used by the community to fit X-ray observations of neutron stars.

Recent observations of thermally emitting isolated neutron stars revealed spectral features that could be interpreted as radiative transitions of He in a magnetized neutron star atmosphere. We present Hartree-Fock calculations of the polarization-dependent photoionization cross sections of the He atom in strong magnetic fields ranging from 10^12 G to 10^14 G. Convenient fitting formulae for the cross sections are given as well as related oscillator strengths for various bound-bound transitions. The effects of finite nucleus mass on the radiative absorption cross sections are examined using perturbation theory.

We review current state of neutron star cooling theory and discuss the prospects to constrain the equation of state, neutrino emission and superfluid properties of neutron star cores by comparing the cooling theory with observations of thermal radiation from isolated neutron stars.

Strong (B >> 10^9 G) and superstrong (B > 10^{14} G) magnetic fields profoundly affect many thermodynamic and kinetic characteristics of dense plasmas in neutron star envelopes. In particular, they produce strongly anisotropic thermal conductivity in the neutron star crust and modify the equation of state and radiative opacities in the atmosphere, which are major ingredients of the cooling theory and spectral atmosphere models. As a result, both the radiation spectrum and the thermal luminosity of a neutron star can be affected by the magnetic field. We briefly review these effects and demonstrate the influence of magnetic field strength on the thermal structure of an isolated neutron star, putting emphasis on the differences brought about by the superstrong fields and high temperatures of magnetars. For the latter objects, it is important to take proper account of a combined effect of the magnetic field on thermal conduction and neutrino emission at densities \rho > 10^{10} g cm^{-3}. We show that the neutrino emission puts a B-dependent upper limit on the effective surface temperature of a cooling neutron star.

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.

We review the theory of electron-conduction opacity, a fundamental ingredient in the computation of low-mass stellar models; shortcomings and limitations of the existing calculations used in stellar evolution are discussed. We then present new determinations of the electron-conduction opacity in stellar conditions for an arbitrary chemical composition, that improve over previous works and, most importantly, cover the whole parameter space relevant to stellar evolution models (i.e., both the regime of partial and high electron degeneracy). A detailed comparison with the currently used tabulations is also performed. The impact of our new opacities on the evolution of low-mass stars is assessed by computing stellar models along both the H- and He-burning evolutionary phases, as well as Main Sequence models of very low-mass stars and white dwarf cooling tracks.

We calculate the transition rates between proton Landau levels due to non-radiative and radiative Coulomb collisions in an electron-proton plasma with strong magnetic field B. Both electron-proton collisions and proton-proton collisions are considered. The roles of the first-order cyclotron absorption and second-order free-free absorption and scattering in determining the line strength and shape as well as the continuum are analysed in detail. We solve the statistical balance equation for the populations of proton Landau levels. For temperatures \sim 10^6-10^7 K, the deviations of the proton populations from LTE are appreciable at density \rho < 0.1 B_{14}^{3.5} g cm^{-3}, where B_{14}=B/(10^{14} G). We present general formulae for the plasma emissivity and absorption coefficents under a wide range of physical conditions. Our results are useful for studying the possibility and the conditions of proton/ion cyclotron line formation in the near vicinity of highly magnetized neutron stars.

RX J1856.5-3754 is one of the brightest nearby isolated neutron stars, and considerable observational resources have been devoted to it. However, current models are unable to satisfactorily explain the data. We show that our latest models of a thin, magnetic, partially ionized hydrogen atmosphere on top of a condensed surface can fit the entire spectrum, from X-rays to optical, of RX J1856.5-3754, within the uncertainties. In our simplest model, the best-fit parameters are an interstellar column density N_H \approx 1x10^20 cm^-2 and an emitting area with R^infty \approx 17 km (assuming a distance to RX J1856.5-3754 of 140 pc), temperature T^infty \approx 4.3x10^5 K, gravitational redshift z_g \sim 0.22, atmospheric hydrogen column y_H \approx 1 g cm^-2, and magnetic field B \approx (3-4)x10^12 G; the values for the temperature and magnetic field indicate an effective average over the surface. We also calculate a more realistic model, which accounts for magnetic field and temperature variations over the neutron star surface as well as general relativistic effects, to determine pulsations; we find there exist viewing geometries that produce pulsations near the currently observed limits. The origin of the thin atmospheres required to fit the data is an important question, and we briefly discuss mechanisms for producing these atmospheres. Our model thus represents the most self-consistent picture to date for explaining all the observations of RX J1856.5-3754.

We study thermal structure and evolution of magnetars as cooling neutron stars with a phenomenological heat source in a spherical internal layer. We explore the location of this layer as well as the heating rate that could explain high observable thermal luminosities of magnetars and would be consistent with the energy budget of neutron stars. We conclude that the heat source should be located in an outer magnetar's crust, at densities rho < 5e11 g/cm^3, and should have the heat intensity of the order of 1e20 erg/s/cm^3. Otherwise the heat energy is mainly emitted by neutrinos and cannot warm up the surface.

We study the equation of state, polarization and radiation properties for nonideal, strongly magnetized plasmas which compose outer envelopes of magnetic neutron stars. Detailed calculations are performed for partially ionized hydrogen atmospheres and for condensed hydrogen or iron surfaces of these stars. This is a companion paper to astro-ph/0511803

We briefly examine the properties of dense plasmas characteristic of the interior of giant planets and the atmospheres of neutron stars. Special attention is devoted to the equation of state of hydrogen and helium at high density and to the effect of magnetic fields on the properties of dense matter.

Equation of state for partially ionized carbon at temperatures T > ~ 10^5 K is calculated in 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 includes the internal partition functions of bound species. The latter 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.

We study the relation between the mean effective surface temperature T_s and the internal temperature T_b for magnetic neutron stars, assuming that the magnetic field near the surface has a presumably small-scale structure. The heavy-element (iron) and light-element (accreted) heat-blanketing envelopes are considered, and the results are compared with the case of a dipole magnetic field. We argue that the difference in the T_b(T_s)-relation for different magnetic configurations is always much smaller than a possible difference caused by variations of the chemical composition in the envelope.

Recent observations show that the thermal X-ray spectra of many isolated neutron stars are featureless and in some cases (e.g., RX J1856.5-3754) well fit by a blackbody. Such a perfect blackbody spectrum is puzzling since radiative transport through typical neutron star atmospheres causes noticeable deviation from blackbody. Previous studies have shown that in a strong magnetic field, the outermost layer of the neutron star may be in a condensed solid or liquid form because of the greatly enhanced cohesive energy of the condensed matter. The critical temperature of condensation increases with the magnetic field strength, and can be as high as 10^6 K (for Fe surface at B \sim 10^{13} G or H surface at B \sim a few times 10^{14} G). Thus the thermal radiation can directly emerge from the degenerate metallic condensed surface, without going through a gaseous atmosphere. Here we calculate the emission properties (spectrum and polarization) of the condensed Fe and H surfaces of magnetic neutron stars in the regimes where such condensation may be possible. For a smooth condensed surface, the overall emission is reduced from the blackbody by less than a factor of 2. The spectrum exhibits modest deviation from blackbody across a wide energy range, and shows mild absorption features associated with the ion cyclotron frequency and the electron plasma frequency in the condensed matter. The roughness of the solid condensate (in the Fe case) tends to decrease the reflectivity of the surface, and make the emission spectrum even closer to blackbody. We discuss the implications of our results for observations of dim, isolated neutron stars and magnetars.

We construct hydrogen atmosphere models for strongly magnetized neutron stars in thermodynamic equilibrium, taking into account partial ionization. The presence of bound states affects the equation of state, absorption coefficients, and polarizability tensor of a strongly magnetized plasma. Therefore the partial ionization influences the polarization vectors and opacities of normal electromagnetic waves, and thus the spectra of outgoing radiation. Here we review a model suitable for the most typical neutron-star atmospheres and focus on the problems that remain to be solved for its extension to other atmospheric parameters.

The impact of nuclear physics theories on cooling of isolated neutron stars is analyzed. Physical properties of neutron star matter important for cooling are reviewed such as composition, the equation of state, superfluidity of various baryon species, neutrino emission mechanisms. Theoretical results are compared with observations of thermal radiation from neutron stars. Current constraints on theoretical models of dense matter, derived from such a comparison, are formulated.

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.

Polarizability tensor of a strongly magnetized plasma and the polarization vectors and opacities of normal electromagnetic waves are studied for the conditions typical of neutron star atmospheres, taking account of partial ionization effects. Vacuum polarization is also included using a new set of fitting formulae that are accurate for wide range of field strengths. The full account of the coupling of the quantum mechanical structure of the atoms to their center-of-mass motion across the magnetic field is shown to be crucial for the correct evaluation of the polarization properties and opacities of the plasma. The self-consistent treatment of the polarizability and absorption coefficients proves to be necessary if the ionization degree of the plasma is low, which can occur in the atmospheres of cool or ultramagnetized neutron stars. Atmosphere models and spectra are presented to illustrate the importance of such self-consistent treatment.

We briefly review recent theoretical studies of the effects of accreted envelopes, magnetic fields and crustal superfluidity on the cooling of neutron stars. These effects are especially important for slowly cooling low-mass neutron stars, where direct Urca process of neutrino emission is forbidden. The effects are useful for interpretation of observations of several isolated middle-aged neutron stars.

Analytical approximations are constructed for binding energies, quantum-mechanical sizes and oscillator strengths of main radiative transitions of hydrogen atoms arbitrarily moving in magnetic fields 10^{12}-10^{13} G. Examples of using the obtained approximations for determination of maximum transverse velocity of an atom and for evaluation of absorption spectra in magnetic neutron star atmospheres are presented.

We calculate the thermal structure and quiescent thermal luminosity of accreting neutron stars (warmed by deep crustal heating in accreted matter) in soft X-ray transients (SXTs). We consider neutron stars with nucleon and hyperon cores and with accreted envelopes. It is assumed that an envelope has an outer helium layer (of variable depth) and deeper layers of heavier elements, either with iron or with much heavier nuclei (of atomic weight A > 100) on the top (Haensel & Zdunik 1990, 2003, astro-ph/0305220). The relation between the internal and surface stellar temperatures is obtained and fitted. The quiescent luminosity of the hottest (low-mass) and coldest (high-mass) neutron stars is calculated, together with the ranges of its possible variations due to variable thickness of the helium layer. The results are compared with observations of SXTs, particularly, containing the coldest (SAX J1808.4-3658) and the hottest (Aql X-1) neutron stars. The observations of SAX J1808.4-3658 in a quiescent state on March 24, 2001 (Campana et al. 2002, astro-ph/0206376) can be explained only if this SXT contains a massive neutron star with a nucleon/hyperon core; a hyperon core with a not too low fraction of electrons is preferable. Future observations may discriminate between the various models of hyperon/nucleon dense matter. The thermal emission of SAX J1808.4-3658 is also sensitive to the models of plasma ionization in the outermost surface layers and can serve for testing such models.

We study the thermal structure of neutron stars with magnetized envelopes composed of accreted material, using updated thermal conductivities of plasmas in quantizing magnetic fields, as well as equation of state and radiative opacities for partially ionized hydrogen in strong magnetic fields. The relation between the internal and local surface temperatures is calculated and fitted by an analytic function of the internal temperature, magnetic field strength, angle between the field lines and the normal to the surface, surface gravity, and the mass of the accreted material. The luminosity of a neutron star with a dipole magnetic field is calculated for various values of the accreted mass, internal temperature, and magnetic field strength. Using these results, we simulate cooling of superfluid neutron stars with magnetized accreted envelopes. We consider slow and fast cooling regimes, paying special attention to very slow cooling of low-mass superfluid neutron stars. In the latter case, the cooling is strongly affected by the combined effect of magnetized accreted envelopes and neutron superfluidity in the stellar crust. Our results are important for interpretation of observations of isolated neutron stars hottest for their age, such as RX J0822-43 and PSR B1055-52.

The equation of state and radiative opacities of partially ionized, strongly magnetized hydrogen plasmas, presented in a previous paper [ApJ 585, 955 (2003), astro-ph/0212062] for the magnetic field strengths 8.e11 G < B < 3.e13 G, are extended to the field strengths 3.e13 G < B < 1.e15 G, relevant for magnetars. The first- and second-order thermodynamic functions and radiative opacities are calculated and tabulated for 5.e5 < T < 4.e7 K in a wide range of densities. We show that bound-free transitions give an important contribution to the opacities in the considered range of B in the outer neutron-star atmosphere layers. Unlike the case of weaker fields, bound-bound transitions are unimportant.

We construct partially ionized hydrogen atmosphere models for magnetized neutron stars in radiative equilibrium with surface fields B=10^12-5 \times 10^14 G and effective temperatures T_eff \sim a few \times 10^5-10^6 K. These models are based on the latest equation of state and opacity results for magnetized, partially ionized hydrogen plasmas that take into account various magnetic and dense medium effects. The atmospheres directly determine the characteristics of thermal emission from isolated neutron stars. For the models with B=10^12-10^13 G, the spectral features due to neutral atoms lie at extreme UV and very soft X-ray energy bands and therefore are difficult to observe. However, the continuum flux is also different from the fully ionized case, especially at lower energies. For the superstrong field models (B\ga 10^14 G), we show that the vacuum polarization effect not only suppresses the proton cyclotron line as shown previously, but also suppresses spectral features due to bound species; therefore spectral lines or features in thermal radiation are more difficult to observe when the neutron star magnetic field is \ga 10^14 G.

The cooling theory of isolated neutron stars is reviewed. The main cooling regulators are discussed, first of all, operation of direct Urca process (or similar processes in exotic phases of dense matter) and superfluidity in stellar interiors. The prospects to constrain gross parameters of supranuclear matter in neutron-star interiors by confronting cooling theory with observations of isolated neutron stars are outlined. A related problem of thermal states of transiently accreting neutron stars with deep crustal heating of accreted matter is discussed in application to soft X-ray transients.

We construct hydrogen atmosphere models for magnetized neutron stars in radiative equilibrium with surface fields B=10^12-5x10^14 G and effective temperatures T_eff a few x 10^5-10^6 K by solving the full radiative transfer equations for both polarization modes in the magnetized hydrogen plasma. The atmospheres directly determine the characteristics of thermal emission from isolated neutron stars. We study the effects of vacuum polarization and bound atoms on the atmosphere structure and spectra. For the lower magnetic field models (B 10^12 G), the spectral features due to neutral atoms lie at extreme UV and very soft X-ray energies and therefore are not likely to be observed. However, the continuum flux is also different from the fully ionized case, especially at lower energies. For the higher magnetic field models, we find that vacuum polarization softens the high energy tail of the thermal spectrum. We show that this depression of continuum flux strongly suppresses not only the proton cyclotron line but also spectral features due to bound species; therefore spectral lines or features in thermal radiation are more difficult to observe when the neutron star magnetic field is > 10^14 G.

We briefly examine the properties of dense plasmas characteristic of the atmospheres of neutron stars and of the interior of massive white dwarfs. These astrophysical bodies are natural laboratories to study respectively the problem of pressure ionization of hydrogen in a strong magnetic field and the crystallization of the quantum one-component-plasma at finite temperature.

We study the thermodynamic properties of a partially ionized hydrogen plasma in strong magnetic fields, B ~ 10^{12}-10^{13} G, typical of neutron stars. The properties of the plasma depend significantly on the quantum-mechanical sizes and binding energies of the atoms, which are strongly modified by thermal motion across the field. We use new fitting formulas for the atomic binding energies and sizes, based on accurate numerical calculations and valid for any state of motion of the atom. In particular, we take into account decentered atomic states, neglected in previous studies of thermodynamics of magnetized plasmas. We also employ analytic fits for the thermodynamic functions of nonideal fully ionized electron-ion Coulomb plasmas. This enables us to construct an analytic model of the free energy. An ionization equilibrium equation is derived, taking into account the strong magnetic field effects and the nonideality effects. This equation is solved by an iteration technique. Ionization degrees, occupancies, and the equation of state are calculated.

Cooling of neutron stars (NSs) with the cores composed of neutrons, protons, and electrons is analyzed. The main cooling regulators are discussed: opening of direct Urca process in a NS central kernel; superfluidity of nucleons in NS interiors; surface layers of light (accreted) elements; strong surface magnetic fields. An emphasis is paid on the cooling scenario with strong $^1$S$_0$ pairing of protons and weak $^3$P$_2$ pairing of neutrons in the NS core, as well as strong $^1$S$_0$ pairing of neutrons in the NS crust. The theory predicts three types of isolated cooling middle-aged NSs with distinctly different properties: low-mass, slowly cooling NSs; medium-mass, moderately cooling NSs; massive, rapidly cooling NSs. The theory is compared with observations of isolated NSs -- pulsars and radio quiet NSs in supernova remnants. The constraints on physical properties of NSs which can be inferred from such a comparison are outlined.

The analytic equation of state of nonideal Coulomb plasmas consisting of pointlike ions immersed in a polarizable electron background (physics/9807042) is improved, and its applicability range is considerably extended. First, the fit of the electron screening contribution in the free energy of the Coulomb liquid is refined at high densities where the electrons are relativistic. Second, we calculate the screening contribution for the Coulomb solid (bcc and fcc) and derive an analytic fitting expression. Third, we propose a simple approximation to the internal and free energy of the liquid one-component plasma of ions, accurate within the numerical errors of the most recent Monte Carlo simulations. We obtain an updated value of the coupling parameter at the solid-liquid phase transition for the one-component plasma: Gamma_m = 175.0 (+/- 0.4).

The low-frequency electric microfield distribution in a Coulomb plasma is calculated for various plasma parameters, from weak to strong Coulomb coupling and from zero to strong electron screening. Two methods of numerical calculations are employed: the adjustable-parameter exponential approximation and the Monte Carlo simulation. The results are represented by analytic fitting formulas suitable for applications.

Phonon frequency moments and thermodynamic functions (electrostatic and vibrational parts of the free energy, internal energy, and heat capacity) are calculated for bcc and fcc Coulomb crystals in the harmonic approximation with a fractional accuracy within one part on 10^5. Temperature dependence of thermodynamic functions is fitted by analytic formulas with an accuracy of a few parts in 10^5. The static-lattice (Madelung) part of the free energy is calculated with an accuracy about 10^{-12}. The Madelung constant and frequency moments of hcp crystals are also computed.

The thermal structure of neutron stars with magnetized envelopes is studied using modern physics input. The relation between the internal (T_i) and local surface temperatures is calculated and fitted by analytic expressions for magnetic field strengths B from 0 to 10^{16} G and arbitrary inclination of the field lines to the surface. The luminosity of a neutron star with dipole magnetic field is calculated and fitted as a function of B, T_i, stellar mass and radius. In addition, we simulate cooling of neutron stars with magnetized envelopes. In particular, we analyse ultramagnetized envelopes of magnetars and also the effects of the magnetic field of the Vela pulsar on the determination of critical temperatures of neutron and proton superfluids in its core.

The quadrupole moment of a hydrogen atom in a magnetic field for field strengths from 0 to 4.414e13 G is calculated by two different methods. The first method is variational, and based on a single trial function. The second method deals with a solution of the Schroedinger equation in the form of a linear combination of Landau orbitals.

The thermal structure of neutron star envelopes is discussed with emphasis on analytic results. Recent progress on the effect of chemical constitution and high magnetic fields on the opacities and the thermal structure is further reviewed in view of the application to pulsar cooling and magnetars.

Matter of subnuclear density in the inner crusts of neutron stars consists of neutron-rich atomic nuclei immersed in strongly degenerate relativistic gas of electrons and strongly nonideal liquid of neutrons. We give a brief account of thermodynamic and kinetic properties of this matter which are greatly affected by Coulomb and nuclear interactions (a companion paper on the outer envelopes is astro-ph/0012316) and show how they can be studied from observations of thermal radiation of young (age < 100 yr) neutron stars (for a discussion of the latter possibility, see astro-ph/0012306).

Outer envelopes of neutron stars consist mostly of fully ionized, strongly coupled Coulomb plasmas characterized by typical densities about 10^4-10^{11} g/cc and temperatures about 10^4-10^9 K. Many neutron stars possess magnetic fields about 10^{11}-10^{14} G. Here we briefly review recent theoretical advances which allow one to calculate thermodynamic functions and electron transport coefficients for such plasmas with an accuracy required for theoretical interpretation of observations.

The internal properties of the neutron star crust can be probed by observing the epoch of thermal relaxation. After the supernova explosion, powerful neutrino emission quickly cools the stellar core, while the crust stays hot. The cooling wave then propagates through the crust, due to its finite thermal conductivity. When the cooling wave reaches the surface (age 10-100 yr), the effective temperature drops sharply from 250 eV to 30 or 100 eV, depending on the cooling model. The crust relaxation time is sensitive to the (poorly known) microscopic properties of matter of subnuclear density, such as the heat capacity, thermal conductivity, and superfluidity of free neutrons. We calculate the cooling models with the new values of the electron thermal conductivity in the inner crust, based on a realistic treatment of the shapes of atomic nuclei. Superfluid effects may shorten the relaxation time by a factor of 4. The comparison of theoretical cooling curves with observations provides a potentially powerful method of studying the properties of the neutron superfluid and highly unusual atomic nuclei in the inner crust.

We discuss experiments and observations aimed at testing the possible space-time variability of fundamental physical constants, predicted by the modern theory. Specifically, we consider two of the dimensionless physical parameters which are important for atomic and molecular physics: the fine-structure constant and the electron-to-proton mass ratio. We review the current status of such experiments and critically analyze recent claims of a detection of the variability of the fine-structure constant on the cosmological time scale. We stress that such a detection remains to be checked by future experiments and observations. The tightest of the firmly established upper limits read that the considered constants could not vary by more than 0.015% on the scale of 10 Gigayears.

One of the topical problems of contemporary physics is a possible variability of the fundamental constants. Here we consider possible variability of two dimensionless constants which are most important for calculation of atomic and molecular spectra (in particular, the X-ray ones): the fine-structure constant \alpha=e^2/\hbar c and the proton-to-electron mass ratio \mu=m_p/m_e. Values of the physical constants in the early epochs are estimated directly from observations of quasars - the most powerful sources of radiation, whose spectra were formed when the Universe was several times younger than now. A critical analysis of the available results leads to the conclusion that present-day data do not reveal any statistically significant evidence for variations of the fundamental constants under study. The most reliable upper limits to possible variation rates at the 95% confidence level, obtained in our work, read: |\dot\alpha/\alpha| < (1.4e-14)/yr, |\dot\mu/\mu| < (1.5e-14)/yr on the average over the last ten billion years.

Practical expressions are derived for evaluation of electrical and thermal conductivities and thermopower of degenerate electrons in the outer envelopes of neutron stars with magnetic fields. All tensor components of the kinetic coefficients are calculated (those related to conduction along and across magnetic field and to the Hall currents). The kinetic coefficients are presented as energy averages of expressions containing energy dependent effective relaxation times of two types, associated either with longitudinal or with transverse currents. The calculation is based on the effective scattering potential proposed in the previous paper (astro-ph/9903127), which describes the electron-ion and electron-phonon scattering, taking into account correlation effects in strongly coupled Coulomb liquid and multi-phonon scattering in Coulomb crystal, respectively. Analytic fitting formulae are devised for the effective relaxation times at arbitrary field strength. Basing on these results, we calculate the transport coefficients at various temperatures, densities, and magnetic fields pertinent to the neutron star envelopes.

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.

We study partially ionized hydrogen plasma at temperatures T = 10^{5.5}-10^{6.5} K and magnetic fields B = 10^{12}-10^{13} G, typical of photospheres of isolated neutron stars. We construct an analytic model of the plasma free energy and derive a generalized Saha equation which is used to obtain the opacities.

Neutrino-pair bremsstrahlung by relativistic degenerate electrons in a neutron-star crust at densities (10^9 - 1.5x10^{14}) g/cm^3 is analyzed. The processes taken into account are neutrino emission due to Coulomb scattering of electrons by atomic nuclei in a Coulomb liquid, and electron-phonon scattering and Bragg diffraction (the static-lattice contribution) in a Coulomb crystal. The static-lattice contribution is calculated including the electron band-structure effects for cubic Coulomb crystals of different types and also for the liquid crystal phases composed of rod- and plate-like nuclei in the neutron-star mantle (at 10^{14} - 1.5x10^{14} g/cm^3). The phonon contribution is evaluated with proper treatment of the multi-phonon processes which removes a jump in the neutrino bremsstrahlung emissivity at the melting point obtained in previous works. Below 10^{13} g/cm^3, the results are rather insensitive to the nuclear form factor, but results for the solid state near the melting point are affected significantly by the Debye-Waller factor and multi-phonon processes. At higher densities, the nuclear form factor becomes more significant. A comparison of the various neutrino generation mechanisms in neutron star crusts shows that electron bremsstrahlung is among the most important ones.

We study hydrogen plasmas at magnetic fields B ~ 10^{12}-10^{13} Gauss, densities ~ 10^{-3}-10^3 g/cc and temperatures T ~ 10^{5.5}-10^{6.5} K, typical of photospheres of middle-aged cooling neutron stars. We construct an analytical free energy model of the partially ionized plasma, including into consideration the decentred atomic states, which arise due to the thermal motion across the strong field. We show that these states, neglected in previous studies, may contribute appreciably into thermodynamics of the outer atmospheric layers at density below 1 g/cc and typical B and T. We take into account Coulomb non-ideality of the ionized component of the plasma affected by intense magnetic field. Ionization degree, occupancies and equation of state are calculated, and their dependences on the temperature, density and magnetic field are studied.

The dynamical structure factor of a Coulomb crystal of ions is calculated at arbitrary temperature below the melting point taking into account multi-phonon processes in the harmonic approximation. In a strongly coupled Coulomb ion liquid, the static structure factor is split into two parts, a Bragg-diffraction-like one, describing incipient long-range order structures, and an inelastic part corresponding to thermal ion density fluctuations. It is assumed that the diffractionlike scattering does not lead to the electron relaxation in the liquid phase. This assumption, together with the inclusion of multi-phonon processes in the crystalline phase, eliminates large discontinuities of the transport coefficients (jumps of the thermal and electric conductivities, as well as shear viscosity, reported previously) at a melting point.

Endeavours of the unification of the four fundamental interactions have resulted in a development of theories having cosmological solutions in which low-energy limits of fundamental physical constants vary with time. The validity of such theoretical models should be checked by comparison of the theoretical predictions with observational and experimental bounds on possible time-dependences of the fundamental constants. Based on high-resolution measurements of quasar spectra, we obtain the following direct limits on the average rate of the cosmological time variation of the fine-structure constant \alpha: |\dot{\alpha}/\alpha| < 1.9 \times 10^{-14} yr^{-1} is the most likely limit, and |\dot{\alpha}/\alpha| < 3.1 \times 10^{-14} yr^{-1} is the most conservative limit. Analogous estimates published previously, as well as other contemporary tests for possible variations of \alpha (those based on the "Oklo phenomenon", on the primordial nucleosynthesis models, and others) are discussed and compared with the present upper limit. We argue that the present result is the most conservative one.

Thermodynamic quantities of Coulomb plasmas consisting of point-like ions immersed in a compressible, polarizable electron background are calculated for ion charges Z=1 to 26 and for a wide domain of plasma parameters ranging from the Debye-Hueckel limit to the crystallization point and from the region of nondegenerate to fully degenerate nonrelativistic or relativistic electrons. The calculations are based on the linear-response theory for the electron-ion interaction, including the local-field corrections in the electronic dielectric function. The thermodynamic quantities are calculated in the framework of the N-body hypernetted-chain equations and fitted by analytic expressions. We present also accurate analytic approximations for the free energy of the ideal electron gas at arbitrary degeneracy and relativity and for the excess free energy of the one-component plasma of ions (OCP) derived from Monte Carlo simulations. The extension to multi-ionic mixtures is discussed within the framework of the linear mixing rule. These formulae provide a completely analytic, accurate description of the thermodynamic quantities of fully ionized electron-ion Coulomb plasmas, a useful tool for various applications from liquid state theory to dense stellar matter.

Multidimensional cosmologies allow for variations of fundamental physical constants over the course of cosmological evolution, and different versions of the theories predict different time dependences. In particular, such variations could manifest themselves as changes of the proton-to-electron mass ratio \mu=m_p/m_e over the period of ~ 10^{10} years since the moment of formation of high-redshift QSO spectra. Here we analyze a new, high-resolution spectrum of the z=2.81080 molecular hydrogen absorption system toward the quasar PKS 0528-250 to derive a new observational constraint to the time-averaged variation rate of the proton-to-electron mass ratio. We find |\dot{\mu} / \mu| < 1.5 \times 10^{-14}/year, which is much tighter than previously measured limits.

Recent calculations indicate that in the outer parts of neutron stars nuclei are rod-like or slab-like, rather than roughly spherical. We consider the elastic properties of these phases, and argue that they behave as liquid crystals, rather than rigid solids. We estimate elastic constants and discuss implications of our results for neutron star behavior.

The relationships between the effective surface (T_e) and internal temperatures of neutron stars (NSs) with and without accreted envelopes are calculated for T_e > 5\times10^4 K using new data on the equation of state and opacities in the outer NS layers. We examine various models of accreted layers (H, He, C, O shells produced by nuclear transformations in accreted matter). We employ new Opacity Library (OPAL) radiative opacities for H, He, and Fe. In the outermost NS layers, we implement the modern OPAL equation of state for Fe, and the Saumon-Chabrier equation of state for H and He. The updated thermal conductivities of degenerate electrons include the Debye-Waller factor for the electron-phonon scattering in solidified matter, while in liquid matter they include the contributions from electron-ion collisions (evaluated with non-Born corrections and with the ion structure factors in responsive electron background) and from the electron-electron collisions. For T_e < 10^{5.5} K, the electron conduction in non-degenerate layers of the envelope becomes important, reducing noticeably the temperature gradient. The accreted matter further decreases this gradient at T_e > 10^5 K. Even a small amount of accreted matter (with mass > 10^{-16} Msun) affects appreciably the NS cooling, leading to higher T_e at the neutrino cooling stage and to lower T_e at the subsequent photon stage.

Analytical approximations are constructed for binding energies of hydrogen atoms arbitrarily moving in magnetic fields 7\times10^{11} G to 2\times10^{13} G, typical of magnetic neutron stars.

Transport properties of degenerate relativistic electrons along quantizing magnetic fields in neutron star crusts are considered. A kinetic equation is derived for the spin polarization density matrix of electrons. Its solution does not depend on the choice of basic electron wave functions unlike previous solutions of the traditional kinetic equation for the distribution function. The density matrix formalism shows that one can always reach high accuracy with the traditional method by a proper choice of the basic functions. Electron Coulomb scattering on ions is considered in liquid matter, and on high-temperature phonons or on charged impurities in solid matter. In the solid regime, the Debye -- Waller reduction of phonon scattering can strongly enhance the longitudinal thermal or electric conductivity. An efficient numerical method is proposed for calculating the transport properties of electron gas at any magnetic field of practical interest.

Practical expressions are derived for a rapid and accurate evaluation of electric and thermal conductivities and thermopower of degenerate relativistic electrons along quantizing magnetic fields in outer neutron star crusts. The electron Coulomb scattering on ions is considered in liquid matter, and on high-temperature phonons or charged impurities in solid matter. A semi-quantitative treatment of low-temperature phonons is also proposed. The transport properties are expressed through the energy dependent effective electron relaxation time, which is calculated using the theoretical formalism of the previous work (astro-ph/9603133). Accurate fitting expressions for the relaxation time are obtained, and an efficient energy averaging procedure proposed. A Fortran code, which calculates the longitudinal transport properties of degenerate electrons in strong magnetic fields for any parameters of dense stellar matter of practical interest, is available from the authors upon request.

The strong magnetic fields (B ~ 10^{12} - 10^{13} G) characteristic of neutron stars make all the properties of an atom strongly dependent on the transverse component K_\perp of its generalized momentum. In particular, the photoionization process is modified substantially: (i) threshold energies are decreased as compared with those for an atom at rest, (ii) cross section values are changed significantly, and (iii) selection rules valid for atoms at rest are violated by the motion so that new photoionization channels become allowed. To calculate the photoionization cross sections, we, for the first time, employ exact numerical treatment of both initial and final atomic states. This enables us to take into account the quasi-bound (autoionizing) atomic states as well as coupling of different ionization channels. We extend the previous consideration, restricted to the so-called centered states corresponding to relatively small values of K_\perp, to arbitrary states of atomic motion. We fold the cross sections with the thermal distribution of atoms over K. For typical temperatures of neutron star atmospheres, the averaged cross sections differ substantially from those of atoms at rest. In particular, the photoionization edges are strongly broadened by the thermal motion of atoms; this "magnetic broadening" exceeds the usual Doppler broadening by orders of magnitude. The decentered states of the atoms give rise to the low-energy component of the photoionization cross section. This new component grows significantly with increasing temperature above 10^{5.5} K and decreasing density below 1 g/cm^3, i.e., for the conditions expected in atmospheres of middle-aged neutron stars.

Contemporary multidimensional cosmological theories predict different variations of fundamental physical constants in course of the cosmological evolution. On the basis of the QSO spectra analysis, we show that the fine-structure constant \alpha=e^2/(\hbar c) and the proton-to-electron mass ratio \mu=m_p/m_e reveal no statistically significant variation over the last 90% of the lifetime of the Universe. At the 2\sigma significance level, the following upper bounds are obtained for the epoch corresponding to the cosmological redshifts z ~ 3 (i.e., ~ 10 Gyr ago): |\Delta\alpha/\alpha| < 0.00016 and |\Delta\mu/\mu| < 0.00022. The corresponding upper limits to the time-average rates of the constant variations are |d\alpha/(\alpha dt)| < 1.6\times 10^{-14} yr^{-1} and |d\mu/(\mu dt)| < 2.2\times10^{-14} yr^{-1}. These limits serve as criteria for selection of those theoretical models which predict \alpha and \mu variation with the cosmological time. In addition, we test a possible anisotropy of the high-redshift fine splitting over the celestial sphere, which might reveal a non-equality of \alpha values in causally disconnected areas of the Universe.

The hydrogen plasma is studied at temperatures T ~ 10^4 - 10^6 K using the free energy minimization method. A simple analytic free energy model is proposed which is accurate at densities up to 1 g/cc and yields convergent internal partition function of atoms. The occupation probability formalism is modified for solving the ionization equilibrium problem. The ionization degree and equation of state are calculated and compared with the results of other models.

High magnetic fields in neutron stars, B ~ 10^{11} - 10^{13} G, substantially modify the properties of atoms and their interaction with radiation. In particular, the photoionization cross section becomes anisotropic and strongly polarization dependent. In a number of previous works based on the adiabatic approximation the conclusion was drawn that the transverse cross section vanishes for frequencies smaller than the electron cyclotron frequency. In other works (which employed a different form of the interaction operator) appreciable finite values were obtained. An adequate interpretation of the neutron star thermal-like radiation requires a resolution of this controversy. In this work the atomic wave functions for both discrete and continuum states are calculated by solving the coupled channel equations allowing the admixture between different Landau levels, which provides much higher accuracy than the adiabatic approximation. This enables to resolve the above contradiction in favour of the finite transverse cross sections. The non-adiabatic treatment of the continuum includes coupling between closed and open channels, which leads to the autoionization of quasi-bound energy levels associated with the electron cyclotron excitations and gives rise to Beutler - Fano resonances of the photoionization cross section. Autoionization widths of these quasi-bound levels are calculated and compared with the radiative widths. The results are important for investigations of the radiation emergent from the surface layers of neutron stars.