Supporting data: Suppressing grid instability and noise in particle-in-cell simulation by smoothing

This archive contains data supporting the research described in the paper, "Suppressing grid instability and noise in particle-in-cell simulation by smoothing," as well as scripts to reproduce the data and analyze the data.

Description of files
--------------------

- finiteDifferenceCauchySqrSmeared3ptPoissonPlasmaDispersion.nb
  A Mathematica notebook that calculates roots of the finite difference
  plasma dispersion.

- pbEsPtclSim1d.py
  A Python script that runs 1D MC-PIC simulations with periodic boundaries.

- figN.zip
  contains a directory figN/ with the data needed for figN,
  as well as a Python script figNplot*.py that will produce pdfs for
  Figure N.  These scripts demonstrate how to extract desired data from
  the raw data.

For fig2, the data the in cauchySqr*.py files output by the
  finiteDifferenceCauchySqrSmeared3ptPoissonPlasmaDispersion.nb
  notebook.  fig2plotGrowthRates.py is the script that reads those
  data files and produces plots.

For fig3-8, the data was produced by pbEsPtclSim1d.py and
  is in sub-directories like
  pbcVdOverOmegapDx0p1debyeOverDx0p02Nx100ppc10000Tp1000dtf0p1stab0Nwaves0amp0/
  The data output is in hdf5 format.
  The pbEsPtclSim1d.py script can be consulted to determine precisely
  what information is stored in each file.
  Roughly:
  - sim.h5 contains parameters describing the simulation -- number of cells,
    Debye length resolution, etc.
  - History.h5 contains multiple datasets representing quantities varying
    in time.  E.g., it contains the total particle energy vs. time.