Theoretical methods thermodynamic property calculations

The common theme in the evolution of methods for property and parameter prediction is the development of equations, either theoretical or empirical, containing quantities that can be calculated from theoretical considerations or experimental data. Mathematical expressions for correlating thermodynamic data may take several forms. 

In the last decade, quantum-chemical investigations have become an integral part of modern chemical research. The appearance of chemistry as a purely experimental discipline has been changed by the development of electronic structure methods that are now widely used. This change became possible because contemporary quantum-chemical programs provide reliable data and important information about structures and reactivities of molecules and solids that complement results of experimental studies. Theoretical methods are now available for compounds of all elements of the periodic table, including heavy metals, as reliable procedures for the calculation of relativistic effects and efficient treatments of many-electron systems have been developed [1, 2] For transition metal (TM) compounds, accurate calculations of thermodynamic properties are of particularly great usefulness due to the sparsity of experimental data. 

Traube s rule accommodates the balance between hydrophobicity and hydro-philicity. It has been extended somewhat and formalized with the development of quantitative methods to estimate the surface area of molecules based on their structures [19, 237]. The molecular surface area approach suggests that the number of water molecules that can be packed around the solute molecule plays an important role in the theoretical calculation of the thermodynamic properties of the solution. Hence, the molecular surface area of the solute is an important parameter in the theory. In compounds other than simple normal alkanes, the functional groups will tend to be more or less polar and thus relatively compatible with the polar water matrix [227,240]. Hence, the total surface area of the molecule can be subdivided into functional group surface area and hydro carbonaceous surface area . These quantities maybe determined for simple compounds as an additive function of constituent groups with subtractions made for the areas where intramolecular contact is made and thus no external surface is presented. 

All these methods have found applications in theoretical considerations of numerous problems more or less directly related to solvent extraction. The MM calculated structures and strain energies of cobalt(III) amino acid complexes have been related to the experimental distribution of isomers, their thermodynamic stability, and some kinetic data connected with transition state energies [15]. The influence of steric strain upon chelate stability, the preference of metal ions for ligands forming five- and six-membered chelate rings, the conformational isomerism of macrocyclic ligands, and the size-match selectivity were analyzed [16] as well as the relation between ligand structures, coordination stereochemistry, and the thermodynamic properties of TM complexes [17]. 

Most of this book concerns the development and application of theoretical thermodynamic models, as these are the basis of the CALPHAD method. However, none of this would be possible without the existence of the computational methods and software which allow these models to be applied in practice. In essence, the issues involved in computational methods are less diverse and mainly revolve around Gibbs energy minimisation. In addition, there are optimiser codes which are used for the thermodynamic assessment of phase equilibria. The essential aim of these codes is to reduce the statistical error between calculated phase equilibria, thermodynamic properties and the equivalent experimentally measured quantities. 

Mary Jo Nye enumerates the topics treated in the Journal as molecular spectroscopy and molecular structures, the quantum mechanical treatment of electronic structure of molecules and crystals and the problem of chemical binding, the kinetics of chemical reactions from the standpoint of basic physical principles, the thermodynamic properties of substances and calculation by statistical mechanical methods, the structure of crystals, and surface phenomena. M.J. Nye, From Chemical Philosophy to Theoretical Chemistry. Dynamics of Matter and Dynamics of Disciplines (Berkeley University of California Press 1993), 254. Many of these were considered by Barriol, as we will see later in this chapter. 

The transport and adsorption properties of hydrocarbons on microporous zeolites have been of practical interest due to the important properties of zeolites as shape-selective adsorbents and catalysts. The system of benzene adsorbed on synthetic faujasite-type zeolites has been thoroughly studied because benzene is an ideal probe molecule and the related role of aromatics in zeolitic catalysts for alkylation and cracking reactions. For instance, its mobility and thermodynamic properties have been studied by conventional diffusion 1-6) and adsorption 7-9) techniques. Moreover, the adsorbate-zeolite interactions and related motion and location of the adsorbate molecules within the zeolite cavities have been investigated by theoretical calculations 10-15) and by various spectroscopic methods such as UV (16, 17), IR 17-23), neutron 24-27), Raman 28), and NMR 29-39). 

