We formulate hydrodynamic equations for nonsuperfluid multicomponent magnetized charged relativistic mixtures, taking into account chemical reactions as well as viscosity, diffusion, thermodiffusion, and thermal conductivity effects. The resulting equations have a rather simple form and can be readily applied, e.g., for studying magnetothermal evolution of neutron stars. We also establish a link between our formalism and the results known in the literature, and express the phenomenological diffusion coefficients through momentum transfer rates which are calculated from microscopic theory.

We consider the thermal conductivity, shear viscosity, and momentum relaxation rates in the nucleon cores of the neutron stars. We study how the choice of the nuclear interaction and the model for three-body forces may affect these transport coefficients calculated within the Brueckner-Hartree-Fock many-body nuclear theory. We find that at relatively large densities the model dependence of the results is substantial. In addition we provide the analytical approximations which allow to incorporate our results in practical simulations.

We investigate a contribution of the Jeans smoothing to the minimal width of Ly$\alpha$ forest lines and discuss how the accounting for this additional broadening affects the inferred parameters of the intergalactic matter equation of state. We estimate a power-law index $\gamma$ of the equation of state, a temperature at the mean density $T_0$ and a hydrogen photoionization rate $\Gamma$ within 4 redshift bins. Furthermore, in each bin we obtain an upper limit on the scale-parameter $f_{\rm J}$, which sets the relation between the Jeans length and the characteristic physical size of the absorber clouds.

We present a new method to infer parameters of the temperature-density relation in the intergalactic medium in the post-reionization epoch at $z\sim 2-4$. This method is based on the analysis of the Ly$\alpha$ absorbers distribution over column densities and Doppler parameters by the model joint probability density function. This approach allows us to measure the power-law index $\gamma$ of the temperature-density relation and a certain combination of the temperature at the mean density $T_0$ and hydrogen photoionization rate $\Gamma$. To estimate $T_0$ and $\Gamma$ separately, we employ measurements of the Ly$\alpha$ forest effective opacity and the model gas probability density function. We show that $\gamma$ tends to be lower than 1.6 and reaches 1.3 at redshift $\sim3$. The inferred temperatures at the mean density are $\sim(2\pm0.5)\times10^4$ K in the studied redshift range. Both these estimates favour HeII reionization at $z\gtrsim3$. We find that the hydrogen photoionization rate is $\sim0.6\times10^{-12}$ s$^{-1}$, which is consistent with previous measurements.

The bright radio pulsar B1727$-$47 with a characteristic age of 80 kyr is among the first pulsars discovered 50 yr ago. Using its regular timing observations and interferometric positions at three epochs, we measured, for the first time, the pulsar proper motion of $151 \pm 19$ mas yr$^{-1}$. At the dispersion measure distance of $\gtrsim 2.7$ kpc, this would suggest a record transverse velocity of the pulsar $\gtrsim 1900$ km s$^{-1}$. However, a backward extrapolation of the pulsar track to its birth epoch points remarkably close to the center of the evolved nearby supernova remnant RCW 114, which suggests genuine association of the two objects. In this case, the pulsar is substantially closer ($\sim 0.6$ kpc) and younger ($\sim 50$ kyr), and its velocity ($\sim400$ km s$^{-1}$) is compatible with the observed pulsar velocity distribution. We also identified two new glitches of the pulsar. We discuss implications of our results on the pulsar and remnant properties.

Understanding signals from neutron stars requires knowledge about the transport inside the star. We review the transport properties and the underlying reaction rates of dense hadronic and quark matter in the crust and the core of neutron stars and point out open problems and future directions.

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

