Fast radio bursts are bright, millisecond-scale radio flashes of yet unknown physical origin. Recently, their extragalactic nature has been demonstrated, and an increasing number of the sources have been found to repeat. Young, highly magnetized, isolated neutron stars - magnetars - have been suggested as the most promising candidates for fast radio burst progenitors owing to their energetics and high X-ray flaring activity. Here we report the detection with the Konus-Wind detector of a hard X-ray event of April 28, 2020, temporarily coincident with a bright, two-peak radio burst with properties remarkably similar to those of fast radio bursts. The source of the radio burst is located in the direction to the Galactic magnetar SGR 1935+2154, which recently entered an active state. We show that a separation between two peaks of the double-peaked X-ray burst is nearly the same as that of the radio peaks, confirming that the X-ray and radio emission most likely have a common origin. Thus, this is the first simultaneous detection of a fast radio burst from a Galactic magnetar and its high-energy counterpart. The total energy emitted in X-rays in this burst is typical of bright short magnetar bursts, but its spiky light curve and an unusual hardness of the energy spectrum distinguish the April 28 event among multiple 'ordinary' flares detected from SGR 1935+2154 previously. This, and a recent non-detection of radio emission from two typical bright soft bursts from the same magnetar, may imply the existence of a special class of magnetar flares capable of producing fast radio bursts.

The giant gamma-ray flares of the Crab nebula discovered by AGILE and Fermi observatories came as a surprise and have challenged the existing models of pulsar wind nebulae. We have carried out an analysis of 10.5 years of Fermi-LAT observations (Aug 2008 -- Feb 2019) and investigated variability of the Crab nebula in the 100-300 MeV range. Besides the flares, we found several month long depressions of the gamma-ray flux and identified several cases of sharp flux drops, where during one week the flux decreased by an order of magnitude with respect to its average value. No statistically significant variations of the nebula flux in the E>10 GeV range were found in the data. We discuss possible implications of the observed gamma-ray flux depressions on the model of synchrotron emission of the Crab nebula.

Non-thermal particles and high-energy radiation can play a role in the dynamical processes in star-forming regions and provide an important piece of the multiwavelength observational picture of their structure and components. Powerful stellar winds and supernovae in compact clusters of massive stars and OB associations are known to be favourable sites of high-energy particle acceleration and sources of non-thermal radiation and neutrinos. Namely, young massive stellar clusters are likely sources of the PeV (petaelectronvolt) regime cosmic rays (CRs). They can also be responsible for the cosmic ray composition, e.g., 22Ne/20Ne anomalous isotopic ratio in CRs. Efficient particle acceleration can be accompanied by super-adiabatic amplification of the fluctuating magnetic fields in the systems converting a part of kinetic power of the winds and supernovae into the magnetic energy through the CR-driven instabilities. The escape and CR propagation in the vicinity of the sources are affected by the non-linear CR feedback. These effects are expected to be important in starburst galaxies, which produce high-energy neutrinos and gamma-rays. We give a brief review of the theoretical models and observational data on high-energy particle acceleration and their radiation in star-forming regions with young stellar population.

In recent years, exciting developments have taken place in the identification of the role of cosmic rays in star-forming environments. Observations from radio to infrared wavelengths and theoretical modelling have shown that low-energy cosmic rays (<1 TeV) play a fundamental role in shaping the chemical richness of the interstellar medium, determining the dynamical evolution of molecular clouds. In this review we summarise in a coherent picture the main results obtained by observations and by theoretical models of propagation and generation of cosmic rays, from the smallest scales of protostars and circumstellar discs, to young stellar clusters, up to Galactic and extragalactic scales. We also discuss the new fields that will be explored in the near future thanks to new generation instruments, such as: CTA, for the $\gamma$-ray emission from high-mass protostars; SKA and precursors, for the synchrotron emission at different scales; and ELT/HIRES, JWST, and ARIEL, for the impact of cosmic rays on exoplanetary atmospheres and habitability.

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.

For most elements, the isotopic ratios seen in cosmic rays (CRs) are similar to those in the solar wind. The most important exception to this is $^{22}$Ne/$^{20}$Ne where the CR value is $\sim 5$ times that of the solar wind. According to most recent models of nucleosynthesis, a large amount of $^{22}\mathrm{Ne}$ is generated in Wolf-Rayet (WR) stars. In the winds of carbon sequence WR stars, i.e., WC stars, the isotopic ratio $^{22}\mathrm{Ne}$/$^{20}\mathrm{Ne}$ can be much larger than in the solar wind. Here, we consider CRs produced by $^{22}\mathrm{Ne}$-enriched WR winds in compact clusters of young massive stars like Westerlund~1. We assume that efficient CR acceleration in clusters occurs over the WR lifetime in an ensemble of shock waves from multiple massive star winds. We estimate the fraction of all Galactic CRs such sources may produce for a given set of parameters.

It has been shown that supernova blast waves interacting with winds from massive stars in compact star clusters may be capable of producing cosmic-ray (CR) protons to above $10^{17}$ eV. We give a brief description of the colliding-shock-flows mechanism and look at generalizations of the diffusion of ~ 100 PeV CRs in the turbulent galactic magnetic field present in the galactic disk. We calculate the temporal evolution of the CR anisotropy from a possible distribution of young compact massive star clusters assuming the sources are intermittent on time scales of a few million years, i.e., comparable to their residence time in the Milky Way. Within the confines of our model, we determine the galactic/extra-galactic fraction of high-energy CRs resulting in anisotropies consistent with observed values. We find that galactic star clusters may contribute a substantial fraction of ~ 100 PeV CRs without producing anisotropies above observed limits.

Relativistic particle acceleration in collisionless shocks of supernova remnants is accompanied by magnetic field amplification from cosmic ray (CR) driven plasma instabilities. Bell's fast CR-current instability is predicted to produce turbulence with a non-zero mean electric field in the shock precursor. We present a Monte Carlo model of Fermi shock acceleration explicitly taking into account an effective mean upstream electric field. Our model is nonlinear and includes the backreaction effects of efficient Fermi acceleration on the shock structure.

All-Sky-ASTROGAM is a gamma-ray observatory operating in a broad energy range, 100 keV to a few hundred MeV, recently proposed as the "Fast" (F) mission of the European Space Agency for a launch in 2028 to an L2 orbit. The scientific payload is composed of a unique gamma-ray imaging monitor for astrophysical transients, with very large field of view (almost 4$\pi$ sr) and optimal sensitivity to detect bright and intermediate flux sources (gamma-ray bursts, active galactic nuclei, X-ray binaries, supernovae and novae) at different timescales ranging from seconds to months. The mission will operate in a maturing gravitational wave and multi-messenger epoch, opening up new and exciting synergies.

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.

This white paper highlights compact object and fundamental physics science opportunities afforded by high-throughput broadband (0.1-60 keV) X-ray polarization observations. X-ray polarimetry gives new observables with geometric information about stellar remnants which are many orders of magnitude too small for direct imaging. The X-ray polarimetric data also reveal details about the emission mechanisms and the structure of the magnetic fields in and around the most extreme objects in the Universe. Whereas the Imaging X-ray Polarimetry Explorer (IXPE) to be launched in 2021 will obtain first results for bright objects, a follow-up mission could be one order of magnitude more sensitive and would be able to use a broader bandpass to perform physics type experiments for representative samples of sources.

Latest precise cosmic-ray (CR) measurements and present gamma-ray observations have started challenging our understanding of CR transport and interaction in the Galaxy. Moreover, because the density of CRs is similar to the density of the magnetic field, gas, and starlight in the interstellar medium (ISM), CRs are expected to affect the ISM dynamics, including the physical and chemical processes that determine transport and star formation. In this context, observations of gamma-ray emission at MeV energies produced by the low-energy CRs are very important and urgent. A telescope covering the energy range between ~0.1 MeV and a few GeV with a sensitivity more than an order of magnitude better than previous instruments would allow for the first time to study in detail the low-energy CRs, providing information on their sources, their spectra throughout the Galaxy, their abundances, transport properties, and their role on the evolution of the Galaxy and star formation. Here we discuss the scientific prospects for studies of CRs, ISM (gas, interstellar photons, and magnetic fields) and associated gamma-ray emissions with such an instrument.

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.

e-ASTROGAM (enhanced ASTROGAM) is a breakthrough Observatory space mission, with a detector composed by a Silicon tracker, a calorimeter, and an anticoincidence system, dedicated to the study of the non-thermal Universe in the photon energy range from 0.3 MeV to 3 GeV - the lower energy limit can be pushed to energies as low as 150 keV for the tracker, and to 30 keV for calorimetric detection. The mission is based on an advanced space-proven detector technology, with unprecedented sensitivity, angular and energy resolution, combined with polarimetric capability. Thanks to its performance in the MeV-GeV domain, substantially improving its predecessors, e-ASTROGAM will open a new window on the non-thermal Universe, making pioneering observations of the most powerful Galactic and extragalactic sources, elucidating the nature of their relativistic outflows and their effects on the surroundings. With a line sensitivity in the MeV energy range one to two orders of magnitude better than previous generation instruments, e-ASTROGAM will determine the origin of key isotopes fundamental for the understanding of supernova explosion and the chemical evolution of our Galaxy. The mission will provide unique data of significant interest to a broad astronomical community, complementary to powerful observatories such as LIGO-Virgo-GEO600-KAGRA, SKA, ALMA, E-ELT, TMT, LSST, JWST, Athena, CTA, IceCube, KM3NeT, and LISA.

e-ASTROGAM is a concept for a breakthrough observatory space mission carrying a gamma-ray telescope dedicated to the study of the non-thermal Universe in the photon energy range from 0.15 MeV to 3 GeV. The lower energy limit can be pushed down to energies as low as 30 keV for gamma-ray burst detection with the calorimeter. The mission is based on an advanced space-proven detector technology, with unprecedented sensitivity, angular and energy resolution, combined with remarkable polarimetric capability. Thanks to its performance in the MeV-GeV domain, substantially improving its predecessors, e-ASTROGAM will open a new window on the non-thermal Universe, making pioneering observations of the most powerful Galactic and extragalactic sources, elucidating the nature of their relativistic outflows and their effects on the surroundings. With a line sensitivity in the MeV energy range one to two orders of magnitude better than previous and current generation instruments, e-ASTROGAM will determine the origin of key isotopes fundamental for the understanding of supernova explosion and the chemical evolution of our Galaxy. The mission will be a major player of the multiwavelength, multimessenger time-domain astronomy of the 2030s, and provide unique data of significant interest to a broad astronomical community, complementary to powerful observatories such as LISA, LIGO, Virgo, KAGRA, the Einstein Telescope and the Cosmic Explorer, IceCube-Gen2 and KM3NeT, SKA, ALMA, JWST, E-ELT, LSST, Athena, and the Cherenkov Telescope Array.

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.

