Electric field generation by the electron beam filamentation instability: Filament size effects

The filamentation instability (FI) of counter-propagating beams of electrons is modelled with a particle-in-cell simulation in one spatial dimension and with a high statistical plasma representation. The simulation direction is orthogonal to the beam velocity vector. Both electron beams have initially equal densities, temperatures and moduli of their nonrelativistic mean velocities. The FI is electromagnetic in this case. A previous study of a small filament demonstrated, that the magnetic pressure gradient force (MPGF) results in a nonlinearly driven electrostatic field. The probably small contribution of the thermal pressure gradient to the force balance implied, that the electrostatic field performed undamped oscillations around a background electric field. Here we consider larger filaments, which reach a stronger electrostatic potential when they saturate. The electron heating is enhanced and electrostatic electron phase space holes form. The competition of several smaller filaments, which grow simultaneously with the large filament, also perturbs the balance between the electrostatic and magnetic fields. The oscillations are damped but the final electric field amplitude is still determined by the MPGF.Comment: 14 pages, 10 plots, accepted for publication in Physica Script

PIC Simulations of the Temperature Anisotropy-Driven Weibel Instability: Analyzing the perpendicular mode

An instability driven by the thermal anisotropy of a single electron species is investigated in a 2D particle-in-cell (PIC) simulation. This instability is the one considered by Weibel and it differs from the beam driven filamentation instability. A comparison of the simulation results with analytic theory provides similar exponential growth rates of the magnetic field during the linear growth phase of the instability. We observe in accordance with previous works the growth of electric fields during the saturation phase of the instability. Some components of this electric field are not accounted for by the linearized theory. A single-fluid-based theory is used to determine the source of this nonlinear electric field. It is demonstrated that the magnetic stress tensor, which vanishes in a 1D geometry, is more important in this 2-dimensional model used here. The electric field grows to an amplitude, which yields a force on the electrons that is comparable to the magnetic one. The peak energy density of each magnetic field component in the simulation plane agrees with previous estimates. Eddy currents develop, which let the amplitude of the third magnetic field component grow, which is not observed in a 1D simulation.Comment: accepted by Plasma Physics and Controlled Fusio

Particle-in-cell simulation of a mildly relativistic collision of an electron-ion plasma carrying a quasi-parallel magnetic field: Electron acceleration and magnetic field amplification at supernova shocks

Plasma processes close to SNR shocks result in the amplification of magnetic fields and in the acceleration of electrons, injecting them into the diffusive acceleration mechanism. The acceleration of electrons and the B field amplification by the collision of two plasma clouds, each consisting of electrons and ions, at a speed of 0.5c is investigated. A quasi-parallel guiding magnetic field, a cloud density ratio of 10 and a plasma temperature of 25 keV are considered. A quasi-planar shock forms at the front of the dense plasma cloud. It is mediated by a circularly left-hand polarized electromagnetic wave with an electric field component along the guiding magnetic field. Its propagation direction is close to that of the guiding field and orthogonal to the collision boundary. It has a low frequency and a wavelength that equals several times the ion inertial length, which would be indicative of a dispersive Alfven wave close to the ion cyclotron resonance frequency of the left-handed mode (ion whistler), provided that the frequency is appropriate. However, it moves with the super-alfvenic plasma collision speed, suggesting that it is an Alfven precursor or a nonlinear MHD wave such as a Short Large-Amplitude Magnetic Structure (SLAMS). The growth of the magnetic amplitude of this wave to values well in excess of those of the quasi-parallel guiding field and of the filamentation modes results in a quasi-perpendicular shock. We present evidence for the instability of this mode to a four wave interaction. The waves developing upstream of the dense cloud give rise to electron acceleration ahead of the collision boundary. Energy equipartition between the ions and the electrons is established at the shock and the electrons are accelerated to relativistic speeds.Comment: 16 pages, 18 figures, Accepted for publication by Astron & Astrophy

Oligomorphic Dynamics for Analyzing the Quantitative Genetics of Adaptive Speciation

