Monte Carlo methods are often employed to numerically integrate kinetic equations, such as the particle-in-cell method for the plasma kinetic equation, but these methods suffer from the introduction of counting noise to the solution. We report on a cautionary tale of counting noise modifying the nonlinear saturation of kinetic instabilities driven by unstable beams of plasma. We find a saturated magnetic field in under-resolved particle-in-cell simulations due to the sampling error in the current density. The noise-induced magnetic field is anomalous, as the magnetic field damps away in continuum kinetic and increased particle count particle-in-cell simulations. This modification of the saturated state has implications for a broad array of astrophysical phenomena beyond the simple plasma system considered here, and it stresses the care that must be taken when using particle methods for kinetic equations.

We present the first 2X2V continuum Vlasov-Maxwell simulations of interpenetrating, unmagnetized plasmas to study the competition between two-stream, Oblique, and filamentation modes in the weakly relativistic regime. We find that after nonlinear saturation of the fastest-growing two-stream and Oblique modes, the effective temperature anisotropy, which drives current filament formation via the secular Weibel instability, has a strong dependence on the internal temperature of the counter-streaming plasmas. The effective temperature anisotropy is significantly more reduced in colder than in hotter plasmas, leading to orders of magnitude lower magnetization for colder plasmas. A strong dependence of the energy conversion efficiency of Weibel-type instabilities on internal beam temperature has implications for determining their contribution to the observed magnetization of many astrophysical and laboratory plasmas.