We give a brief review of the origin and acceleration of cosmic rays (CRs), emphasizing the production of CRs at different stages of supernova evolution by the first-order Fermi shock acceleration mechanism. We suggest that supernovae with trans-relativistic outflows, despite being rather rare, may accelerate CRs to energies above 10$^{18}$ eV over the first year of their evolution. Supernovae in young compact clusters of massive stars, and interaction powered superluminous supernovae, may accelerate CRs well above the PeV regime. We discuss the acceleration of the bulk of the galactic CRs in isolated supernova remnants and re-acceleration of escaped CRs by the multiple shocks present in superbubbles produced by associations of OB stars. The effects of magnetic field amplification by CR driven instabilities, as well as superdiffusive CR transport, are discussed for nonthermal radiation produced by nonlinear shocks of all speeds including trans-relativistic ones.

We discuss the physics case for and the concept of a medium-scale axion helioscope with sensitivities in the axion-photon coupling a few times better than CERN Axion Solar Telescope (CAST). Search for an axion-like particle with these couplings is motivated by several persistent astrophysical anomalies. We present early conceptual design, existing infrastructure, projected sensitivity and timeline of such a helioscope (Troitsk Axion Solar Telescope Experiment, TASTE) to be constructed in the Institute for Nuclear Research, Troitsk, Russia. The proposed instrument may be also used for the search of dark-matter halo axions.

We discuss a specific population of galactic PeVatrons which may be the main source of the galactic cosmic-ray (CR) component well above PeV energies. Supernovae in compact clusters of massive stars are proposed as powerful sources of CRs, neutrinos, and gamma-ray emission. Numerical simulations of non-linear Fermi acceleration at converging shock flows have revealed that these accelerators can provide very hard spectra of protons up to $10^{16}-10^{17}$ eV which is well above the "knee" in the all-particle CR spectrum at about $3\times10^{15}$ eV. We suggest that known supernova remnants interacting with stellar winds in the compact clusters of young massive stars Westerlund I and Cl*1806-20 can be associated with the sources of the TeV gamma-ray emission detected by H.E.S.S. and may be responsible for a fraction of the high-energy neutrinos detected with the IceCube observatory. A recent CR composition measurement with the LOFAR array has revealed a light-mass component possibly dominating the all-particle spectrum at energies around $10^{17}$ eV. Such a strong light component (mainly protons and helium) may require specific galactic CR sources such as supernovae interacting with compact clusters of massive stars in addition to isolated supernova remnants.