Ecological interaction, including competition for resources, often causes frequency-dependent disruptive selection, which, when accompanied by reproductive isolation, may act as driving forces of adaptive speciation. While adaptive dynamics models have added new perspectives to our understanding of the ecological dimensions of speciation processes, it remains an open question how best to incorporate and analyze genetic detail in such models. Conventional approaches, based on quantitative genetics theory, typically assume a unimodal character distribution and examine how its moments change over time. Such approximations inevitably fail when a character distribution becomes multimodal. Here, we propose a new approximation, oligomorphic dynamics, to the quantitative genetics of populations that include several morphs and that therefore exhibit multiple peaks in their character distribution. To this end we first decompose the character distribution into a sum of unimodal distributions corresponding to individual morphs. Characterizing these morphs by their frequency (fraction of individuals belonging to each morph), position (mean character of each morph), and width (standard deviation of each morph) we then derive the coupled eco-evolutionary dynamics of morphs through a double Taylor expansion. We show that the demographic, convergence, and evolutionary stability of a populations character distribution correspond, respectively, to the asymptotic stability of frequencies, positions, and widths under the oligomorphic dynamics introduced here. As first applications of oligomorphic dynamics theory we analytically derive the effects (a) of the strength of disruptive selection on waiting times until speciation, (b) of mutation on conditions for speciation, and (c) of the fourth moments of competition kernels on patterns of speciation

Resource heterogeneity can facilitate cooperation

Although social structure is known to promote cooperation, by locally exposing selfish agents to their own deeds, studies to date assumed that all agents have access to the same level of resources. This is clearly unrealistic. Here we find that cooperation can be maintained when some agents have access to more resources than others. Cooperation can then emerge even in populations in which the temptation to defect is so strong that players would act fully selfishly if their resources were distributed uniformly. Resource heterogeneity can thus be crucial for the emergence and maintenance of cooperation. We also show that resource heterogeneity can hinder cooperation once the temptation to defect is significantly lowered. In all cases, the level of cooperation can be maximized by managing resource heterogeneity

The filamentation instability driven by warm electron beams: Statistics and electric field generation

The filamentation instability of counterpropagating symmetric beams of electrons is examined with 1D and 2D particle-in-cell (PIC) simulations, which are oriented orthogonally to the beam velocity vector. The beams are uniform, warm and their relative speed is mildly relativistic. The dynamics of the filaments is examined in 2D and it is confirmed that their characteristic size increases linearly in time. Currents orthogonal to the beam velocity vector are driven through the magnetic and electric fields in the simulation plane. The fields are tied to the filament boundaries and the scale size of the flow-aligned and the perpendicular currents are thus equal. It is confirmed that the electrostatic and the magnetic forces are equally important, when the filamentation instability saturates in 1D. Their balance is apparently the saturation mechanism of the filamentation instability for our initial conditions. The electric force is relatively weaker but not negligible in the 2D simulation, where the electron temperature is set higher to reduce the computational cost. The magnetic pressure gradient is the principal source of the electrostatic field, when and after the instability saturates in the 1D simulation and in the 2D simulation.Comment: 10 pages, 6 figures, accepted by the Plasma Physics and Controlled Fusion (Special Issue EPS 2009

Insights from unifying modern approximations to infections on networks

Networks are increasingly central to modern science owing to their ability to conceptualize multiple interacting components of a complex system. As a specific example of this, understanding the implications of contact network structure for the transmission of infectious diseases remains a key issue in epidemiology. Three broad approaches to this problem exist: explicit simulation; derivation of exact results for special networks; and dynamical approximations. This paper focuses on the last of these approaches, and makes two main contributions. Firstly, formal mathematical links are demonstrated between several prima facie unrelated dynamical approximations. And secondly, these links are used to derive two novel dynamical models for network epidemiology, which are compared against explicit stochastic simulation. The success of these new models provides improved understanding about the interaction of network structure and transmission dynamics

How large can the electron to proton mass ratio be in Particle-In-Cell simulations of unstable systems?

Particle-in-cell (PIC) simulations are widely used as a tool to investigate instabilities that develop between a collisionless plasma and beams of charged particles. However, even on contemporary supercomputers, it is not always possible to resolve the ion dynamics in more than one spatial dimension with such simulations. The ion mass is thus reduced below 1836 electron masses, which can affect the plasma dynamics during the initial exponential growth phase of the instability and during the subsequent nonlinear saturation. The goal of this article is to assess how far the electron to ion mass ratio can be increased, without changing qualitatively the physics. It is first demonstrated that there can be no exact similarity law, which balances a change of the mass ratio with that of another plasma parameter, leaving the physics unchanged. Restricting then the analysis to the linear phase, a criterion allowing to define a maximum ratio is explicated in terms of the hierarchy of the linear unstable modes. The criterion is applied to the case of a relativistic electron beam crossing an unmagnetized electron-ion plasma.Comment: To appear in Physics of Plasma