The results of hybrid simulation of low-beta supercritical quasi-parallel shocks in high metallicity plasma are presented. The structure of upstream and downstream turbulence is addressed and velocities of the corresponding scattering centers are derived. It is shown that independently of their chemical composition the shocks experience self-reformation process. However, the period of self-reformation as well as the wave spectrum is greatly affected by the presence of substantial admixture of weakly charged heavy ions. Also the downstream magnetic field amplification is stronger for the high metallicity case.

We consider acceleration of leptons up to GeV-TeV energies in the bow shock wind nebula of PSR J0437-4715 and their subsequent diffusion through the interstellar magnetic fields. The leptons accelerated at the pulsar wind termination surface are injected into re-acceleration in colliding shock flows. Modelled spectra of synchrotron emission from the accelerated electrons and positrons are consistent with the far-ultraviolet and X-ray observations of the nebula carried out with the Hubble Space Telescope and Chandra X-ray Observatory. These observations are employed to constrain the absolute fluxes of relativistic leptons, which are escaping from the nebula and eventually reaching the Solar System after energy-dependent diffusion through the local interstellar medium accompanied by synchrotron and Compton losses. It is shown that accelerated leptons from the nebula of PSR J0437-4715 can be responsible both for the enhancement of the positron fraction above a few GeV detected by PAMELA and AMS-02 spectrometers and for the TeV range lepton fluxes observed with H.E.S.S., VERITAS, Fermi, CALET, and DAMPE.

Galaxy clusters grow by gas accretion, mostly from mergers of substructures, which release powerful shock waves into cosmic plasmas and convert a fraction of kinetic energy into thermal energy, amplification of magnetic fields and into the acceleration of energetic particles. The modeling of the radio signature of cosmic shocks, combined with the lack of detected gamma-rays from cosmic ray (CR) protons, poses challenges to our understanding of how cosmic rays get accelerated and stored in the intracluster medium. Here we review the injection of CRs by cosmic shocks of different strengths, combining the detailed "microscopic" view of collisionless processes governing the creation of non-thermal distributions of electrons and protons in cluster shocks (based on analytic theory and particle-in-cell simulations), with the "macroscopic" view of the large-scale distribution of cosmic rays, suggested by modern cosmological simulations. We discuss time dependent non-linear kinetic models of particle acceleration by multiple internal shocks with large scale compressible motions of plasma. The models produce soft CR spectra containing a noticeable energy density in the super-thermal protons of energies below a few GeV which are difficult to constrain by Fermi observations. We consider the effect of plasma composition on CR injection and super-thermal particle population in the hot intracluster matter which can be constrained by fine high resolution X-ray spectroscopy of Fe ions.

A second-order accurate divergence-conserving hybrid particle-in-cell code Maximus has been developed for microscopic modeling of collisionless plasmas. The main specifics of the code include a constrained transport algorithm for exact conservation of magnetic field divergence, a Boris-type particle pusher, a weighted particle momentum deposit on the cells of the 3D spatial grid, an ability to model multispecies plasmas and an adaptive time step. The code is efficiently parallelized for running on supercomputers by means of the message passing interface (MPI) technology; an analysis of parallelization efficiency and overall resource intensity is presented. A Maximus simulation of the shocked flow in the Solar wind is shown to agree well with the observations of the Ion Release Module (IRM) aboard the Active Magnetospheric Particle Tracer Explorers interplanetary mission.