e-ASTROGAM (`enhanced ASTROGAM') is a breakthrough Observatory mission dedicated to the study of the non-thermal Universe in the photon energy range from 0.3 MeV to 3 GeV. The mission is based on an advanced space-proven detector technology, with unprecedented sensitivity, angular and energy resolution, combined with polarimetric capability. In the largely unexplored MeV-GeV domain, e-ASTROGAM will open a new window on the non-thermal Universe, making pioneering observations of the most powerful Galactic and extragalactic sources, elucidating the nature of their relativistic outflows and their effects on Galactic ecosystems. With a line sensitivity in the MeV energy range one to two orders of magnitude better than previous generation instruments, will determine the origin of key isotopes fundamental for the understanding of supernova explosion and the chemical evolution of our Galaxy. The mission will provide unique data of significant interest to a broad astronomical community, complementary to powerful observatories such as LIGO-Virgo-GEO600-KAGRA, SKA, ALMA, E-ELT, TMT, LSST, JWST, Athena, CTA, IceCube, KM3NeT, and the promise of eLISA. Keywords: High-energy gamma-ray astronomy, High-energy astrophysics, Nuclear Astrophysics, Compton and Pair creation telescope, Gamma-ray bursts, Active Galactic Nuclei, Jets, Outflows, Multiwavelength observations of the Universe, Counterparts of gravitational waves, Fermi, Dark Matter, Nucleosynthesis, Early Universe, Supernovae, Cosmic Rays, Cosmic antimatter.

Relativistic astrophysical collisionless shocks represent outstanding dissipation agents of the huge power of relativistic outflows produced by accreting black holes, core collapsed supernovae and other objects into multi-messenger radiation (cosmic rays, neutrinos, electromagnetic radiation). This article provides a theoretical discussion of the fundamental physical ingredients of these extreme phenomena. In the context of weakly magnetized shocks, in particular, it is shown how the filamentation type instabilities, which develop in the precursor of pair dominated or electron-ion shocks, provide the seeds for the scattering of high energy particles as well as the agent which preheats and slows down the incoming precursor plasma. This analytical discussion is completed with a mesoscopic, non-linear model of particle acceleration in relativistic shocks based on Monte Carlo techniques. This Monte Carlo model uses a semi-phenomenological description of particle scattering which allows it to calculate the back-reaction of accelerated particles on the shock structure on length and momentum scales which are currently beyond the range of microscopic particle-in-cell (PIC) simulations.

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)

Fast collisionless shocks in cosmic plasmas convert their kinetic energy flow into the hot downstream thermal plasma with a substantial fraction of energy going into a broad spectrum of superthermal charged particles and magnetic fluctuations. The superthermal particles can penetrate into the shock upstream region producing an extended shock precursor. The cold upstream plasma flow is decelerated by the force provided by the superthermal particle pressure gradient. In high Mach number collisionless shocks, efficient particle acceleration is likely coupled with turbulent magnetic field amplification (MFA) generated by the anisotropic distribution of accelerated particles. This anisotropy is determined by the fast particle transport making the problem strongly nonlinear and multi-scale. Here, we present a nonlinear Monte Carlo model of collisionless shock structure with super-diffusive propagation of high-energy Fermi accelerated particles coupled to particle acceleration and MFA which affords a consistent description of strong shocks. A distinctive feature of the Monte Carlo technique is that it includes the full angular anisotropy of the particle distribution at all precursor positions. The model reveals that the super-diffusive transport of energetic particles (i.e., Levy-walk propagation) generates a strong quadruple anisotropy in the precursor particle distribution. The resultant pressure anisotropy of the high-energy particles produces a non-resonant mirror-type instability which amplifies compressible wave modes with wavelengths longer than the gyroradii of the highest energy protons produced by the shock.

The e-ASTROGAM is a gamma-ray space mission to be proposed as the M5 Medium-size mission of the European Space Agency. It is dedicated to the observation of the Universe with unprecedented sensitivity in the energy range 0.2 - 100 MeV, extending up to GeV energies, together with a groundbreaking polarization capability. It is designed to substantially improve the COMPTEL and Fermi sensitivities in the MeV-GeV energy range and to open new windows of opportunity for astrophysical and fundamental physics space research. e-ASTROGAM will operate as an open astronomical observatory, with a core science focused on (1) the activity from extreme particle accelerators, including gamma-ray bursts and active galactic nuclei and the link of jet astrophysics to the new astronomy of gravitational waves, neutrinos, ultra-high energy cosmic rays, (2) the high-energy mysteries of the Galactic center and inner Galaxy, including the activity of the supermassive black hole, the Fermi Bubbles, the origin of the Galactic positrons, and the search for dark matter signatures in a new energy window; (3) nucleosynthesis and chemical evolution, including the life cycle of elements produced by supernovae in the Milky Way and the Local Group of galaxies. e-ASTROGAM will be ideal for the study of high-energy sources in general, including pulsars and pulsar wind nebulae, accreting neutron stars and black holes, novae, supernova remnants, and magnetars. And it will also provide important contributions to solar and terrestrial physics. The e-ASTROGAM telescope is optimized for the simultaneous detection of Compton and pair-producing gamma-ray events over a large spectral band. It is based on a very high technology readiness level for all subsystems and includes many innovative features for the detectors and associated electronics.

We report on six new Chandra observations of the Geminga pulsar wind nebula (PWN). The PWN consists of three distinct elongated structures - two $\approx 0.2 d_{250}$ pc long lateral tails and a segmented axial tail of $\approx 0.05 d_{250}$ pc length, where $d_{250}=d/(250 {\rm pc})$. The photon indices of the power law spectra of the lateral tails, $\Gamma \approx 1$, are significantly harder than those of the pulsar ($\Gamma \approx 1.5$) and the axial tail ($\Gamma \approx 1.6$). There is no significant diffuse X-ray emission between the lateral tails -- the ratio of the X-ray surface brightness between the south tail and this sky area is at least 12. The lateral tails apparently connect directly to the pulsar and show indication of moving footpoints. The axial tail comprises time-variable emission blobs. However, there is no evidence for constant or decelerated outward motion of these blobs. Different physical models are consistent with the observed morphology and spectra of the Geminga PWN. In one scenario, the lateral tails could represent an azimuthally asymmetric shell whose hard emission is caused by the Fermi acceleration mechanism of colliding winds. In another scenario, the lateral tails could be luminous, bent polar outflows, while the blobs in the axial tail could represent a crushed torus. In a resemblance to planetary magnetotails, the blobs of the axial tail might also represent short-lived plasmoids which are formed by magnetic field reconnection in the relativistic plasma of the pulsar wind tail.

We report on Chandra X-ray Observatory (CXO) observations of the pulsar wind nebula (PWN) associated with PSR B0355+54 (eight observations with a 395 ks total exposure, performed over an 8 month period). We investigated the spatial and spectral properties of the emission coincident with the pulsar, compact nebula (CN), and extended tail. We find that the CN morphology can be interpreted in a way that suggests a small angle between the pulsar spin axis and our line-of-sight, as inferred from the radio data. On larger scales, emission from the 7' (2 pc) tail is clearly seen. We also found hints of two faint extensions nearly orthogonal to the direction of the pulsar's proper motion. The spectrum extracted at the pulsar position can be described with an absorbed power-law + blackbody model. The nonthermal component can be attributed to magnetospheric emission, while the thermal component can be attributed to emission from either a hot spot (e.g., a polar cap) or the entire neutron star surface. Surprisingly, the spectrum of the tail shows only a slight hint of cooling with increasing distance from the pulsar. This implies either a low magnetic field with fast flow speed, or particle re-acceleration within the tail. We estimate physical properties of the PWN and compare the morphologies of the CN and the extended tail with those of other bow shock PWNe observed with long CXO exposures.

The dynamics, energetics and evolution of pulsar wind nebulae (PWNe) and supernova remnants (SNRs), are strongly affected by their magnetic field strength and distribution. They are usually strong, extended, sources of non-thermal X-ray radiation, producing intrinsically polarised radiation. The energetic wind around pulsars produces a highly-magnetised, structured flow, often displaying a jet and a torus and different features (i.e. wisps, knots). This magnetic-dominant wind evolves as it moves away from the pulsar magnetosphere to the surrounding large-scale nebula, becoming kinetic-dominant. Basic aspects such how this conversion is produced, or how the jets and torus are formed, as well as the level of turbulence in the nebula are still unknown. Likewise, the processes ruling the acceleration of particles in shell-like SNRs up to 1e15 eV, including the amplification of the magnetic field, are not clear yet. Imaging polarimetry in this regard is crucial to localise the regions of shock acceleration and to measure the strength and the orientation of the magnetic field at these emission sites. X-ray polarimetry with the X-ray Imaging Polarimetry Explorer (XIPE) will allow the understanding of the magnetic field structure and intensity on different regions in SNRs and PWNe, helping to unveil long-standing questions such as i.e. acceleration of cosmic rays in SNRs or magnetic-to-kinetic energy transfer. SNRs and PWNe also represent the largest population of Galactic very-high energy gamma-ray sources, therefore the study of their magnetic distribution with XIPE will provide fundamental ingredients on the investigation of those sources at very high energies. We will discuss the physics case related to SNRs and PWNe and the expectations of the XIPE observations of some of the most prominent SNRs and PWNe.

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.

Collisionless shocks, that is shocks mediated by electromagnetic processes, are customary in space physics and in astrophysics. They are to be found in a great variety of objects and environments: magnetospheric and heliospheric shocks, supernova remnants, pulsar winds and their nebul\ae, active galactic nuclei, gamma-ray bursts and clusters of galaxies shock waves. Collisionless shock microphysics enters at different stages of shock formation, shock dynamics and particle energization and/or acceleration. It turns out that the shock phenomenon is a multi-scale non-linear problem in time and space. It is complexified by the impact due to high-energy cosmic rays in astrophysical environments. This review adresses the physics of shock formation, shock dynamics and particle acceleration based on a close examination of available multi-wavelength or in-situ observations, analytical and numerical developments. A particular emphasize is made on the different instabilities triggered during the shock formation and in association with particle acceleration processes with regards to the properties of the background upstream medium. It appears that among the most important parameters the background magnetic field through the magnetization and its obliquity is the dominant one. The shock velocity that can reach relativistic speeds has also a strong impact over the development of the micro-instabilities and the fate of particle acceleration. Recent developments of laboratory shock experiments has started to bring some new insights in the physics of space plasma and astrophysical shock waves. A special section is dedicated to new laser plasma experiments probing shock physics

We include a general form for the scattering mean free path in a nonlinear Monte Carlo model of relativistic shock formation and Fermi acceleration. Particle-in-cell (PIC) simulations, as well as analytic work, suggest that relativistic shocks tend to produce short-scale, self-generated magnetic turbulence that leads to a scattering mean free path (mfp) with a stronger momentum dependence than the mfp ~ p dependence for Bohm diffusion. In unmagnetized shocks, this turbulence is strong enough to dominate the background magnetic field so the shock can be treated as parallel regardless of the initial magnetic field orientation, making application to gamma-ray bursts (GRBs), pulsar winds, Type Ibc supernovae, and extra-galactic radio sources more straightforward and realistic. In addition to changing the scale of the shock precursor, we show that, when nonlinear effects from efficient Fermi acceleration are taken into account, the momentum dependence of the mfp has an important influence on the efficiency of cosmic-ray production as well as the accelerated particle spectral shape. These effects are absent in nonrelativistic shocks and do not appear in relativistic shock models unless nonlinear effects are self-consistently described. We show, for limited examples, how the changes in Fermi acceleration translate to changes in the intensity and spectral shape of gamma-ray emission from proton-proton interactions and pion-decay radiation.

Clusters of galaxies are the largest gravitationally bounded structures in the Universe dominated by dark matter. We review the observational appearance and physical models of plasma structures in clusters of galaxies. Bubbles of relativistic plasma which are inflated by supermassive black holes of AGNs, cooling and heating of the gas, large scale plasma shocks, cold fronts, non-thermal halos and relics are observed in clusters. These constituents are reflecting both the formation history and the dynamical properties of clusters of galaxies. We discuss X-ray spectroscopy as a tool to study the metal enrichment in clusters and fine spectroscopy of Fe X-ray lines as a powerful diagnostics of both the turbulent plasma motions and the energetics of the non-thermal electron populations. The knowledge of the complex dynamical and feedback processes is necessary to understand the energy and matter balance as well as to constrain the role of the non-thermal components of clusters.

Starforming factories in galaxies produce compact clusters and loose associations of young massive stars. Fast radiation-driven winds and supernovae input their huge kinetic power into the interstellar medium in the form of highly supersonic and superalfvenic outflows. Apart from gas heating, collisionless relaxation of fast plasma outflows results in fluctuating magnetic fields and energetic particles. The energetic particles comprise a long-lived component which may contain a sizeable fraction of the kinetic energy released by the winds and supernova ejecta and thus modify the magnetohydrodynamic flows in the systems. We present a concise review of observational data and models of nonthermal emission from starburst galaxies, superbubbles, and compact clusters of massive stars. Efficient mechanisms of particle acceleration and amplification of fluctuating magnetic fields with a wide dynamical range in starburst regions are discussed. Sources of cosmic rays, neutrinos and multi-wavelength nonthermal emission associated with starburst regions including potential galactic "PeVatrons" are reviewed in the global galactic ecology context.

Starburst regions with multiple powerful winds of young massive stars and supernova remnants are favorable sites for high-energy cosmic ray acceleration. A supernova shock colliding with a fast wind from a compact cluster of young stars allows the acceleration of protons to energies well above the standard limits of diffusive shock acceleration in an isolated SN. The proton spectrum in such a wind-supernova PeV accelerator is hard with a large flux in the high-energy-end of the spectrum producing copious gamma-rays and neutrinos in inelastic nuclear collisions. We argue that SN shocks in the Westerlund 1 cluster in the Milky Way may accelerate protons to about 40 PeV. Once accelerated, these CRs will diffuse into surrounding dense clouds and produce neutrinos with fluxes sufficient to explain a fraction of the events detected by IceCube Observatory from the inner Galaxy.

We have modeled the simultaneous first-order Fermi shock acceleration of protons, electrons, and helium nuclei by relativistic shocks. By parameterizing the particle diffusion, our steady-state Monte Carlo simulation allows us to follow particles from particle injection at nonthermal thermal energies to above PeV energies, including the nonlinear smoothing of the shock structure due to cosmic-ray (CR) backpressure. We observe the mass-to-charge (A/Z) enhancement effect believed to occur in efficient Fermi acceleration in non-relativistic shocks and we parameterize the transfer of ion energy to electrons seen in particle-in-cell (PIC) simulations. For a given set of environmental and model parameters, the Monte Carlo simulation determines the absolute normalization of the particle distributions and the resulting synchrotron, inverse-Compton, and pion-decay emission in a largely self-consistent manner. The simulation is flexible and can be readily used with a wide range of parameters typical of gamma-ray burst (GRB) afterglows. We describe some preliminary results for photon emission from shocks of different Lorentz factors and outline how the Monte Carlo simulation can be generalized and coupled to hydrodynamic simulations of GRB blast waves. We assume Bohm diffusion for simplicity but emphasize that the nonlinear effects we describe stem mainly from an extended shock precursor where higher energy particles diffuse further upstream. Quantitative differences will occur with different diffusion models, particularly for the maximum CR energy and photon emission, but these nonlinear effects should be qualitatively similar as long as the scattering mean free path is an increasing function of momentum.

Astrophysical shocks are often collisionless shocks. An open question about collisionless shocks is whether electrons and ions each establish their own post-shock temperature, or whether they quickly equilibrate in the shock region. Here we provide simple relations for the minimal amount of equilibration to expect. The basic assumption is that the enthalpy-flux of the electrons is conserved separately, but that all particle species should undergo the same density jump across the the shock. This assumption results in an analytic treatment of electron-ion equilibration that agrees with observations of collisionless shocks: at low Mach numbers ($<2$) the electrons and ions are close to equilibration, whereas for Mach numbers above $M \sim 60$ the electron-ion temperature ratio scales with the particle masses $T_e/T_i = m_e/m_i$. In between these two extremes the electron-ion temperature ratio scales as $T_e/T_i \propto 1/M_s^2$. This relation also hold if adiabatic compression of the electrons is taken into account. For magnetised plasmas the compression is governed by the magnetosonic Mach number, whereas the electron-ion temperatures are governed by the sonic Mach number. The derived equations are in agreement with observational data at low Mach numbers, but for supernova remnants the relation requires that the inferred Mach numbers for the observations are over- estimated, perhaps as a result of upstream heating in the cosmic-ray precursor. In addition to predicting a minimal electron/ion temperature ratio, we also heuristically incorporate ion-electron heat exchange at the shock, quantified with a dimensionless parameter ${\xi}$. Comparing the model to existing observations in the solar system and supernova remnants suggests that the data are best described by ${\xi} \sim 5$ percent. (Abridged abstract.)

We obtained six observations of PSR J1741-2054 using the $Chandra$ ACIS-S detector totaling $\sim$300 ks. By registering this new epoch of observations to an archival observation taken 3.2 years earlier using X-ray point sources in the field of view, we have measured the pulsar proper motion at $\mu =109 \pm 10 {\rm mas\ yr}^{-1}$ in a direction consistent with the symmetry axis of the observed H$\alpha$ nebula. We investigated the inferred past trajectory of the pulsar but find no compelling association with OB associations in which the progenitor may have originated. We confirm previous measurements of the pulsar spectrum as an absorbed power law with photon index $\Gamma$=2.68$\pm$0.04, plus a blackbody with an emission radius of (4.5$^{+3.2}_{-2.5})d_{0.38}$ km, for a DM-estimated distance of $0.38d_{0.38}$ kpc and a temperature of $61.7\pm3.0$ eV. Emission from the compact nebula is well described by an absorbed power law model with a photon index of $\Gamma$ = 1.67$\pm$0.06, while the diffuse emission seen as a trail extending northeast of the pulsar shows no evidence of synchrotron cooling. We also applied image deconvolution techniques to search for small-scale structures in the immediate vicinity of the pulsar, but found no conclusive evidence for such structures.

The giant molecular clouds (MCs) found in the Milky Way and similar galaxies play a crucial role in the evolution of these systems. The supernova explosions that mark the death of massive stars in these regions often lead to interactions between the supernova remnants (SNRs) and the clouds. These interactions have a profound effect on our understanding of SNRs. Shocks in SNRs should be capable of accelerating particles to cosmic ray (CR) energies with efficiencies high enough to power Galactic CRs. X-ray and gamma-ray studies have established the presence of relativistic electrons and protons is some SNRs and provided strong evidence for diffusive shock acceleration as the primary acceleration mechanism, including strongly amplified magnetic fields, temperature and ionization effects on the shock-heated plasmas, and modifications to the dynamical evolution of some systems. Because protons dominate the overall energetics of the CRs, it is crucial to understand this hadronic component even though electrons are much more efficient radiators and it can be difficult to identify the hadronic component. However, near MCs the densities are sufficiently high to allow the gamma-ray emission to be dominated by protons. Thus, these interaction sites provide some of our best opportunities to constrain the overall energetics of these particle accelerators. Here we summarize some key properties of interactions between SNRs and MCs, with an emphasis on recent X-ray and gamma-ray studies that are providing important constraints on our understanding of cosmic rays in our Galaxy.

We present a nonlinear Monte Carlo model of efficient diffusive shock acceleration (DSA) where the magnetic turbulence responsible for particle diffusion is calculated self-consistently from the resonant cosmic-ray (CR) streaming instability, together with non-resonant short- and long-wavelength CR-current-driven instabilities. We include the backpressure from CRs interacting with the strongly amplified magnetic turbulence which decelerates and heats the super-alfvenic flow in the extended shock precursor. Uniquely, in our plane-parallel, steady-state, multi-scale model, the full range of particles, from thermal (~eV) injected at the viscous subshock, to the escape of the highest energy CRs (~PeV) from the shock precursor, are calculated consistently with the shock structure, precursor heating, magnetic field amplification (MFA), and scattering center drift relative to the background plasma. In addition, we show how the cascade of turbulence to shorter wavelengths influences the total shock compression, the downstream proton temperature, the magnetic fluctuation spectra, and accelerated particle spectra. A parameter survey is included where we vary shock parameters, the mode of magnetic turbulence generation, and turbulence cascading. From our survey results, we obtain scaling relations for the maximum particle momentum and amplified magnetic field as functions of shock speed, ambient density, and shock size.

The large-scale structure of the Universe formed from initially small perturbations in the cosmic density field, leading to galaxy clusters with up to 10^15 Msun at the present day. Here, we review the formation of structures in the Universe, considering the first primordial galaxies and the most massive galaxy clusters as extreme cases of structure formation where fundamental processes such as gravity, turbulence, cooling and feedback are particularly relevant. The first non-linear objects in the Universe formed in dark matter halos with 10^5-10^8 Msun at redshifts 10-30, leading to the first stars and massive black holes. At later stages, larger scales became non-linear, leading to the formation of galaxy clusters, the most massive objects in the Universe. We describe here their formation via gravitational processes, including the self-similar scaling relations, as well as the observed deviations from such self-similarity and the related non-gravitational physics (cooling, stellar feedback, AGN). While on intermediate cluster scales the self-similar model is in good agreement with the observations, deviations from such self-similarity are apparent in the core regions, where numerical simulations do not reproduce the current observational results. The latter indicates that the interaction of different feedback processes may not be correctly accounted for in current simulations. Both in the most massive clusters of galaxies as well as during the formation of the first objects in the Universe, turbulent structures and shock waves appear to be common, suggesting them to be ubiquitous in the non-linear regime.

Aims. Observations show nebulae around some massive stars but not around others. If observed, their chemical composition is far from homogeneous. Our goal is to put these observational features into the context of the evolution of massive stars and their circumstellar medium (CSM) and, more generally, to quantify the role of massive stars for the chemical and dynamical evolution of the ISM. Methods. Using the A-MAZE code, we perform 2d-axisymmetric hydrodynamical simulations of the evolution of the CSM, shaped by stellar winds, for a whole grid of massive stellar models from 15 to 120 Msun and following the stellar evolution from the zero-age main-sequence to the time of supernova explosion. In addition to the usual quantities, we also follow five chemical species: H, He, C, N, and O. Results. We show how various quantities evolve as a function of time: size of the bubble, position of the wind termination shock, chemical composition of the bubble, etc. The chemical composition of the bubble changes considerably compared to the initial composition, particularly during the red-supergiant (RSG) and Wolf-Rayet (WR) phases. In some extreme cases, the inner region of the bubble can be completely depleted in hydrogen and nitrogen, and is mainly composed of carbon, helium and oxygen. We argue why the bubble typically expands at a lower rate than predicted by self-similarity theory. In particular, the size of the bubble is very sensitive to the density of the ISM, decreasing by a factor of around 2.5 for each additional dex in ISM density. The bubble size also decreases with the metallicity of the central star, as low-metallicity stars have weaker winds. Our models qualitatively fit the observations of WR ejecta nebulae.

We present results from a Monte Carlo simulation of a parallel collisionless shock undergoing particle acceleration. Our simulation, which contains parameterized scattering and a particular thermal leakage injection model, calculates the feedback between accelerated particles ahead of the shock, which influence the shock precursor and "smooth" the shock, and thermal particle injection. We show that there is a transition between nonrelativistic shocks, where the acceleration efficiency can be extremely high and the nonlinear compression ratio can be substantially greater than the Rankine-Hugoniot value, and fully relativistic shocks, where diffusive shock acceleration is less efficient and the compression ratio remains at the Rankine-Hugoniot value. This transition occurs in the trans-relativistic regime and, for the particular parameters we use, occurs around a shock Lorentz factor = 1.5. We also find that nonlinear shock smoothing dramatically reduces the acceleration efficiency presumed to occur with large-angle scattering in ultra-relativistic shocks. Our ability to seamlessly treat the transition from ultra-relativistic to trans-relativistic to nonrelativistic shocks may be important for evolving relativistic systems, such as gamma-ray bursts and type Ibc supernovae. We expect a substantial evolution of shock accelerated spectra during this transition from soft early on to much harder when the blast-wave shock becomes nonrelativistic.

Energetic nonthermal particles (cosmic rays, CRs) are accelerated in supernova remnants, relativistic jets and other astrophysical objects. The CR energy density is typically comparable with that of the thermal components and magnetic fields. In this review we discuss mechanisms of magnetic field amplification due to instabilities induced by CRs. We derive CR kinetic and magnetohydrodynamic equations that govern cosmic plasma systems comprising the thermal background plasma, comic rays and fluctuating magnetic fields to study CR-driven instabilities. Both resonant and non-resonant instabilities are reviewed, including the Bell short-wavelength instability, and the firehose instability. Special attention is paid to the longwavelength instabilities driven by the CR current and pressure gradient. The helicity production by the CR current-driven instabilities is discussed in connection with the dynamo mechanisms of cosmic magnetic field amplification.

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.

We discuss processes in galactic cosmic ray (GCR) acceleration sites - supernova remnants, compact associations of young massive stars, and superbubbles. Mechanisms of efficient conversion of the mechanical power of the outflows driven by supernova shocks and fast stellar winds of young stars into magnetic fields and relativistic particles are discussed. The high efficiency of particle acceleration in the sources implies the importance of nonlinear feedback effects in a symbiotic relationship where the magnetic turbulence required to accelerate the CRs is created by the accelerated CRs themselves. Non-thermal emission produced by relativistic particles (both those confined in and those that escape from the cosmic accelerators) can be used to constrain the basic physical models of the GCR sources. High resolution X-ray synchrotron imaging, combined with GeV-TeV gamma ray spectra, is a powerful tool to probe the maximum energies of accelerated particles. Future MeV regime spectroscopy will provide unique information on the composition of accelerated particles.

Powerful stellar winds and supernova explosions with intense energy release in the form of strong shock waves can convert a sizeable part of the kinetic energy release into energetic particles. The starforming regions are argued as a favorable site of energetic particle acceleration and could be efficient sources of nonthermal emission. We present here a non-linear time-dependent model of particle acceleration in the vicinity of two closely approaching fast magnetohydrodynamic (MHD) shocks. Such MHD flows are expected to occur in rich young stellar cluster where a supernova is exploding in the vicinity of a strong stellar wind of a nearby massive star. We find that the spectrum of the high energy particles accelerated at the stage of two closely approaching shocks can be harder than that formed at a forward shock of an isolated supernova remnant. The presented method can be applied to model particle acceleration in a variety of systems with colliding MHD flows.

We evaluate the current status of supernova remnants as the sources of Galactic cosmic rays. We summarize observations of supernova remnants, covering the whole electromagnetic spectrum and describe what these obser- vations tell us about the acceleration processes by high Mach number shock fronts. We discuss the shock modification by cosmic rays, the shape and maximum energy of the cosmic-ray spectrum and the total energy budget of cosmic rays in and surrounding supernova remnants. Additionally, we discuss problems with supernova remnants as main sources of Galactic cosmic rays, as well as alternative sources.

In this review we confront the current theoretical understanding of particle acceleration at relativistic outflows with recent observational results on various source classes thought to involve such outflows, e.g. gamma-ray bursts, active galactic nuclei, and pulsar wind nebulae. We highlight the possible contributions of these sources to ultra-high-energy cosmic rays.

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.

Diffusive shock acceleration is the theory of particle acceleration through multiple shock crossings. In order for this process to proceed at a rate that can be reconciled with observations of high-energy electrons in the vicinity of the shock, and for cosmic rays protons to be accelerated to energies up to observed galactic values, significant magnetic field amplification is required. In this review we will discuss various theories on how magnetic field amplification can proceed in the presence of a cosmic ray population. On both short and long scales, cosmic ray streaming can induce instabilities that act to amplify the magnetic field. Developments in this area that have occurred over the past decade are the main focus of this paper.

Synchrotron radiation of ultra-relativistic particles accelerated in a pulsar wind nebula may dominate its spectrum up to gamma-ray energies. Because of the short cooling time of the gamma-ray emitting electrons, the gamma-ray emission zone is in the immediate vicinity of the acceleration site. The particle acceleration likely occurs at the termination shock of the relativistic striped wind, where multiple forced magnetic field reconnections provide strong magnetic fluctuations facilitating Fermi acceleration processes. The acceleration mechanisms imply the presence of stochastic magnetic fields in the particle acceleration region, which cause stochastic variability of the synchrotron emission. This variability is particularly strong in the steep gamma-ray tail of the spectrum, where modest fluctuations of the magnetic field lead to strong flares of spectral flux. In particular, stochastic variations of magnetic field, which may lead to quasi-cyclic gamma-ray flares, can be produced by the relativistic cyclotron ion instability at the termination shock. Our model calculations of the spectral and temporal evolution of synchrotron emission in the spectral cut-off regime demonstrate that the intermittent magnetic field concentrations dominate the gamma-ray emission from highest energy electrons and provide fast, strong variability even for a quasi-steady distribution of radiating particles. The simulated light curves and spectra can explain the very strong gamma-ray flares observed in the Crab nebula and the lack of strong variations at other wavelengths. The model predicts high polarization in the flare phase, which can be tested with future polarimetry observations.

We briefly discuss models of energetic particle acceleration by supernova shock in active starforming regions at different stages of their evolution. Strong shocks may strongly amplify magnetic fields due to cosmic ray driven instabilities. We discuss the magnetic field amplification emphasizing the role of the long-wavelength instabilities. Supernova shock propagating in the vicinity of a powerful stellar wind in a young stellar cluster is argued to increase the maximal CR energies at a given evolution stage of supernova remnant (SNR) and can convert a sizeable fraction of the kinetic energy release into energetic particles.

We present a spherically symmetric, core-collapse model of SNR RX J1713.7-3946 that includes a hydrodynamic simulation of the remnant evolution coupled to the efficient production of cosmic rays (CRs) by nonlinear diffusive shock acceleration (DSA). High-energy CRs that escape from the forward shock (FS) are propagated in surrounding dense material that simulates either a swept-up, pre-supernova shell or a nearby molecular cloud. The continuum emission from trapped and escaping CRs, along with the thermal X-ray emission from the shocked heated ISM behind the FS, integrated over the remnant, is compared against broadband observations. Our results show conclusively that, overall, the GeV-TeV emission is dominated by inverse-Compton from CR electrons if the supernova is isolated regardless of its type, i.e., not interacting with a >>100 Msun shell or cloud. If the SNR is interacting with a much larger mass >10^4 Msun, pion-decay from the escaping CRs may dominate the TeV emission, although a precise fit at high energy will depend on the still uncertain details of how the highest energy CRs are accelerated by, and escape from, the FS. Based on morphological and other constraints, we consider the 10^4 Msun pion-decay scenario highly unlikely for SNR RX J1713.7-3946 regardless of the details of CR escape. Importantly, even though CR electrons dominate the GeV-TeV emission, the efficient production of CR ions is an essential part of our leptonic model.

It is only now, with low-frequency radio telescopes, long exposures with high-resolution X-ray satellites and gamma-ray telescopes, that we are beginning to learn about the physics in the periphery of galaxy clusters. In the coming years, Sunyaev-Zeldovich telescopes are going to deliver further great insights into the plasma physics of these special regions in the Universe. The last years have already shown tremendous progress with detections of shocks, estimates of magnetic field strengths and constraints on the particle acceleration efficiency. X-ray observations have revealed shock fronts in cluster outskirts which have allowed inferences about the microphysical structure of shocks fronts in such extreme environments. The best indications for magnetic fields and relativistic particles in cluster outskirts come from observations of so-called radio relics, which are megaparsec-sized regions of radio emission from the edges of galaxy clusters. As these are difficult to detect due to their low surface brightness, only few of these objects are known. But they have provided unprecedented evidence for the acceleration of relativistic particles at shock fronts and the existence of muG strength fields as far out as the virial radius of clusters. In this review we summarise the observational and theoretical state of our knowledge of magnetic fields, relativistic particles and shocks in cluster outskirts.

We review recent progress on collisionless relativistic shocks. Kinetic instability theory is briefed including its predictions and limitations. The main focus is on numerical experiments in (i) pair and (ii) electron-nucleon plasmas. The main results are: (i) confirmation of shock evolution in non-magnetised relativistic plasma in 3D due to either the lepton-Weibel instability or the ion-Weibel instability; (ii) sensitive dependence on upstream magnetisation ; (iii) the sensitive dependence of particle dynamics on the upstream magnetic inclination angle $\thetabn$, where particles of $\thetabn>34^\circ$ cannot escape upstream, leading to the distinction between `sub-luminal' and `super-luminal' shocks; (iv) particles in ultra-relativistic shocks can hardly overturn the shock and escape to upstream; they may oscillate around the shock ramp for a long time, so to speak `surfing it' and thereby becoming accelerated by a kind of SDA; (v) these particles form a power law tail on the downstream distribution; their limitations are pointed out; (vi) recently developed methods permit the calculation of the radiation spectra emitted by the downstream high-energy particles; (vii) the Weibel-generated downstream magnetic fields form large amplitude vortices which could be advected by the downstream flow to large distances from the shock and possibly contribute to an extended strong field region; (viii) if cosmic rays are included, Bell-like modes can generate upstream magnetic turbulence at short and, by diffusive re-coupling, also long wavelengths in nearly parallel magnetic field shocks; (ix) advection of such large-amplitude waves should cause periodic reformation of the quasi-parallel shock and eject large amplitude magnetic field vortices downstream where they contribute to turbulence and to maintaining an extended region of large magnetic fields.

High-resolution Chandra observations of Tycho's SNR have revealed several sets of quasi-steady, high-emissivity, nearly-parallel X-ray stripes in some localized regions of the SNR. These stripes are most likely the result of cosmic-ray (CR) generated magnetic turbulence at the SNR blast wave. However, for the amazingly regular pattern of these stripes to appear requires the simultaneous action of a number of shock-plasma phenomena and is not predicted by most models of magnetic field amplification. A consistent explanation of these stripes yields information on the complex nonlinear plasma processes connecting efficient CR acceleration and magnetic field fluctuations in strong collisionless shocks. The nonlinear diffusive shock acceleration (NL-DSA) model described here, which includes magnetic field amplification from a cosmic-ray current driven instability, does predict stripes consistent with the synchrotron emission observations of Tycho's SNR. We argue that the local ambient mean magnetic field geometry determines the orientation of the stripes and therefore it can be reconstructed with the high resolution X-ray imaging. The estimated maximum energy of the CR protons responsible for the stripes is $\sim 10^{15}$\,eV. Furthermore, the model predicts that a specific X-ray polarization pattern, with a polarized fraction ~ 50%, accompanies the stripes, which can be tested with future X-ray polarimeter missions.

We propose to perform a continuously scanning all-sky survey from 200 keV to 80 MeV achieving a sensitivity which is better by a factor of 40 or more compared to the previous missions in this energy range. The Gamma-Ray Imaging, Polarimetry and Spectroscopy (GRIPS) mission addresses fundamental questions in ESA's Cosmic Vision plan. Among the major themes of the strategic plan, GRIPS has its focus on the evolving, violent Universe, exploring a unique energy window. We propose to investigate $\gamma$-ray bursts and blazars, the mechanisms behind supernova explosions, nucleosynthesis and spallation, the enigmatic origin of positrons in our Galaxy, and the nature of radiation processes and particle acceleration in extreme cosmic sources including pulsars and magnetars. The natural energy scale for these non-thermal processes is of the order of MeV. Although they can be partially and indirectly studied using other methods, only the proposed GRIPS measurements will provide direct access to their primary photons. GRIPS will be a driver for the study of transient sources in the era of neutrino and gravitational wave observatories such as IceCUBE and LISA, establishing a new type of diagnostics in relativistic and nuclear astrophysics. This will support extrapolations to investigate star formation, galaxy evolution, and black hole formation at high redshifts.

We review here some magnetic phenomena in astrophysical particle accelerators associated with collisionless shocks in supernova remnants, radio galaxies and clusters of galaxies. A specific feature is that the accelerated particles can play an important role in magnetic field evolution in the objects. We discuss a number of CR-driven, magnetic field amplification processes that are likely to operate when diffusive shock acceleration (DSA) becomes efficient and nonlinear. The turbulent magnetic fields produced by these processes determine the maximum energies of accelerated particles and result in specific features in the observed photon radiation of the sources. Equally important, magnetic field amplification by the CR currents and pressure anisotropies may affect the shocked gas temperatures and compression, both in the shock precursor and in the downstream flow, if the shock is an efficient CR accelerator. Strong fluctuations of the magnetic field on scales above the radiation formation length in the shock vicinity result in intermittent structures observable in synchrotron emission images. Resonant and non-resonant CR streaming instabilities in the shock precursor can generate mesoscale magnetic fields with scale-sizes comparable to supernova remnants and even superbubbles. This opens the possibility that magnetic fields in the earliest galaxies were produced by the first generation Population III supernova remnants and by clustered supernovae in star forming regions.

We briefly review sources of cosmic rays, their composition and spectra as well as their propagation in the galactic and extragalactic magnetic fields, both regular and fluctuating. A special attention is paid to the recent results of the X-ray and gamma-ray observations that shed light on the origin of the galactic cosmic rays and the challenging results of Pierre Auger Observatory on the ultra high energy cosmic rays. The perspectives of both high energy astrophysics and cosmic-ray astronomy to identify the sources of ultra high energy cosmic rays, the mechanisms of particle acceleration, to measure the intergalactic radiation fields and to reveal the structure of magnetic fields of very different scales are outlined.

We present a model of gamma-ray emission from core-collapse supernovae originating from the explosions of massive young stars. The fast forward shock of the supernova remnant (SNR) can accelerate particles by diffusive shock acceleration (DSA) in a cavern blown by a strong, pre-supernova stellar wind. As a fundamental part of nonlinear DSA, some fraction of the accelerated particles escape the shock and interact with a surrounding massive dense shell producing hard photon emission. To calculate this emission, we have developed a new Monte Carlo technique for propagating the cosmic rays (CRs) produced by the forward shock of the SNR, into the dense, external material. This technique is incorporated in a hydrodynamic model of an evolving SNR which includes the nonlinear feedback of CRs on the SNR evolution, the production of escaping CRs along with those that remain trapped within the remnant, and the broad-band emission of radiation from trapped and escaping CRs. While our combined CR-hydro-escape model is quite general and applies to both core collapse and thermonuclear supernovae, the parameters we choose for our discussion here are more typical of SNRs from very massive stars whose emission spectra differ somewhat from those produced by lower mass progenitors directly interacting with a molecular cloud.

In order for diffusive shock acceleration (DSA) to accelerate particles to high energies, the energetic particles must be able to interact with magnetic turbulence over a broad wavelength range. The weakly anisotropic distribution of accelerated particles, i.e., cosmic rays (CRs), is believed capable of producing this turbulence in a symbiotic relationship where the magnetic turbulence required to accelerate the CRs is created by the accelerated CRs themselves. In efficient DSA, this wave-particle interaction can be strongly nonlinear where CRs modify the plasma flow and the specific mechanisms of magnetic field amplification. Resonant interactions have long been known to amplify magnetic fluctuations on the scale of the CR gyroradius, and Bell (2004) showed that the CR current can efficiently amplify magnetic fluctuations with scales smaller than the CR gyroradius. Here, we show with a multi-scale, quasi-linear analysis that the presence of turbulence with scales shorter than the CR gyroradius enhances the growth of modes with scales longer than the gyroradius, at least for particular polarizations. We use a mean-field approach to average the equation of motion and the induction equation over the ensemble of magnetic field oscillations accounting for the anisotropy of relativistic particles on the background plasma. We derive the response of the magnetized CR current on magnetic field fluctuations and show that, in the presence of short-scale, Bell-type turbulence, long wavelength modes are amplified. The polarization, helicity, and angular dependence of the growth rates are calculated for obliquely propagating modes for wavelengths both below and above the CR mean free path. The long wavelength growth rates we estimate for typical supernova remnant parameters are sufficiently fast to suggest a fundamental increase in the maximum CR energy a given shock can produce.

The first gamma-ray line originating from outside the solar system that was ever detected is the 511 keV emission from positron annihilation in the Galaxy. Despite 30 years of intense theoretical and observational investigation, the main sources of positrons have not been identified up to now. Observations in the 1990's with OSSE/CGRO showed that the emission is strongly concentrated towards the Galactic bulge. In the 2000's, the SPI instrument aboard ESA's INTEGRAL gamma-ray observatory allowed scientists to measure that emission across the entire Galaxy, revealing that the bulge/disk luminosity ratio is larger than observed in any other wavelength. This mapping prompted a number of novel explanations, including rather "exotic ones (e.g. dark matter annihilation). However, conventional astrophysical sources, like type Ia supernovae, microquasars or X-ray binaries, are still plausible candidates for a large fraction of the observed total 511 keV emission of the bulge. A closer study of the subject reveals new layers of complexity, since positrons may propagate far away from their production sites, making it difficult to infer the underlying source distribution from the observed map of 511 keV emission. However, contrary to the rather well understood propagation of high energy (>GeV) particles of Galactic cosmic rays, understanding the propagation of low energy (~MeV) positrons in the turbulent, magnetized interstellar medium, still remains a formidable challenge. We review the spectral and imaging properties of the observed 511 keV emission and we critically discuss candidate positron sources and models of positron propagation in the Galaxy.

We model strong forward shocks in young supernova remnants with efficient particle acceleration where a nonresonant instability driven by the cosmic ray current amplifies magnetic turbulence in the shock precursor. Particle injection, magnetic field amplification (MFA) and the nonlinear feedback of particles and fields on the bulk flow are derived consistently. The shock structure depends critically on the efficiency of turbulence cascading. If cascading is suppressed, MFA is strong, the shock precursor is stratified, and the turbulence spectrum contains several discrete peaks. These peaks, as well as the amount of MFA, should influence synchrotron X-rays, allowing observational tests of cascading and other assumptions intrinsic to the nonlinear model of nonresonant wave growth.

Synchrotron X-ray emission components were recently detected in many young supernova remnants (SNRs). There is even an emerging class - SN1006, RXJ1713.72-3946, Vela Jr, and others - that is dominated by non-thermal emission in X-rays, also probably of synchrotron origin. Such emission results from electrons/positrons accelerated well above TeV energies in the spectral cut-off regime. In the case of diffusive shock acceleration, which is the most promising acceleration mechanism in SNRs, very strong magnetic fluctuations with amplitudes well above the mean magnetic field must be present. Starting from such a fluctuating field, we have simulated images of polarized X-ray emission of SNR shells and show that these are highly clumpy with high polarizations up to 50%. Another distinct characteristic of this emission is the strong intermittency, resulting from the fluctuating field amplifications. The details of this "twinkling" polarized X-ray emission of SNRs depend strongly on the magnetic-field fluctuation spectra, providing a potentially sensitive diagnostic tool. We demonstrate that the predicted characteristics can be studied with instruments that are currently being considered. These can give unique information on magnetic-field characteristics and high-energy particle acceleration in SNRs.

Thermal X-ray spectra of clusters of galaxies and other sources are commonly calculated assuming Maxwellian electron distributions. There are situations where this approximation is not valid, for instance near interfaces of hot and cold gas and near shocks. The presence of non-thermal electrons affects the X-ray spectrum. To study the role of these electrons in clusters and other environments, an efficient algorithm to calculate the X-ray spectra is needed. We approximate an arbitrary electron distribution by the sum of Maxwellian components. The decomposition is done using either a genetic algorithm or an analytical approximation. The X-ray spectrum is then evaluated using a linear combination of those Maxwellian components. Our method is fast and leads to an accurate evaluation of the spectrum. The use of Maxwellian components allows to use the standard collisional rates that are available in plasma codes such as SPEX. We give an example of a spectrum for the supra-thermal electron distribution behind a shock in a cluster of galaxies. The relative intensities of the satellite lines in such a spectrum are sensitive to the presence of the supra-thermal electrons. These lines can only be investigated with high spectral resolution. We show that the instruments on future missions like Astro-H and IXO will be able to demonstrate the presence or absence of these supra-thermal electrons.

We propose a nonlinear self-consistent model of the turbulent non-resonant particle acceleration in solar flares. We simulate temporal evolution of the spectra of charged particles accelerated by strong long-wavelength MHD turbulence taking into account back reaction of the accelerated particles on the turbulence. The main finding is that the nonlinear coupling of accelerated particles and MHD turbulence result in prominent evolution of the spectra of accelerated particles, which can be either soft-hard-soft or soft-hard-harder depending on the particle injection efficiency. Such evolution patterns are widely observed in hard X-ray and gamma-ray emission from solar flares.

Non-thermal X-ray emission in some supernova remnants originates from synchrotron radiation of ultra-relativistic particles in turbulent magnetic fields. We address the effect of a random magnetic field on synchrotron emission images and spectra. A random magnetic field is simulated to construct synchrotron emission maps of a source with a steady distribution of ultra-relativistic electrons. Non-steady localized structures (dots, clumps and filaments), in which the magnetic field reaches exceptionally high values, typically arise in the random field sample. These magnetic field concentrations dominate the synchrotron emission (integrated along the line of sight) from the highest energy electrons in the cut-off regime of the distribution, resulting in an evolving, intermittent, clumpy appearance. The simulated structures resemble those observed in X-ray images of some young supernova remnants. The lifetime of X-ray clumps can be short enough to be consistent with that observed even in the case of a steady particle distribution. The efficiency of synchrotron radiation from the cut-off regime in the electron spectrum is strongly enhanced in a turbulent field compared to emission from a uniform field of the same magnitude.

This is a report on the findings of the SNR/cosmic-ray working group for the white paper on the status and future of ground-based gamma-ray astronomy. The white paper is an APS commissioned document, and the overall version has also been released and can be found on astro-ph. This detailed section of the white paper discusses the status of past and current attempts to observe shell-type supernova remnants and diffuse emission from cosmic rays at GeV-TeV energies. We concentrate on the potential of future ground-based gamma-ray experiments to study the acceleration of relativistic charged particles which is one of the main unsolved, yet fundamental, problems in modern astrophysics. The acceleration of particles relies on interactions between energetic particles and magnetic turbulence. In the case of SNRs we can perform spatially resolved studies in systems with known geometry, and the plasma physics deduced from these observations will help us to understand other systems where rapid particle acceleration is believed to occur and where observations as detailed as those of SNRs are not possible.

The highly amplified magnetic fields suggested by observations of some supernova remnant (SNR) shells are most likely an intrinsic part of efficient particle acceleration by shocks. This strong turbulence, which may result from cosmic ray driven instabilities, both resonant and non-resonant, in the shock precursor, is certain to play a critical role in self-consistent, nonlinear models of strong, cosmic ray modified shocks. Here we present a Monte Carlo model of nonlinear diffusive shock acceleration (DSA) accounting for magnetic field amplification through resonant instabilities induced by accelerated particles, and including the effects of dissipation of turbulence upstream of a shock and the subsequent precursor plasma heating. Feedback effects between the plasma heating due to turbulence dissipation and particle injection are strong, adding to the nonlinear nature of efficient DSA. Describing the turbulence damping in a parameterized way, we reach two important results: first, for conditions typical of supernova remnant shocks, even a small amount of dissipated turbulence energy (~10%) is sufficient to significantly heat the precursor plasma, and second, the heating upstream of the shock leads to an increase in the injection of thermal particles at the subshock by a factor of several. In our results, the response of the non-linear shock structure to the boost in particle injection prevented the efficiency of particle acceleration and magnetic field amplification from increasing. We argue, however, that more advanced (possibly, non-resonant) models of turbulence generation and dissipation may lead to a scenario in which particle injection boost due to turbulence dissipation results in more efficient acceleration and even stronger amplified magnetic fields than without the dissipation.

Superbubbles (SBs) are amongst the greatest injectors of energy into the Galaxy, and have been proposed to be the acceleration site of Galactic cosmic rays. They are thought to be powered by the fast stellar winds and powerful supernova explosions of massive stars in dense stellar clusters and associations. Observations of the SB 'DEM L192' in the neighboring Large Magellenic Cloud (LMC) galaxy show that it contains only about one-third the energy injected by its constituent stars via fast stellar winds and supernovae. It is not yet understood where the excess energy is going, thus, the so-called 'energy crisis'. We show here that it is very likely that a significant fraction of the unaccounted for energy is being taken up in accelerating cosmic rays, thus bolstering the argument for the SB origin of cosmic rays.

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 review the possible mechanisms for production of non-thermal electrons which are responsible for non-thermal radiation in clusters of galaxies. Our primary focus is on non-thermal Bremsstrahlung and inverse Compton scattering, that produce hard X-ray emission. We briefly review acceleration mechanisms and point out that in most astrophysical situations, and in particular for the intracluster medium, shocks, turbulence and plasma waves play a crucial role. We consider two scenarios for production of non-thermal radiation. The first is hard X-ray emission due to non-thermal Bremsstrahlung by nonrelativistic particles. Non-thermal tails are produced by accelerating electrons from the background plasma with an initial Maxwellian distribution. However, these tails are accompanied by significant heating and they are present for a short time of <10^6 yr, which is also the time that the tail will be thermalised. Such non-thermal tails, even if possible, can only explain the hard X-ray but not the radio emission which needs GeV or higher energy electrons. For these and for production of hard X-rays by the inverse Compton model, we need the second scenario where there is injection and subsequent acceleration of relativistic electrons. It is shown that a steady state situation, for example arising from secondary electrons produced from cosmic ray proton scattering by background protons, will most likely lead to flatter than required electron spectra or it requires a short escape time of the electrons from the cluster. An episodic injection of relativistic electrons, presumably from galaxies or AGN, and/or episodic generation of turbulence and shocks by mergers can result in an electron spectrum consistent with observations but for only a short period of less than one billion years.

We briefly review capabilities and requirements for future instrumentation in UV- and X-ray astronomy that can contribute to advancing our understanding of the diffuse, highly ionised intergalactic medium.

Non-thermal components are key ingredients for understanding clusters of galaxies. In the hierarchical model of structure formation, shocks and large-scale turbulence are unavoidable in the cluster formation processes. Understanding the amplification and evolution of the magnetic field in galaxy clusters is necessary for modelling both the heat transport and the dissipative processes in the hot intra-cluster plasma. The acceleration, transport and interactions of non-thermal energetic particles are essential for modelling the observed emissions. Therefore, the inclusion of the non-thermal components will be mandatory for simulating accurately the global dynamical processes in clusters. In this review, we summarise the results obtained with the simulations of the formation of galaxy clusters which address the issues of shocks, magnetic field, cosmic ray particles and turbulence.

Modern cosmological observations allow us to study in great detail the evolution and history of the large scale structure hierarchy. The fundamental problem of accurate constraints on the cosmological parameters, within a given cosmological model, requires precise modelling of the observed structure. In this paper we briefly review the current most effective techniques of large scale structure simulations, emphasising both their advantages and shortcomings. Starting with basics of the direct N-body simulations appropriate to modelling cold dark matter evolution, we then discuss the direct-sum technique GRAPE, particle-mesh (PM) and hybrid methods, combining the PM and the tree algorithms. Simulations of baryonic matter in the Universe often use hydrodynamic codes based on both particle methods that discretise mass, and grid-based methods. We briefly describe Eulerian grid methods, and also some variants of Lagrangian smoothed particle hydrodynamics (SPH) methods.

In this paper we review the possible radiation mechanisms for the observed non-thermal emission in clusters of galaxies, with a primary focus on the radio and hard X-ray emission. We show that the difficulty with the non-thermal, non-relativistic Bremsstrahlung model for the hard X-ray emission, first pointed out by Petrosian (2001) using a cold target approximation, is somewhat alleviated when one treats the problem more exactly by including the fact that the background plasma particle energies are on average a factor of 10 below the energy of the non-thermal particles. This increases the lifetime of the non-thermal particles, and as a result decreases the extreme energy requirement, but at most by a factor of three. We then review the synchrotron and so-called inverse Compton emission by relativistic electrons, which when compared with observations can constrain the value of the magnetic field and energy of relativistic electrons. This model requires a low value of the magnetic field which is far from the equipartition value. We briefly review the possibilities of gamma-ray emission and prospects for GLAST observations. We also present a toy model of the non-thermal electron spectra that are produced by the acceleration mechanisms discussed in an accompanying paper.

The Warm-Hot Intergalactic Medium (WHIM) is thought to contribute about 40-50 % to the baryonic budget at the present evolution stage of the universe. The observed large scale structure is likely to be due to gravitational growth of density fluctuations in the post-inflation era. The evolving cosmic web is governed by non-linear gravitational growth of the initially weak density fluctuations in the dark energy dominated cosmology. Non-linear structure formation, accretion and merging processes, star forming and AGN activity produce gas shocks in the WHIM. Shock waves are converting a fraction of the gravitation power to thermal and non-thermal emission of baryonic/leptonic matter. They provide the most likely way to power the luminous matter in the WHIM. The plasma shocks in the WHIM are expected to be collisionless. Collisionless shocks produce a highly non-equilibrium state with anisotropic temperatures and a large differences in ion and electron temperatures. We discuss the ion and electron heating by the collisionless shocks and then review the plasma processes responsible for the Coulomb equilibration and collisional ionisation equilibrium of oxygen ions in the WHIM. MHD-turbulence produced by the strong collisionless shocks could provide a sizeable non-thermal contribution to the observed Doppler parameter of the UV line spectra of the WHIM.

Large-scale structure formation, accretion and merging processes, AGN activity produce cosmological gas shocks. The shocks convert a fraction of the energy of gravitationally accelerated flows to internal energy of the gas. Being the main gas-heating agent, cosmological shocks could amplify magnetic fields and accelerate energetic particles via the multi-fluid plasma relaxation processes. We first discuss the basic properties of standard single-fluid shocks. Cosmological plasma shocks are expected to be collisionless. We then review the plasma processes responsible for the microscopic structure of collisionless shocks. A tiny fraction of the particles crossing the shock is injected into the non-thermal energetic component that could get a substantial part of the ram pressure power dissipated at the shock. The energetic particles penetrate deep into the shock upstream producing an extended shock precursor. Scaling relations for postshock ion temperature and entropy as functions of shock velocity in strong collisionless multi-fluid shocks are discussed. We show that the multi-fluid nature of collisionless shocks results in excessive gas compression, energetic particle acceleration, precursor gas heating, magnetic field amplification and non-thermal emission. Multi-fluid shocks provide a reduced gas entropy production and could also modify the observable thermodynamic scaling relations for clusters of galaxies.

We review observations of extended regions of radio emission in clusters; these include diffuse emission in `relics', and the large central regions commonly referred to as `halos'. The spectral observations, as well as Faraday rotation measurements of background and cluster radio sources, provide the main evidence for large-scale intracluster magnetic fields and significant densities of relativistic electrons. Implications from these observations on acceleration mechanisms of these electrons are reviewed, including turbulent and shock acceleration, and also the origin of some of the electrons in collisions of relativistic protons by ambient protons in the (thermal) gas. Improved knowledge of non-thermal phenomena in clusters requires more extensive and detailed radio measurements; we briefly review prospects for future observations.

Recent observations of high energy (> 20 keV) X-ray emission in a few clusters of galaxies broaden our knowledge of physical phenomena in the intracluster space. This emission is likely to be nonthermal, probably resulting from Compton scattering of relativistic electrons by the cosmic microwave background (CMB) radiation. Direct evidence for the presence of relativistic electrons in some 50 clusters comes from measurements of extended radio emission in their central regions. We briefly review the main results from observations of extended regions of radio emission, and Faraday rotation measurements of background and cluster radio sources. The main focus of the review are searches for nonthermal X-ray emission conducted with past and currently operating satellites, which yielded appreciable evidence for nonthermal emission components in the spectra of a few clusters. This evidence is clearly not unequivocal, due to substantial observational and systematic uncertainties, in addition to virtually complete lack of spatial information. If indeed the emission has its origin in Compton scattering of relativistic electrons by the CMB, then the mean magnetic field strength and density of relativistic electrons in the cluster can be directly determined. Knowledge of these basic nonthermal quantities is valuable for the detailed description of processes in intracluster gas and for the origin of magnetic fields.

We present the work of an international team at the International Space Science Institute (ISSI) in Bern that worked together to review the current observational and theoretical status of the non-virialised X-ray emission components in clusters of galaxies. The subject is important for the study of large-scale hierarchical structure formation and to shed light on the "missing baryon" problem. The topics of the team work include thermal emission and absorption from the warm-hot intergalactic medium, non-thermal X-ray emission in clusters of galaxies, physical processes and chemical enrichment of this medium and clusters of galaxies, and the relationship between all these processes. One of the main goals of the team is to write and discuss a series of review papers on this subject. These reviews are intended as introductory text and reference for scientists wishing to work actively in this field. The team consists of sixteen experts in observations, theory and numerical simulations.

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 introduce a Monte Carlo model of nonlinear diffusive shock acceleration allowing for the generation of large-amplitude magnetic turbulence. The model is the first to include strong wave generation, efficient particle acceleration to relativistic energies in nonrelativistic shocks, and thermal particle injection in an internally self-consistent manner. We find that the upstream magnetic field can be amplified by large factors and show that this amplification depends strongly on the ambient Alfven Mach number. We also show that in the nonlinear model large increases in the magnetic field do not necessarily translate into a large increase in the maximum particle momentum a particular shock can produce, a consequence of high momentum particles diffusing in the shock precursor where the large amplified field converges to the low ambient value. To deal with the field growth rate in the regime of strong fluctuations, we extend to strong turbulence a parameterization that is consistent with the resonant quasi-linear growth rate in the weak turbulence limit. We believe our parameterization spans the maximum and minimum range of the fluctuation growth and, within these limits, we show that the nonlinear shock structure, acceleration efficiency, and thermal particle injection rates depend strongly on the yet to be determined details of wave growth in strongly turbulent fields. The most direct application of our results will be to estimate magnetic fields amplified by strong cosmic-ray modified shocks in supernova remnants.

We report here X-ray imaging spectroscopy observations of the northeastern shell of the supernova remnant RCW 86 with Chandra and XMM-Newton. Along this part of the shell the dominant X-ray radiation mechanism changes from thermal to synchrotron emission. We argue that both the presence of X-ray synchrotron radiation and the width of the synchrotron emitting region suggest a locally higher shock velocity of V_s = 2700 km/s and a magnetic field of B = 24+/-5 microGauss. Moreover, we also show that a simple power law cosmic ray electron spectrum with an exponential cut-off cannot explain the broad band synchrotron emission. Instead a concave electron spectrum is needed, as predicted by non-linear shock acceleration models. Finally, we show that the derived shock velocity strengthens the case that RCW 86 is the remnant of SN 185.

We report the detection of both the 67.9 and 78.4 keV 44Sc gamma-ray lines in Cassiopeia A with the INTEGRAL IBIS/ISGRI instrument. Besides the robustness provided by spectro-imaging observations, the main improvements compared to previous measurements are a clear separation of the two 44Sc lines together with an improved significance of the detection of the hard X-ray continuum up to 100 keV. These allow us to refine the determination of the 44Ti yield and to constrain the nature of the nonthermal continuum emission. By combining COMPTEL, BeppoSAX/PDS and ISGRI measurements, we find a line flux of (2.5 +/- 0.3)*10(-5) cm(-2) s(-1) leading to a synthesized 44Ti mass of 1.6 (+0.6-0.3)*10(-4) solar mass. This high value suggests that Cas A is peculiar in comparison to other young supernova remnants, from which so far no line emission from 44Ti decay has been unambiguously detected.

Using a new long X-ray observation of the cluster of galaxies Sersic 159-03 with XMM-Newton, we derive radial temperature and abundance profiles using single- and multi-temperature models. The fits to the EPIC and RGS spectra prefer multi-temperature models especially in the core. The radial profiles of oxygen and iron measured with EPIC/RGS and the line profiles in RGS suggest that there is a dip in the O/Fe ratio in the centre of the cluster compared to its immediate surroundings. A possible explanation for the large scale metallicity distribution is that SNIa and SNII products are released in the ICM through ram-pressure stripping of in-falling galaxies. This causes a peaked metallicity distribution. In addition, SNIa in the central cD galaxy enrich mainly the centre of the cluster with iron. This excess of SNIa products is consistent with the low O/Fe ratio we detect in the centre of the cluster. We fit the abundances we obtain with yields from SNIa, SNII and Population-III stars to derive the clusters chemical evolution. We find that the measured abundance pattern does not require a Population-III star contribution. The relative contribution of the number of SNIa with respect to the total number of SNe which enrich the ICM is about 25-50%. Furthermore, we discuss the possible presence of a non-thermal component in the EPIC spectra. A potential source of this non-thermal emission can be inverse-Compton scattering between Cosmic Microwave Background (CMB) photons and relativistic electrons, which are accelerated in bow shocks associated with ram-pressure stripping of in-falling galaxies.

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.

We report INTEGRAL/IBIS gamma-ray and VLA radio observations of G0.570-0.018, a diffuse X-ray source recently discovered by ASCA and Chandra in the Galactic center region. Based on its spectrum and morphology, G0.570-0.018 has been proposed to be a very young supernova remnant. In this scenario, the presence of gamma-ray lines coming from the short-lived radioactive nucleus 44Ti as well as synchrotron radio continuum emission are expected. The first could provide informations on nucleosynthesis environments in the interior of exploding stars, the latter could probe the interaction between the supernova blast wave and the circumstellar/interstellar matter. We have not detected 44Ti lines nor any conspicuous radio feature associated with this source down to the achieved sensitivities. From the derived upper limits we set constraints on the nature of G0.570-0.018.

A Chandra observation of the isolated hard X-ray source XMMU J061804.3+222732, located in the region of apparent interaction of the supernova remnant IC 443 with a molecular cloud, resolved the complex structure of the source in a few bright clumps embedded in an extended emission of a ~ 30 arcsec size. The X-ray spectra of the clumps and the extended emission are dominated by a hard power-law component with a photon index of 1.2--1.4. In addition, we see some indications of an optically thin thermal plasma of a ~ 0.3 keV temperature. The observed X-ray morphology and spectra are consistent with those expected for an isolated supernova ejecta fragment interacting with a dense ambient medium. A possible alternative interpretation is a pulsar wind nebula associated with either IC 443 or another SNR, G189.6+3.3.

The nonthermal phenomena in clusters of galaxies are considered in the context of the hierarchical model of cosmic structure formation by accretion and merging of the dark matter (DM) substructures.Accretion and merging processes produce large-scale gas shocks. The plasma shocks are expected to be collisionless. In the course of cluster's aggregation, the shocks, being the main gas-heating agent, generate turbulent magnetic fields and accelerate energetic particles via collisionless multi-fluid plasma relaxation processes. The intracluster gas heating and entropy production rate by a collisionless shock may differ significantly from that in a single-fluid collisional shock. Simple scaling relations for postshock ion temperature and entropy as functions of shock velocity in strong collisionless multi-fluid shocks are presented. We show that the multi-fluid nature of the collisionless shocks results in high gas compression, reduced entropy production and modified sigma_v-T, M-T and L_x-T scalings. The scaling indexes estimated for a simple model of a strong accretion multi-fluid shock are generally consistent with observations. Soft X-ray and extreme ultraviolet photons dominate the emission of strong accretion shock precursors that appear as large-scale filaments. Magnetic fields, turbulence and energetic particles constitute the nonthermal components contributing into the pressure balance, energy transport and emission of clusters. Nonthermal emission of energetic particles could be a test to constrain the cluster properties.

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.

Observations indicate that most massive stars in the Galaxy appear in groups, called OB associations, where their strong wind activity generates large structures known as superbubbles, inside which the subsequent supernovae (SNe) explode, in tight space and time correlation. Acknowledging this fact, we investigate four main questions: 1) does the clustering of massive stars and SN explosions influence the particle acceleration process usually associated with SNe, and induce collective effects which would not manifest around isolated supernova remnants?; 2) does it make a difference for the general phenomenology of Galactic Cosmic Rays (GCRs), notably for their energy spectrum and composition?; 3) can this help alleviate some of the problems encountered within the standard GCR source model?; and 4) Is the link between superbubbles and energetic particles supported by observational data, and can it be further tested and constrained? We argue for a positive answer to all these questions. Theoretical, phenomenological and observational aspects are treated in separate papers. Here, we discuss the interaction of massive stellar winds and SN shocks inside superbubbles and indicate how this leads to specific acceleration effects. We also show that due to the high SN explosion rate and low diffusion coefficient, low-energy particles experience repeated shock acceleration inside superbubbles.

