Matrix equations, thermodynamics

It is important to recognize that the small subset of matrix equations introduced in the main text (typically, restricted to real matrix elements) will be found sufficient to exploit the geometrical simplicity that underlies equilibrium thermodynamics. Nevertheless, it is useful to introduce the thermodynamic vector geometry in the broader framework of matrix theory and Dirac notation that is broadly applicable to the advanced thermodynamic topics of Chapters 11-13, as well as to many other areas of modem physical chemistry research. 

Physicist P. A. M. Dirac suggested an inspired notation for the Hilbert space of quantum mechanics [essentially, the Euclidean space of (9.20a, b) for / — oo, which introduces some subtleties not required for the finite-dimensional thermodynamic geometry]. Dirac s notation applies equally well to matrix equations [such as (9.7)-(9.19)] and to differential equations [such as Schrodinger s equation] that relate operators (mathematical objects that change functions or vectors of the space) and wavefunctions in quantum theory. Dirac s notation shows explicitly that the disparate-looking matrix mechanical vs. wave mechanical representations of quantum theory are actually equivalent, by exhibiting them in unified symbols that are free of the extraneous details of a particular mathematical representation. Dirac s notation can also help us to recognize such commonality in alternative mathematical representations of equilibrium thermodynamics. 

The extension of ideal phase analysis of the Maxwell-Stefan equations to nonideal liquid mixtures requires the sufficiently accurate estimation of composition-dependent mutual diffusion coefficients and the matrix of thermodynamic factors. However, experimental data on mutual diffusion coefficients are rare, and prediction methods are satisfactory only for certain types of liquid mixtures. The thermodynamic factor may be calculated from activity coefficient models such as NRTL or UNIQUAC, which have adjustable parameters estimated from experimental phase equilibrium data. The group contribution method of UNIFAC may also be helpful, as it has a readily available parameter table consisting of mam7 species. If, however, reliable data are not available, then the averaged values of the generalized Maxwell-Stefan diffusion coefficients and the matrix of thermodynamic factors are calculated at some mean composition between x0i and xzi. Hence, the matrix of zero flux mass transfer coefficients [k ] is estimated by... 

The subscripts 0 and 5 serve as reminders that the appropriate compositions, x q and respectively, have to be used in the defining equations for (Eqs. 2.1.20 and 2.1.21). The subscript av on [F] emphasizes the assumption that the matrix of thermodynamic factors evaluated at the average mole fraction and is assumed constant along the diffusion path. 

In Sec. 3 our presentation is focused on the most important results obtained by different authors in the framework of the rephca Ornstein-Zernike (ROZ) integral equations and by simulations of simple fluids in microporous matrices. For illustrative purposes, we discuss some original results obtained recently in our laboratory. Those allow us to show the application of the ROZ equations to the structure and thermodynamics of fluids adsorbed in disordered porous media. In particular, we present a solution of the ROZ equations for a hard sphere mixture that is highly asymmetric by size, adsorbed in a matrix of hard spheres. This example is relevant in describing the structure of colloidal dispersions in a disordered microporous medium. On the other hand, we present some of the results for the adsorption of a hard sphere fluid in a disordered medium of spherical permeable membranes. The theory developed for the description of this model agrees well with computer simulation data. Finally, in this section we demonstrate the applications of the ROZ theory and present simulation data for adsorption of a hard sphere fluid in a matrix of short chain molecules. This example serves to show the relevance of the theory of Wertheim to chemical association for a set of problems focused on adsorption of fluids and mixtures in disordered microporous matrices prepared by polymerization of species. 

Here Jta(x) denotes the a-th component of the stationary vector x of the Markov chain with transition matrix Q whose elements depend on the monomer mixture composition in microreactor x according to formula (8). To have the set of Eq. (24) closed it is necessary to determine the dependence of x on X in the thermodynamic equilibrium, i.e. to solve the problem of equilibrium partitioning of monomers between microreactors and their environment. This thermodynamic problem has been solved within the framework of the mean-field Flory approximation [48] for copolymerization of any number of monomers and solvents. The dependencies xa=Fa(X)(a=l,...,m) found there in combination with Eqs. (24) constitute a closed set of dynamic equations whose solution permits the determination of the evolution of the composition of macroradical X(Z) with the growth of its length Z, as well as the corresponding change in the monomer mixture composition in the microreactor. 

Perhaps the best starting point in a review of the nonequilibrium field, and certainly the work that most directly influenced the present theory, is Onsager s celebrated 1931 paper on the reciprocal relations [10]. This showed that the symmetry of the linear hydrodynamic transport matrix was a consequence of the time reversibility of Hamilton s equations of motion. This is an early example of the overlap between macroscopic thermodynamics and microscopic statistical mechanics. The consequences of time reversibility play an essential role in the present nonequilibrium theory, and in various fluctuation and work theorems to be discussed shortly. 

A rather general method of the calculation of the tunneling taking account of the dissipation was given in Ref. 82. The cases of rather strong dissipation were considered in Refs. 81 and 82, where it was assumed that a thermodynamical equilibrium in the initial potential well exists. The case of extremely weak friction has been considered using the equations for the density matrix in Ref. 83. A quantum analogue of the Focker-Planck equation for the adiabatic and nonadiabatic processes in condensed media was obtained in Refs. 105 and 106. 

When applying the matrix method to very large values of m, it is only the largest root of the characteristic equation that is significant for calculating the thermodynamic properties of the system. For the particular matrix defined in Eq. (7.1.9), the characteristic equation is... 

Kirkaldy J.S., Weichert D., and Haq Z.U. (1963) Diffusion in multicomponent metallic systems, VI some thermodynamic properties of the D matrix and the corresponding solutions of the diffusion equations. Can. f. Phys. 41, 2166-2173. 

It will be useful practice for the physical chemistry student to rewrite other matrix and vector equations of Sidebar 9.1 in Dirac notation, both for future applications to quantum theory as well as the intended present application to equilibrium thermodynamics. 

Taken together, (11.22) and (11.23) lead to various thermodynamic identities between measured response functions, as will be illustrated below. Equation (11.23) shows that the inverse metric matrix M-1 plays a role for conjugate vectors R/) that is highly analogous to the role played by M itself for the intensive vectors R,). In view of this far-reaching relationship, we can define the conjugate metric M,... 

Exercise. In particular this situation obtains if the equation refers to a closed thermodynamic system. The corresponding matrix S oo) = is known from equilib-... 

Solve the Liouville equation for the density matrix p(0 at time t, given that the system is initially in thermodynamic equilibrium. 

The frequency matrix Qy and the memory function matrix Ty, in the relaxation equation are equivalent to the Liouville operator matrix Ly and the Uy matrix, respectively. The later two matrices were introduced by Kadanoff and Swift [37] (see Section V). Thus the frequency matrix can be identified with the static variables (the wavenumber-dependent thermodynamic quantities) associated with the nondissipative part, and the memory kernel matrix can be identified with the transport coefficients associated with the dissipative part. 

The composition boundary values entering into Eqs. (All) represent external values for Eqs. (A10). With some further assumptions concerning the diffusion and reaction terms, this allows an analytical solution of the boundary-value problem [Eqs. (A10) and (All)] in a closed matrix form (see Refs. 58 and 135). On the other hand, the boundary values need to be determined from the total system of equations describing the process. The bulk values in both phases are found from the balance relations, Eqs. (Al) and (A2). The interfacial liquid-phase concentrations xj are related to the relevant concentrations of the second fluid phase, y , by the thermodynamic equilibrium relationships and by the continuity condition for the molar fluxes at the interface (57,135). 

---