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.

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)

Pulsars traveling at supersonic speeds are often accompanied by cometary bow shocks seen in Halpha. We report on the first detection of a pulsar bow shock in the far-ultraviolet (FUV). We detected it in FUV images of the nearest millisecond pulsar J0437-4715 obtained with the Hubble Space Telescope. The images reveal a bow-like structure positionally coincident with part of the previously detected Halpha bow shock, with an apex at 10" ahead of the moving pulsar. Its FUV luminosity, L(1250-2000 A) ~ 5x10^28 erg/s, exceeds the Halpha luminosity from the same area by a factor of 10. The FUV emission could be produced by the shocked ISM matter or, less likely, by relativistic pulsar wind electrons confined by strong magnetic field fluctuations in the bow shock. In addition, in the FUV images we found a puzzling extended (~3" in size) structure overlapping with the limb of the bow shock. If related to the bow shock, it could be produced by an inhomogeneity in the ambient medium or an instability in the bow shock. We also report on a previously undetected X-ray emission extending for about 5" ahead of the pulsar, possibly a pulsar wind nebula created by shocked pulsar wind, with a luminosity L(0.5-8 keV) ~ 3x10^28 erg/s.

In this review we discuss some observational aspects and theoretical models of astrophysical collisionless shocks in partly ionized plasma with the presence of non-thermal components. A specific feature of fast strong collisionless shocks is their ability to accelerate energetic particles that can modify the shock upstream flow and form the shock precursors. We discuss the effects of energetic particle acceleration and associated magnetic field amplification and decay in the extended shock precursors on the line and continuum multi-wavelength emission spectra of the shocks. Both Balmer-type and radiative astrophysical shocks are discussed in connection to supernova remnants interacting with partially neutral clouds. Quantitative models described in the review predict a number of observable line-like emission features that can be used to reveal the physical state of the matter in the shock precursors and the character of nonthermal processes in the shocks. Implications of recent progress of gamma-ray observations of supernova remnants in molecular clouds are highlighted.

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 nature of the extended hard X-ray source XMMU J061804.3+222732 and its surroundings is investigated using XMM-Newton, Chandra, and Spitzer observations. This source is located in an interaction region of the IC 443 supernova remnant with a neighboring molecular cloud. The X-ray emission consists of a number of bright clumps embedded in an extended structured non-thermal X-ray nebula larger than 30" in size. Some clumps show evidence for line emission at ~1.9 keV and ~3.7 keV at the 99% confidence level. Large-scale diffuse radio emission of IC 443 passes over the source region, with an enhancement near the source. An IR source of about 14" x 7" size is prominent in the 24 um, 70 um, and 2.2 um bands, adjacent to a putative Si K-shell X-ray line emission region. The observed IR/X-ray morphology and spectra are consistent with those expected for J/C-type shocks of different velocities driven by fragmented supernova ejecta colliding with the dense medium of a molecular cloud. The IR emission of the source detected by Spitzer can be attributed to both continuum emission from an HII region created by the ejecta fragment and line emission excited by shocks. This source region in IC 443 may be an example of a rather numerous population of hard X-ray/IR sources created by supernova explosions in the dense environment of star-forming regions. Alternative Galactic and extragalactic interpretations of the observed source are also discussed.