I consider the thermal conductivity and shear viscosity of leptons (electrons and muons) in the nucleon NS cores where protons are in the superconducting state. I restrict the consideration to the case of not too high temperatures $T\lesssim 0.35T_{\mathrm{c}p}$, where $T_{\mathrm{c}p}$ is the critical temperature of the proton pairing. In this case, lepton collisions with protons can be neglected. Charged lepton collision frequencies are mainly determined by the transverse plasmon exchange and are mediated by the character of the transverse plasma screening. In our previous works [Shternin \& Yakovlev, Phys. Rev. D {\bf 75} 103004 (2007); {\bf 78} 063006 (2008)] the superconducting proton contribution to the transverse screening was considered in the Pippard limit $\Delta \ll \hbar q v_{\mathrm{F}p}$, where $\Delta$ is the proton pairing gap, $v_{\mathrm{F}p}$ is the proton Fermi velocity, and $\hbar q$ is the typical transferred momentum in collisions. However, for large critical temperatures (large $\Delta$) and relatively small densities (small $q$) the Pippard limit may become invalid. In the present study I show that this is indeed the case and that the older calculations severely underestimated the screening in a certain range of the parameters appropriate to the neutron star cores. As a consequence, the kinetic coefficients at $T\ll T_{\mathrm{c}p}$ are found to be smaller than in previous calculations.

Now 50 years since the existence of the neutron star crust was proposed, we review the current understanding of the nuclear physics of the outer layers of accreting neutron stars. Nuclei produced during nuclear burning replace the nascent composition of the neutron star ocean and crust. Non-equilibrium nuclear reactions driven by compression alter the outer thermal structure and chemical composition, leaving observable imprints on astronomical phenomena. As observations of bursting neutron stars and cooling neutron stars have increased, the recent volume of astronomical data allows new insights into the microphysics of the neutron star interior and the possibility to test nuclear physics input in model calculations. Despite numerous advances in our understanding of neutron star interiors and observed neutron star phenomena, many challenges remain in the astrophysics theory of accreting neutron stars, the nuclear theory of neutron-rich nuclei, and the reach and precision of terrestrial nuclear physics experiments.

Quasar spectra provide a unique opportunity to investigate the intergalactic medium at high redshifts. The spectral analysis of the Ly$\alpha$ forest allows us to constrain the effective equation of state of the intergalactic medium and its changes during the evolution of the Universe. Based on the Voigt profile fitting of Ly$\alpha$ forest lines in 50 high-resolution QSO spectra obtained at Keck telescope we present new measurements of the power-law index $\gamma$ of temperature-density relation in IGM for six redshift bins in the range $z=1.6-3.7$. We find that the IGM state is close to isothermal ($\gamma\approx 1$) at $z\sim 3$ which may indicate that HeII reionization occurred at this redshift.

X-ray observations of transiently accreting neutron stars during quiescence provide information about the structure of neutron star crusts and the properties of dense matter. Interpretation of the observational data requires an understanding of the nuclear reactions that heat and cool the crust during accretion, and define its nonequilibrium composition. We identify here in detail the typical nuclear reaction sequences down to a depth in the inner crust where the mass density is 2E12 g/cm^3 using a full nuclear reaction network for a range of initial compositions. The reaction sequences differ substantially from previous work. We find a robust reduction of crust impurity at the transition to the inner crust regardless of initial composition, though shell effects can delay the formation of a pure crust somewhat to densities beyond 2E12 g/cm^3. This naturally explains the small inner crust impurity inferred from observations of a broad range of systems. The exception are initial compositions with A >= 102 nuclei, where the inner crust remains impure with an impurity parameter of Qimp~20 due to the N = 82 shell closure. In agreement with previous work we find that nuclear heating is relatively robust and independent of initial composition, while cooling via nuclear Urca cycles in the outer crust depends strongly on initial composition. This work forms a basis for future studies of the sensitivity of crust models to nuclear physics and provides profiles of composition for realistic crust models.

Spectroscopy of the Ly$\alpha$ forest in quasar spectra proved to be a useful tool for probing the intergalactic gas. We developed the automatic program for Voigt profile fitting of Ly$\alpha$ forest lines. We run this code on 9 high resolution ($\sim 50000$) quasars spectra with a high signal-to-noise ratio ($\sim 50 -100$) from the Keck telescope archive and obtained the sample of single well-fitted Ly$\alpha$ lines. Fitting the joint 2d distribution of column density and Doppler parameter from this sample by physically reasonable model we estimate a power law index $\gamma$ of the equation of state of the intergalactic medium in the redshift range $z\sim 2-3$. We found that our measurement is in an agreement with measurements by other groups obtained with Voigt profile fitting technique.