Theoretical calculations share with gas-phase kinetic and thermodynamic measurements the common aim of the understanding of the intrinsic reactivity properties of heteroaromatic compounds. The purpose of this subsection is to consider the predictive value of theoretical methods insofar as ionic substitution reactions on simple heteroaromatics are concerned. The topic under discussion is inherently limited by the wide range of interest in the understanding of the principles of these processes in solution. It is exactly in this field that an appropriate amount of data concerning gas-phase structural and reactivity properties of heteroaromatic compounds is at present available from modern experimental techniques that can be tested against theoretical predictions. 

Excess chemical potentials As early as 2002 Lynden-Bell et al. published an investigation about the chemical potential of water and organic solutes in [C,mim] [Cl] [9], The authors stated ... the chemical potential is the most important thermodynamic property of a solute in solution because it determines the solubility and chemical reactivity of a solute. Within this seminal article the authors determined the excess chemical potential (pf) by means of theoretical methods. The excess chemical potentials pA of a series of molecules dissolved in the IL [C,mim][Cl] were calculated by a sequence of transformations [9], It is defined by... 

The experimental and theoretical methods of measuring or calculating the kinetic and thermodynamic properties of reactions of thermal electrons are important to physical chemistry. The ECD method will be compared to other methods. The CURES-EC procedure can be applied to other energetic quantities. We present a method to consolidate diverse data into pseudo-two-dimensional potential energy curves for these reactions [9, 10]. 

The efficiency may depend on the two components selected as the keys because, although the rigorous calculations can closely match measured key component compositions, other components may be off due to inaccuracies in the theoretical column model. Also, since the results of the column calculations depend on the thermodynamic properties prediction methods, the efficiency may also depend on the selection of these methods. 

In the present paper we review recent advances in the symmetry-adapted perturbation theory calculations of interaction potentials and interaction-induced properties. We will give a brief description of the theoretical methods needed on the route from the intermolecular potential and properties to rovibrational spectra and collision-induced Raman spectra. We also discuss applications of the interaction potentials and interaction-induced polarizabilities to compute (thermodynamic and dielectric) second virial coefficients. Finally, we illustrate these theoretical approaches on several examples from our own work. 

Perturbation theory has been applied to anharmonic calculations of spectroscopy from ab initio potentials in a large number of studies [19-25,115-121]. In nearly all cases so far, second-order perturbation theory was employed. The representation of the anharmonic potential generally used in these studies is a polynomial in the normal modes, most often a quartic force field. A code implementing this vibrational method was recently incorporated by V. Barone in gaussian [24]. Calculations were carried out for relatively large molecules, such as pyrrole and furan [25], uracil and thiouracil [118], and azabenzenes [119]. We note that in addition to spectroscopy, the ab initio perturbation theoretic algorithms were also applied to the calculation of thermodynamic properties... 

To date, theoretical studies of the thermodynamic behavior of water clusters have been limited to model potentials [41-45,47,84], although there have been studies that examined the temperature dependence of selected isomers [67,73] and others that have used DFT-based MD simulations to optimize structures and to calculate vibrational spectra [7]. In the absence of either experimental or ab initio data on the thermodynamic properties of small water clusters, it is difficult to assess the reliability of the simulations carried out with various model potentials. Although simulations of the thermodynamic properties of small water clusters are feasible with DFT methods, it has been found that DFT calculations with commonly employed functionals such as Becke3LYP [86,87] incorrectly order various isomers of small water clusters [88]. This appears to be due to the inability of current density functionals to describe long-range dispersion interactions [89,90]. For this reason, it is preferable that finite temperature simulations of water clusters be carried out using an appropriate wavefunction-based electronic stmcture method. 

Other theoretical methods. In addition to the procedure described above, three other theoretical methods have been developed for calculating fractionation factors involving minerals. The first is based on computer simulation of crystal structures and first principles prediction of their thermodynamic properties (Patel et al. 1991 Dove et al. 

---