Continuous plasma coherent emission is maintained by repetitive Langmuir collapse driven by the nonlinear evolution of a strong electron two-stream instability. The Langmuir waves are modulated by solitary waves in the linear stage, and by electrostatic whistler waves in the nonlinear stage. Modulational instability leads to Langmuir collapse and electron heating that fills in cavitons. The high pressure is released via excitation of a short wavelength ion acoustic mode that is damped by electrons and that re-excites small-scale Langmuir waves---this process closes a feedback loop that maintains the continuous coherent emission.

The unexpected energy spectrum of the positron/electron ratio is interpreted astrophysically, with a possible exception of the 100-300 GeV range. The data indicate that this ratio, after a decline between $0.5-8$ GeV, rises steadily with a trend towards saturation at 200-400GeV. These observations (except for the trend) appear to be in conflict with the diffusive shock acceleration (DSA) mechanism, operating in a \emph{single} supernova remnant (SNR) shock. We argue that $e^{+}/e^{-}$ ratio can still be explained by the DSA if positrons are accelerated in a \emph{subset} of SNR shocks which: (i) propagate in clumpy gas media, and (ii) are modified by accelerated CR \emph{protons}. The protons penetrate into the dense gas clumps upstream to produce positrons and, \emph{charge the clumps positively}. The induced electric field expels positrons into the upstream plasma where they are shock-accelerated. Since the shock is modified, these positrons develop a harder spectrum than that of the CR electrons accelerated in other SNRs. Mixing these populations explains the increase in the $e^{+}/e^{-}$ ratio at $E>8$ GeV. It decreases at $E<8$ GeV because of a subshock weakening which also results from the shock modification. Contrary to the expelled positrons, most of the antiprotons, electrons, and heavier nuclei, are left unaccelerated inside the clumps. Scenarios for the 100-300 GeV AMS-02 fraction exceeding the model prediction, including, but not limited to, possible dark matter contribution, are also discussed.

Cosmic rays (CR), constrained by scattering on magnetic irregularities, are believed to propagate diffusively. But a well-known defect of diffusive approximation whereby some of the particles propagate unrealistically fast has directed interest towards an alternative CR transport model based on the "telegraph" equation. However, its derivations often lack rigor and transparency leading to inconsistent results. We apply the classic Chapman-Enskog method to the CR transport problem. We show that no "telegraph" (second order time derivative) term emerges in any order of a proper asymptotic expansion with systematically eliminated short time scales. Nevertheless, this term may formally be\emph{ converted} from the \emph{fourth} order hyper-diffusive term of the expansion. But, both the telegraph and hyperdiffusive terms may only be important for a short relaxation period associated with either strong pitch-angle anisotropy or spatial inhomogeneity of the initial CR distribution. Beyond this period the system evolves diffusively in both cases. The term conversion, that makes the telegraph and Chapman-Enskog approaches reasonably equivalent, is possible only after this relaxation period. During this period, the telegraph solution is argued to be unphysical. Unlike the hyperdiffusion correction, it is not uniformly valid and introduces implausible singular components to the solution. These dominate the solution during the relaxation period. As they are shown not to be inherent in the underlying scattering problem, we argue that the telegraph term is involuntarily acquired in an asymptotic reduction of the problem.

Supernova remnants (SNRs), as the major contributors to the galactic cosmic rays (CR), are believed to maintain an average CR spectrum by diffusive shock acceleration (DSA) regardless of the way they release CRs into the interstellar medium (ISM). However, the interaction of the CRs with nearby gas clouds crucially depends on the release mechanism. We call into question two aspects of a popular paradigm of the CR injection into the ISM, according to which they passively and isotropically diffuse in the prescribed magnetic fluctuations as test particles. First, we treat the escaping CR and the Alfven waves excited by them on an equal footing. Second, we adopt field aligned CR escape outside the source, where the waves become weak. An exact analytic self-similar solution for a CR "cloud" released by a dimmed accelerator strongly deviates from the test-particle result. The CR diffusion coefficient $D_{NL}$ is strongly suppressed compared to its background ISM value $D_{ISM}$: $D_{NL}\sim D_{ISM}\exp(-\Pi)<< D_{ISM}$ for sufficiently high field-line-integrated CR partial pressure, $\Pi$. When $\Pi>>1$, the CRs drive Alfven waves efficiently enough to build a transport barrier that strongly reduces the leakage. The solution has a spectral break at $p=p_{br}$, where $p_{br}$ satisfies the following equation $D_{NL}(p_{br})\simeq z^{2}/t$.

The proof of cosmic ray (CR) origin in supernova remnants (SNR) must hinge on full consistency of the CR acceleration theory with the observations; direct proof is impossible because of the orbit stochasticity of CR particles. Recent observations of a number of galactic SNR strongly support the SNR-CR connection in general and the Fermi mechanism of CR acceleration, in particular. However, many SNR expand into weakly ionized dense gases, and so a significant revision of the mechanism is required to fit the data. We argue that strong ion-neutral collisions in the remnant surrounding lead to the steepening of the energy spectrum of accelerated particles by \emph{exactly one power}. The spectral break is caused by a partial evanescence of Alfven waves that confine particles to the accelerator. The gamma-ray spectrum generated in collisions of the accelerated protons with the ambient gas is also calculated. Using the recent Fermi spacecraft observation of the SNR W44 as an example, we demonstrate that the parent proton spectrum is a classical test particle power law $\propto E^{-2}$, steepening to $E^{-3}$ at $E_{br}\approx7GeV$.