Long Chandra observations of 17 arcmin x 17 arcmin field in the Galactic Center (GC) region reported by Muno et. al. (2003) have revealed a population of about 2,000 hard faint X-ray sources limited by luminosities Lx > 10(31} erg/s. We show that fast moving knots (FMKs) of supernovae ejecta could comprise a sizeable fraction of the sources. Each supernova event can produce hundreds of FMKs. The presence of ~ 3 supernova remnants of age ~ 1,000 years in the GC region could provide the required number of the FMKs. Simulated X-ray spectra of the FMKs contain both thermal and nonthermal components. The nonthermal spectra are power-laws with hard photon indexes 0 < Gamma < 1.5 and prominent lines of some metals (e.g. Fe). The logN-logS distribution and the hardness ratios of the FMKs in the GC region are consistent with the Chandra findings. Spatial distribution of the FMKs in the GC region should appear as an ensemble of coherent shell and jet-like structures of a few arcmin-scale.

We present \xmm observations of hard X-ray emission from the field of IC443, a supernova remnant interacting with a molecular cloud. The hard emission from the field is dominated by 12 isolated sources having 2--10 keV flux $\ga 10^{-14}$ $\enf$. Only a fraction of the sources are expected to be extragalactic or stars on statistical grounds, while the others may be associated with the remnant. We have analyzed near-infrared K band and also DSS optical data for all of the detected sources, finding that six X-ray sources are located in a relatively small $15^\prime\times 15^\prime$ region where there is strong 2.2 $\mu {\rm m}$ infrared emission, indicating interaction with a molecular cloud. The source 1SAX J0618.0+2227, the brightest in this region (excluding the plerion), is resolved with \xmm into two sources, one of which is extended and has a hard power law spectrum photon index $\lsim$ 1.5) and shows some indications of spectral line signatures (e.g. Si), while the other is point-like and has a featureless spectrum of steeper photon index $\sim$ 2.2. Possible interpretations of some of the discrete sources in terms of interaction between the SNR and the molecular cloud are discussed.