We report preliminary results of the analysis of the proper motion of the bright radio pulsar B1727$-$47. Using archival Parkes timing data, as well as original and archival ATCA interferometry observations, we, for the first time, constrain the pulsar proper motion at the level of 148$\pm$11 mas yr$^{-1}$. The backward extrapolation of the proper motion vector to the pulsar birth epoch points at the center of the Galactic supernova remnant RCW 114 suggesting the genuine association between the two objects. We discuss the implications of the association and argue that the distance to the system is less than $1$ kpc. This value is at least two times lower than the dispersion measure distance estimates. This suggests that the existing Galaxy electron density models are incomplete in the direction to the pulsar.

Thermal conductivity and shear viscosity of npe$\mu$ matter in non-superfluid neutron star cores are considered in the framework of Brueckner-Hartree-Fock many-body theory. We extend our previous work (Shternin, Baldo and Haensel, 2013) by analysing different nucleon-nucleon potentials and different three-body forces. We find that the use of different potentials leads up to one order of magnitude variations in the values of the nucleon contribution to transport coefficients. The nucleon contribution dominates the thermal conductivity, but for all considered models the shear viscosity is dominated by leptons.

We present the analysis of the archival Suzaku and Swift X-ray observations of the young $\gamma$-ray pulsar J1932+1916 field. The data revealed a point-like object at the $\gamma$-ray position of the pulsar and diffuse X-ray emission around it. Spectra of the point-like source and diffuse emission are well-described by absorbed power-law models with spectral parameters typical for pulsar plus pulsar wind nebula systems. Therefore we suggest that Suzaku and Swift detected the X-ray counterpart of PSR J1932+1916. Assuming this interpretation, we constrain the distance to the pulsar in the range of 2-6 kpc. We also suggest possible association of the pulsar with the nearby supernova remnant G54.4-0.3 and discuss its implications for the pulsar proper motion, age and distance.

Calvera is an unusual isolated neutron star with pure thermal X-ray spectrum typical for central compact objects in supernova remnants. On the other hand, its rotation period and spin-down rate are typical for ordinary rotation-powered pulsars. It was discovered and studied in X-rays and not yet detected in other spectral domains. We present deep optical imaging of the Calvera field obtained with the Gran Telescopio Canarias in $g'$ and $i'$ bands. Within $\approx 1^{\prime\prime}$ vicinity of Calvera, we detected two point-like objects invisible at previous shallow observations. However, accurate astrometry showed that none of them can be identified with the pulsar. We put new upper limits on its optical brightness of $g' > 27.87$ and $i' > 26.84$. We also reanalyzed all available archival X-ray data on Calvera. Comparison of the Calvera thermal emission parameters and upper limits on optical and non-thermal X-ray emission with respective data on rotation-powered pulsars shows that Calvera might belong to the class of ordinary middle-aged pulsars, if we assume that its distance is in the range of $1.5-5$ kpc.

We report the results of XMM-Newton observations of the Galactic mixed-morphology supernova remnant G350.0$-$2.0. Diffuse thermal X-ray emission fills the north-western part of the remnant surrounded by radio shell-like structures. We did not detect any X-ray counterpart of the latter structures, but found several bright blobs within the diffuse emission. The X-ray spectrum of the most part of the remnant can be described by a collisionally-ionized plasma model VAPEC with solar abundances and a temperature of $\approx 0.8$ keV. The solar abundances of plasma indicate that the X-ray emission comes from the shocked interstellar material. The overabundance of Fe was found in some of the bright blobs. We also analysed the brightest point-like X-ray source 1RXS J172653.4$-$382157 projected on the extended emission. Its spectrum is well described by the two-temperature optically thin thermal plasma model MEKAL typical for cataclysmic variable stars. The cataclysmic variable source nature is supported by the presence of a faint ($g\approx21$) optical source with non-stellar spectral energy distribution at the X-ray position of 1RXS J172653.4$-$382157. It was detected with the XMM-Newton optical/UV monitor in the $U$ filter and was also found in the archival H$\alpha$ and optical/near-infrared broadband sky survey images. On the other hand, the X-ray spectrum is also described by the power law plus thermal component model typical for a rotation powered pulsar. Therefore, the pulsar interpretation of the source cannot be excluded. For this source, we derived the upper limit for the pulsed fraction of 27 per cent.