The much-anticipated proof of cosmic ray (CR) acceleration in supernova remnants (SNR) must hinge on full consistency of acceleration theory with the observations; direct proof is impossible because of the orbit scrambling of CR particles. The PAMELA orbital telescope revealed deviation between helium and proton CR spectra deemed inconsistent with the theory, since the latter does not differentiate between elements of ultrarelativistic rigidity. By considering an initial (injection-) phase of the diffusive shock acceleration (DSA), where elemental similarity does not apply, we demonstrate that the spectral difference is, in fact, a unique signature of the DSA. Collisionless plasma SNR shocks inject more He2+ relative to protons when they are stronger and so produce harder helium spectra. The injection bias is due to Alfven waves driven by the more abundant protons, so the He2+ ions are harder to trap by these waves because of the larger gyroradii. By fitting the p/He ratio to the PAMELA data, we bolster the DSA-case for resolving the century-old mystery of CR origin.

In shock precursors populated by accelerated cosmic rays (CR), the CR return current instability is believed to significantly enhance the pre-shock perturbations of magnetic field. We have obtained fully-nonlinear exact ideal MHD solutions supported by the CR return current. The solutions occur as localized spikes of circularly polarized Alfven envelopes (solitons, or breathers). As the conventional (undriven) solitons, the obtained magnetic spikes propagate at a speed $C$ proportional to their amplitude, $C=C_{A}B_{{\rm max}}/\sqrt{2}B_{0}$. The sufficiently strong solitons run thus ahead of the main shock and stand in the precursor, being supported by the return current. This property of the nonlinear solutions is strikingly different from the linear theory that predicts non-propagating (that is, convected downstream) circularly polarized waves. The nonlinear solutions may come either in isolated pulses (solitons) or in soliton-trains (cnoidal waves). The morphological similarity of such quasi-periodic soliton chains with recently observed X-ray stripes in Tycho supernova remnant (SNR) is briefly discussed. The magnetic field amplification determined by the suggested saturation process is obtained as a function of decreasing SNR blast wave velocity during its evolution from the ejecta-dominated to the Sedov-Taylor stage.

A mechanism for proton acceleration to ~10^21eV is suggested. It may operate in accretion flows onto thin dark matter filaments of cosmic structure formation. The flow compresses the ambient magnetic field to strongly increase and align it with the filament. Particles begin the acceleration by the ExB drift with the accretion flow. The energy gain in the drift regime is limited by the conservation of the adiabatic invariant p_perp^2/B. Upon approaching the filament, the drift turns into the gyro-motion around the filament so that the particle moves parallel to the azimuthal electric field. In this 'betatron' regime the acceleration speeds up to rapidly reach the electrodynamic limit $cp_{max}=eBR$ for an accelerator with magnetic field $B$ and the orbit radius $R$ (Larmor radius). The periodic orbit becomes unstable and the particle slings out of the filament to the region of a weak (uncompressed) magnetic field, which terminates the acceleration. The mechanism requires pre-acceleration that is likely to occur in structure formation shocks upstream or nearby the filament accretion flow. Previous studies identify such shocks as efficient proton accelerators to a firm upper limit ~10^19.5 eV placed by the catastrophic photo-pion losses. The present mechanism combines explosive energy gain in its final (betatron) phase with prompt particle release from the region of strong magnetic field. It is this combination that allows protons to overcome both the photo-pion and the synchrotron-Compton losses and therefore attain energy 10^21 eV. A requirement on accelerator to reach a given E_max placed by the accelerator energy dissipation \propto E_{max}^{2}/Z_0 due to the finite vacuum impedance Z_0 is circumvented by the cyclic operation of the accelerator.

Recent observations of the supernova remnant W44 by the \emph{Fermi} spacecraft observatory strongly support the idea that the bulk of galactic cosmic rays is accelerated in such remnants by a Fermi mechanism, also known as diffusive shock acceleration. However, the W44 expands into weakly ionized dense gas, and so a significant revision of the mechanism is required. In this paper we provide the necessary modifications and demonstrate that strong ion-neutral collisions in the remnant surrounding lead to the steepening of the energy spectrum of accelerated particles by \emph{exactly one power}. The spectral break is caused by Alfven wave evanescence leading to the fractional particle losses. The gamma-ray spectrum generated in collisions of the accelerated protons with the ambient gas is also calculated and successfully fitted to the Fermi Observatory data. The parent proton spectrum is best represented by a classical test particle power law E^-2, steepening to E^-3 at E_br~7GeV due to deteriorated particle confinement.

