Monte Carlo Sampling of Negative-temperature Plasma States
A Monte Carlo procedure is used to generate N-particle configurations compatible with two-temperature canonical equilibria in two dimensions, with particular attention to nonlinear plasma gyrokinetics. An unusual feature of the problem is the importance of a nontrivial probability density function R0(PHI), the probability of realizing a set {Phi} of Fourier amplitudes associated with an ensemble of uniformly distributed, independent particles. This quantity arises because the equilibrium distribution is specified in terms of {Phi}, whereas the sampling procedure naturally produces particles states gamma; {Phi} and gamma are related via a gyrokinetic Poisson equation, highly nonlinear in its dependence on gamma. Expansion and asymptotic methods are used to calculate R0(PHI) analytically; excellent agreement is found between the large-N asymptotic result and a direct numerical calculation. The algorithm is tested by successfully generating a variety of states of both positive and negative temperature, including ones in which either the longest- or shortest-wavelength modes are excited to relatively very large amplitudes.
- Research Organization:
- Princeton Plasma Physics Lab. (PPPL), Princeton, NJ (United States)
- Sponsoring Organization:
- USDOE Office of Science (US)
- DOE Contract Number:
- AC02-76CH03073
- OSTI ID:
- 808375
- Report Number(s):
- PPPL-3729; TRN: US0301994
- Resource Relation:
- Other Information: PBD: 19 Jul 2002
- Country of Publication:
- United States
- Language:
- English
Similar Records
Nonlinear Gyrokinetics: A Powerful Tool for the Description of Microturbulence in Magnetized Plasmas
Accurate and precise determination of critical properties from Gibbs ensemble Monte Carlo simulations