Probability density function theories

One of the simplest chemical reactions involving a barrier, H2 + H —> [H—H—H] —> II + H2, has been investigated in some detail in a number of publications. The theoretical description of this hydrogen abstraction sequence turns out to be quite involved for post-Hartree-Fock methods and is anything but a trivial task for density functional theory approaches. Table 13-7 shows results reported by Johnson et al., 1994, and Csonka and Johnson, 1998, for computed classical barrier heights (without consideration of zero-point vibrational corrections or tunneling effects) obtained with various methods. The CCSD(T) result of 9.9 kcal/mol is probably very accurate and serves as a reference (the experimental barrier, which of course includes zero-point energy contributions, amounts to 9.7 kcal/mol). 

Since the early days of quantum mechanics, the wave function theory has proven to be very successful in describing many different quantum processes and phenomena. However, in many problems of quantum chemistry and solid-state physics, where the dimensionality of the systems studied is relatively high, ab initio calculations of the structure of atoms, molecules, clusters, and crystals, and their interactions are very often prohibitive. Hence, alternative formulations based on the direct use of the probability density, gathered under what is generally known as the density matrix theory [1], were also developed since the very beginning of the new mechanics. The independent electron approximation or Thomas-Fermi model, and the Hartree and Hartree-Fock approaches are former statistical models developed in that direction [2]. These models can be considered direct predecessors of the more recent density functional theory (DFT) [3], whose principles were established by Hohenberg,... 

A theoretical framework based on the one-point, one-time joint probability density function (PDF) is developed. It is shown that all commonly employed models for turbulent reacting flows can be formulated in terms of the joint PDF of the chemical species and enthalpy. Models based on direct closures for the chemical source term as well as transported PDF methods, are covered in detail. An introduction to the theory of turbulence and turbulent scalar transport is provided for completeness. 

In order to compare various reacting-flow models, it is necessary to present them all in the same conceptual framework. In this book, a statistical approach based on the one-point, one-time joint probability density function (PDF) has been chosen as the common theoretical framework. A similar approach can be taken to describe turbulent flows (Pope 2000). This choice was made due to the fact that nearly all CFD models currently in use for turbulent reacting flows can be expressed in terms of quantities derived from a joint PDF (e.g., low-order moments, conditional moments, conditional PDF, etc.). Ample introductory material on PDF methods is provided for readers unfamiliar with the subject area. Additional discussion on the application of PDF methods in turbulence can be found in Pope (2000). Some previous exposure to engineering statistics or elementary probability theory should suffice for understanding most of the material presented in this book. 

The importance of N-representability for pair-density functional theory was not fully appreciated probably because most research on pair-density theories has been performed by people from the density functional theory community, and there is no W-representability problem in conventional density functional theory. Perhaps this also explains why most work on the pair density has been performed in the first-quantized spatial representation (p2(xi,X2) = r2(xi,X2 xi,X2)) instead of the second-quantized orbital representation... 

The probability density function W(0) at r = 0 and the directional correlation factors for homologous PACA sequences with x = 2 to 7 units are evaluated. The influences of these factors on the cyclization equilibria constants Kx are determined. Agreement of theory with experimental results for x = 3 - 6 is within limits set by uncertainties in the calculations combined with experimental errors, i.e., within about 15 %. 

Here, we present an approach for the description of such anomalous transport processes that is based on the continuous-time random walk theory for a power-law waiting time distribution w(t) but which can be used to find the probability density function of the random walker in the presence of an external force field, or in phase space. This framework is fractional dynamics, and we show how the traditional kinetic equations can be generalized and solved within this approach. 

---