Both the acceleration of cosmic rays (CR) in supernova remnant shocks and their subsequent propagation through the random magnetic field of the Galaxy deem to result in an almost isotropic CR spectrum. Yet the MILAGRO TeV observatory discovered a sharp ($\sim10^{\circ})$ arrival anisotropy of CR nuclei. We suggest a mechanism for producing a weak and narrow CR beam which operates en route to the observer. The key assumption is that CRs are scattered by a strongly anisotropic Alfven wave spectrum formed by the turbulent cascade across the local field direction. The strongest pitch-angle scattering occurs for particles moving almost precisely along the field line. Partly because this direction is also the direction of minimum of the large scale CR angular distribution, the enhanced scattering results in a weak but narrow particle excess. The width, the fractional excess and the maximum momentum of the beam are calculated from a systematic transport theory depending on a single scale $l$ which can be associated with the longest Alfven wave, efficiently scattering the beam. The best match to all the three characteristics of the beam is achieved at $l\sim1$pc. The distance to a possible source of the beam is estimated to be within a few 100pc. Possible approaches to determination of the scale $l$ from the characteristics of the source are discussed. Alternative scenarios of drawing the beam from the galactic CR background are considered. The beam related large scale anisotropic CR component is found to be energy independent which is also consistent with the observations.

Strong astrophysical shocks, diffusively accelerating cosmic rays (CR) ought to develop CR precursors. The length of such precursor $L_{p}$ is believed to be set by the ratio of the CR mean free path $\lambda$ to the shock speed, i.e., $L_{p}\sim c\lambda/V_{sh}\sim cr_{g}/V_{sh}$, which is formally independent of the CR pressure $P_{c}$. However, the X-ray observations of supernova remnant shocks suggest that the precursor scale may be significantly shorter than $L_{p}$ which would question the above estimate unless the magnetic field is strongly amplified and the gyroradius $r_{g}$ is strongly reduced over a short (unresolved) spatial scale. We argue that while the CR pressure builds up ahead of the shock, the acceleration enters into a strongly nonlinear phase in which an acoustic instability, driven by the CR pressure gradient, dominates other instabilities (at least in the case of low $\beta$ plasma). In this regime the precursor steepens into a strongly nonlinear front whose size scales with \emph{the CR pressure}as $L_{f}\sim L_{p}\cdot(L_{s}/L_{p})^{2}(P_{c}/P_{g})^{2}$, where $L_{s}$ is the scale of the developed acoustic turbulence, and $P_{c}/P_{g}$ is the ratio of CR to gas pressure. Since $L_{s}\ll L_{p}$, the precursor scale reduction may be strong in the case of even a moderate gas heating by the CRs through the acoustic and (possibly also) the other instabilities driven by the CRs.

Many stormy events in astrophysics occur due to the sudden magnetic energy release. This is possible if a magnetic configuration abruptly changes its topology, an event usually referred to as magnetic reconnection. It is known that pure Ohmic decay is inefficient, occurring during cosmological times (due to the huge characteristic scales $L$). It is recognized that the presence of current sheets speeds up the process, but still insufficiently$^{1,2,3,4,5}$. We show that, in highly compressible and substantially gravitational media, the reconnection is fast enough to account for stormy events. Thus, highly compressible situations offer exiting opportunities in explanations of violent events, although full-scale compressible and gravitational simulations proved to be quite challenging.

Numerical simulations of strongly compressible MHD corresponding to a stellar atmosphere with substantial gravity and near force-free magnetic fields show that the current sheet collapses (its width decreasing substantially). As a result, the reconnection rate increases dramatically.

Special features of magnetohydrodynamic waves linear dynamics in smooth shear flows are studied. Quantitative asymptotic and numerical analysis are performed for wide range of system parameters when basic flow has constant shear of velocity and uniform magnetic field is parallel to the basic flow. The special features consist of magnetohydrodynamic wave mutual transformation and over-reflection phenomena. The transformation takes place for arbitrary shear rates and involves all magnetohydrodynamic wave modes. While the over-reflection occurs only for slow magnetosonic and Alfv\'en waves at high shear rates. Studied phenomena should be decisive in the elaboration of the self-sustaining model of magnetohydrodynamic turbulence in the shear flows.

Cosmic rays (CRs) are thought to be accelerated in SNRs. The most favorable situation for proving that the main, hadronic CR component is accelerated there is when CRs interact with dense gases, such as molecular clouds (MC) which surround the SN shock. Here, a new mechanism of spectrum formation in partially ionized gases near SNRs is proposed. Using an analytic model of nonlinear diffusive shock acceleration, we calculate the spectra of protons and estimate the resulting $\gamma$-ray emission occurring when the SNR shock approaches a MC. We show that the spectrum develops a break in the TeV range and that its GeV component is suppressed. These modifications to the standard theory occur because of the proximity of the partially ionized MC-gas and because of the physics of particle and Alfven wave propagation inside the gas. Possible applications of the new spectra to the recent CANGAROO and HESS observations of the SNR RX J1713.7-3946 are discussed.