Molecular system thermodynamical quantities

Conformational free energy simulations are being widely used in modeling of complex molecular systems [1]. Recent examples of applications include study of torsions in n-butane [2] and peptide sidechains [3, 4], as well as aggregation of methane [5] and a helix bundle protein in water [6]. Calculating free energy differences between molecular states is valuable because they are observable thermodynamic quantities, related to equilibrium constants and... 

Ire boundary element method of Kashin is similar in spirit to the polarisable continuum model, lut the surface of the cavity is taken to be the molecular surface of the solute [Kashin and lamboodiri 1987 Kashin 1990]. This cavity surface is divided into small boimdary elements, he solute is modelled as a set of atoms with point polarisabilities. The electric field induces 1 dipole proportional to its polarisability. The electric field at an atom has contributions from lipoles on other atoms in the molecule, from polarisation charges on the boundary, and where appropriate) from the charges of electrolytes in the solution. The charge density is issumed to be constant within each boundary element but is not reduced to a single )oint as in the PCM model. A set of linear equations can be set up to describe the electrostatic nteractions within the system. The solutions to these equations give the boundary element harge distribution and the induced dipoles, from which thermodynamic quantities can be letermined. 

From the standpoint of thermodynamics, the essential quantity which governs the course of a chemical reaction is the chemical potential of the system or, in the case of interest, the free energy storage within the molecular coil. This quantity, unfortunately, is difficult to evaluate in non-steady flow. At modest extension ratios (X < 4), the free energy storage of a freely-jointed bead-spring chain is... 

The net free energy is, in the end, the thermodynamic quantity that dictates molecular behavior. However, to understand why the free energy profile for a system looks as it does, it is valuable to also determine the potential and entropic components of the net free energy ... 

The molecular modelling of systems consisting of many molecules is the field of statistical mechanics, sometimes called statistical thermodynamics [28,29], Basically, the idea is to go from a molecular model to partition functions, and then, from these, to predict thermodynamic observables and dynamic and structural quantities. As in classical thermodynamics, in statistical mechanics it is essential to define which state variables are fixed and which quantities are allowed to fluctuate, i.e. it is essential to specify the macroscopic boundary conditions. In the present context, there are a few types of molecular systems of interest, which are linked to so-called ensembles. 

Equation (4.70) is a starting point in the determination of diffusivities in liquid metal alloys, but in most real systems, experimental values are difficult to obtain to confirm theoretical expressions, and pair potentials and molecular interactions that exist in liquid alloys are not sufficiently quantified. Even semiempirical approaches do not fare well when applied to liquid alloy systems. There have been some attempts to correlate diffusivities with thermodynamic quantities such as partial molar enthalpy and free energy of solution, but their application has been limited to only a few systems. 

In discussing quantitatively the molecular recognition phenomena in supramolecu-lar chemistry, it is indispensable to determine the thermodynamic parameters for each (supra)molecular interaction of interest. In principle, these parameters should be determined by means of calorimetry. However, in spite of its long history as an established methodology,29,30 the calorimetry does not appear to be the first choice to determine the thermodynamic quantities for various supramolecular systems. This is probably because a relatively large sample amount, sophisticated equipment, and some experience are required in the precise calorimetric measurement.29,30... 

There have been several approaches to the expression of thermodynamic quantities of solutions. Scatchard published a series of papers (see 1801 and previous papers, especially 1802) based on the classical approach via chemical potentizil. Barker (130, 128, 131) applied the theory of conformal solutions (1254) to some H bonding systems after modifying it to allow for dipole attractions or, more generally, molecular orientations. The curves are similar in both cases. 

Perhaps one of the greatest successes of the molecular dynamics (MD) method is its ability both to predict macroscopically observable properties of systems, such as thermodynamic quantities, structural properties, and time correlation functions, and to allow modeling of the microscopic motions of individual atoms. From modeling, one can infer detailed mechanisms of structural transformations, diffusion processes, and even chemical reactions (using, for example, the method of ab initio molecular dynamics).Such information is extremely difficult, if not impossible, to obtain experimentally, especially when detailed behavior of a local defect is sought. The variety of different experimental conditions that can be mimicked in an MD simulation, such as... 

All the methods we have thus far considered use Monte Carlo or molecular dynamics simulations. An alternative approach to get thermodynamic quantities in an effective Hamiltonian+discrete model is based on the use of the integral equation RISM method. The RISM-SCF method proposed by Ten-no and coworkers (Ten-no et al., 1993, 1994 Kawata et al., 1995) combines a QM description of the solute with a RISM description of the whole system in a way which deserves attention. 

Because changes in the thermodynamic quantity entropy may be understood in terms of changes in molecular disorder, we can often predict the sign of The following illustrations emphasize several common types of processes that result in predictable entropy changes for the system. 

According to the Lorentz-Lorenz equation (4.3.21) for the molar refraction at optical frequencies, Y is directly proportional to the molecular polarizability p. The Koppel-Palm equation has also been applied to the analysis of solvent effects on thermodynamic quantities related to the solvation of electrolytes [48, 49]. In the case of the systems considered in table 4.11, addition of the parameter X to the linear equation describing the solvent effect improves the quality of the fit to the experimental data, especially in the case of alkali metal halide electrolytes involving larger ions. The parameter Y is not important for these systems but does assist in the interpretation of other thermodynamic quantities which are solvent dependent [48, 49]. Addition of these parameters to the analysis is only possible when the solvent-dependent phenomenon has been studied in a large number of solvents. 

With Eq. (1.88), we conclude our discussion of phenomenological thermodynamics of confined fluids. In Chapter 2, we shall turn to an interpretation of the various thermodynamic quantities introduced above in terms of interactions between the microscopic constituents forming the system at a molecular level of description (i.e., atoms and molecules). 

It is well appreciated that thermodynamic and kinetic parameters are difficult to compute for organometallic molecular systems (see, e.g.. Refs 12-14 and Chap. 4 by Frenking in the present volume). In particular, such quantities cannot be predicted within an independent-particle, single-determinant Hartree-Fock type of approach electron correlation must be included in the computational methods applied to achieve reliable and accurate results. In this work, we examine the performance of three first-principles methods, generally acknowledged by the abbreviations BLYP, B3LYP, and MP2. The first two are methods based on density functional theory (DFT) (15) the latter is an ab initio, molecular orbital (MO)-based method (16). 

---