A Self-Consistent Marginally Stable State for Parallel Ion Cyclotron Waves

We derive an equation whose solutions describe self-consistent states of marginal stability for a proton-electron plasma interacting with parallel-propagating ion cyclotron waves. Ion cyclotron waves propagating through this marginally stable plasma will neither grow nor damp. The dispersion relation of these waves, {\omega} (k), smoothly rises from the usual MHD behavior at small |k| to reach {\omega} = {\Omega}p as k \rightarrow \pm\infty. The proton distribution function has constant phase-space density along the characteristic resonant surfaces defined by this dispersion relation. Our equation contains a free function describing the variation of the proton phase-space density across these surfaces. Taking this free function to be a simple "box function", we obtain specific solutions of the marginally stable state for a range of proton parallel betas. The phase speeds of these waves are larger than those given by the cold plasma dispersion relation, and the characteristic surfaces are more sharply peaked in the v\bot direction. The threshold anisotropy for generation of ion cyclotron waves is also larger than that given by estimates which assume bi-Maxwellian proton distributions.Comment: in press in Physics of Plasma

Three Dimensional Structure and Energy Balance of a Coronal Mass Ejection

The Ultraviolet Coronagraph Spectrometer (UVCS) observed Doppler shifted material of a partial Halo Coronal Mass Ejection (CME) on December 13 2001. The observed ratio of [O V]/O V] is a reliable density diagnostic important for assessing the state of the plasma. Earlier UVCS observations of CMEs found evidence that the ejected plasma is heated long after the eruption. We have investigated the heating rates, which represent a significant fraction of the CME energy budget. The parameterized heating and radiative and adiabatic cooling have been used to evaluate the temperature evolution of the CME material with a time dependent ionization state model. The functional form of a flux rope model for interplanetary magnetic clouds was also used to parameterize the heating. We find that continuous heating is required to match the UVCS observations. To match the O VI-bright knots, a higher heating rate is required such that the heating energy is greater than the kinetic energy. The temperatures for the knots bright in Ly$\alpha$ and C III emission indicate that smaller heating rates are required for those regions. In the context of the flux rope model, about 75% of the magnetic energy must go into heat in order to match the O VI observations. We derive tighter constraints on the heating than earlier analyses, and we show that thermal conduction with the Spitzer conductivity is not sufficient to account for the heating at large heights.Comment: 40 pages, 16 figures, accepted for publication in ApJ For associated mpeg file, please see https://www.cora.nwra.com/~jylee/mpg/f5.mp

Suppression of energetic electron transport in flares by double layers

During flares and coronal mass ejections, energetic electrons from coronal sources typically have very long lifetimes compared to the transit times across the systems, suggesting confinement in the source region. Particle-in-cell simulations are carried out to explore the mechanisms of energetic electron transport from the corona to the chromosphere and possible confinement. We set up an initial system of pre-accelerated hot electrons in contact with ambient cold electrons along the local magnetic field, and let it evolve over time. Suppression of transport by a nonlinear, highly localized electrostatic electric field (in the form of a double layer) is observed after a short phase of free-streaming by hot electrons. The double layer (DL) emerges at the contact of the two electron populations. It is driven by an ion-electron streaming instability due to the drift of the back-streaming return current electrons interacting with the ions. The DL grows over time and supports a significant drop in temperature and hence reduces heat flux between the two regions that is sustained for the duration of the simulation. This study shows transport suppression begins when the energetic electrons start to propagate away from a coronal acceleration site. It also implies confinement of energetic electrons with kinetic energies less than the electrostatic energy of the DL for the DL lifetime, which is much longer than the electron transit time through the source region

Self-consistent nonlinear kinetic simulations of the anomalous Doppler instability of suprathermal electrons in plasmas

Suprathermal tails in the distributions of electron velocities parallel to the magnetic field are found in many areas of plasma physics, from magnetic confinement fusion to solar system plasmas. Parallel electron kinetic energy can be transferred into plasma waves and perpendicular gyration energy of particles through the anomalous Doppler instability (ADI), provided that energetic electrons with parallel velocities v ≥ (ω + Ωce )/k are present; here Ωce denotes electron cyclotron frequency, ω the wave angular frequency and k the component of wavenumber parallel to the magnetic field. This phenomenon is widely observed in tokamak plasmas. Here we present the first fully self-consistent relativistic particle-in-cell simulations of the ADI, spanning the linear and nonlinear regimes of the ADI. We test the robustness of the analytical theory in the linear regime and follow the ADI through to the steady state. By directly evaluating the parallel and perpendicular dynamical contributions to j · E in the simulations, we follow the energy transfer between the excited waves and the bulk and tail electron populations for the first time. We find that the ratio Ωce /(ωpe + Ωce ) of energy transfer between parallel and perpendicular, obtained from linear analysis, does not apply when damping is fully included, when we find it to be ωpe /(ωpe + Ωce ); here ωpe denotes the electron plasma frequency. We also find that the ADI can arise beyond the previously expected range of plasma parameters, in particular when Ωce > ωpe . The simulations also exhibit a spectral feature which may correspond to observations of suprathermal narrowband emission at ωpe detected from low density tokamak plasmas

Dispersion Relations for Bernstein Waves in a Relativistic Pair Plasma

