A possibility of DM being multicomponent has a strong implication on resolving decades-long known cosmological problems on small scale. In addition to elastic scattering, the model allows for inelastic interactions, which can be characterized by a 'velocity kick' parameter. The simplest 2cDM model with cross section $0.01\lesssim\sigma/m<1\textrm{ cm}^{2}{ \rm g}^{-1}$ and the kick velocity $V_{k}\simeq 100\textrm{ km s}^{-1}$ has been shown to robustly resolve the missing satellites, core-cusp, and too-big-to-fail problems in $N$-body cosmological simulations tested on MW-like haloes of a virial mass $\sim5 \times 10^{11}$ M$_{\odot}$ (Paper I $\&$ II). With the aim of further constraining the parameter space available for the 2cDM model, we extend our analysis to dwarf and galaxy cluster haloes with their virial mass of $\sim 10^7 - 10^8$ and $\sim 10^{13} - 10^{14}$ M$_{\odot}$, respectively. We find $\sigma_{0} / m \gtrsim 0.1 \textrm{ cm}^{2}{\rm g}^{-1}$ is preferentially disfavored for both dwarfs and galaxy cluster haloes in comparison with observations, while $\sigma_{0} / m = 0.001 \textrm{ cm}^{2}{\rm g}^{-1}$ causes little perceptible difference from that of the CDM counterpart for most of the cross section's velocity dependence studied in this work. Our main result is that within the reasonable set of parameters the 2cDM model can successfully explain the observational trends seen in dwarf galaxy and galaxy cluster haloes and the model leaves us an open window for other possible alternative DM models.

This is a science white paper submitted to the Astro-2020 and Plasma-2020 Decadal Surveys. The paper describes the present status and emerging opportunities in Extreme Plasma Astrophysics -- a study of astrophysically-relevant plasma processes taking place under extreme conditions that necessitate taking into account relativistic, radiation, and QED effects.

Multicomponent dark matter with self-interactions, which allows for inter-conversions of species into one another, is a promising paradigm that is known to successfully and simultaneously resolve major problems of the conventional $\Lambda$CDM cosmology at galactic and sub-galactic scales. In this paper, we present $N$-body simulations of the simplest two-component (2cDM) model aimed at studying the distribution of dark matter halos with masses $M\lesssim10^{12}M_\odot$. In particular, we investigate how the maximum circular velocity function of the halos is affected by the velocity dependence of the self-interaction cross-sections, $\sigma(v)\propto v^a$, and compare them with available observational data. The results demonstrate that the 2cDM paradigm with the range of self-interaction cross-section per particle mass (evaluated at $v=100$ km s$^{-1}$) of $0.01\lesssim \sigma_0/m\lesssim 1 $ cm$^2$g$^{-1}$ and the mass degeneracy $\Delta m/m\sim 10^{-7}-10^{-8}$ is robustly resolving the substructure and too-big-to-fail problems by suppressing the substructure having small maximum circular velocities, $V_{\rm max}\lesssim100$ km s$^{-1}$. We also discuss the disagreement between the radial distribution of dwarfs in a host halo observed in the Local Group and simulated with CDM. This can be considered as one more small-scale problem of CDM. We demonstrate that such a disagreement is alleviated in 2cDM. Finally, the computed matter power-spectra of the 2cDM structure indicate the model's consistency with the existing Ly-$\alpha$ forest constraints.

The multicomponent dark matter model with self-scattering and inter-conversions of species into one another is an alternative dark matter paradigm that is capable of resolving the long-standing problems of $\Lambda$CDM cosmology at small scales. In this paper, we have studied in detail the properties of dark matter halos with $M \sim 4-5 \times10^{11} M_{\odot}$ obtained in $N$-body cosmological simulations with the simplest two-component (2cDM) model. A large set of velocity-dependent cross-section prescriptions for elastic scattering and mass conversions, $\sigma_s(v)\propto v^{a_s}$ and $\sigma_c(v)\propto v^{a_c}$, has been explored and the results were compared with observational data. The results demonstrate that self-interactions with the cross-section per particle mass evaluated at $v=100$ km s$^{-1}$ being in the range of $0.01\lesssim \sigma_0/m\lesssim 1$ cm$^2$g$^{-1}$ robustly suppress central cusps, thus resolving the core-cusp problem. The core radii are controlled by the values of $\sigma_0/m$ and the DM cross-section's velocity-dependent power-law indices $(a_s,a_c)$, but are largely insensitive to the species' mass degeneracy. These values are in full agreement with those resolving the substructure and too-big-to-fail problems. We have also studied the evolution of halos in the 2cDM model with cosmic time.

Existing theoretical and observational constraints on the abundance of magnetic monopoles are limited. Here we demonstrate that an ensemble of monopoles forms a plasma whose properties are well determined and whose collective effects place new tight constraints on the cosmological abundance of monopoles. In particular, the existence of micro-Gauss magnetic fields in galaxy clusters and radio relics implies that the scales of these structures are below the Debye screening length, thus setting an upper limit on the cosmological density parameter of monopoles, $\Omega_M\lesssim3\times10^{-4}$, which precludes them from being the dark matter. Future detection of Gpc-scale coherent magnetic fields could improve this limit by a few orders of magnitude. In addition, we predict the existence of magnetic Langmuir waves and turbulence which may appear on the sky as "zebra patterns" of an alternating magnetic field with ${\bf k\cdot B}\not=0$. We also show that magnetic monopole Langmuir turbulence excited near the accretion shock of galaxy clusters may be an efficient mechanism for generating the observed intracluster magnetic fields.

We quantitatively study the stationary, axisymmetric, force-free magnetospheres of spinning (Kerr) black holes (BHs) and the conditions needed for relativistic jets to be powered by the Blandford-Znajek mechanism. These jets could be from active galactic nuclei, blazars, quasars, micro-quasars, radio active galaxies, and other systems that host Kerr BHs. The structure of the magnetosphere determines how the BH energy is extracted, e.g., via Blandford-Znajek mechanism, which converts the BH rotational energy into Poynting flux. The key assumption is the force-free condition, which requires the presence of plasma with the density being above the Goldreich-Julian density. Unlike neutron stars, which in principle can supply electrons from the surface, BH cannot supply plasma at all. The plasma must be generated \em{in situ} via an electron-positron cascade, presumably in the gap region. Here we study varying conditions that provide a sufficient amount of plasma for the Blandford-Znajek mechanism to work effectively.

Astrophysical and high-energy-density laboratory plasmas often have large-amplitude, sub-Larmor-scale electromagnetic fluctuations excited by various kinetic-streaming or anisotropy-driven instabilities. The Weibel (or the filamentation) instability is particularly important because it can rapidly generate strong magnetic fields, even in the absence of seed fields. Particles propagating in collisionless plasmas with such small-scale magnetic fields undergo stochastic deflections similar to Coulomb collisions, with the magnetic pitch-angle diffusion coefficient representing the effective "collision" frequency. We show that this effect of the plasma "quasi-collisionality" can strongly affect the growth rate and evolution of the Weibel instability in the deeply nonlinear regime. This result is especially important for understanding cosmic-ray-driven turbulence in an upstream region of a collisionless shock of a gamma-ray burst or a supernova. We demonstrate that the quasi-collisions caused by the fields generated in the upstream suppress the instability slightly but can never shut it down completely. This confirms the assumptions made in the self-similar model of the collisionless foreshock.

Understanding of the observed structure in the universe can be reached only in the theoretical framework of dark matter. N-body simulations are indispensable for the analysis of the formation and evolution of the dark matter web. Two primary fields - density and velocity fields - are used in most of studies. However dark matter provides two additional fields which are unique for collisionless media only. These are the multi- stream field in Eulerian space and flip-flop field in Lagrangian space. The flip-flop field represents the number of sign reversals of an elementary volume of each collisionless fluid element. This field can be estimated by counting the sign reversals of the Jacobian at each particle at every time step of the simulation. The Jacobian is evaluated by numerical differentiation of the Lagrangian submanifold, i.e., the three-dimensional dark matter sheet in the six-dimensional space formed by three Lagrangian and three Eulerian coordinates. We present the results of the statistical study of the evolution of the flip-flop field from z = 50 to the present time z = 0. A number of statistical characteristics show that the pattern of the flip-flop field remains remarkably stable from z = 30 to the present time. As a result the flip-flop field evaluated at z = 0 stores a wealth of information about the dynamical history of the dark matter web. In particular one of the most intriguing properties of the flip-flop is a unique capability to preserve the information about the merging history of dark matter haloes.

High-amplitude, chaotic/turbulent electromagnetic fluctuations are ubiquitous in high-energy-density laboratory and astrophysical plasmas, where they can be excited by various kinetic-streaming and/or anisotropy-driven instabilities, such as the Weibel instability. These fields typically exist on "sub-Larmor scales" -- scales smaller than the electron Larmor radius. Electrons moving through such magnetic fields undergo small-angle stochastic deflections of their pitch-angles, thus establishing diffusive transport on long time-scales. We show that this behavior, under certain conditions, is equivalent to Coulomb collisions in collisional plasmas. The magnetic pitch-angle diffusion coefficient, which acts as an effective "collision" frequency, may be substantial in these, otherwise, collisionless environments. We show that this effect, colloquially referred to as the plasma "quasicollisionality", may radically alter the expected radiative transport properties of candidate plasmas. We argue that the modified magneto-optic effects in these plasmas provide an attractive, novel radiative diagnostic tool for the exploration and characterization of small-scale magnetic turbulence, as well as affect inertial confinement fusion and other laser-plasma experiments.

Cosmic ray propagation is diffusive because of pitch angle scattering by waves. We demonstrate that if the high-amplitude magnetohydrodynamic turbulence with $\tilde B/\langle B\rangle \sim 1$ is present on top of the mean field gradient, the diffusion becomes asymmetric. As an example, we consider the vertical transport of cosmic rays in our Galaxy propagating away from a point-like source. We solve this diffusion problem analytically using a one-dimensional Markov chain analysis. We obtained that the cosmic ray density markedly differs from the standard diffusion prediction and has a sizable effect on their distribution throughout the galaxy. The equation for the continuous limit is also derived, which shows limitations of the convection-diffusion equation.

