Equilibrium properties, calculation

Two generally accepted models for the vapor phase were discussed in Chapter 3 and one particular model for the liquid phase (UNIQUAC) was discussed in Chapter 4. Unfortunately, these, and all other presently available models, are only approximate when used to calculate equilibrium properties of dense fluid mixtures. Therefore, any such model must contain a number of adjustable parameters, which can only be obtained from experimental measurements. The predictions of the model may be sensitive to the values selected for model parameters, and the data available may contain significant measurement errors. Thus, it is of major importance that serious consideration be given to the proper treatment of experimental measurements for mixtures to obtain the most appropriate values for parameters in models such as UNIQUAC. 

The probability fimction P r ) cannot be factored into contributions from individual particles, since they are coupled by their interactions. However, integration over the coordinates of all but a few particles leads to reduced probability fiinctions containing the infomration necessary to calculate the equilibrium properties of the system. 

The most conunon choice for a reference system is one with hard cores (e.g. hard spheres or hard spheroidal particles) whose equilibrium properties are necessarily independent of temperature. Although exact results are lacking in tluee dimensions, excellent approximations for the free energy and pair correlation fiinctions of hard spheres are now available to make the calculations feasible. 

Critical Temperature The critical temperature of a compound is the temperature above which a hquid phase cannot be formed, no matter what the pressure on the system. The critical temperature is important in determining the phase boundaries of any compound and is a required input parameter for most phase equilibrium thermal property or volumetric property calculations using analytic equations of state or the theorem of corresponding states. Critical temperatures are predicted by various empirical methods according to the type of compound or mixture being considered. 

In a recent paper [11] this approach has been generalized to deal with reactions at surfaces, notably dissociation of molecules. A lattice gas model is employed for homonuclear molecules with both atoms and molecules present on the surface, also accounting for lateral interactions between all species. In a series of model calculations equilibrium properties, such as heats of adsorption, are discussed, and the role of dissociation disequilibrium on the time evolution of an adsorbate during temperature-programmed desorption is examined. This approach is adaptable to more complicated systems, provided the individual species remain in local equilibrium, allowing of course for dissociation and reaction disequilibria. 

Eqs. (1,4,5) show that to determine the equilibrium properties of an adsorbate and also the adsorption-desorption and dissociation kinetics under quasi-equilibrium conditions we need to calculate the chemical potential as a function of coverage and temperature. We illustrate this by considering a single-component adsorbate. The case of dissociative equilibrium with both atoms and molecules present on the surface has recently been given elsewhere [11]. 

In order to calculate an equilibrium property A of the system, we have to evaluate the integral... 

In statistical mechanics the properties of a system in equilibrium are calculated from the partition function, which depending on the choice for the ensemble considered involves a sum over different states of the system. In the very popular canonical ensemble, that implies a constant number of particles N, volume V, and temperature T conditions, the quasiclassical partition function Q is... 

It should be kept in mind that an objective function which does not require any phase equilibrium calculations during each minimization step is the basis for a robust and efficient estimation method. The development of implicit objective functions is based on the phase equilibrium criteria (Englezos et al. 1990a). Finally, it should be noted that one important underlying assumption in applying ML estimation is that the model is capable of representing the data without any systematic deviation. Cubic equations of state compute equilibrium properties of fluid mixtures with a variable degree of success and hence the ML method should be used with caution. 

Equilibrium properties are surprisingly accurately predicted by molecular-level SCF calculations. MC simulations help us to understand why the SCF theory works so well for these densely packed layers. In effect, the high density screens the correlations for chain packing and chain conformation effects to such a large extent that the properties of a single chain in an external field are rather accurate. Cooperative fluctuations, such as undulations, are not included in the SCF approach. Also, undulations cannot easily develop in an MD box. To see undulations, one needs to perform molecularly realistic simulations on very large membrane systems, which are extremely expensive in terms of computation time. 

The theoretical prediction of these properties for branched molecules has to take into account the peculiar aspects of these chains. It is possible to obtain these properties as the low gradient Hmits of non-equilibrium averages, calculated from dynamic models. The basic approach to the dynamics of flexible chains is given by the Rouse or the Rouse-Zimm theories [12,13,15,21]. How-ever,both the friction coefficient and the intrinsic viscosity can also be evaluated from equilibrium averages that involve the forces acting on each one of the units. This description is known as the Kirkwood-Riseman (KR) theory [15,71 ]. Thus, the translational friction coefficient, fl, relates the force applied to the center of masses of the molecule and its velocity... 

