Search

Nonlinear evolution equation for Alfven waves, propagating in streaming plasmas with nonuniform densities and inhomogeneous magnetic fields, is obtained by using the reductive perturbation technique.

A fully nonlinear, time-asymptotic theory of nonlinear Landau damping and resonant particle trapping in finite-amplitude waves is presented. The virial theorem and the conservation of the parallel adiabatic invariant are used to determine the time-asymptotic distribution function. The effect of trapped particles on the nonlinear wave dynamics is highly nontrivial and forces a significant departure from the conventional models of finite-amplitude waves.

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., 6 pages, 4 figures

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.

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 give a theoretical summary of 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.

The radial structure of tokamak profiles determined by anomalous transport is elucidated by studying the statistical mechanics of a sand pile automaton for which the toppling conditions depend on local gradient, alone. In this representation, the sand pile dynamics is Markovian, and spatial profiles may be obtained from calculated expectation values of the local gradient. The Markovian structure of the dynamics is exploited to analytically derive a local gradient probability distribution function from a generalized kinetic equation. For homogeneous, weak noise, the calculated expectation value of the gradient is well below the marginally stable state. In the over‐driven limit (i.e., strong homogeneous noise), a region of super‐critical gradient is shown to form near the bottom of the pile. For the case of localized noise, the mean self‐organized profile is always sub‐critical. These results are consistent with numerical studies of simple automata. Their relevance to and implications for tokamak confinement are discussed.

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$)., Comment: 6 pages, 6 figures, presented at Dynamical processes in space plasmas II, Isradinamic 2012, in press, ANGEO. arXiv admin note: text overlap with arXiv:1303.2569

In this article existing key points detection and matching methods are observed. The new wavelet transformation based key point detection algorithm is proposed and the descriptor creation is implemented.

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., Comment: 16 pages, 8 figures. Accepted for publication in JCAP

The ideal magnetohydrodynamic ballooning stability of a poloidally rotating plasma in a sheared radial electric field is investigated for an axisymmetric toroidal configuration. An analytical criterion for the ballooning instability of such a system is derived. Electric field shearing is the dominant effect. The sign of the electric field gradient enters the stability criterion via the angular momentum gradient a la Taylor–Couette flow. Thus, Er’

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., 18 pages, 2 figures

Plasma instabilities excited in collisionless shocks are responsible for particle acceleration, generation of magnetic fields , and associated radiation. We have investigated the particle acceleration and shock structure associated with an unmagnetized relativistic jet propagating into an unmagnetized plasma. Cold jet electrons are thermalized and slowed while the ambient electrons are swept up to create a partially developed hydrodynamic-like shock structure. The shock structure depends on the composition of the jet and ambient plasma (electron-positron or electron-ions). Strong electromagnetic fields are generated in the reverse , jet shock and provide an emission site. These magnetic fields contribute to the electron's transverse deflection behind the shock. We have calculated, self-consistently, the radiation from electrons accelerated in the turbulent magnetic fields. We found that the synthetic spectra depend on the Lorentz factor of the jet, its thermal temperature and strength of the generated magnetic fields. The detailed properties of the radiation are important for understanding the complex time evolution and/or spectral structure in gamma-ray bursts, relativistic jet shocks, and supernova remnants

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., Comment: 8 pages, 1 figure

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., 17 pages, 4 figures, submitted to Phys. Plasmas

Plasma filamentation is often encountered in collisionless shocks and inertial confinement fusion. We develop a general analytical description of the two-dimensional relativistic filamentary equilibrium and derive the conditions for existence of potential-free equilibria. A pseudopotential equation for the vector-potential is constructed for cold and relativistic Maxwellian distributions. The role of counter-streaming is explained. We present single current sheet and periodic current sheet solutions, and analyze the equilibria with electric potential. These solutions can be used to study linear and nonlinear evolution of the relativistic filamentation instability.