We analysed Chandra observations of the bright Fermi pulsar J0633+0632 and found evidence of an absorption feature in its spectrum at $804^{+42}_{-26}$ eV (the errors here and below are at 90% confidence) with equivalent width of $63^{+47}_{-36}$ eV. In addition, we analysed in detail the X-ray spectral continuum taking into account correlations between the interstellar absorption and the distance to the source. We confirm early findings by Ray et al. (2011) that the spectrum contains non-thermal and thermal components. The latter is equally well described by the blackbody and magnetised atmosphere models and can be attributed to the emission from the bulk of the stellar surface in both cases. The distance to the pulsar is constrained in a range of 1--4 kpc from the spectral fits. We infer the blackbody surface temperature of $108^{+22}_{-14}$ eV, while for the atmosphere model, the temperature, as seen by a distant observer, is $53^{+12}_{-7}$ eV. In the latter case J0633+0632 is one of the coldest middle-aged isolated neutron stars with measured temperatures. Finally, it powers an extended pulsar wind nebula whose shape suggests a high pulsar proper motion. Looking backwards the direction of the presumed proper motion we found a likely birthplace of the pulsar -- the Rosette nebula, a 50-Myr-old active star-forming region located at about 1$.\!\!^\circ$5 from the pulsar. If true, this constrains the distance to the pulsar in the range of 1.2--1.8 kpc.

PSR J1357$-$6429 is a young and energetic radio pulsar detected in X-rays and $\gamma$-rays. It powers a compact pulsar wind nebula with a jet visible in X-rays and a large scale plerion detected in X-ray and TeV ranges. Previous multiwavelength studies suggested that the pulsar has a significant proper motion of about 180 mas yr$^{-1}$ implying an extremely high transverse velocity of about 2000 km s$^{-1}$. In order to verify that, we performed radio-interferometric observations of PSR J1357$-$6429 with the the Australia Telescope Compact Array (ATCA) in the 2.1 GHz band. We detected the pulsar with a mean flux density of $212\pm5$ $\mu$Jy and obtained the most accurate pulsar position, RA = 13:57:02.525(14) and Dec = $-$64:29:29.89(15). Using the new and archival ATCA data, we did not find any proper motion and estimated its 90 per cent upper limit $\mu < 106$ mas yr$^{-1}$. The pulsar shows a highly polarised single pulse, as it was earlier observed at 1.4 GHz. Spectral analysis revealed a shallow spectral index $\alpha_{\nu}$ = $0.5 \pm 0.1$. Based on our new radio position of the pulsar, we disclaim its optical counterpart candidate reported before.

PSR J2021+3651 is a 17 kyr old rotation powered pulsar detected in the radio, X-rays, and $\gamma$-rays. It powers a torus-like pulsar wind nebula with jets, dubbed the Dragonfly, which is very similar to that of the Vela pulsar. The Dragonfly is likely associated with the extended TeV source VER J2019+368 and extended radio emission. We conducted first deep optical observations with the GTC in the Sloan $r'$ band to search for optical counterparts of the pulsar and its nebula. No counterparts were detected down to $r'\gtrsim27.2$ and $\gtrsim24.8$ for the point-like pulsar and the compact X-ray nebula, respectively. We also reanalyzed Chandra archival X-ray data taking into account an interstellar extinction--distance relation, constructed by us for the Dragonfly line of sight using the red-clump stars as standard candles. This allowed us to constrain the distance to the pulsar, $D=1.8^{+1.7}_{-1.4}$ kpc at 90% confidence. It is much smaller than the dispersion measure distance of $\sim$12 kpc but compatible with a $\gamma$-ray "pseudo-distance" of 1 kpc. Based on that and the optical upper limits, we conclude that PSR J2021+3651, similar to the Vela pulsar, is a very inefficient nonthermal emitter in the optical and X-rays, while its $\gamma$-ray efficiency is consistent with an average efficiency for $\gamma$-pulsars of similar age. Our optical flux upper limit for the pulsar is consistent with the long-wavelength extrapolation of its X-ray spectrum while the nebula flux upper limit does not constrain the respective extrapolation.

