Impact of Resonant Magnetic Perturbations on Zonal Modes, Drift-Wave Turbulence and the L-H Transition Threshold

We study the effects of Resonant Magnetic Perturbations (RMPs) on turbulence, flows and confinement in the framework of resistive drift-wave turbulence. This work was motivated, in parts, by experiments reported at the IAEA 2010 conference [Y. Xu {\it et al}, Nucl. Fusion \textbf{51}, 062030] which showed a decrease of long-range correlations during the application of RMPs. We derive and apply a zero-dimensional predator-prey model coupling the Drift-Wave Zonal Mode system [M. Leconte and P.H. Diamond, Phys. Plasmas \textbf{19}, 055903] to the evolution of mean quantities. This model has both density gradient drive and RMP amplitude as control parameters and predicts a novel type of transport bifurcation in the presence of RMPs. This model allows a description of the full L-H transition evolution with RMPs, including the mean sheared flow evolution. The key results are: i) The L-I and I-H power thresholds \emph{both} increase with RMP amplitude |\bx|, the relative increase of the L-I threshold scales as \Delta P_{\rm LI} \propto |\bx|^2 \nu_*^{-2} \gyro^{-2}, where $\nu_*$ is edge collisionality and \gyro is the sound gyroradius. ii) RMPs are predicted to \emph{decrease} the hysteresis between the forward and back-transition. iii) Taking into account the mean density evolution, the density profile - sustained by the particle source - has an increased turbulent diffusion compared with the reference case without RMPs which provides one possible explanation for the \emph{density pump-out} effect.Comment: 30 pages, IAEA-based articl

Plus Charge Prevalence in Cosmic Rays: Room for Dark Matter in the Positron Spectrum

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

Breaking Kelvin: Circulation conservation and vortex breakup in MHD at low Magnetic Prandtl Number

In this paper we examine the role of weak magnetic fields in breaking Kelvin's circulation theorem and in vortex breakup in two-dimensional magnetohydrodynamics for the physically important case of a low magnetic Prandtl number (low $Pm$) fluid. We consider three canonical inviscid solutions for the purely hydrodynamical problem, namely a Gaussian vortex, a circular vortex patch and an elliptical vortex patch. We examine how magnetic fields lead to an initial loss of circulation $\Gamma$ and attempt to derive scaling laws for the loss of circulation as a function of field strength and diffusion as measured by two non-dimensional parameters. We show that for all cases the loss of circulation depends on the integrated effects of the Lorentz force, with the patch cases leading to significantly greater circulation loss. For the case of the elliptical vortex the loss of circulation depends on the total area swept out by the rotating vortex and so this leads to more efficient circulation loss than for a circular vortex.Comment: 21 pages, 12 figure

Modern theory of Fermi acceleration: a new challenge to plasma physics

One of the main features of astrophysical shocks is their ability to accelerate particles to extremely high energies. The leading acceleration mechanism, the diffusive shock acceleration is reviewed. It is demonstrated that its efficiency critically depends on the injection of thermal plasma into acceleration which takes place at the subshock of the collisionless shock structure that, in turn, can be significantly smoothed by energetic particles. Furthermore, their inhomogeneous distribution provides free energy for MHD turbulence regulating the subshock strength and injection rate. Moreover, the MHD turbulence confines particles to the shock front controlling their maximum energy and bootstrapping acceleration. Therefore, the study of the MHD turbulence in a compressive plasma flow near a shock is a key to understanding of the entire process. The calculation of the injection rate became part of the collisionless shock theory. It is argued that the further progress in diffusive shock acceleration theory is impossible without a significant advance in these two areas of plasma physics.Comment: 12 pages, 4 figures, invited talk at APS/ICPP, Quebec 2000, to appear in Phys. of Plasma

Critical self-organization of astrophysical shocks

There are two distinct regimes of the first order Fermi acceleration at shocks. The first is a linear (test particle) regime in which most of the shock energy goes into thermal and bulk motion of the plasma. The second is an efficient regime when it goes into accelerated particles. Although the transition region between them is narrow, we identify the factors that drive the system to a {\it self-organized critical state} between those two. Using an analytic solution, we determine this critical state and calculate the spectra and maximum energy of accelerated particles.Comment: To appear in ApJL, Sec.3 extensively rewritten, 4 pages, Latex, emulateapj.sty, eps

Polarization morphology of SiO masers in the circumstellar envelope of the AGB star R Cassiopeiae

Silicon monoxide maser emission has been detected in the circumstellar envelopes of many evolved stars in various vibrationally-excited rotational transitions. It is considered a good tracer of the wind dynamics close to the photosphere of the star. We have investigated the polarization morphology in the circumstellar envelope of an AGB star, R Cas. We mapped the linear and circular polarization of SiO masers in the v=1, J=1-0 transition. The linear polarization is typically a few tens of percent while the circular polarization is a few percent. The fractional polarization tends to be higher for emission of lower total intensity. We found that, in some isolated features the fractional linear polarization appears to exceed 100%. We found the Faraday rotation is not negligible but is ~15 deg., which could produce small scale structure in polarized emission whilst total intensity is smoother and partly resolved out. The polarization angles vary considerably from feature to feature but there is a tendency to favour the directions parallel or perpendicular to the radial direction with respect to the star. In some features, the polarization angle abruptly flips 90 deg. We found that our data are in the regime where the model of Goldreich et al (1973) can be applied and the polarization angle flip is caused when the magnetic field is at close to 55 deg. to the line of sight. The polarization angle configuration is consistent with a radial magnetic field although other configurations are not excluded.Comment: 14 pages, 15 figures. Accepted for publication in MNRA

Fluid Models for Kinetic Effects on Coherent Nonlinear Alfven Waves. II. Numerical Solutions

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.Comment: 38 pages (including 23 figures and 1 table