We show that fast moving isolated fragments of a supernova ejecta composed of heavy elements should be sources of K_alpha X-ray line emission of the SN nuclear-processed products. Supersonic motion of the knots in the intercloud medium will result in a bow-shock/knot-shock structure creation. Fast nonthermal particles accelerated by Fermi mechanism in the MHD collisionless shocks diffuse through a cold metallic knot, producing the X-ray emission. We modeled the X-ray emission from a fast moving knot of a mass M ~ 10^{-3} \Msun, containing about 10^{-4} \Msun of any metal impurities like Si, S, Ar, Ca, Fe. The fast electron distribution was simulated using the kinetic description. We accounted for nonlinear effects of shock modification by the nonthermal particles pressure. The K_alpha line emission is most prominent for the knots propagating through dense molecular clouds. The bow shock should be a radiative wave with a prominent infrared and optical emission. In that case the X-ray line spectrum of an ejecta fragment is dominated by the low ionization states of the ions with the metal line luminosities reaching L_x \gsim 10^{31} erg/s. High resolution XMM and Chandra observations are able to detect the line emission from the knots at distances up to a few kpcs. The knots propagating through tenuous interstellar matter are of much fainter surface brightness but long-lived. The line spectra with higher ionization states of the ions are expected in that case. Compact dense knots could be opaque for some X-ray lines and that is important for their abundances interpretation. The ensemble of unresolved knots of galactic supernovae can contribute substantially to the iron line emission observed from the Galactic Center region and the Galactic ridge.