We consider models of cooling neutron stars with nucleon cores which possess moderately strong triplet-state superfluidity of neutrons. When the internal temperature drops below the maximum of the critical temperature over the core, $T_{\rm C}$, this superfluidity sets in. It produces a neutrino outburst due to Cooper pairing of neutrons which greatly accelerates the cooling. We show that the cooling of the star with internal temperature $T$ within $0.6\,T_{\rm C} \lesssim T \leq T_{\rm C}$ is described by analytic self-similar relations. A measurement of the effective surface temperature of the star and its decline, supplemented by assumptions on star's mass, radius and composition of heat-blanketing envelope, allows one to construct a family of cooling models parametrized by the value of $T_{\rm C}$. Each model reconstructs cooling history of the star including its neutrino emission level before neutron superfluidity onset and the intensity of Cooper pairing neutrinos. The results are applied to interpret the observations of the neutron star in the Cassiopeia A supernova remnant.

We present results of the spectral analysis of the X-ray emission from the middle-aged Fermi pulsar J1741$-$2054 using all Chandra archival data collected in 2010 and 2013. We confirm early findings by Romani et al. 2010 that the pulsar spectrum contains a thermal emission component. The component is best described by the blackbody model with the temperature $\approx$ 60 eV and the emitting area radius $\approx$ 17 $D_{\rm kpc}$ km. The thermal emission likely originates from the entire surface of the cooling neutron star if the distance to the pulsar is $\approx$ 0.8 kpc. The latter is supported by a large absorbing column density inferred from the X-ray fit and empirical optical extinction--distance relations along the pulsar line of sight. The neutron star surface temperature and characteristic age make it similar to the well studied middle-aged pulsar B1055$-$52. Like the latter PSR J1741$-$2054 is hotter than the standard cooling scenario predicts.

A middle-aged radio-quiet pulsar J0357+3205 was discovered in gamma-rays with $Fermi$ and later in X-rays with $Chandra$ and $XMM$-$Newton$ observatories. It produces an unusual thermally-emitting pulsar wind nebula observed in X-rays. Deep optical observations were obtained to search for the pulsar optical counterpart and its nebula using the Gran Telescopio Canarias (GTC). The direct imaging mode in the Sloan $g'$ band was used. Archival X-ray data were reanalysed and compared with the optical data. No pulsar optical counterpart was detected down to $g'\geq~28_{\cdotp}^{\text{m}}1$. No pulsar nebula was either identified in the optical. We confirm early results that the X-ray spectrum of the pulsar consists of a nonthermal power-law component of the pulsar magnetospheric origin dominating at high energies and a soft thermal component from the neutron star surface. Using magnetised partially ionised hydrogen atmosphere models in X-ray spectral fits we found that the thermal component can come from entire surface of the cooling neutron star with a temperature of 36$^{+9}_{-6}$ eV, making it one of the coldest among cooling neutron stars known. The surface temperature agrees with the standard neutron star cooling scenario. The optical upper limit does not put any additional constraints on the thermal component, however it implies a strong spectral break for the nonthermal component between the optical and X-rays as is observed in other middle-aged pulsars. The thermal emission from the entire surface of the neutron star likely dominates over the nonthermal emission in the UV range. Observations of the PSR J0357+3205 in this range are promising to put more stringent constraints on its thermal properties.

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

