Statistical mechanics configurational integral

The Metropolis MC [41] was originally developed as a method suited to electronic computers for calculating statistical mechanical configurational integrals. Since the MC sampling is a Markovian process, if we introduce a time scale t that actually labels the order of subsequent configurations X, the dynamic evolution of the probability distribution function P(X, t) is governed by the master equation... 

Given an approximation to V( 0s ) acceptable for the purposes at hand, one can proceed to compute equilibrium, i.e., statistically mechanically averaged, values

for properties P( 0s ) of interest using standard procedures which weight each conformation of the carbohydrate molecule by the Boltzmann factor of V( 0s ) normalized by the configuration integral given in eqn. (9). 

The isotherm model of Schirmer et al. (T.) for sorption in molecular sieves is based on statistical thermodynamics in which the configuration integrals describing the sorption behaviour are extracted from the available data. The model does not presuppose any specific kind of sorption mechanism. The multi-component form of this isotherm derived by Loughlin and Roberts (8 ) is also not limited to any particular sorption mechanism,... 

In many cases in statistical mechanics, we are not interested in the configurational part of the partition function itself, but in averages of the type in Eq. (1.1), where the ratio between integrals is involved. Metropolis et a1. [1] showed that it is possible to devise an efficient Monte Carlo scheme to sample such a ratio even when we do not know the probability density P(q) in configuration space ... 

The topic arises from the following sequence of aspects of entropy when entropy is introduced on a thermodynamic basis the issue is the motion of heat (Jaynes, 1988), and the assessment involves calorimetry an entropy change is evaluated. When entropy is formalized with the classical view of statistical thermodynamics, the entropy is found by evaluating a configurational integral (Bennett, 1976). But a macroscopic physical system at a particular thermodynamic state has a particular entropy, a state function, and the whole description of the physical system shouldn t involve more than a mechanical trajectory for the system in a stationary, equilibrium condition. How are these different concepts compatible ... 

Another way to view MD simulation is as a technique to probe the atomic positions and momenta that are available to a molecular system under certain conditions. In other words, MD is a statistical mechanics method that can be used to obtain a set of configurations distributed according to a certain statistical ensemble. The natural ensemble for MD simulation is the microcanonical ensemble, where the total energy E, volume V, and amount of particles N (NVE) are constant. Modifications of the integration algorithm also allow for the sampling of other ensembles, such as the canonical ensemble (NVT) with constant temperature... 

The path integral technique was first proposed by Feynmann (Feynmann Hibbs, 1965). The purpose of this technique was to deal with questions in quantum mechanics. It has been applied to the study of the statistical mechanics of polymer systems (Kreed, 1972 Doi Edwards, 1986) and liquid crystalline polymers as well (Jahnig, 1981 Warner et al, 1985 Wang Warner, 1986). The path integrals relate the configurations of a polymer chain to the paths of a particle when the particle is undergoing Brownian or diffusive motion. 

The method of Monte Carlo integrations over configuration space seems to be a feasible approach to statistical mechanical problems as yet not analytically soluble. For the computing time of a few hours with presently available electronic computers, it seems possible to obtain pressure for a given volume and temperature to an accuracy of a few percent." ... 

Here we review well-known principles of quantum statistical mechanics as necessary to develop a path-integral representation of the partition function. The equations of quantum statistical mechanics are, like so many equations, easy to write down and difficult to implement (at least, for interesting systems). Our purpose here is not to solve these equations but rather to write them down as integrals over configuration space. These integrals can be seen to have a form that is isomorphic to the discretized path-integral representation of the kernel developed in the previous section. 

This is not to say that the subject is not without significant unanswered theoretical questions. How does one justify applying statistical mechanics to a system of so few particles as the electrons in a molecule How docs one refine the maximum hardness principle to cover the prediction of equilibrium nuclear configurations in a molecule Isolated atoms have integral numbers of electrons, but in a molecule an atom can bear a noninlegral number of electrons how then does one best describe the process of molecule formation from constituent atoms Where is the simple model of chemical bonding itself, which for example explains why the covalent radius of an atom can be determined from the electron density of the isolated atom in the simple way we have proposed And so on. 

To evaluate this configurational integral, we make use of the maximum term method of statistical mechanics, i.e., in the limit N co, V co, but N/V = p = constant, can be approximated by the largest term t in the sum. The problem therefore reduces to one of maximizing / with respect to the occupation numbers. Using Stirling s approximation and introducing the mole fraction jc = N ol)/N of the various components. 

J. E. Mayer used linear graphs as notational abbreviations for the integrands in expansions of the configurational integral and distribution functions for various model systems. The use of graphs in standard combinatorial theory,however, indicated how graphs could be utilized to solve the combinatorial problems involved in treating many-body statistical mechanical... 

U will typically be the sum of potential energies between pairs of particles although, of course, many-body forces or (with minor modification) external forces can be included. Other examples of mechanical quantities would include the molecular distribution functions and the pressure. It would be nice to estimate as well integrals like the configuration integral Q itself since this would give the statistical thermodynamic quantities such as the entropy and free energy. However, this represents a more difficult Monte Carlo (MC) problem, for reasons that will shortly be clear. Some unconventional approaches to this problem are discussed in Chapter 5 of this volume. 

For many-particle systems in classical statistical mechanics, the key numerical problem is the solution of the configurational integral that appears in the averages for a property Q (Eqs. 1, 2). 

The thermodynamic properties for a system of N molecules (or N atoms) can be rigorously accounted for using statistical mechanics. Monte Carlo simulation methods provide the foundation for numerically simulating the configurational integral shown in Eq. (B13) that arise from the statistical mechanics treatment. 

D. A. McQuarrie, Statistical Mechanics, Harper Row, New York, 1976. [This well-known book provides an extensive treatment of the statistical thermodynamics of gases, liquids, and sohds. Chapter 13 provides a comprehensive description of configurational integral equations. Also see the discussion in J. M. Ziman, Models of Disorder The Theoretical Physics of Homogeneously Disordered Systems, Cambridge University Press, Cambridge 1979.]... 

---