Prompt GRB emission can carry information about local conditions in a jet and provide a link toward understanding of a GRB progenitor. Conventional interpretation of the prompt emission, regardless of the jet origin, is either synchrotron, perhaps with the (self)‐comptonized component, or thermal. Whereas some GRBs are consistent with this interpretation, many are at odds with it. Particularly, the ubiquitous spectral variability and correlations still awaits its physical interpretation. In view of the recent discovery of the generation of small‐scale magnetic fields during reconnection in pair plasmas and the claimed non‐discovery of synchrotron‐violating spectra in PIC simulations of relativistic shocks the reassessment of physical models of GRB jets and its radiation becomes imperative. Here we argue that the conventional assumption of the isotropy (in the co‐moving frame) of the produced radiation pattern is too simplistic to explain the wealth of observational data. We show that the anisotropic spectra fit observations much better. The anisotropy of radiation can be due to the ambient (progenitor) magnetic field, it can also be induced by the shock dynamics and geometry, or reflect the anisotropy of the post‐shock turbulence. Interestingly, different emission mechanisms (e.g., synchrotron, small‐pitch‐angle, jitter, thermal) can dominate in different environments (e.g., magnetic, baryonic or leptonic jets) and result in different spectro‐temporal characteristics. We show how the corresponding spectral correlations can distinguish between the jet origin and composition or as least rule our certain models.

This is the publisher's version, also available electronically from http://onlinelibrary.wiley.com.

International audience; The filamentation instability of counterstreaming beams is a nonresonant hydrodynamic-type instability whose growth rate is a smooth function of the wavelength (scale). As a result, perturbations with all unstable wavelengths develop, and the growth saturates due to the saturation of available current. For a given scale, the magnetic field at saturation is proportional to the scale. As a result, the instability develops in a nearly linear regime, where the unstable modes stop growing as soon as the saturation of the corresponding wavelength is reached. At each moment there exists a dominant scale of the magnetic field which is the scale that reached saturation at this particular time. The smaller scales do not disappear and can be easily distinguished in the current structure. The overall growth of the instability stops when the loss of the streaming ion energy because of deceleration is comparable to the initial ion energy.

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., Subm. to proceedings of the Annual International Astrophysics Conference (AIAC-8), Hawaii, 2009

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., Comment: 6 pages; submitted to proc. of Huntsville 2008 symposium on GRBs