Pulsars with high spin-down power produce relativistic winds radiating a fraction of the power in the range from radio to gamma-rays in the pulsar wind nebulae (PWNe). The rest of the power is dissipated in the interactions of the PWNe with the interstellar medium (ISM). Some of the PWNe are moving relative to the ISM with supersonic speeds producing bow shocks. In this case, the ultrarelativistic particles accelerated at the termination surface of the pulsar wind may be reaccelerated in the converging flow system formed by the outflow from the wind termination shock and the inflow from the bow shock. An outcome of this reacceleration is the creation of particle distributions with hard spectra, such as required to explain the observed synchrotron spectra with photon indices Gamma <~ 1.5. The presence of this hard component is specific to PWNe with bow shocks (BSPWNe). The accelerated particles may end up containing a substantial fraction of the shock ram pressure. For typical ISM and pulsar parameters, the positrons released by these systems are numerous enough to contribute a substantial fraction of the positrons detected as galactic cosmic ray particles above few tens and up to several hundred GeV. The escape of ultrarelativistic particles from a BSPWN and its appearance in the far-UV and X-ray bands is determined by the directions of the interstellar magnetic field, the velocity of the astrosphere and the pulsar rotation axis. In this respect we review the observed appearance and multiwavelength spectra of three different types of BSPWNe: PSR J0437-4715, the Guitar and Lighthouse nebulae, and Vela-like objects. We argue that high resolution imaging of such objects provides unique information on pulsar winds and on the ISM. We discuss the interpretation of imaging observations and estimate the BSPWN contribution to the positron flux observed at the Earth. (ABRIDGED)

We have studied the fascinating dynamics of the nearby Vela pulsar's nebula in a campaign comprising eleven 40ks observations with Chandra X-ray Observatory (CXO). The deepest yet images revealed the shape, structure, and motion of the 2-arcminute-long pulsar jet. We find that the jet's shape and dynamics are remarkably consistent with that of a steadily turning helix projected on the sky. We discuss possible implications of our results, including free precession of the neutron star and MHD instability scenarios.

Recent X-ray observations of the supernova remnant IC443 interacting with molecular clouds have shown the presence of a new population of hard X-ray sources related to the remnant itself, which has been interpreted in terms of fast ejecta fragment propagating inside the dense environment. Prompted by these studies, we have obtained a deep {\sl XMM-Newton} observation of the supernova remnant (SNR) Kes 69, which also shows signs of shock-cloud interaction. We report on the detection of 18 hard X-ray sources in the field of Kes 69, a significant excess of the expected galactic source population in the field, spatially correlated with CO emission from the cloud in the remnant environment. The spectra of 3 of the 18 sources can be described as hard power laws with photon index <2 plus line emission associated to K-shell transitions. We discuss the two most promising scenarios for the interpretation of the sources, namely fast ejecta fragments (as in IC443) and cataclysmic variables. While most of the observational evidences are consistent with the former interpretation, we cannot rule out the latter.

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

We discuss pycnonuclear burning of highly exotic atomic nuclei in deep crusts of neutron stars, at densities up to 1e13 g/cc. As an application, we consider pycnonuclear burning of matter accreted on a neutron star in a soft X-ray transient (SXT, a compact binary containing a neutron star and a low-mass companion). The energy released in this burning, while the matter sinks into the stellar crust under the weight of newly accreted material, is sufficient to warm up the star and initiate neutrino emission in its core. The surface thermal radiation of the star in quiescent states becomes dependent of poorly known equation of state (EOS) of supranuclear matter in the stellar core, which gives a method to explore this EOS. Four qualitatively different model EOSs are tested against observations of SXTs. They imply different levels of the enhancement of neutrino emission in massive neutron stars by (1) the direct Urca process in nucleon/hyperon matter; (2) pion condensates; (3) kaon condensates; (4) Cooper pairing of neutrons in nucleon matter with the forbidden direct Urca process. A low level of the thermal quiescent emission of two SXTs, SAX J1808.4-3658 and Cen X-4, contradicts model (4). Observations of SXTs test the same physics of dense matter as observations of thermal radiation from cooling isolated neutron stars, but the data on SXTs are currently more conclusive.

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.

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 thermal states of transiently accreting neutron stars (with mean accretion rates $\dot{M} \sim 10^{-14}-10^{-9}$ M$_\odot$ yr$^{-1}$) determined by the deep crustal heating of accreted matter sinking into stellar interiors. We formalize a direct correspondence of this problem to the problem of cooling neutron stars. Using a simple toy model we analyze the most important factors which affect the thermal states of accreting stars: a strong superfluidity in the cores of low-mass stars and a fast neutrino emission (in nucleon, pion-condensed, kaon-condensed, or quark phases of dense matter) in the cores of high-mass stars. We briefly compare the results with the observations of soft X-ray transients in quiescence. If the upper limit on the quiescent thermal luminosity of the neutron star in SAX J1808.4-3658 (Campana et al. 2002) is associated with the deep crustal heating, it favors the model of nucleon neutron-star cores with switched-on direct Urca process.

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 outline the phenomenon of deep crustal heating in transiently accreting neutron stars. It is produced by nuclear transformations (mostly, by pycnonuclear reactions) in accreted matter while this matter sinks to densities rho > 10^{10} g/cc under the weight of freshly accreted material. We consider then thermal states of transiently accreting neutron stars (with mean mass accretion rates \dot{M}=(10^{-14}-10^{-9}) M_\odot/yr) determined by deep crustal heating. In a simplified fashion we study how the thermal flux emergent from such stars depends on the properties of superdense matter in stellar interiors. We analyze the most important regulators of the thermal flux: strong superfluidity in the cores of low-mass stars and fast neutrino emission (in nucleon, pion-condensed, kaon-condensed, or quark phases of dense matter) in the cores of high-mass stars. We compare the results with observations of soft X-ray transients in quiescent states.

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

