  Collisionless plasmas exhibit nonthermal particle distributions after being
energized; as a consequence, they enter a state of low Boltzmann-Gibbs (BG)
entropy relative to the thermal state. The Vlasov equations predict that in a
collisionless plasma with closed boundaries, BG entropy is formally conserved,
along with an infinite set of Casimir invariants; this provides a seemingly
strong constraint that may explain how plasmas maintain low entropy.
Nevertheless, it is commonly believed that entropy production is enabled by
phase mixing or nonlinear entropy cascades. The question of whether such
anomalous entropy production occurs, and of how to characterize it
quantitatively, is a fundamental problem in plasma physics. We construct a new
theoretical framework for characterizing entropy production (in a generalized
sense) based on ideally conserved dimensional quantities derived from the
Casimir invariants; these are referred to as ``Casimir momenta'' and they
generalize BG entropy. The growth of Casimir momenta relative to the average
particle momentum indicates entropy production. We apply this framework to
quantify entropy production in particle-in-cell simulations of laminar flows
and turbulent flows driven in relativistic plasma, where efficient nonthermal
particle acceleration is enabled. We demonstrate that a large amount of
anomalous entropy is produced by turbulence despite nonthermal features. These
results have implications for reduced modeling of nonthermal particle
acceleration and for diagnosing irreversible dissipation in collisionless
plasmas such as the solar wind and Earth's magnetosphere. Dimensional
representations of generalized entropy analogous to Casimir momenta may be
useful for other problems in statistical physics.