\xmm ~observations of the X-ray feature 1SAX J0617.1+2221 in the IC443 supernova remnant are reported.We resolve the structure of the nebula into a compact core with a hard spectrum of photon index $\gamma= 1.63^{+0.11}_{-0.10}$ in the 2--10 keV energy range. The nebula also has an extended ($\sim 8\arcmin \times 5\arcmin$) X-ray halo, much larger than the radio emission extension. The photon index softens, following a linear scaling with distance from the centroid, similar to other known X-ray plerions. The index range is compatible with synchrotron burn-off models. All the observational evidence points toward a confirmation of the plerionic nature of the nebula, as recently suggested by a \ch observation, but with characteristics more similar to "non Crab-like" plerions. We discuss the implications on the synchrotron nebula magnetic field if the $>100$ MeV emission reported by {\it CGRO EGRET} is produced by the synchrotron emission. We also constrain the thermal emission of the central object, arguing that the surface temperature should be around 0.1 keV, although other possible fits cannot be excluded on the base of the \xmm data

This report discusses the possibilities for observational tests of current models and ideas regarding the origin of the Galactic Cosmic Rays using facilities which are already available.

We report BeppoSAX observations of spatially resolved hard X-ray emission from IC 443, a supernova remnant interacting with a molecular cloud. The emission is shown to come from two localized features spatially correlated with bright molecular emission regions. Both hard X-ray features have soft X-ray counterparts, in one case shifted by ~ 2 arcmin toward the remnant interior. The spectra of X-ray emission from both isolated features have photon index ~2 in the MECS regime. The emission detected from the remnant with PDS detector extends up to 100 keV. We discuss the observed properties of the hard X-ray features in relation to non-thermal emission from shocked molecular clumps and pulsar wind nebula.