A fully relativistic treatment of Bernstein waves in an electron-positron pair plasma has remained too formidable a task owing to the very complex nature of the problem. In this article, we perform contour integration of the dielectric response function and numerically compute the dispersion curves for a uniform, magnetized, relativistic electron-positron pair plasma. The behavior of the dispersion solution for several cases with different plasma temperatures is highlighted. In particular, we find two wave modes that exist only for large wavelengths and frequencies similar to the cyclotron frequency in a moderately relativistic pair plasma. The results presented here have important implications for the study of those objects where a hot magnetized electron-positron plasma plays a fundamental role in generating the observed radiation.Comment: 8 pages, 8 figures, Accepted for publication by Phys. Rev. E with minor change

Heavy-Fermion Instability in Double-Degenerate Plasmas

In this work we study the propagations of normal frequency modes for quantum hydrodynamic (QHD) waves in the linear limit and introduce a new kind of instability in a double-degenerate plasma. Three different regimes, namely, low, intermediate and high magnetic field strengths are considered which span the applicability of the work to a wide variety of environments. Distinct behavior is observed for different regimes, for instance, in the laboratory-scale field regime no frequency-mode instability occurs unlike those of intermediate and high magnetic-field strength regimes. It is also found that the instability of this kind is due to the heavy-fermions which appear below a critical effective-mass parameter ($\mu_{cr}=\sqrt{3}$) and that the responses of the two (lower and upper frequency) modes to fractional effective-mass change in different effective-mass parameter ranges (below and above the critical value) are quite opposite to each other. It is shown that, the heavy-fermion instability due to extremely high magnetic field such as that encountered for a neutron-star crust can lead to confinement of stable propagations in both lower and upper frequency modes to the magnetic poles. Current study can have important implications for linear wave dynamics in both laboratory and astrophysical environments possessing high magnetic fields

A Coronal Hole's Effects on CME Shock Morphology in the Inner Heliosphere

We use STEREO imagery to study the morphology of a shock driven by a fast coronal mass ejection (CME) launched from the Sun on 2011 March 7. The source region of the CME is located just to the east of a coronal hole. The CME ejecta is deflected away from the hole, in contrast with the shock, which readily expands into the fast outflow from the coronal hole. The result is a CME with ejecta not well centered within the shock surrounding it. The shock shape inferred from the imaging is compared with in situ data at 1 AU, where the shock is observed near Earth by the Wind spacecraft, and at STEREO-A. Shock normals computed from the in situ data are consistent with the shock morphology inferred from imaging.Comment: to appear in The Astrophysical Journa

Nonlinear Electron Oscillations in a Viscous and Resistive Plasma

New non-linear, spatially periodic, long wavelength electrostatic modes of an electron fluid oscillating against a motionless ion fluid (Langmuir waves) are given, with viscous and resistive effects included. The cold plasma approximation is adopted, which requires the wavelength to be sufficiently large. The pertinent requirement valid for large amplitude waves is determined. The general non-linear solution of the continuity and momentum transfer equations for the electron fluid along with Poisson's equation is obtained in simple parametric form. It is shown that in all typical hydrogen plasmas, the influence of plasma resistivity on the modes in question is negligible. Within the limitations of the solution found, the non-linear time evolution of any (periodic) initial electron number density profile n_e(x, t=0) can be determined (examples). For the modes in question, an idealized model of a strictly cold and collisionless plasma is shown to be applicable to any real plasma, provided that the wavelength lambda >> lambda_{min}(n_0,T_e), where n_0 = const and T_e are the equilibrium values of the electron number density and electron temperature. Within this idealized model, the minimum of the initial electron density n_e(x_{min}, t=0) must be larger than half its equilibrium value, n_0/2. Otherwise, the corresponding maximum n_e(x_{max},t=tau_p/2), obtained after half a period of the plasma oscillation blows up. Relaxation of this restriction on n_e(x, t=0) as one decreases lambda, due to the increase of the electron viscosity effects, is examined in detail. Strong plasma viscosity is shown to change considerably the density profile during the time evolution, e.g., by splitting the largest maximum in two.Comment: 16 one column pages, 11 figures, Abstract and Sec. I, extended, Sec. VIII modified, Phys. Rev. E in pres

Attracted Diffusion-Limited Aggregation

In this paper, we present results of extensive Monte Carlo simulations of diffusion-limited aggregation (DLA) with a seed placed on an attractive plane as a simple model in connection with the electrical double layers. We compute the fractal dimension of the aggregated patterns as a function of the attraction strength \alpha. For the patterns grown in both two and three dimensions, the fractal dimension shows a significant dependence on the attraction strength for small values of \alpha, and approaches to that of the ordinary two-dimensional (2D) DLA in the limit of large \alpha. For non-attracting case with \alpha=1, our results in three dimensions reproduce the patterns of 3D ordinary DLA, while in two dimensions our model leads to formation of a compact cluster with dimension two. For intermediate \alpha, the 3D clusters have quasi-2D structure with a fractal dimension very close to that of the ordinary 2D-DLA. This allows one to control morphology of a growing cluster by tuning a single external parameter \alpha.Comment: 6 pages, 6 figures, to appear in Phys. Rev. E (2012