Steady State Fluctuations

9.1 Introduction

Nonequilibrium steady states are fascinating systems to study. Although there are many parallels between these states and equilibrium states, a convincing theoretical description of steady states, particularly far from equilibrium, has yet to be found. Close to equilibrium, linear response theory and linear irreversible thermodynamics provide a relatively complete treatment, (§2.1 - 2.3). However, in systems where local thermodynamic equilibrium has broken down, and thermodynamic properties are not the same local functions of thermodynamic state variables that they are at equilibrium, our understanding is very primitive indeed.

In §7.3 we gave a statistical mechanical description of thermostatted, nonequilibrium steady states far from equilibrium - the transient time correlation function (TTCF) and Kawasaki formalisms. The Transient Time Correlation Function is the nonlinear analogue of the Green-Kubo correlation functions. For linear transport processes the Green-Kubo relations play a role which is analogous to that of the partition function at equilibrium. Like the partition function, Green-Kubo relations are highly nontrivial to evaluate. They do however provide an exact starting point from which one can derive exact interrelations between thermodynamic quantities. The Green-Kubo relations also provide a basis for approximate theoretical treatments as well as being used directly in equilibrium molecular dynamics simulations.

The TTCF and Kawasaki expressions may be used as nonlinear, nonequilibrium partition functions. For example if a particular derivative commutes with the thermostatted, field-dependent propagator then one can formally differentiate the TTCF and Kawasaki expressions for steady state phase averages, yielding fluctuation expressions for the so-called derived properties. The key point in such derivations is that the particular derivative should commute with the relevant propagators. If this is not so one cannot derive tractable or useful results.

In order to constrain thermodynamic variables two basic feedback mechanisms can be employed: the integral feedback mechanism employed for example in the Nose-Hoover thermostat, (§5.2) and the differential mechanism employed in the Gaussian thermostat. A third mechanism, the proportional mechanism has not found much use either in simulations or in theory because it necessarily employs irreversible equations of motion.

In this chapter we will derive fluctuation expressions for the derivatives of steady state phase averages. We will derive expressions for derivatives with respect to temperature, pressure and the mean value of the dissipative flux. Applying these derivatives in turn to averages of the internal energy, the volume and the thermodynamic driving force yields expressions for the specific heats, the compressibility and the inverse Burnett coefficients respectively. In order to ensure the commutivity of the respective derivatives and propagators, we will employ the Gaussian feedback mechanism exclusively. Corresponding derivations using Nose-Hoover feedback are presently unknown.

Rather than giving a general but necessarily formal derivation of the fluctuation formulae, we will instead concentrate on two specific systems: planar Couette flow and colour conductivity. By concentrating on specific systems we hope to make the discussion more concrete and simultaneously illustrate particular applications of the theory of nonequilibrium steady states discussed in Chapter 7.