A statistical physics of stationary and metastable states

We present a generalization of Gibbs statistical mechanics designed to describe a general class of stationary and metastable equilibrium states. It is assumed that the physical system maximizes the entropy functional S subject to the standard conditions plus an extra conserved constraint function F, imposed to force the system to remain in the metastable configuration. After requiring additivity for two quasi-independent subsystems, and the commutation of the new constraint with the density matrix ρ, it is argued that F should be a homogeneous function of ρ, at least for systems in which the spectrum is sufficiently dense to be considered as continuous. Therefore, surprisingly, the analytic form of F turns out to be of the kind F(p_i) = p_i^q, where the p_i are the eigenvalues of the density matrix and q is a real number to be determined. Thus, the discussion identifies the physical relevance of Lagrange multiplier constraints of the Tsallis kind and their q parameter, as enforced by the additivity of the constraint F which fixes the metastable state. An approximate analytic solution for the probability density is found for q close to unity. The procedure is applied to describe the results from the plasma experiment of Huang and Driscoll. For small and medium values of the radial distance, the measured density is predicted with a precision similar to that achieved by minimal enstrophy and Tsallis procedures. Also, the particle density is predicted at all the radial positions. Thus, the discussion gives a solution to the conceptual difficulties of the two above mentioned approaches as applied to this problem, which both predict a non-analytic abrupt vanishing of the density above a critical radial distance.