The thermal evolution of young neutron stars (NSs) reflects the neutrino emission properties of their cores. Heinke et al. (2010) measured a 3.6+/-0.6% decay in the surface temperature of the Cassiopeia A (Cas A) NS between 2000 and 2009, using archival data from the Chandra X-ray Observatory ACIS-S detector in Graded mode. Page et al. (2011) and Shternin et al. (2011) attributed this decay to enhanced neutrino emission from a superfluid neutron transition in the core. Here we test this decline, combining analysis of the Cas A NS using all Chandra X-ray detectors and modes (HRC-S, HRC-I, ACIS-I, ACIS-S in Faint mode, and ACIS-S in Graded mode) and adding a 2012 May ACIS-S Graded mode observation, using the most current calibrations (CALDB 4.5.5.1). We measure the temperature changes from each detector separately and test for systematic effects due to the nearby filaments of the supernova remnant. We find a 0.92%-2.0% decay over 10 years in the effective temperature, inferred from HRC-S data, depending on the choice of source and background extraction regions, with a best-fit decay of 1.0+/-0.7%. In comparison, the ACIS-S Graded data indicate a temperature decay of 3.1%-5.0% over 10 years, with a best-fit decay of 3.5+/-0.4%. Shallower observations using the other detectors yield temperature decays of 2.6+/-1.9% (ACIS-I), 2.1+/-1.0% (HRC-I), and 2.1+/-1.9% (ACIS-S Faint mode) over 10 years. Our best estimate indicates a decline of 2.9+/-0.9 (stat) +1.6/-0.3 (sys) % over 10 years. The complexity of the bright and varying supernova remnant background makes a definitive interpretation of archival Cas A Chandra observations difficult. A temperature decline of 1-3.5% over 10 years would indicate extraordinarily fast cooling of the NS that can be regulated by superfluidity of nucleons in the stellar core.

We consider the astrophysical reaction rates for radiative neutron capture reactions ($n,\gamma$) in the crust of a neutron star. The presence of degenerate neutrons at high densities (mainly in the inner crust) can drastically affect the reaction rates. Standard rates assuming a Maxwell-Boltzmann distribution for neutrons can underestimate the rates by several orders of magnitude. We derive simple analytical expressions for reaction rates at a variety of conditions with account for neutron degeneracy. We also discuss the plasma effects on the outgoing radiative transition channel in neutron radiative capture reactions and show that these effects can also increase the reaction rates by a few orders of magnitude. In addition, using detailed balance, we analyze the effects of neutron degeneracy and plasma physics on reverse ($\gamma,n$) photodisintegration. We discuss the dependence of the reaction rates on temperature and neutron chemical potential and outline the efficiency of these reactions in the neutron star crust.

According to recent results of Ho & Heinke (2009) and Heinke & Ho (2010), the Cassiopeia A supernova remnant contains a young neutron star which has carbon atmosphere and shows noticeable decline of the effective surface temperature. We report a new (November 2010) Chandra observation which confirms the previously reported decline rate. The decline is naturally explained if neutrons have recently become superfluid (in triplet-state) in the NS core, producing a splash of neutrino emission due to Cooper pair formation (CPF) process that currently accelerates the cooling. This scenario puts stringent constraints on poorly known properties of NS cores: on density dependence of the temperature $T_\mathrm{cn}(\rho)$ for the onset of neutron superfluidity [$T_\mathrm{cn}(\rho)$ should have a wide peak with maximum $\approx (7-9)\times 10^8$ K], on the reduction factor $q$ of CPF process by collective effects in superfluid matter ($q > 0.4$), and on the intensity of neutrino emission before the onset of neutron superfluidity (30--100 times weaker than the standard modified Urca process). This is serious evidence for nucleon superfluidity in NS cores that comes from observations of cooling NSs.

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.

We study plasma effects on radiative transitions (e.g., decay of excited states of atoms or atomic nuclei) in a dense plasma at the transition frequencies $\omega \lesssim \omega_p$ (where $\omega_p$ is the electron plasma frequency). The decay goes through four channels -- the emission of real transverse and longitudinal plasmons as well as the emission of virtual transverse and longitudinal plasmons with subsequent absorption of such plasmons by the plasma. The emission of real plasmons dies out at $\omega \leq \omega_p$, but the processes with virtual plasmons strongly enhance the radiative decay. Applications of these results to radiative processes in white dwarf cores and neutron star envelopes are discussed.