An equilibrium-flash calculation (using the same equations as in case A above) is made at each point in time to find the vapor and liquid flow rates and properties immediately after the pressure letdown valve (the variables with the primes F , F l, y], x j,.. . shown in Fig. 3.8). These two streams are then fed into the vapor and liquid phases. The equations describing the two phases will be similar to Eqs. (3.40) to (3.42) and (3.44) to (3.46) with the addition of (1) a multi-component vapor-liquid equilibrium equation to calculate Pi and (2) NC — 1 component continuity equations for each phase. Controller equations relating 1 to Fi and P to F complete the model. 

We will consider dipolar interaction in zero field so that the total Hamiltonian is given by the sum of the anisotropy and dipolar energies = E -TEi. By restricting the calculation of thermal equilibrium properties to the case 1. we can use thermodynamical perturbation theory [27,28] to expand the Boltzmann distribution in powers of This leads to an expression of the form [23]... 

Molecular modeling, like all other technical disciplines, has its own jargon. Much of this is described in Appendix B (Common Terms and Acronyms), and only one aspect will be addressed here. This concerns specification of theoretical model used for property calculation together with theoretical model used for equilibrium (or transition-state) geometry calculation. 

This is a question of considerable practical importance, given that optimization of equilibrium geometry can easily require one or two orders of magnitude more computation than an energy (property) calculation at a single geometry (see Chapter 11). Rephrased, the question might read ... 

In terms of layout, it might be preferable from a historic sense to start with quantum theories and then develop classical theories as an approximation to the more rigorous formulation. However, I think it is more pedagogically straightforward (and far easier on the student) to begin with classical models, which are in the widest use by experimentalists and tend to feel very intuitive to the modern chemist, and move from there to increasingly more complex theories. In that same vein, early emphasis will be on single-molecule (gas-phase) calculations followed by a discussion of extensions to include condensed-phase effects. While the book focuses primarily on the calculation of equilibrium properties, excited states and reaction dynamics arc dealt with as advanced subjects in later chapters. 

In recent years we have undertaken a systematic investigation of the volumes and heat capacities of transfer of alkali halides and tetraalkylammonium bromides from water to mixed aqueous solvents (1-6). These properties are important because, when combined with enthalpies and free energies, they can be used to calculate the temperature and pressure dependences of various equilibrium properties of electrolytes in mixed solvents. Since the properties of electrolytes in mixed aqueous solvents are closely related to the corresponding properties of the nonelectrolyte in an electrolyte solution, infor-... 

These calculated intracrystalline diffusion coefficients are particularly appropriate for comparison with those determined from pulsed field gradient (PFG) NMR experiments. Time-independent equilibrium properties such as adsorbate conformations are also readily accessible. The classical nature of the simulations allows a particle s trajectory to be followed, and from this it is possible to determine all kinds of information, such as how often a particle diffuses through a certain region. 

The present state in the theory of time-dependent processes in liquids is the following. We know which correlation functions determine the results of certain physical measurements. We also know certain general properties of these correlation functions. However, because of the mathematical complexities of the V-body problem, the direct calculation of the fulltime dependence of these functions is, in general, an extremely difficult affair. This is analogous to the theory of equilibrium properties of liquids. That is, in equilibrium statistical mechanics the equilibrium properties of a system can be found if certain multidimensional integrals involving the system s partition function are evaluated. However, the exact evaluation of these integrals is usually extremely difficult especially for liquids. 

Molecular Dynamics and Monte Carlo Simulations. At the heart of the method of molecular dynamics is a simulation model consisting of potential energy functions, or force fields. Molecular dynamics calculations represent a deterministic method, ie, one based on the assumption that atoms move according to laws of Newtonian mechanics. Molecular dynamics simulations can be performed for short time-periods, eg, 50—100 picoseconds, to examine localized very high frequency motions, such as bond length distortions, or, over much longer periods of time, eg, 500—2000 ps, in order to derive equilibrium properties. It is worthwhile to summarize what properties researchers can expect to evaluate by performing molecular simulations ... 

The (liquid 4- liquid) equilibria diagram for (cyclohexane + methanol) was taken from D. C. Jones and S. Amstell, The Critical Solution Temperature of the System Methyl Alcohol-Cyclohexane as a Means of Detecting and Estimating Water in Methyl Alcohol , J. Chem. Soc., 1930, 1316-1323 (1930). The G results were calculated from the (vapor 4- liquid) results of K. Strubl, V. Svoboda, R. Holub, and J. Pick, Liquid-Vapour Equilibrium. XIV. Isothermal Equilibrium and Calculation of Excess Functions in the Systems Methanol -Cyclohexane and Cyclohexane-Propanol , Collect. Czech. Chem. Commun., 35, 3004-3019 (1970). The results are from M. Dai and J.-P.Chao, Studies on Thermodynamic Properties of Binary Systems Containing Alcohols. II. Excess Enthalpies of C to C5 Normal Alcohols + 1,4-Dioxane , Fluid Phase Equilib., 23, 321-326 (1985). 

---