Bulk viscosity of neutron star cores containing hyperons is studied taking into account non-equilibrium weak process $n+n \rightleftharpoons p+\Sigma^-$. Rapid growth of the bulk viscosity within the neutron star core associated with switching on new reactions (modified Urca process, direct Urca process, hyperon reactions) is analyzed. The suppression of the bulk viscosity by superfluidity of baryons is considered and found out to be very important.

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

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

We study the heat capacity and neutrino emission reactions (direct and modified Urca processes, nucleon-nucleon bremsstrahlung, Cooper pairing of nucleons) in matter of supranuclear density of the neutron star cores with superfluid neutrons and protons. Various superfluidity types are analysed (singlet-state pairing and two types of triplet-state pairing, without and with nodes of the gap at a nucleon Fermi surface). The results are used for cooling simulations of isolated neutron stars. Both, the standard cooling and the cooling enhanced by the direct Urca process, are strongly affected by nucleon superfluidity. Comparison of cooling theory of isolated neutron stars with observations of their thermal radiation may give stringent constraints on the critical temperatures of the neutron and proton superfluidities in the neutron star cores.

Cooling of neutron stars (NSs) with superfluid cores is simulated taking into account neutrino emission produced by Cooper pairing of nucleons. The critical temperatures of neutron and proton superfluidities, $T_{cn}$ and $T_{cp}$, are assumed to be constant over the NS core, and treated as free parameters. They are constrained using the surface temperatures $T_s$ of isolated NSs (RX J0822-43, PSR 1055-52, 1E 1207-52, Vela, Geminga, PSR 0656+14, RX J0002+62), obtained by interpretation of observed thermal radiation either with black body spectrum or with hydrogen atmosphere models.

The neutrino energy emission rate due to formation of Cooper pairs of neutrons and protons in the superfluid cores of neutron stars is studied. The cases of singlet-state pairing with isotropic superfluid gap and triplet-state pairing with anisotropic gap are analysed. The neutrino emission due to singlet-state pairing of protons is found to be greatly suppressed with respect to the cases of singlet- and triplet-state pairings of neutrons. The neutrino emission due to pairing of neutrons is shown to be very important in the superfluid neutron-star cores with the standard neutrino luminosity and with the luminosity enhanced by the direct Urca process. It can greatly accelerate both, standard and enhanced, cooling of neutron stars with superfluid cores. This enables one to interpret the data on surface temperatures of six neutron stars, obtained by fitting the observed spectra with the hydrogen atmosphere models, by the standard cooling with moderate nucleon superfluidity.

Calculations are performed of the cooling of neutron stars with standard and enhanced neutrino energy losses in the presence of neutron and proton superfluidities in the stellar cores. The effects of superfluidity on the heat capacity and the neutrino luminosity produced by the direct and modified Urca processes and by the neutrino bremsstrahlung emission in nucleon-nucleon collisions are taken into account. The constraints are analyzed on the critical temperatures $\tn$ and $\tp$ of the transition of neutrons and protons to the superfluid state, which can be obtained from a comparison of observational data on thermal radiation from neutron stars with theoretical cooling curves for the standard and enhanced neutrino energy losses. Possible ranges of $\tn$ and $\tp$ for the pulsar Geminga are discussed.