We calculate the shear viscosity $\eta = \eta_{e\mu}+\eta_{n}$ in a neutron star core composed of nucleons, electrons and muons ($\eta_{e\mu}$ being the electron-muon viscosity, mediated by collisions of electrons and muons with charged particles, and $\eta_{n}$ the neutron viscosity, mediated by neutron-neutron and neutron-proton collisions). Deriving $\eta_{e\mu}$, we take into account the Landau damping in collisions of electrons and muons with charged particles via the exchange of transverse plasmons. It lowers $\eta_{e\mu}$ and leads to the non-standard temperature behavior $\eta_{e\mu}\propto T^{-5/3}$. The viscosity $\eta_{n}$ is calculated taking into account that in-medium effects modify nucleon effective masses in dense matter. Both viscosities, $\eta_{e\mu}$ and $\eta_{n}$, can be important, and both are calculated including the effects of proton superfluidity. They are presented in the form valid for any equation of state of nucleon dense matter. We analyze the density and temperature dependence of $\eta$ for different equations of state in neutron star cores, and compare $\eta$ with the bulk viscosity in the core and with the shear viscosity in the crust.

We calculate the partial electron shear viscosity $\eta_{ee}$ limited by electron-electron collisions in a strongly degenerate electron gas taking into account the Landau damping of transverse plasmons. The Landau damping strongly suppresses $\eta_{ee}$ in the domain of ultrarelativistic degenerate electrons and modifies its %asymptotic temperature behavior. The efficiency of the electron shear viscosity in the cores of white dwarfs and envelopes of neutron stars is analyzed.

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 calculate the thermal conductivity of electrons and muons kappa_{e-mu} produced owing to electromagnetic interactions of charged particles in neutron star cores and show that these interactions are dominated by the exchange of transverse plasmons (via the Landau damping of these plasmons in nonsuperconducting matter and via a specific plasma screening in the presence of proton superconductivity). For normal protons, the Landau damping strongly reduces kappa_{e-mu} and makes it temperature independent. Proton superconductivity suppresses the reduction and restores the Fermi-liquid behavior kappa_{e-mu} ~ 1/T. Comparing with the thermal conductivity of neutrons kappa_n, we obtain kappa_{e-mu}> kappa_n for T>2 GK in normal matter and for any T in superconducting matter with proton critical temperatures T_c>3e9 K. The results are described by simple analytic formulae.

We calculate the thermal conductivity of electrons produced by electron-electron Coulomb scattering in a strongly degenerate electron gas taking into account the Landau damping of transverse plasmons. The Landau damping strongly reduces this conductivity in the domain of ultrarelativistic electrons at temperatures below the electron plasma temperature. In the inner crust of a neutron star at temperatures T < 1e7 K this thermal conductivity completely dominates over the electron conductivity due to electron-ion (electron-phonon) scattering and becomes competitive with the the electron conductivity due to scattering of electrons by impurity ions.

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 calculate the polarization of radiation from thick accretion disks with vertically averaged global magnetic field. The polarization arises as a result of the radiation scattering by free electrons in magnetized plasma of a disk. We consider as a basic effect the Faraday rotation of polarization plane along the photon propagation in a magnetized disk. The various models of optically thick accretion disk with a vertically averaged magnetic field are considered. The main goal of this paper is to obtain simple asymptotic formulae for the polarization of radiation when the Faraday rotation angle >> 1 at the Thomson optical length ~1. The results of our calculation allows us to estimate the magnetic field magnitude near the marginally stable orbit region of a black hole via the data of polarimetric observations including the expected future X-ray polarimetric observations. The wavelength dependence of polarization is strongly dependent on various models of an accretion disks and thus allows to choose a real model from the polarization data.

The numerical solution of the Milne problem for semi-infinite plane-parallel magnetized electron atmosphere is obtained. It is assumed that magnetic field is directed along the normal to the atmosphere. The angular dependence, the polarization degree and positional angle of outgoing radiation are presented in the tables for various values of the Faraday rotation parameter and the degree of absorption q=0, 0.2 and 0.4. We assume that magnetic field B~1000000G when one can neglect the effects of circular dichroism and take into account only the Faraday rotation effect.