1,184 research outputs found
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 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 () 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
- …