We found a very likely counterpart to the recently discovered hard X-ray source IGR J2018+4043 in the multi-wavelength observations of the source field. The source, originally discovered in the 20-40 keV band, is now confidently detected also in the 40-80 keV band, with a flux of (1.4 +/- 0.4) x 10(-11) erg cm(-2) s(-1). A 5 ks Swift observation of the IGR J2018+4043 field revealed a hard point-like source with the observed 0.5-10 keV flux of 3.4(+0.7)(-0.8) x 10(-12) erg cm(-2) s(-1) (90% confidence level) at alpha = 20h18m38.55s, delta = +40d41m00.4s (with a 4.2" uncertainty). The combined Swift-INTEGRAL spectrum can be described by an absorbed power-law model with photon index gamma = 1.3 +/- 0.2 and N_H = 6.1(+3.2)(-2.2) x 10(22) cm(-2). In archival optical and infrared data we found a slightly extended and highly absorbed object at the Swift source position. There is also an extended VLA 1.4 GHz source peaked at a beam-width distance from the optical and X-ray positions. The observed morphology and multiwavelength spectra of IGR J2018+4043 are consistent with those expected for an obscured accreting object, i.e. an AGN or a Galactic X-ray binary. The identification suggests possible connection of IGR J2018+4043 to the bright gamma-ray source GEV J2020+4023 (3EG J2020+4017) detected by COS B and CGRO EGRET in the gamma-Cygni SNR field.

We report the detection of hard X-ray emission from the field of the star-forming region NGC 6334 with the the International Gamma-Ray Astrophysics Laboratory INTEGRAL. The JEM-X monitor and ISGRI imager aboard INTEGRAL and Chandra ACIS imager were used to construct 3-80 keV images and spectra of NGC 6334. The 3-10 keV and 10-35 keV images made with JEM-X show a complex structure of extended emission from NGC 6334. The ISGRI source detected in the energy ranges 20-40 keV and 40-80 keV coincides with the NGC 6334 ridge. The 20-60 keV flux from the source is (1.8+-0.37)*10(-11) erg cm(-2) s(-1). Spectral analysis of the source revealed a hard power-law component with a photon index about 1. The observed X-ray fluxes are in agreement with extrapolations of X-ray imaging observations of NGC 6334 by Chandra ACIS and ASCA GIS. The X-ray data are consistent with two very different physical models. A probable scenario is emission from a heavily absorbed, compact and hard Chandra source that is associated with the AGN candidate radio source NGC 6334B. Another possible model is the extended Chandra source of nonthermal emission from NGC 6334 that can also account for the hard X-ray emission observed by INTEGRAL. The origin of the emission in this scenario is due to electron acceleration in energetic outflows from massive early type stars. The possibility of emission from a young supernova remnant, as suggested by earlier infrared observations of NGC 6334, is constrained by the non-detection of 44Ti lines.

Non-stationary column accretion onto a surface of a magnetized neutron star is studied with a numerical code based on modified first-order Godunov method with splitting. Formation and evolution of shocks in the column is modeled for accretion rates ranging from 10^{15} g/s to 10^{16} g/s and surface magnetic fields ranging from 5*10^{11} G to 10^{13} G. Non-stationary solutions with plasma deceleration at collisionless oscillating shocks are found. The kinetic energy of the accreting flow efficiently transforms into a cyclotron radiation field. Collisionless stopping of the flow allows a substantial part of accreting CNO nuclei to avoid spallation and reach the surface. The nuclei survival fraction depends on the surface magnetic field, being higher at lower magnetic fields.

Spatially resolved images of the galactic supernova remnant G78.2+2.1 (gamma-Cygni) in hard X-ray energy bands from 25 keV to 120 keV are obtained with the IBIS-ISGRI imager aboard the International Gamma-Ray Astrophysics Laboratory INTEGRAL. The images are dominated by localized clumps of about ten arcmin in size. The flux of the most prominent North-Western (NW) clump is (1.7 +/- 0.4) 10^{-11} erg/cm^2/s in the 25-40 keV band. The observed X-ray fluxes are in agreement with extrapolations of soft X-ray imaging observations of gamma-Cygni by ASCA GIS and spatially unresolved RXTE PCA data. The positions of the hard X-ray clumps correlate with bright patches of optical line emission, possibly indicating the presence of radiative shock waves in a shocked cloud. The observed spatial structure and spectra are consistent with model predictions of hard X-ray emission from nonthermal electrons accelerated by a radiative shock in a supernova interacting with an interstellar cloud, but the powerful stellar wind of the O9V star HD 193322 is a plausible candidate for the NW source as well.