Temperature and pressure estimators next up previous
Next: Energy fluctuations in the Up: No Title Previous: No Title

Temperature and pressure estimators

From the classical virial theorem

displaymath29

we arrived at the equipartition theorem:

displaymath33

where tex2html_wrap_inline420 are the N Cartesian momenta of the N particles in a system. This says that the microscopic function of the N momenta that corresponds to the temperature, a macroscopic observable of the system, is given by

displaymath39

The ensemble average of K can be related directly to the temperature

displaymath45

tex2html_wrap_inline430 is known as an estimator (a term taken over from the Monte Carlo literature) for the temperature. An estimator is some function of the phase space coordinates, i.e., a function of microscopic states, whose ensemble average gives rise to a physical observable.

An estimator for the pressure can be derived as well, starting from the basic thermodynamic relation:

displaymath53

with

displaymath59

The volume dependence of the partition function is contained in the limits of integration, since the range of integration for the coordinates is determined by the size of the physical container. For example, if the system is confined within a cubic box of volume tex2html_wrap_inline432 , with L the length of a side, then the range of each q integration will be from 0 to L. If a change of variables is made to tex2html_wrap_inline440 , then the range of each s integration will be from 0 to 1. The coordinates tex2html_wrap_inline444 are known as scaled coordinates. For containers of a more general shape, a more general transformation is

displaymath66

In order to preserve the phase space volume element, however, we need to ensure that the transformation is a canonical one. Thus, the corresponding momentum transformation is

displaymath69

With this coordinate/momentum transformation, the phase space volume element transforms as

displaymath72

Thus, the volume element remains the same as required. With this transformation, the Hamiltonian becomes

displaymath78

and the canonical partition function becomes

displaymath89

Thus, the pressure can now be calculated by explicit differentiation with respect to the volume, V:

eqnarray99

Thus, the pressure estimator is

displaymath134

and the pressure is given by

displaymath141

For periodic systems, such as solids and currently used models of liquids, an absolute Cartesian coordinate tex2html_wrap_inline448 is ill-defined. Thus, the virial part of the pressure estimator tex2html_wrap_inline450 must be rewritten in a form appropriate for periodic systems. This can be done by recognizing that the force tex2html_wrap_inline452 is obtained as a sum of contributions tex2html_wrap_inline454 , which is the force on particle i due to particle j. Then, the classical virial becomes

eqnarray144

where tex2html_wrap_inline460 is now a relative coordinate. tex2html_wrap_inline460 must be computed consistent with periodic boundary conditions, i.e., the relative coordinate is defined with respect to the closest periodic image of particle j with respect to particle i. This gives rise to surface contributions, which lead to a nonzero pressure, as expected.


next up previous
Next: Energy fluctuations in the Up: No Title Previous: No Title

Mark Tuckerman
Wed Feb 7 17:32:09 EST 2001