Thermodynamic and spectral properties of a photon gas in ${e}^{+}{e}^{\ensuremath{-}}\ensuremath{\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 ultrarelativistic plasma, the photon spectrum turns out to be universal with the temperature normalized plasma frequency cutoff as a fundamental constant independent of plasma parameters. Such a universality does not hold in the nonrelativistic plasma.

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., 15 pages, submitted to ApJ

[1] 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 normally to the galactic plane, thereby inducing an observed terrestrial periodicity in biodiversity. There are a number of plausible mechanisms by which cosmic rays might affect terrestrial biodiversity. Here we investigate one of these mechanisms, the consequent ionization and dissociation in the atmosphere, resulting in changes in atmospheric chemistry, which culminate in the depletion of ozone and a resulting increase in the dangerous solar UVB flux on the ground. We estimate the enhancement of cosmic-ray intensity for a range of reasonable parameters of the galactic wind and galactic magnetic field, and use these to compute steady-state atmospheric effects. At the lower end of this range, we find that the effects are far too small to be of serious consequence. At the upper end of this range, the level of ozone depletion approaches that currently experienced due to anthropogenic effects such as accumulated chlorofluorocarbons, i.e., ∼2.1% global average loss of ozone column density. We discuss some of the possible effects. While much smaller intensity than the atmospheric effects of a nearby galactic gamma-ray burst, the duration of the effects would be about 106 times greater. Current ozone depletion is a documented stress on the biosphere; it is not clear whether its consequences would be severe if of extended duration. We conclude that, for estimates at the upper end of the reasonable range of the cosmic-ray variability over geologic time, the mechanism of atmospheric ozone depletion may provide a small additional stress, enhancing the impact of other events. However, in order to account for large fluctuations in biodiversity correlated with cosmic-ray flux, other mechanisms should be investigated.

We review the models of Gamma‐Ray Bursts (GRBs) with relativistic baryon‐dominated ejecta: the Weibel‐jitter model and the vorticity model. The former may also be relevant to the magnetic‐dominated outflows. We discuss the origin of magnetic fields, mechanisms of electron heating/acceleration, nature and properties of the emitted radiation and the relevance of computer modeling to realistic GRBs. We discuss observational effects and predictions that, hopefully, will allow us to discriminate between the models.

Radiation from GRBs, either in the prompt phase or seen as flares or an afterglow, is believed to be produced at shocks by accelerated electrons in magnetic fields. However, there are alternatives, at least for the prompt and flares, namely the magnetically dominated and Poynting flux driven outflows, where no shocks presumably form at all. The outflow properties are intimately linked to the central engine/progenitor; hence, to understand the nature of the outflow is of great importance. Here we review the present status of the “standard model of GRBs”—the model with internal/external shocks. We briefly discuss two alternative models of the magnetic field generation at shocks: the model based on the Weibel instability and the recently proposed one based on the generation of vorticity at a shock. We also address the electron energization/acceleration mechanisms (including non‐Fermi) in both models. The observational predictions for the spectral and temporal dynamics of the produced radiation have been elab...

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., 16 pages, 1 figure

The origin of the micro-Gauss magnetic fields in galaxy clusters is one of the outstanding problem of modern cosmology. We suggest that these fields could have been produced via the Weibel instability operating at shocks during the Large-Scale Structure formation and at accretion shocks in Galaxy Clusters. We have performed three-dimensional particle-in-cell simulations of the nonrelativistic Weibel instability in an electron-proton plasma (with the ion-to-electron mass ratio of 100), 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 discuss whether and how our results may change with the change the mass ratio to the realistic value of 1836. We stress that even if the Weibel-generated small-scale magnetic fields decay with time, they can serve as seed fields that can be further amplified and inverse-cascade to larger scales by turbulent motions of post-shock MHD turbulence.

The existence of a critical gradient is an important feature of many transport models. [M. Kotschenreuther et al., Phys. Plasmas 2, 2381 (1995)]. However, fundamental differences exist in the dynamics near marginal stability, depending on whether the transport phenomena are controlled by strict (linear) marginal stability or by a self-organized criticality. One of the most striking differences is in the stiffness of the profiles. In this paper, a sand pile model is used to gain some basic understanding of the stiffness of the profile under different dynamics.

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., Comment: 4 pages, 2 color figures

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., 30 pages, to be published in ApJ. Replaced for conformity with published version, including correction of minor typos and updated references

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., Comment: 4 pages, 2 figs; ApJL accepted

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., Comment: 4 pages, Proceedings of "Swift-05" meeting, Washington, DC

X‐rays are produced throughout the heliosphere as a consequence of charge transfer collisions between heavy solar wind ions and neutral atoms. After such a collision the solar wind ion is left in a highly excited state and emits extreme ultraviolet and soft X‐ray photons. In the outer heliosphere, solar wind charge exchange X‐ray emission is mainly due to charge exchange with neutral interstellar hydrogen. We have combined MHD simulations with a comprehensive charge exchange computation code. We trace the full evolution of solar wind ions along stream line in order to produce three‐dimensional emissivities and, subsequently, two‐dimensional X‐ray brightness and spectral maps of the heliosphere as would be observed from the outside. The model treats both the collisionally thin and the collisionally thick regimes. This model can be a diagnostic tool for studying stellar wind properties of nearby Sun‐like stars.

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., Comment: Minor revisions, includes more discussion of reaction rates; 1 new figure. Accepted for publication in ApJ. 68 pages, 21 figures, 12 tables

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., 6 pages, 5 figures. Review talk at GRB-03 conference (Santa Fe)

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., 14 pages, 1 figure. accepted for publication in ApJ

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.., Comment: 9 pages, 4 figures (jpg), accepted for publication in ApJLett

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., 7 pages, 1 figure. Submitted to ApJ

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., Comment: 28 pages, 12 figures, in press, ApJ, September 10, 2014

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., 4 pages, 2 figures. Accepted for publication in ApJ Letters

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., Comment: 5 pages, 2 figures. To appear in proceedings of the Gamma 2001 symposium

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$, 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., Comment: 3 pages with 2 eps figures, aipproc.sty. To appear in Proceedings of the 20th Texas Symposium on Relativistic Astrophysics, Austin, Texas, 2000, edited by J. Craig Wheeler and Hugo Martel (American Institute of Physics)

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., Comment: 5 pages, 4 figures. To appear in proceedings of the Gamma 2001 symposium

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., Comment: 3 pages with 1 eps figure, aipproc.cls. To appear in Proceedings of the 20th Texas Symposium on Relativistic Astrophysics, Austin, Texas, 2000, edited by J. Craig Wheeler and Hugo Martel (American Institute of Physics)

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., Comment: 7 pages, 2 eps figures. Accepted to ApJL

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), 5 pages, 2 figures, emilateapj.sty. Submitted to ApJL

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., 3 pages with 2 eps figures, aipproc.sty. To appear in Proceedings of the 20th Texas Symposium on Relativistic Astrophysics, Austin, Texas, 2000, edited by J. Craig Wheeler and Hugo Martel (American Institute of Physics)

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., 12 pages, 2 figures included. Submitted to ApJ

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, 13 pages, 7 figures, uses emulateapj.sty. Submitted to ApJ