We study the effect of random matter density fluctuations in the sun on resonant neutrino conversion in matter by solving numerically the evolution equation for the neutrino system including the full effect of the random matter density fluctuations of given amplitude and correlation length. In order to establish the possible effect on the MSW solutions to the solar neutrino problem we perform a global analysis of all the existing observables including the measured total rates as well as the Super-Kamiokande measurement on the time dependence of the event rates during the day and night and the recoil electron energy spectrum. We find the effects of random noise to be larger for small mixing angles and they are mostly important for correlation lengths in the range few 100 km $\lesssim L_0\lesssim$ few 1000 km. They can be understood as due to a parametric resonance occuring when the phase acquired by the oscillating neutrino state on one fluctuation length L_0 is a multiple of 2$\pi$. We find that this resonant parametric condition is mainly achieved for low energy neutrinos such as the pp-neutrinos and therefore its effect is mostly seen on the total event rates while the other Super-Kamiokande observables are very marginally sensitive to the presence of noise due to the higher energy threshold.

In evolved supernova remnants (SNRs) interacting with molecular clouds, such as IC 443, W44, and 3C391, a highly inhomogeneous structure consisting of a forward shock of moderate Mach number, a cooling layer, a dense radiative shell and an interior region filled with hot tenuous plasma is expected. We present a kinetic model of nonthermal electron injection, acceleration and propagation in that environment and find that these SNRs are efficient electron accelerators and sources of hard X- and gamma-ray emission. The energy spectrum of the nonthermal electrons is shaped by the joint action of first and second order Fermi acceleration in a turbulent plasma with substantial Coulomb losses. Bremsstrahlung, synchrotron, and inverse Compton radiation of the nonthermal electrons produce multiwavelength photon spectra in quantitative agreement with the radio and the hard emission observed by ASCA and EGRET from IC 443. We distinguish interclump shock wave emission from molecular clump shock wave emission accounting for a complex structure of molecular cloud. Spatially resolved X- and gamma- ray spectra from the supernova remnants IC 443, W44, and 3C391 as might be observed with BeppoSAX, Chandra XRO, XMM, INTEGRAL and GLAST would distinguish the contribution of the energetic lepton component to the gamma-rays observed by EGRET.

In energetic nonthermal sources such as gamma-ray bursts, AGN or galactic jet sources, etc., one expects both relativistic and transrelativistic shocks acompanied by violent motions of moderately relativistic plasma. We present general considerations indicating that these sites are electron and positron accelerators leading to a modified power law spectrum. The electron (or $e^\pm$) energy index is very hard, $\propto \gamma^{-1}$ or flatter up to a comoving frame break energy $\gamma_\ast$, and becomes steeper above that. In the example of gamma-ray bursts the Lorentz factor reaches $\gamma_\ast\sim 10^3$ for $e^{\pm}$ accelerated by the internal shock ensemble on subhydrodynamical time scales. For pairs accelerated on hydrodynamical timescales in the external shocks similarly hard spectra are obtained, and the break Lorentz factor can be as high as $\gamma_\star \siml 10^5$. Radiation from the nonthermal electrons produces photon spectra with shape and characteristic energies in qualitative agreement with observed generic gamma-ray burst and blazar spectra. The scenario described here provides a plausible way to solve one of the crucial problems of nonthermal high energy sources, namely the efficient transfer of energy from the proton flow to an apropriate nonthermal lepton component.