The nature of dark matter is unknown. A number of dark matter candidates are quantum flavor-mixed particles but this property has never been accounted for in cosmology. Here we explore this possibility from the first principles via extensive $N$-body cosmological simulations and demonstrate that the two-component dark matter model agrees with observational data at all scales. Substantial reduction of substructure and flattening of density profiles in the centers of dark matter halos found in simulations can simultaneously resolve several outstanding puzzles of modern cosmology. The model shares the "why now?" fine-tuning caveat pertinent to all self-interacting models. Predictions for direct and indirect detection dark matter experiments are made.

We have investigated particle acceleration and emission from shocks and shear flows associated with an unmagnetized relativistic jet plasma propagating into an unmagnetized ambient plasma. Strong electro-magnetic fields are generated in the jet shock via the filamentation (Weibel) instability. Shock field strength and structure depend on plasma composition (($e^{\pm}$ or $e^-$- $p^+$ plasmas) and Lorentz factor. In the velocity shear between jet and ambient plasmas, strong AC ($e^{\pm}$ plasmas) or DC ($e^-$- $p^+$ plasmas) magnetic fields are generated via the kinetic Kelvin-Helmholtz instability (kKHI), and the magnetic field structure also depends on the jet Lorentz factor. We have calculated, self-consistently, the radiation from electrons accelerated in shock generated magnetic fields. The spectra depend on the jet's initial Lorentz factor and temperature via the resulting particle acceleration and magnetic field generation. Our ongoing "Global" jet simulations containing shocks and velocity shears will provide us with the ability to calculate and model the complex time evolution and/or spectral structure observed from gamma-ray bursts, AGN jets, and supernova remnants.

A new computational paradigm for the analysis of substructure of the Cosmic Web in cosmological cold dark matter simulations is proposed. We introduce a new data-field --- the flip-flop field ---which carries wealth of information about the history and dynamics of the structure formation in the universe. The flip-flop field is an ordered data set in Lagrangian space representing the number of turns inside out sign reversals of an elementary volume of each collisionless fluid element represented by a computational particle in a N-body simulation. This field is computed using the Lagrangian submanifold, i.e. the three-dimensional dark matter sheet in the six-dimensional space formed by three Lagrangian and three Eulerian coordinates of the simulation particles. It is demonstrated that the very rich substructure of dark matter haloes and the void regions can be reliably and unambiguously recovered from the flip-flop field.

We have investigated magnetic field generation in velocity shears via the kinetic Kelvin-Helmholtz instability (kKHI) using a relativistic plasma jet core and stationary plasma sheath. Our three-dimensional particle-in-cell simulations consider plasma jet cores with Lorentz factors of 1.5, 5, and 15 for both electron-proton and electron-positron plasmas. For electron-proton plasmas we find generation of strong large-scale DC currents and magnetic fields which extend over the entire shear-surface and reach thicknesses of a few tens of electron skin depths. For electron-positron plasmas we find generation of alternating currents and magnetic fields. Jet and sheath plasmas are accelerated across the shear surface in the strong magnetic fields generated by the kKHI. The mixing of jet and sheath plasmas generates transverse structure similar to that produced by the Weibel instability.

Evolution of a system of interacting non-relativistic quantum flavor-mixed particles is considered both theoretically and numerically. It was shown that collisions of mixed particles not only scatter them elastically, but can also change their mass eigenstates thus affecting particles' flavor composition and kinetic energy. The mass eigenstate conversions and elastic scattering are related but different processes, hence the conversion $S$-matrix elements can be arbitrarily large even when the elastic scattering $S$-matrix elements vanish. The conversions are efficient when the mass eigenstates are well-separated in space but suppressed if their wave-packets overlap; the suppression is most severe for mass-degenerate eigenstates in flat space-time. The mass eigenstate conversions can lead to an interesting process, called `quantum evaporation,' in which mixed particles, initially confined deep inside a gravitational potential well and scattering only off each other, can escape from it without extra energy supply leaving nothing behind inside the potential at $t\to \infty$. Implications for the cosmic neutrino background and the two-component dark matter model are discussed and a prediction for the direct detection dark matter experiments is made.

We present extensive radio and X-ray observations of SN\, 2012au, the energetic radio luminous supernova of type Ib that may be a link between subsets of hydrogen-poor superluminous and normal core-collapse supernovae. The observations closely follow models of synchrotron emission from shock heated circum-burst medium that has a wind density profile ($\rho \propto r^{-2}$). We infer a sub-relativistic velocity for the shock wave $v \approx 0.2\,c$ and a radius of $r \approx 1.4 \times 10^{16} \rm cm$ at 25 days after the estimated date of explosion. For a constant wind velocity of 1000 km/s we determine the constant mass loss rate of the progenitor to be $\dot{M} = 3.6 \times 10^{-6} \rm M_{\odot} yr^{-1}$, consistent with the estimates from X-ray observations. We estimate the total internal energy of the radio emitting material to be $E \approx 10^{47} \rm erg$, which is intermediate to SN\,1998bw and SN\,2002ap. Evolution of the radio light curves of SN\,2012au is consistent with interaction with a smoothly distributed circum-burst medium and absence of stellar shells ejected from previous outbursts out to $r \approx 10^{17} \rm cm$ from the supernova site. Based on this we conclude that the evolution of the SN\,2012au progenitor star was relatively quiet during the final years preceding explosion. We find that the bright radio emission from SN2012au was not dissimilar from other core collapse supernovae despite it's extraordinary optical properties. We speculate that it was the nature of the explosion that led to the unusual demise of the SN2012au progenitor star.

We have investigated generation of magnetic fields associated with velocity shear between an unmagnetized relativistic jet and an unmagnetized sheath plasma. We have examined the strong magnetic fields generated by kinetic shear (Kelvin-Helmholtz) instabilities. Compared to the previous studies using counter-streaming performed by Alves et al. (2012), the structure of KKHI of our jet-sheath configuration is slightly different even for the global evolution of the strong transverse magnetic field. In our simulations the major components of growing modes are the electric field $E_{\rm z}$ and the magnetic field $B_{\rm y}$. After the $B_{\rm y}$ component is excited, an induced electric field $E_{\rm x}$ becomes significant. However, other field components remain small. We find that the structure and growth rate of KKHI with mass ratios $m_{\rm i}/m_{\rm e} = 1836$ and $m_{\rm i}/m_{\rm e} = 20$ are similar. In our simulations saturation in the nonlinear stage is not as clear as in counter-streaming cases. The growth rate for a mildly-relativistic jet case ($\gamma_{\rm j} = 1.5$) is larger than for a relativistic jet case ($\gamma_{\rm j} = 15$).

The relation of particle transport of relativistic particles in plasmas with high-amplitude isotropic sub-Larmor-scale magnetic turbulence to the spectra of radiation simultaneously produced by these particles is investigated both analytically and numerically. We have found that in the asymptotic regime of very small particle deflections the pitch angle diffusion coefficient is directly related to the spectrum of the emitted radiation. Moreover, this spectrum provides much information about the statistical properties of the underlying magnetic turbulence. The transition from small- to large-scale jitter to synchrotron radiation regimes as a function of turbulence properties has also been explored. These results can readily be used to diagnose laboratory and astrophysical plasmas.

We have investigated particle acceleration and shock structure associated with an unmagnetized relativistic jet propagating into an unmagnetized plasma. Strong magnetic fields generated in the trailing jet shock lead to transverse deflection and acceleration of the electrons. We have self-consistently calculated the radiation from the electrons accelerated in the turbulent magnetic fields. We find that the synthetic spectra depend on the bulk Lorentz factor of the jet, the jet temperature, and the strength of the magnetic fields generated in the shock. We have also begun study of electron acceleration in the strong magnetic fields generated by kinetic shear (Kelvin-Helmholtz) instabilities. Our calculated spectra should lead to a better understanding of the complex time evolution and/or spectral structure from gamma-ray bursts, relativistic jets, and supernova remnants.

Merging binaries of compact relativistic objects (neutron stars and black holes) are thought to be progenitors of short gamma-ray bursts and sources of gravitational waves, hence their study is of great importance for astrophysics. Because of the strong magnetic field of one or both binary members and high orbital frequencies, these binaries are strong sources of energy in the form of Poynting flux (e.g., magnetic-field-dominated outflows, relativistic leptonic winds, electromagnetic and plasma waves). The steady injection of energy by the binary forms a bubble (or a cavity) filled with matter with the relativistic equation of state, which pushes on the surrounding plasma and can drive a shock wave in it. Unlike the Sedov-von Neumann-Taylor blast wave solution for a point-like explosion, the shock wave here is continuously driven by the ever-increasing pressure inside the bubble. We calculate from the first principles the dynamics and evolution of the bubble and the shock surrounding it and predict that such systems can be observed as radio sources a few hours before and after the merger. At much later times, the shock is expected to settle onto the Sedov-von Neumann-Taylor solution, thus resembling an explosion.

Bubbles in the interstellar medium are produced by astrophysical sources, which continuously or explosively deposit large amount of energy into the ambient medium. These expanding bubbles can drive shocks in front of them, which dynamics is markedly different from the widely used Sedov-von Neumann-Taylor blast wave solution. Here we present the theory of a bubble-driven shock and show how its properties and evolution are determined by the temporal history of the source energy output, generally referred to as the source luminosity law, $L(t)$. In particular, we find the analytical solutions for a driven shock in two cases: the self-similar scaling $L\propto (t/t_s)^p$ law (with $p$ and $t_s$ being constants) and the finite activity time case, $L\propto (1-t/t_s)^{-p}$. The latter with $p>0$ describes a finite-time-singular behavior, which is relevant to a wide variety of systems with explosive-type energy release. For both luminosity laws, we derived the conditions needed for the driven shock to exist and predict the shock observational signatures. Our results can be relevant to stellar systems with strong winds, merging neutron star/magnetar/black hole systems, and massive stars evolving to supernovae explosions.

We demonstrate that a stable particle with flavor mixing, confined in a gravitational potential can gradually and irreversibly escape -- or "evaporate" -- from it. This effect is due to mass eigenstate conversions which occur in interactions (scattering) of mass states with other particles even when the energy exchange between them is vanishing. The evaporation and conversion are quantum effects not related to flavor oscillations, particle decay, quantum tunneling or other well-known processes. Apart from their profound academic interest, these effects should have tremendous implications for cosmology, e.g., (1) the cosmic neutrino background distortion is predicted and (2) the softening of central cusps in dark matter halos and smearing out or destruction of dwarf halos were suggested.

lasma instabilities excited in collisionless shocks are responsible for particle acceleration. We have investigated the particle acceleration and shock structure associated with an unmagnetized relativistic electron-positron jet propagating into an unmagnetized electron-positron plasma. Cold jet electrons are thermalized and slowed while the ambient electrons are swept up to create a partially developed hydrodynamic-like shock structure. In the leading shock, electron density increases by a factor of about 3.5 in the simulation frame. Strong electromagnetic fields are generated in the trailing shock and provide an emission site. These magnetic fields contribute to the electron's transverse deflection behind the shock. Our initial results of a jet-ambient interaction with anti-parallel magnetic fields show pile-up of magnetic fields at the colliding shock, which may lead to reconnection and associated particle acceleration. We will investigate the radiation in transient stage as a possible generation mechanism of precursors of prompt emission. In our simulations we calculate the radiation from electrons in the shock region. The detailed properties of this radiation are important for understanding the complex time evolution and spectral structure in gamma-ray bursts, relativistic jets, and supernova remnants.

Spontaneous rapid growth of strong magnetic fields is rather ubiquitous in high-energy density environments ranging from astrophysical sources (e.g., gamma-ray bursts and relativistic shocks), to reconnection, to laser-plasma interaction laboratory experiments, where they are produced by kinetic streaming instabilities of the Weibel type. Relativistic electrons propagating through these sub-Larmor-scale magnetic fields radiate in the jitter regime, in which the anisotropy of the magnetic fields and the particle distribution have a strong effect on the produced radiation. Here we develop the general theory of jitter radiation, which includes (i) anisotropic magnetic fields and electron velocity distributions, (ii) the effects of trapped electrons and (iii) extends the description to large deflection angles of radiating particles thus establishing a cross-over between the classical jitter and synchrotron regimes. Our results are in remarkable agreement with the radiation spectra obtained from particle-in-cell simulations of the classical Weibel instability. Particularly interesting is the onset of the field growth, when the transient hard synchrotron-violating spectra are common as a result of the dominant role of the trapped population. This effect can serve as a distinct observational signature of the violent field growth in astrophysical sources and lab experiments. It is also interesting that a system with small-scale fields tends to evolve toward the small-angle jitter regime, which can, under certain conditions, dominate the overall emission of a source.

Radiative diagnostics of high-energy density plasmas is addressed in this paper. We propose that the radiation produced by energetic particles in small-scale magnetic field turbulence, which can occur in laser-plasma experiments, collisionless shocks, and during magnetic reconnection, can be used to deduce some properties of the turbulent magnetic field. Particles propagating through such turbulence encounter locally strong magnetic fields, but over lengths much shorter than a particle gyroradius. Consequently, the particle is accelerated but not deviated substantially from a straight line path. We develop the general jitter radiation solutions for this case and show that the resulting radiation is directly dependent upon the spectral distribution of the magnetic field through which the particle propagates. We demonstrate the power of this approach in considering the radiation produced by particles moving through a region in which a (Weibel-like) filamentation instability grows magnetic fields randomly oriented in a plane transverse to counterstreaming particle populations. We calculate the spectrum as would be seen from the original particle population and as could be seen by using a quasi-monoenergetic electron beam to probe the turbulent region at various angles to the filamentation axis.

Using our new 3-D relativistic electromagnetic particle (REMP) code parallelized with MPI, we investigated long-term particle acceleration associated with a relativistic electron-positron jet propagating in an unmagnetized ambient electron-positron plasma. The simulations were performed using a much longer simulation system than our previous simulations in order to investigate the full nonlinear stage of the Weibel instability and its particle acceleration mechanism. Cold jet electrons are thermalized and ambient electrons are accelerated in the resulting shocks. Acceleration of ambient electrons leads to a maximum ambient electron density three times larger than the original value. Behind the bow shock in the jet shock strong electromagnetic fields are generated. These fields may lead to time dependent afterglow emission. We calculated radiation from electrons propagating in a uniform parallel magnetic field to verify the technique. We also used the new technique to calculate emission from electrons based on simulations with a small system. We obtained spectra which are consistent with those generated from electrons propagating in turbulent magnetic fields with red noise. This turbulent magnetic field is similar to the magnetic field generated at an early nonlinear stage of the Weibel instability. A fully developed shock within a larger system generates a jitter/synchrotron spectrum.

Radiation from many astrophysical sources, e.g. gamma-ray bursts and active galactic nuclei, is believed to arise from relativistically shocked collisionless plasmas. Such sources often exhibit highly transient spectra evolving rapidly, compared with source lifetimes. Radiation emitted from these sources is typically associated with non-linear plasma physics, complex field topologies and non-thermal particle distributions. In such circumstances a standard synchrotron paradigm may fail to produce accurate conclusions regarding the underlying physics. Simulating spectral emission and spectral evolution numerically in various relativistic shock scenarios is then the only viable method to determine the detailed physical origin of the emitted spectra. In this Letter we present synthetic radiation spectra representing the early stage development of the filamentation (streaming) instability of an initially unmagnetized plasma, which is relevant for both collisionless shock formation and reconnection dynamics in relativistic astrophysical outflows, as well as for laboratory astrophysics experiments. Results were obtained using a highly efficient "in situ" diagnostics method, based on detailed particle-in-cell modeling of collisionless plasmas. The synthetic spectra obtained here are compared with those predicted by a semi-analytical model for jitter radiation from the filamentation instability, the latter including self-consistent generated field topologies and particle distributions obtained from the simulations reported upon here. Spectra exhibit dependence on the presence - or absence - of an inert plasma constituent, when comparing baryonic plasmas (i.e. containing protons) with pair plasmas. The results also illustrate that considerable care should be taken when using lower-dimensional models to obtain information about the astrophysical phenomena generating observed spectra.

In this paper we consider a stable particle with flavor mixing. We demonstrate that incoherent conversion of heavy mass eigenstates into light ones and vice versa can occur, as a result of elastic scattering. This effect is nontrivial for non-relativistic particles, for which the standard flavor oscillation ceases rapidly due to incoherence. We also prove that if a heavy state is bound in a gravitational potential and a light state is unbound, the mass-state conversion can lead to gradual "evaporation" of the mixed particle from the potential. A number of implications, ranging from the cosmic neutrino background distortions to scenarios of cold dark matter evaporation from halos, are addressed.

It has been suggested that galactic shock asymmetry induced by our galaxy's infall toward the Virgo Cluster may be a source of periodicity in cosmic ray exposure as the solar system oscillates perpendicular to the galactic plane. Here we investigate a mechanism by which cosmic rays might affect terrestrial biodiversity, ionization and dissociation in the atmosphere, resulting in depletion of ozone and a resulting increase in the dangerous solar UVB flux on the ground, with an improved ionization background computation averaged over a massive ensemble (about 7 x 10^5) shower simulations. We study minimal and full exposure to the postulated extragalactic background. The atmospheric effects are greater than with our earlier, simplified ionization model. At the lower end of the range effects are too small to be of serious consequence. At the upper end of the range, ~6 % global average loss of ozone column density exceeds that currently experienced due to effects such as accumulated chlorofluorocarbons. The intensity is less than a nearby supernova or galactic gamma-ray burst, but the duration would be about 10^6 times longer. Present UVB enhancement from current ozone depletion ~3% is a documented stress on the biosphere, but a depletion of the magnitude found at the upper end of our range would double the global average UVB flux. For estimates at the upper end of the range of the cosmic ray variability over geologic time, the mechanism of atmospheric ozone depletion may provide a major biological stress, which could easily bring about major loss of biodiversity. Future high energy astrophysical observations will resolve the question of whether such depletion is likely.

Jitter radiation is produced by relativistic electrons moving in turbulent small-scale magnetic fields such as those produced by streaming Weibel-type instabilities at collisionless shocks in weakly magnetized media. Here we present a comprehensive study of the dependence of the jitter radiation spectra on the properties of, in general, anisotropic magnetic turbulence. We have obtained that the radiation spectra do reflect, to some extent, properties of the magnetic field spatial distribution, yet the radiation field is anisotropic and sensitive to the viewing direction with respect to the field anisotropy direction. We explore the parameter space of the magnetic field distribution and its effect on the radiation spectrum. Some important results include: the presence of the harder-than-synchrotron segment below the peak frequency at some viewing angles, the presence of the high-frequency power-law tail even for a monoenergetic distribution of electrons, the dependence of the peak frequency on the field correlation length rather than the field strength, the strong correlation of the spectral parameters with the viewing angle. In general, we have found that even relatively minor changes in the magnetic field properties can produce very significant effects upon the jitter radiation spectra. We consider these results to be important for accurate interpretation of prompt gamma-ray burst spectra and possibly other sources.

Using our new 3-D relativistic electromagnetic particle (REMP) code parallelized with MPI, we investigated long-term particle acceleration associated with a relativistic electron-positron jet propagating in an unmagnetized ambient electron-positron plasma. We have also performed simulations with electron-ion jets. The simulations were performed using a much longer simulation system than our previous simulations in order to investigate the full nonlinear stage of the Weibel instability for electron-positron jets and its particle acceleration mechanism. Cold jet electrons are thermalized and ambient electrons are accelerated in the resulting shocks for both cases. Acceleration of ambient electrons leads to a maximum ambient electron density three times larger than the original value for pair plasmas. Behind the bow shock in the jet shock strong electromagnetic fields are generated. These fields may lead to time dependent afterglow emission. We calculated radiation from electrons propagating in a uniform parallel magnetic field to verify the technique. We also used the new technique to calculate emission from electrons based on simulations with a small system with two different cases for Lorentz factors (15 and 100). We obtained spectra which are consistent with those generated from electrons propagating in turbulent magnetic fields with red noise. This turbulent magnetic field is similar to the magnetic field generated at an early nonlinear stage of the Weibel instability.

Relativistic shocks are usually thought to occur in violent astrophysical explosions. These collisionless shocks are mediated by a plasma kinetic streaming instability, often loosely referred to as the Weibel instability, which generates strong magnetic fields "from scratch" very efficiently. In this review paper we discuss the shock micro-physics and present a recent model of "pre-conditioning" of an initially unmagnetized upstream region via the cosmic-ray-driven Weibel-type instability.

The origin of rapid spectral variability and certain spectral correlations of the prompt gamma-ray burst emission remains an intriguing question. The recently proposed theoretical model of the prompt emission is build upon unique spectral properties of jitter radiation -- the radiation from small-scale magnetic fields generated at a site of strong energy release (e.g., a relativistic collisionless shock in baryonic or pair-dominated ejecta, or a reconnection site in a magnetically-dominated outflow). Here we present the results of implementation of the model. We show that anisotropy of the jitter radiation pattern and relativistic shell kinematics altogether produce effects commonly observed in time-resolved spectra of the prompt emission, e.g., the softening of the spectrum below the peak energy within individual pulses in the prompt light-curve, the so-called "tracking" behavior (correlation of the observed flux with other spectral parameters), the emergence of hard, synchrotron-violating spectra at the beginning of individual spikes. Several observational predictions of the model are discussed.

We discuss how the electrons in relativistic GRB shocks can reach near-equipartition in energy with the protons. We emphasize the non-Fermi origin of such acceleration. We argue that the dynamics of the electrons in the foreshock region and at the shock front plays an important role. We also demonstrate that PIC simulations can directly probe this physics in the regimes relevant to GRBs.

Weibel instability turns out to be the a ubiquitous phenomenon in High-Energy Density environments, ranging from astrophysical sources, e.g., gamma-ray bursts, to laboratory experiments involving laser-produced plasmas. Relativistic particles (electrons) radiate in the Weibel-produced magnetic fields in the Jitter regime. Conventionally, in this regime, the particle deflections are considered to be smaller than the relativistic beaming angle of 1/$\gamma$ ($\gamma$ being the Lorentz factor of an emitting particle) and the particle distribution is assumed to be isotropic. This is a relatively idealized situation as far as lab experiments are concerned. We relax the assumption of the isotropy of radiating particle distribution and present the extension of the jitter theory amenable for comparisons with experimental data.

Plasma instabilities (e.g., Buneman, Weibel and other two-stream instabilities) excited in collisionless shocks are responsible for particle (electron, positron, and ion) acceleration. Using a new 3-D relativistic particle-in-cell code, we have investigated the particle acceleration and shock structure associated with an unmagnetized relativistic electron-positron jet propagating into an unmagnetized electron-positron plasma. The simulation has been performed using a long simulation system in order to study the nonlinear stages of the Weibel instability, the particle acceleration mechanism, and the shock structure. Cold jet electrons are thermalized and slowed while the ambient electrons are swept up to create a partially developed hydrodynamic (HD) like shock structure. In the leading shock, electron density increases by a factor of 3.5 in the simulation frame. Strong electromagnetic fields are generated in the trailing shock and provide an emission site. We discuss the possible implication of our simulation results within the AGN and GRB context.

This is a white paper submitted to the Stars and Stellar Evolution (SSE) Science Frontier Panel (SFP) of the NRC's Astronomy and Astrophysics 2010 Decadal Survey. The white paper is endorsed by the American Physical Society's (APS) Topical Group on Plasma Astrophysics (GPAP).

Using our new 3-D relativistic electromagnetic particle (REMP) code parallelized with MPI, we have investigated long-term particle acceleration associated with an relativistic electron-positron jet propagating in an unmagnetized ambient electron-positron plasma. The simulations have been performed using a much longer simulation system than our previous simulations in order to investigate the full nonlinear stage of the Weibel instability and its particle acceleration mechanism. Cold jet electrons are thermalized and ambient electrons are accelerated in the resulting shocks. The acceleration of ambient electrons leads to a maximum ambient electron density three times larger than the original value. Behind the bow shock in the jet shock strong electromagnetic fields are generated. These fields may lead to the afterglow emission. We have calculated the time evolution of the spectrum from two electrons propagating in a uniform parallel magnetic field to verify the technique.

Radiation from GRBs in the prompt phase, flares and an afterglow is thought to be produced by accelerated electrons in magnetic fields. Such emission may be produced at collisionless shocks of baryonic outflows or at reconnection sites (at least for the prompt and flares) of the magnetically dominated (Poynting flux driven) outflows, where no shocks presumably form at all. An astonishing recent discovery is that during reconnection strong small-scale magnetic fields are produced via the Weibel instability, very much like they are produced at relativistic shocks. The relevant physics has been successfully and extensively studied with the PIC simulations in 2D and, to some extent, in 3D for the past few years. We discuss how these simulations predict the existence of MeV-range synchrotron/jitter emission in some GRBs, which can be observed with Fermi. Recent results on modeling of the spectral variability and spectral correlations of the GRB prompt emission in the Weibel-jitter paradigm applicable to both baryonic and magnetic-dominated outflows is reviewed with the emphasis on observational predictions.

Cosmic rays accelerated by a shock form a streaming distribution of outgoing particles in the foreshock region. If the ambient fields are negligible compared to the shock and cosmic ray energetics, a stronger magnetic field can be generated in the shock upstream via the streaming (Weibel-type) instability. Here we develop a self-similar model of the foreshock region and calculate its structure, e.g., the magnetic field strength, its coherence scale, etc., as a function of the distance from the shock. Our model indicates that the entire foreshock region of thickness $\sim R/(2\Gamma_{\rm sh}^2)$, being comparable to the shock radius in the late afterglow phase when $\Gamma_{\rm sh}\sim1$, can be populated with large-scale and rather strong magnetic fields (of sub-gauss strengths with the coherence length of order $10^{17} {\rm cm}$) compared to the typical interstellar medium magnetic fields. The presence of such fields in the foreshock region is important for high efficiency of Fermi acceleration at the shock. Radiation from accelerated electrons in the foreshock fields can constitute a separate emission region radiating in the UV/optical through radio band, depending on time and shock parameters. We also speculate that these fields being eventually transported into the shock downstream can greatly increase radiative efficiency of a gamma-ray burst afterglow shock.

We address the question of whether numerical particle-in-cell (PIC) simulations and laboratory laser-plasma experiments can (or will be able to, in the near future) model realistic gamma-ray burst (GRB) shocks. For this, we compare the radiative cooling time, t_cool, of relativistic electrons in the shock magnetic fields to the microscopic dynamical time of collisionless relativistic shocks -- the inverse plasma frequency of protons, omega_pp^{-1}. We obtain that for t_cool*omega_pp^{-1}\lesssim ~few hundred, the electrons cool efficiently at or near the shock jump and are capable of emitiing away a large fraction of the shock energy. Such shocks are well-resolved in existing PIC simulations; therefore, the microscopic structure can be studied in detail. Since most of the emission in such shocks would be coming from the vicinity of the shock, the spectral power of the emitted radiation can be directly obtained from finite-length simulations and compared with observational data. Such radiative shocks correspond to the internal baryon-dominated GRB shocks for the conventional range of ejecta parameters. Fermi acceleration of electrons in such shocks is limited by electron cooling, hence the emitted spectrum should be lacking a non-thermal tail, whereas its peak likely falls in the multi-MeV range. Incidentally, the conditions in internal shocks are almost identical to those in laser-produced plasmas; thus, such GRB-like plasmas can be created and studied in laboratory experiments using the presently available Petawatt-scale laser facilities. An analysis of the external shocks shows that only the highly relativistic shocks, corresponding to the extremely early afterglow phase, can have efficient electron cooling in the shock transition. We emphasize the importance of radiative PIC simulations for further studies.

Nonthermal radiation observed from astrophysical systems containing relativistic jets and shocks, e.g., gamma-ray bursts (GRBs), active galactic nuclei (AGNs), and microquasars commonly exhibit power-law emission spectra. Recent PIC simulations of relativistic electron-ion (or electron-positron) jets injected into a stationary medium show that particle acceleration occurs within the downstream jet. In collisionless, relativistic shocks, particle (electron, positron, and ion) acceleration is due to plasma waves and their associated instabilities (e.g., the Weibel (filamentation) instability) created in the shock region. The simulations show that the Weibel instability is responsible for generating and amplifying highly non-uniform, small-scale magnetic fields. These fields contribute to the electron's transverse deflection behind the jet head. The resulting "jitter" radiation from deflected electrons has different properties compared to synchrotron radiation, which assumes a uniform magnetic field. Jitter radiation may be important for understanding the complex time evolution and/or spectra in gamma-ray bursts, relativistic jets in general, and supernova remnants.

We consider radiation emitted by the jitter mechanism in a Blandford-McKee self-similar blastwave. We assume the magnetic field configuration throughout the whole blastwave meets the condition for the emission of jitter radiation and we compute the ensuing images, light curves and spectra. The calculations are performed for both a uniform and a wind environment. We compare our jitter results to synchrotron results. We show that jitter radiation produces slightly different spectra than synchrotron, in particular between the self-absorption and the peak frequency, where the jitter spectrum is flat, while the synchrotron spectrum grows as \nu^{1/3}. The spectral difference is reflected in the early decay slope of the light curves. We conclude that jitter and synchrotron afterglows can be distinguished from each other with good quality observations. However, it is unlikely that the difference can explain the peculiar behavior of several recent observations, such as flat X-ray slopes and uncorrelated optical and X-ray behavior.

Nonthermal radiation observed from astrophysical systems containing relativistic jets and shocks, e.g., gamma-ray bursts (GRBs), active galactic nuclei (AGNs), and Galactic microquasar systems usually have power-law emission spectra. Recent PIC simulations of relativistic electron-ion (electron-positron) jets injected into a stationary medium show that particle acceleration occurs within the downstream jet. In the presence of relativistic jets, instabilities such as the Buneman instability, other two-streaming instability, and the Weibel (filamentation) instability create collisionless shocks, which are responsible for particle (electron, positron, and ion) acceleration. The simulation results show that the Weibel instability is responsible for generating and amplifying highly nonuniform, small-scale magnetic fields. These magnetic fields contribute to the electron's transverse deflection behind the jet head. The ``jitter'' radiation from deflected electrons in small-scale magnetic fields has different properties than synchrotron radiation which is calculated in a uniform magnetic field. This jitter radiation, a case of diffusive synchrotron radiation, may be important to understand the complex time evolution and/or spectral structure in gamma-ray bursts, relativistic jets, and supernova remnants.

Nonthermal radiation observed from astrophysical systems containing relativistic jets and shocks e.g. gamma-ray bursts (GRBs) active galactic nuclei (AGNs) and microquasars commonly exhibit power-law emission spectra. Recent PIC simulations of relativistic electron-ion (or electron-positron) jets injected into a stationary medium show that particle acceleration occurs within the downstream jet. In collisionless relativistic shocks particle (electron, positron and ion) acceleration is due to plasma waves and their associated instabilities (e.g. the Weibel (filamentation) instability) created in the shock region. The simulations show that the Weibel instability is responsible for generating and amplifying highly non-uniform small-scale magnetic fields. These fields contribute to the electron's transverse deflection behind the jet head. The resulting ``jitter'' radiation from deflected electrons has different properties compared to synchrotron radiation which assumes a uniform magnetic field. Jitter radiation may be important for understanding the complex time evolution and/or spectra in gamma-ray bursts, relativistic jets in general and supernova remnants.

Relativistic electrons moving into a highly tangled magnetic field emit jitter radiation. We present a detailed computation of the jitter radiation spectrum, including self-absorption, for electrons inside Weibel-like shock generated magnetic fields. We apply our results to the case of the prompt and afterglow emission of gamma-ray bursts. We show that jitter emission can reproduce most of the observed features with some important differences with respect to standard synchrotron, especially in the frequency range between the self-absorption and the peak frequency. We discuss the similarities and differences between jitter and synchrotron and discuss experiments that can disentangle the two mechanisms.

Resent observations of a number of galaxy clusters using the Sunyaev-Zel'dovich effect indicate that about 1/3 of baryonic mass is missing from the hot intracluster medium (ICM), which is significantly larger than the fraction of stars and cool gas, which account for only about 10%. Here we address the question whether the remaining $22\pm10%$ can be accounted for by thermal evaporation of gas from clusters. We have found that evaporation can occur only from the cluster ``surface'', $r\sim r_{\rm vir}$, and not from it's interior. We evaluated particle diffusion through the magnetized ICM for several scenarios of ISM turbulence and found that diffusivity is suppressed by at least a factor of 100 or more, compared to the Spitzer value. Thus, only particles from radii $r\ga0.9r_{\rm vir}$ can evaporate. Diffusion of particles from inside the cluster, $r\la0.9r_{\rm vir}$, takes longer than the Hubble time. This lowers the cluster-averaged fraction of the evaporated hot gas to few percent or less. However, if the missing hot component {\it is indeed} due to evaporation, this strongly constrains the magnetic field structure in the cluster envelope, namely either (i) the gas is completely unmagnetized ($B\le10^{-21}$ gauss) in the cluster halo or (ii) the magnetic fields in the ICM are rather homogeneous and non-turbulent.

Recent work has revealed a 62 (+/-) 3-million-year cycle in the fossil diversity in the past 542 My, however no plausible mechanism has been found. We propose that the cycle may be caused by modulation of cosmic ray (CR) flux by the Solar system vertical oscillation (64 My period) in the galaxy, the galactic north-south anisotropy of CR production in the galactic halo/wind/termination shock (due to the galactic motion toward the Virgo cluster), and the shielding by galactic magnetic fields. We revisit the mechanism of CR propagation and show that CR flux can vary by a factor of about 4.6 and reach a maximum at north-most displacement of the Sun. The very high statistical significance of (i) the phase agreement between Solar north-ward excursions and the diversity minima and (ii) the correlation of the magnitude of diversity drops with CR amplitudes through all cycles provide solid support for our model. Various observational predictions which can be used to confirm or falsify our hypothesis are presented.

The standard model of GRB afterglows assumes that the radiation observed as a delayed emission is of synchrotron origin, which requires the shock magnetic field to be relatively homogeneous on small scales. An alternative mechanism -- jitter radiation, which traditionally has been applied to the prompt emission -- substitutes synchrotron when the magnetic field is tangled on a microscopic scale. Such fields are produced at relativistic shocks by the Weibel instability. Here we explore the possibility that small-scale fields populate afterglow shocks. We derive the spectrum of jitter radiation under the afterglow conditions. We also derive the afterglow lightcurves for the ISM and Wind profiles of the ambient density. Jitter self-absorption is calculated here for the first time. We find that jitter radiation can produce afterglows similar to synchrotron-generated ones, but with some important differences. We compare the predictions of the two emission mechanisms. By fitting observational data to the synchrotron and jitter afterglow lightcurves, it can be possible to discriminate between the small-scale vs large-scale magnetic field models in afterglow shocks.

Gamma-Ray Bursts (GRBs) directed at Earth from within a few kpc may have damaged the biosphere, primarily though changes in atmospheric chemistry which admit greatly increased Solar UV. However, GRBs are highly variable in spectrum and duration. Recent observations indicate that short (~0.1 s) burst GRBs, which have harder spectra, may be sufficiently abundant at low redshift that they may offer an additional significant effect. A much longer timescale is associated with shock breakout luminosity observed in the soft X-ray (~10^3 s) and UV (~10^5 s) emission, and radioactive decay gamma-ray line radiation emitted during the light curve phase of supernovae (~10^7 s). Here we generalize our atmospheric computations to include a broad range of peak photon energies and investigate the effect of burst duration while holding total fluence and other parameters constant. The results can be used to estimate the probable impact of various kinds of ionizing events (such as short GRBs, X-ray flashes, supernovae) upon the terrestrial atmosphere. We find that the ultimate intensity of atmospheric effects varies only slightly with burst duration from 10^-1 s to 10^8 s. Therefore, the effect of many astrophysical events causing atmospheric ionization can be approximated without including time development. Detailed modeling requires specification of the season and latitude of the event. Harder photon spectra produce greater atmospheric effects for spectra with peaks up to about 20 MeV, because of greater penetration into the stratosphere.

The shock model of gamma-ray bursts (GRBs) contains two equipartition parameters: the magnetic energy density and the kinetic energy density of the electrons relative to the total energy density of the shock, "epsilon_B" and "epsilon_e", respectively. These are free parameters within the model. Whereas the Weibel shock theory and numerical simulations fix "epsilon_B" at the level of ~few times(10^{-3}...10^{-4}), no understanding of "epsilon_e" exists so far. Here we demonstrate that it inevitably follows from the theory that "epsilon_e"~(epsilon_B)^(1/2). The GRB afteglow data fully agree with this theoretical prediction. Our result explains why the electrons are close to equipartition in GRBs. The "epsilon_e"-"epsilon_B" relation can potentially be used to reduce the number of free parameters in afterglow models.

We demonstrate that the rapid spectral variability of prompt GRBs is an inherent property of radiation emitted from shock-generated, highly anisotropic small-scale magnetic fields. We interpret the hard-to-soft evolution and the correlation of the soft index $\alpha$ with the photon flux observed in GRBs as a combined effect of temporal variation of the shock viewing angle and relativistic aberration of an individual thin, instantaneously illuminated shell. The model predicts that about a quarter of time-resolved spectra should have hard spectra, violating the synchrotron $\alpha=-2/3$ limit. The model also naturally explains why the peak of the distribution of $\alpha$ is at $\alpha\sim-1$. The presence of a low-energy break in the jitter spectrum at oblique angles also explains the appearance of a soft X-ray component in some GRBs and their paucity. We emphasize that our theory is based solely on the first principles and contains no ad hoc (phenomenological) assumptions.

The origin of the micro-Gauss magnetic fields in galaxy clusters is one of the outstanding problem of modern cosmology. We have performed three-dimensional particle-in-cell simulations of the nonrelativistic Weibel instability in an electron-proton plasma, in conditions typical of cosmological shocks. These simulations indicate that cluster fields could have been produced by shocks propagating through the intergalactic medium during the formation of large-scale structure or by shocks within the cluster. The strengths of the shock-generated fields range from tens of nano-Gauss in the intercluster medium to a few micro-Gauss inside galaxy clusters.

We develop the theory of jitter radiation from GRB shocks containing small-scale magnetic fields and propagating at an angle with respect to the line of sight. We demonstrate that the spectra vary considerably: the low-energy photon index, $\alpha$, ranges from 0 to -1 as the apparent viewing angle goes from 0 to $\pi/2$. Thus, we interpret the hard-to-soft evolution and the correlation of $\alpha$ with the photon flux observed in GRBs as a combined effect of temporal variation of the viewing angle and relativistic aberration of an individual thin, instantaneously illuminated shell. The model predicts that about a quarter of time-resolved spectra should have hard spectra, violating the synchrotron $\alpha=-2/3$ line of death. The model also naturally explains why the peak of the distribution of $\alpha$ is at $\alpha\approx-1$. The presence of a low-energy break in the jitter spectrum at oblique angles also explains the appearance of a soft X-ray component in some GRBs and a relatively small number of them. We emphasize that our theory is based solely on the first principles and contains no {\it ad hoc} (phenomenological) assumptions.

Self-gravitating systems with nonlocal, long-range interactions are described by nonextensive statistics. Recently, Leubner demonstrated that the nonextensivity parameter $\kappa$ should be negative for self-gravitating, pressureless systems, such as dark matter halos. The equation for the spherically symmetric nonextensive dark matter halos has also been derived. Here we demonstrate that this equation is identical to the classical Lane-Emden equation describing the structure of self-gravitating polytropic spheres. This establishes an intimate connection between self-gravitating polytropes and nonextensive thermostatistics. Moreover, based on this fact and observational data, we put a stronger constraint on $\kappa$, namely $\kappa<-3.4$.

Gamma-Ray Bursts (GRBs) are likely to have made a number of significant impacts on the Earth during the last billion years. We have used a two-dimensional atmospheric model to investigate the effects on the Earth's atmosphere of GRBs delivering a range of fluences, at various latitudes, at the equinoxes and solstices, and at different times of day. We have estimated DNA damage levels caused by increased solar UVB radiation, reduction in solar visible light due to $\mathrm{NO_2}$ opacity, and deposition of nitrates through rainout of $\mathrm{HNO_3}$. For the ``typical'' nearest burst in the last billion years, we find globally averaged ozone depletion up to 38%. Localized depletion reaches as much as 74%. Significant global depletion (at least 10%) persists up to about 7 years after the burst. Our results depend strongly on time of year and latitude over which the burst occurs. We find DNA damage of up to 16 times the normal annual global average, well above lethal levels for simple life forms such as phytoplankton. The greatest damage occurs at low to mid latitudes. We find reductions in visible sunlight of a few percent, primarily in the polar regions. Nitrate deposition similar to or slightly greater than that currently caused by lightning is also observed, lasting several years. We discuss how these results support the hypothesis that the Late Ordovician mass extinction may have been initiated by a GRB.

We strongly disagree with the criticism by G. Fleishman, which has been published on the astro-ph arXiv (astro-ph/0502245), of our paper on jitter radiation from GRB shocks (Medvedev 2000, ApJ, 540, 704). In this note, we present the rebuttal of all critical points raised, demonstrating that our original paper (i) contains no errors, and (ii) correctly describes prompt GRB emission spectra.

A substantial fraction of prompt GRB spectra have soft spectral indexes exceeding the maximum allowed by the synchrotron model alpha_{max}=-2/3. Some spectra also exhibit very sharp break at E_p, inconsistent with the smooth synchrotron spectra. These facts pose a serious problem for the ``optically thin synchrotron'' interpretation of the prompt emission. We review various models suggested in order to resolve this puzzle.

We investigate the long-time evolution of magnetic fields generated by the two-stream instability at ultra- and sub-relativistic astrophysical collisionless shocks. Based on 3D PIC simulation results, we introduce a 2D toy model of interacting current filaments. Within the framework of this model, we demonstrate that the field correlation scale in the region far downstream the shock grows nearly as the light crossing time, lambda(t) ~ ct, thus making the diffusive field dissipation inefficient. The obtained theoretical scaling is tested using numerical PIC simulations. This result extends our understanding of the structure of collisionless shocks in gamma-ray bursts and other astrophysical objects.

(Abridged) A new type of self-similar hot viscous radiative accretion flow onto a rapidly spinning neutron star has recently been discovered. This ``hot brake'' flow forms in the two-temperature zone (close to a central object), but at a sufficiently low accretion rate and a high spin it may extend in the radial direction beyond ~300 Schwarzchild radii into a one-temperature zone. When the spin of the star is small enough, the flow transforms smoothly to an advection-dominated accretion flow. All gas parameters (density, angular velocity, temperature, luminosity, angular momentum flux) except for the radial velocity are independent of the mass accretion rate. The radiative efficiency may be arbitrarily large as M-dot -> 0. The gas angular momentum is transported outward under most conditions, hence the central star is nearly always spun-down. The flow is convectively stable. We also find that themal conduction in the flow is strong enough to make the flow thermally stable. The very fact that the density, temperature and angular velocity of the gas at any radius in the hot brake flow are completely independent of the outer and inner (except for the star spin) boundary conditions implies that the flow cannot be smoothly matched to a general external medium as well as to general conditions on the star surface. We demonstrate that there are two extra self-similar solutions: one bridges the gap between the original solution and the external medium, and another represents a boundary layer between the bulk of the flow and the star surface, in which the gas temperature rapidly drops while the density builds up. Finally, we briefly discuss that a hot brake flow may form around other rapidly spinning compact objects: white dwarfs and black holes.

We consider hot accretion onto a rapidly spinning neutron star (or any other compact object with a surface). A radiative hot settling flow has been discovered at low accretion rates in the early work by Medvedev & Narayan (2001) and analytical solution has been presented. It was shown later that this flow can match external medium smoothly, thus enforcing its physical feasibility. Here we complete the study of the global structure of such hot accretion by presenting the analytical solution for the boundary later, which forms between the bulk of the flow and the stellar surface. We confirm our results via full numerical solution of height-integrated two-temperature hydrodynamic equations.

Gamma-ray bursts (hereafter GRB) produce a flux of radiation detectable across the observable Universe, and at least some of them are associated with galaxies. A GRB within our own Ggalaxy could do considerable damage to the Earth's biosphere; rate estimates suggest that a dangerously near GRB should occur on average two or more times per billion years. At least five times in the history of life, the Earth experienced mass extinctions that eliminated a large percentage of the biota. Many possible causes have been documented, and GRB may also have contributed. The late Ordovician mass extinction approximately 440 million years ago may be at least partly the result of a GRB. A special feature of GRB in terms of terrestrial effects is a nearly impulsive energy input of order 10 s. Due to expected severe depletion of the ozone layer, intense solar ultraviolet radiation would result from a nearby GRB, and some of the patterns of extinction and survivorship at this time may be attributable to elevated levels of UV radiation reaching the Earth. In addition a GRB could trigger the global cooling which occurs at the end of the Ordovician period that follows an interval of relatively warm climate. Intense rapid cooling and glaciation at that time, previously identified as the probable cause of this mass extinction, may have resulted from a GRB.

Charged plasma particles form an anisotropic counter-streaming distribution at the front of a shock. In the near- and ultra-relativistic regimes, Weibel (or two-stream) instability produces near-equipartition, chaotic, small-scale magnetic fields. The fields introduce effective collisions and "thermalize" the plasma particles via pitch-angle scattering. The properties of jitter radiation emitted by accelerated electrons from these small-scale fields are markedly different from synchrotron spectra. Here we theoretical summary and the results of recent numerical 3D PIC plasma kinetic simulations. The relation of the obtained results to the theory of gamma-ray bursts is outlined.

We report Chandra observations of the dwarf nova WX Hyi in quiescence. The X-ray spectrum displays strong and narrow emission lines of N, O, Mg, Ne, Si, S and Fe. The various ionization states implied by the lines suggest that the emission is produced within a flow spanning a wide temperature range, from T ~ 10^6 K to T >~ 10^8 K. Line diagnostics indicate that most of the radiation originates from a very dense region, with n ~ 10^{13}-10^{14} cm^{-3}. The Chandra data allow the first tests of specific models proposed in the literature for the X-ray emission in quiescent dwarf novae. We have computed the spectra for a set of models ranging from hot boundary layers, to hot settling flows solutions, to X-ray emitting coronae. WX Hyi differs from other dwarf novae observed at minimum in having much stronger low temperature lines, which prove difficult to fit with existing models, and possibly a very strong, broad O VII line, perhaps produced in a wind moving at a few x 10^3 km/s. The accretion rate inferred from the X-rays is lower than the value inferred from the UV. The presence of high-velocity mass ejection could account for this discrepancy while at the same time explaining the presence of the broad O VII line. If this interpretation is correct, it would provide the first detection of a wind from a dwarf nova in quiescence.

We present the first three-dimensional fully kinetic electromagnetic relativistic particle-in-cell simulations of the collision of two interpenetrating plasma shells. The highly accurate plasma-kinetic "particle-in-cell" (with the total of $10^8$ particles) parallel code OSIRIS has been used. Our simulations show: (i) the generation of long-lived near-equipartition (electro)magnetic fields, (ii) non-thermal particle acceleration, and (iii) short-scale to long-scale magnetic field evolution, in the collision region. Our results provide new insights into the magnetic field generation and particle acceleration in relativistic and sub-relativistic colliding streams of particles, which are present in gamma-ray bursters, supernova remnants, relativistic jets, pulsar winds, etc..

Medvedev & Narayan have described a hot accretion flow onto a spinning neutron star in which the gas viscously brakes the spin of the star. Their self-similar solution has the surprising property that the density, temperature and angular velocity of the gas at any radius are completely independent of the outer boundary conditions. Hence, the solution cannot be matched to a general external medium. We resolve this paradoxical situation by showing that there is a second self-similar solution which bridges the gap between the original solution and the external medium. This new solution has an extra degree of freedom which permits it to match general outer boundary conditions. We confirm the main features of the analytical results with a full numerical solution.

We have calculated lower limits for the global effective thermal conductivity in a sample of 165 Abell clusters. We assumed that cluster X-ray luminosity is compensated by a conductive heat flux which we then compare with an upper limit to the temperature gradient inferred from the cluster temperature and radius. This gives a lower limit to the thermal conductivity and therefore to the relative suppression from the Spitzer conductivity, kappa_{Sp}. Not a single cluster requires super-Spitzer values of thermal conduction to balance the observed X-ray luminosity. The suppression coefficient f=kappa/kappa_{Sp} is clustered in a range 10^{-2} < f < 0.4. A weak dependence of $f$ versus $z$ is observed over $ 0 < z <0.41$. Possible biases and/or selection effects are discussed.

The energetics of electromagnetic acceleration of ultra-high-energy cosmic rays (UHECRs) is constrained both by confinement of a particle within an acceleration site and by radiative energy losses of the particle in the confining magnetic fields. We demonstrate that the detection of ~ 3 x 10^{20} eV events is inconsistent with the hypothesis that compact cosmic accelerators with high magnetic fields can be the sources of UHECRs. This rules out the most popular candidates, namely spinning neutron stars, active galactic nuclei (AGNs), and gamma-ray burst blast waves. Galaxy clusters and, perhaps, AGN radio lobes remain the only possible (although not very strong) candidates for UHECR acceleration sites. Our analysis places no limit on linear accelerators. With the data from the future Auger experiment one should be able to answer whether a conventional theory works or some new physics is required to explain the origin of UHECRs.

We study the structure and properties of hot MHD accretion onto a Kerr black hole. In such a system, the hole is magnetically coupled to the inflowing gas and exerts a torque onto the accretion flow. A hot settling flow can form around the hole and transport the angular momentum outward, to the outer edge of the flow. Unlike other hot flows, such as advection- and convection-dominated flows and inflow-outflow solutions (ADAFs, CDAFs, and ADIOS), the properties of the hot settling flow are determined by the spin of the central black hole, but are insensitive to the mass accretion rate. Therefore, it may be possible to identify rapidly spinning BHs simply from their broad-band spectra. Observationally, the hot settling flow around a Kerr hole is somewhat similar to other hot flows in that they all have hard, power-law spectra and relatively low luminosities. Thus, most black hole candidates in the low/hard and, perhaps, intermediate X-ray state may potentially accrete via the hot settling flow. However, a settling flow will be somewhat more luminous than ADAFs/CDAFs/ADIOS, will exhibit high variability in X-rays, and may have relativistic jets. This suggests that galactic microquasars and active galactic nuclei may be powered by hot settling flows. We identify several galactic X-ray sources as the best candidates.

We present dynamically consistent solutions for hot accretion onto unmagnetized, rotating white dwarfs (WDs) in five quiescent dwarf novae. The measured WD rotation rates (and other system parameters) in RX And, SS Cyg, U Gem, VW Hyi and WZ Sge imply spindown of the WD by an extended hot flow emitting most of its X-rays in the vicinity of the stellar surface. In general, energy advection is absent and the flow is stable to convection and hydrodynamical outflows. In rapidly rotating systems, the X-ray luminosity provides only an upper limit on the quiescent accretion rate because of substantial stellar spindown luminosity. We suggest that the presence of hot flows in quiescent dwarf novae may limit the long-term WD rotation rates to significantly sub-Keplerian values.

We estimate the thermal conductivity of a weakly collisional magnetized plasma with chaotic magnetic field fluctuations. When the fluctuation spectrum extends over two or more decades in wave-vector, we find that thermal conduction is very efficient; the conduction coefficient is only a factor ~5 below the classical Spitzer estimate. We suggest that conduction could play a significant role in cooling flows in clusters of galaxies.

Recently, Medvedev & Narayan (2001) discovered a new type of accretion flow, a hot settling flow around a rapidly rotating neutron star. The flow is cooling-dominated and energetically similar to the Shapiro, Lightman, & Eardley (1976, SLE) solution. Since the SLE solution is known to be thermally unstable, one might suspect that the new solution would also be unstable. However, due to the very high temperature of the accreting gas, thermal conduction is very strong and could suppress the thermal instability. We analyze the role of thermal conduction in both the hot settling flow and the SLE solution. In the hot settling flow collisions are very rare. Therefore, thermal transport occurs via free streaming of electrons along tangled magnetic field lines. We find that conduction is strong enough to make the flow marginally stable. In contrast, in the cooler SLE solution, conduction is via collisional, Spitzer-type transport. In this case, conduction is weaker, and we find that the SLE solution is thermally unstable even in the presence of conduction.

We present a theory of ultrarelativistic collisionless shocks based on the relativistic kinetic two-stream instability. We demonstrate that the shock front is unstable to the generation of small-scale, randomly tangled magnetic fields. These fields are strong enough to scatter the energetic incoming (in the shock frame) protons and electrons over pitch angle and, therefore, to convert their kinetic energy of bulk motion into heat with very high efficiency. This validates the use of MHD approximation and the shock jump conditions in particular. The effective collisions are also necessary for the diffusive Fermi acceleration of electrons to operate and produce an observed power-law. Finally, these strong (sub-equipartition) magnetic fields are also required for the efficient synchrotron-type radiation emission from the shocks. The predicted magnetic fields have an impact on polarization properties of the observed radiation (e.g., a linear polarization from a jet-like ejecta and polarization scintillations in radio for a spherical one) and on its spectrum. We present an analytical theory of jitter radiation, which is emitted when the magnetic field is correlated on scales smaller then the gyration (Larmor) radius of the accelerated electrons. A composite jitter+synchrotron model of GRB $\gamma$-ray emission from internal shocks is capable of resolving many puzzles of GRB spectra, such as the violation of the ``line of death'', sharp spectral breaks, and multiple spectral components seen in some bursts (good examples are GRB910503, GRB910402, etc.). We stress that simultaneous detection of both spectral components opens a way to a precise diagnostics of the conditions in GRB shocks. We also discuss the relation of our results to other systems, such as internal shocks in blazars, radio lobes, and supernova shocks.

We present analytical and numerical solutions which describe a hot, viscous, two-temperature accretion flow onto a rotating neutron star or any other rotating compact star with a surface. We assume Coulomb coupling between the protons and electrons, and free-free cooling from the electrons. Outside a thin boundary layer, where the accretion flow meets the star, we show that there is an extended settling region which is well-described by two self-similar solutions: (i) a two-temperature solution which is valid in an inner zone $r\le10^{2.5}$ ($r$ is in Schwarzchild units), and (ii) a one-temperature solution at larger radii. In both zones, $\rho\propto r^{-2}, \Omega\propto r^{-3/2}, v\propto r^0, T_p\propto r^{-1}$; in the two-temperature zone, $T_e\propto r^{-1/2}$. The luminosity of the settling zone arises from the rotational energy of the star as the star is braked by viscosity. Hence the luminosity and the flow parameters (density, temperature, angular velocity) are independent of $\dot M$. The settling solution described here is not advection-dominated, and is thus different from the self-similar ADAF found around black holes. When the spin of the star is small enough, however, the present solution transforms smoothly to a (settling) ADAF. We carried out a stability analysis of the settling flow. The flow is convectively and viscously stable and is unlikely to have strong winds or outflows. Unlike another cooling-dominated system --- the SLE disk, --- the settling flow is thermally stable provided that thermal conduction is taken into account. This strong saturated-like thermoconduction does not change the structure of the flow.

The self-interacting dark matter (SIDM) model with flavor mixing (astro-ph/0010616) was proposed to resolve problems of the CDM model on small scales by keeping attractive features of both SIDM and annihilating dark matter, and simultaneously avoid their drawbacks. A dark particle produced in a flavor eigenstate will separate into two mass eigenstates because they propagate with different velocities and, in a gravitational filed, along different geodesics, see Fig. 1. Thus, in the flavor-mixed SIDM, dark halos are made of heavy eigenstates, whereas light eigenstates may leave the halo. Collisions (elastic scattering) of mass states results in eigenstate conversion, see Fig. 2, which leads to the gradual decrease of the halo mass in high-density cusps. On the other hand, in the early Universe, one may expect a problem of over-production of light (hot) particles over heavy (cold) when the temperature of dark matter falls below the mass of the heavy component. We show how this problem is avoided.

We present analytical and numerical solutions which describe a hot, viscous, two-temperature accretion flow onto a neutron star or any other compact star with a surface. We assume Coulomb coupling between the protons and electrons, and free-free cooling from the electrons. Outside a thin boundary layer, where the accretion flow meets the star, we show that there is an extended settling region which is well-described by two self-similar solutions: (1) a two-temperature solution which is valid in an inner zone $r\le10^{2.5}$ ($r$ is in Schwarzchild units), and (2) a one-temperature solution at larger radii. In both zones, $\rho\propto r^{-2}, \Omega\propto r^{-3/2}, v\propto r^0,\ T_p\propto r^{-1}$; in the two-temperature zone, $T_e\propto r^{-1/2}$. The luminosity of the settling zone arises from the rotational energy of the star as the star is braked by viscosity; hence the luminosity is independent of $\dot M$. The settling solution is convectively and viscously stable and is unlikely to have strong winds or outflows. The flow is thermally unstable, but the instability may be stabilized by thermal conduction. The settling solution described here is not advection-dominated, and is thus different from the self-similar ADAF found around black holes. When the spin of the star is small enough, however, the present solution transforms smoothly to a (settling) ADAF.

We present a theory of generation of strong (sub-equipartition) magnetic fields in relativistic collisionless GRB shocks. These fields produced by the kinetic two-stream instability are tangled on very small spatial scales. This has a clear signature in the otherwise synchrotron(-self-Compton) $\gamma$-ray spectrum. Second, we present an analytical theory of jitter radiation, which is emitted when the correlation length of the magnetic field is smaller then the gyration (Larmor) radius of the accelerated electrons. We demonstrate that the spectral power $P(\nu)$ for pure jitter radiation is well-described by a sharply broken power-law: $P(\nu)\propto\nu^1$ for $\nu<\nu_j$ and $P(\nu)\propto\nu^{-(p-1)/2}$ for $\nu>\nu_j$, where $p$ is the electron power-law index and $\nu_j$ is the jitter break, which is independent of the magnetic field strength and depends on the shock energetics and kinematics. Here we mostly focus on the first problem. The radiation theory and comparison with observations will be discussed in the forthcoming publications.

The crisis of the cold dark matter and problems of the self-interacting dark matter models is resolved by postulating flavor mixing of dark matter particles. Flavor-mixed particles segregate in the gravitational field to form dark halos composed of heavy mass eigenstates. Since these particles are mixed in the interaction basis, elastic collisions convert some of heavy eigenstates into light ones which leave dense central regions of the halo. This annihilation-like process will soften dense central cusps of halos. The proposed model accumulates most of the attractive features of self-interacting and annihilating dark matter models, but does not suffer from their severe drawbacks. This model is natural; it does not require fine tuning.

We investigate the structure of self-gravitating polytropic stellar systems. We present a method which allows to obtain approximate analytical solutions, $\psi_{n+\epsilon}({\bf x})$, of the nonlinear Poisson equation with the polytropic index $n+\epsilon$, given the solution $\psi_n({\bf x})$ with the polytropic index n, for any positive or negative $\epsilon$ such that $|\epsilon|\ll1$. Application of this method to the spherically symmetric stellar polytropes with $n\simeq5$ yields the solutions which describe spatially bound systems if n<5 and the formation of a second core if n>5. A heuristic approximate expressions for the radial profiles are also presented. Due to the duality between stellar and gas polytropes, our results are valid for gaseous, self-gravitating polytropic systems (e.g., molecular clouds) with the index $\gamma\simeq 6/5$. Stability of such systems and observational consequences for both stellar and gaseous systems are discussed.

We consider hot, two-temperature, viscous accretion onto a rotating, unmagnetized neutron star. We assume Coulomb coupling betweenthe protons and electrons, and free-free cooling from the electrons. We show that the accretion flow has an extended settling region which can be described by means of two analytical self-similar solutions: a two-temperature solution which is valid in an inner zone, r<10^{2.5}, where r is the radius in Schwarzchild units; and a one-temperature solution which is valid in an outer zone, r>10^{2.5}. In both zones the density varies as \rho ~ r^{-2} and the angular velocity as \Omega ~ r^{-3/2}. We solve the flow equations numerically and confirm that the analytical solutions are accurate. The settling flow radiates the energy dissipated by viscosity; so it is not advection-dominated. Except for the radial velocity, all other gas properties - density, angular velocity, temperature, luminosity, angular momentum flux - are independent of the mass accretion rate; these quantities do depend sensitively on the spin of the neutron star. The angular momentum flux is outward under most conditions; therefore, the central star is nearly always spun-down. The luminosity of the settling zone arises from the rotational energy that is released as the star is braked by viscosity, and the contribution from gravity is small; hence the radiative efficiency can be arbitrarily large at low $\mdot$. For reasonable values of the gas adiabatic index, the Bernoulli parameter is negative; therefore, a strong outflow or wind is not expected. The flow is convectively stable, but may be thermally unstable. The described solution is not advection-dominated; however, when the spin of the star is small enough, it transforms smoothly to an advection-dominated branch of solution.

The crisis of both cold and collisional dark matter (DM) models is resolved by postulating flavor mixing of DM particles. Flavor-mixed particles segregate in the gravitational field to form dark halos composed of heavy mass eigenstates. Since these particles are mixed in the interaction basis, elastic collisions convert some of heavy eigenstates into light ones which leave the halo. This annihilation-like process will soften dense central cusps of halos. The proposed model accumulates most of the attractive features of self-interacting and annihilating DM models, but does not suffer from their severe drawbacks. This model is natural does not require fine tuning.

Particle heating in advection-dominated accretion flows (ADAFs) by nonlinear MHD (Alfvenic) turbulence is investigated. Such turbulence with highly-fluctuating magnetic fields, $\tilde B\sim B_0$, is believed to be naturally produced by the magnetic shearing instability near the nonlinear saturation. It is shown that the energy is dissipated in the parallel cascade, which occurs due to nonlinear compressibility of high-amplitude turbulence, and predominantly heats protons, but not electrons. The conservative limit on the electron--to--proton heating fraction is $\delta\lesssim{\rm few}\times10^{-2}$.

Abridged.-- We demonstrate that the radiation emitted by ultrarelativistic electrons in highly nonuniform, small-scale magnetic fields is different from synchrotron radiation if the electron's transverse deflections in these fields are much smaller than the beaming angle. A quantitative analytical theory of this radiation, which we refer to as jitter radiation, is developed. It is shown that the emergent spectrum is determined by statistical properties of the magnetic field. As an example,we then use the model of a magnetic field in internal shocks of GRBs. The spectral power distribution of radiation produced by the power-law electrons is well described by a sharply broken power-law with indices 1 and -(p-1)/2 and the jitter break frequency is independent of the field strength but depends on the electron density in the ejecta. Since large-scale fields may also be present in the ejecta, we construct a two-component, jitter+synchrotron spectral model of the prompt $\gamma$-ray emission. Quite surprisingly, this model seems to be readily capable of explaining several properties of time-resolved spectra of some GRBs, such as (i) the violation of the constraint on the low-energy spectral index called the synchrotron ``line of death'', (ii) the sharp spectral break at the peak frequency, inconsistent with the broad synchrotron bump, (iii) the evidence for two spectral sub-components, and (iv) possible existence of emission features called ``GRB lines''. We believe these facts strongly support both the existence of small-scale magnetic fields and the proposed radiation mechanism from GRB shocks. As an example, we use the composite model to analyze GRB 910503 which has two spectral peaks.

We demostrate that the measurement of fluctuations of polarization due to the galactic interstellar scintillations may be used to study the structure of the radiation field at compact radio sources. We develop a mathematical formalism and demonstrate it on a simple analytical model in which the scale of the polarization variation through the source is comparable to the source size. The predicted amplitude of modulation of the polarized radiation flux is ~20% x (pi_s) x (m_sc), where (pi_s) is the characteristic degree of polarization of radiation at the source and (m_sc) is the typical modulation index due to scattering, i.e., (m_sc)~1 for diffractive scintillations and (m_sc)<1 for refractive scintillations.

All axisymmetric self-similar equilibria of self-gravitating, rotating, isothermal systems are identified by solving the nonlinear Poisson equation analytically. There are two families of equilibria: (1) Cylindrically symmetric solutions in which the density varies with cylindrical radius as R^(-alpha), with 0 <= alpha <= 2. (2) Axially symmetric solutions in which the density varies as f(theta)/r^2, where `r' is the spherical radius and `theta' is the co-latitude. The singular isothermal sphere is a special case of the latter class with f(theta)=constant. The axially symmetric equilibrium configurations form a two-parameter family of solutions and include equilibria which are surprisingly asymmetric with respect to the equatorial plane. The asymmetric equilibria are, however, not force-free at the singular points r=0, infinity, and their relevance to real systems is unclear. For each hydrodynamic equilibrium, we determine the phase-space distribution of the collisionless analog.

We show that the relativistic two-stream instability can naturally generate strong magnetic fields with 1e-5 - 1e-1 of the equipartition energy density, in the collisionless shocks of Gamma-Ray-Burst (GRB) sources. The generated fields are parallel to the shock front and fluctuate on the very short scale of the plasma skin depth. The synchrotron radiation emitted from the limb-brightened source image is linearly polarized in the radial direction relative to the source center. Although the net polarization vanishes under circular symmetry, GRB sources should exhibit polarization scintillations as their radio afterglow radiation gets scattered by the Galactic interstellar medium. Detection of polarization scintillations could therefore test the above mechanism for magnetic field generation.

Thermodynamic and spectral properties of a photon gas in $e^+e^-\gamma$ plasmas are studied. The effect of a finite effective mass of a photon, associated with the plasma frequency cutoff, is self-consistently included. In the ultra-relativistic plasma, the photon spectrum turns out to be universal with the temperature normalized plasma frequency cutoff being a fundamental constant independent of plasma parameters. Such a universality does not hold in the non-relativistic plasma.

The magnetic energy of nonlinear Alfven waves in compressible plasmas may be ponderomotively coupled only to ion-acoustic quasi-modes which modulate the wave phase velocity and cause wave-front steepening. In the collisionless plasma with $\beta\not=0$, the dynamics of nonlinear Alfven wave is also affected by the resonant particle-wave interactions. Upon relatively rapid evolution (compared to the particle bounce time), the quasi-stationary wave structures, identical to the so called (Alfvenic) Rotational Discontinuities, form, the emergence and dynamics of which has not been previously understood. Collisionless (Landau) dissipation of nonlinear Alfven waves is also a plausible and natural mechanism of the solar wind heating. Considering a strong, compressible, Alfvenic turbulence as an ensemble of randomly interacting Alfvenic discontinuities and nonlinear waves, it is shown that there exist two distinct phases of turbulence. What phase realizes depends on whether this collisionless damping is strong enough to provide adequate energy sink at all scales and, thus, to support a steady-state cascade of the wave energy. In long-time asymptotics, however, the particle distribution function is affected by the wave magnetic fields. In this regime of nonlinear Landau damping, resonant particles are trapped in the quasi-stationary Alfvenic discontinuities, giving rise to a formation of a plateau on the distribution function and quenching collisionless damping. Using the virial theorem for trapped particles, it is analytically demonstrated that their effect on the nonlinear dynamics of such discontinuities is non-trivial and forces a significant departure of the theory from the conventional paradigm.

The spatial structure, fluctuations as well as all state probabilities of self-organized (steady) states of cellular automata can be found (almost) exactly and {\em explicitly} from their Markovian dynamics. The method is shown on an example of a natural sand pile model with a gradient threshold.

It is shown that photon--photino inter-conversions (if exist) may result in a detectable CMBR spectrum distortion which amplitude depends on photino properties, such as its mass. An upper bound on the distortion parameter determined from the recent COBE-FIRAS data, hence, sets a lower bound on the photino mass, $m_{\tilde\gamma} \gtrsim 300 eV$.

A fully nonlinear, time-asymptotic theory of resonant particle trapping in large-amplitude quasi-parallel Alfven waves is presented. The effect of trapped particles on the nonlinear dynamics of quasi-stationary Alfvenic discontinuities and coherent Alfven waves is highly non-trivial and forces to a significant departure of the theory from the conventional DNLS and KNLS equation models. The virial theorem is used to determine the time-asymptotic distribution function.

A simple model collisionless, dissipative, compressible MHD (Alfvenic) turbulence in a magnetized system is investigated. In contrast to more familiar paradigms of turbulence, dissipation arises from Landau damping, enters via nonlinearity, and is distributed over all scales. The theory predicts that two different regimes or phases of turbulence are possible, depending on the ratio of steepening to damping coefficient (m_1/m_2). For strong damping (|m_1/m_2|<1), a regime of smooth, hydrodynamic turbulence is predicted. For |m_1/m_2|>1, steady state turbulence does not exist in the hydrodynamic limit. Rather, spikey, small scale structure is predicted.

The nonlinear dynamics of collisionless Alfven trains, including resonant particle effects is studied using the kinetic nonlinear Schroedinger (KNLS) equation model. Numerical solutions of the KNLS reveal the dynamics of Alfven waves to be sensitive to the sense of polarization as well as the angle of propagation with respect to the ambient magnetic field. The combined effects of both wave nonlinearity and Landau damping result in the evolutionary formation of stationaryOA S- and arc-polarized directional and rotational discontinuities. These waveforms are freqently observed in the interplanetary plasma.

Collisionless regime kinetic models for coherent nonlinear Alfven wave dynamics are studied using fluid moment equations with an approximate closure anzatz. Resonant particle effects are modelled by incorporating an additional term representing dissipation akin to parallel heat conduction. Unlike collisional dissipation, parallel heat conduction is presented by an integral operator. The modified derivative nonlinear Schrodinger equation thus has a spatially nonlocal nonlinear term describing the long-time evolution of the envelope of parallel-propagating Alfven waves, as well. Coefficients in the nonlinear terms are free of the 1/(1-beta) singularity usually encountered in previous analyses, and have very a simple form which clarifies the physical processes governing the large amplitude Alfvenic nonlinear dynamics. The nonlinearity appears via coupling of an Alfvenic mode to a kinetic ion-acoustic mode. Damping of the nonlinear Alfven wave appears via strong Landau damping of the ion-acoustic wave when the electron-to-ion temperature ratio is close to unity. For a (slightly) obliquely propagating wave, there are finite Larmor radius corrections in the dynamical equation. This effect depends on the angle of wave propagation relative to B_0 and vanishes for the limit of strictly parallel propagation. Explicit magnetic perturbation envelope equations amenable to further analysis and numerical solution are obtained. Implications of these models for collisionless shock dynamics are discussed.

The influence of various kinetic effects (e.g. Landau damping, diffusive and collisional dissipation, and finite Larmor radius terms) on the nonlinear evolution of finite amplitude Alfvenic wave trains in a finite-beta environment is systematically investigated using a novel, kinetic nonlinear Schrodinger (KNLS) equation. The dynamics of Alfven waves is sensitive to the sense of polarization as well as the angle of propagation with respect to the ambient magnetic field. Numerical solution for the case with Landau damping reveals the formation of dissipative structures, which are quasi-stationary, S-polarized directional (and rotational) discontinuities which self-organize from parallel propagating, linearly polarized waves. Parallel propagating circularly polarized packets evolve to a few circularly polarized Alfven harmonics on large scales. Stationary arc-polarized rotational discontinuities form from obliquely propagating waves. Collisional dissipation, even if weak, introduces enhanced wave damping when beta is very close to unity. Cyclotron motion effects on resonant particle interactions introduce cyclotron resonance into the nonlinear Alfven wave dynamics.

The saturation of ion cyclotron Alfven turbulence excited by beam particles is investigated using resonance broadening theory. The stochastic scattering which decorrelates particles, includes both random acceleration by electric fields and a turbulent magnetic mirroring effect. Turbulent mirroring is shown to yield non-Gaussian corrections to the orbits even if the random electric and magnetic fields are Gaussian. The predicted steady-state turbulence level exhibits a peaked anglular distribution, with a maximum near Theta ~ 60 degrees.