Nonequilibrium Equations

VER in liquid O 2 is far too slow to be studied directly by nonequilibrium simulations. The force-correlation function, equation (C3.5.2), was computed from an equilibrium simulation of rigid O2. The VER rate constant given in equation (C3.5.3) is proportional to the Fourier transfonn of the force-correlation function at the Oj frequency. Fiowever, there are two significant practical difficulties. First, the Fourier transfonn, denoted [Pg.3041]

Jarzynski, 1997a] Jarzynski, C. Equilibrium free-energy differences from nonequilibrium measurements A master equation approach. Phys. Rev. E. 56 (1997a) 5018-5035... 

The nonequilibrium-model equations for the stage in Fig. 13-56 are as follows in residual form, where i = component (i = 1 to C), j = stage number (j = 1 to N), and V = a stage in another column that supplies an interlink. 

The jump conditions must be satisfied by a steady compression wave, but cannot be used by themselves to predict the behavior of a specific material under shock loading. For that, another equation is needed to independently relate pressure (more generally, the normal stress) to the density (or strain). This equation is a property of the material itself, and every material has its own unique description. When the material behind the shock wave is a uniform, equilibrium state, the equation that is used is the material s thermodynamic equation of state. A more general expression, which can include time-dependent and nonequilibrium behavior, is called the constitutive equation. 

There are basically two different computer simulation techniques known as molecular dynamics (MD) and Monte Carlo (MC) simulation. In MD molecular trajectories are computed by solving an equation of motion for equilibrium or nonequilibrium situations. Since the MD time scale is a physical one, this method permits investigations of time-dependent phenomena like, for example, transport processes [25,61-63]. In MC, on the other hand, trajectories are generated by a (biased) random walk in configuration space and, therefore, do not per se permit investigations of processes on a physical time scale (with the dynamics of spin lattices as an exception [64]). However, MC has the advantage that it can easily be applied to virtually all statistical-physical ensembles, which is of particular interest in the context of this chapter. On account of limitations of space and because excellent texts exist for the MD method [25,61-63,65], the present discussion will be restricted to the MC technique with particular emphasis on mixed stress-strain ensembles. 

In one dimension the truncation of the equations of motion has been worked out in detail [59]. This has allowed an accurate examination of the role of diffusion in desorption, and implications for the Arrhenius analysis in nonequilibrium situations. The largest deviations from the desorption kinetics of a mobile adsorbate obviously occur for an immobile adsorbate... 

Molecular dynamics, in contrast to MC simulations, is a typical model in which hydrodynamic effects are incorporated in the behavior of polymer solutions and may be properly accounted for. In the so-called nonequilibrium molecular dynamics method [54], Newton s equations of a (classical) many-particle problem are iteratively solved whereby quantities of both macroscopic and microscopic interest are expressed in terms of the configurational quantities such as the space coordinates or velocities of all particles. In addition, shear flow may be imposed by the homogeneous shear flow algorithm of Evans [56]. 

A method is outlined by which it is possible to calculate exactly the behavior of several hundred interacting classical particles. The study of this many-body problem is carried out by an electronic computer which solves numerically the simultaneous equations of motion. The limitations of this numerical scheme are enumerated and the important steps in making the program efficient on the computer are indicated. The applicability of this method to the solution of many problems in both equilibrium and nonequilibrium statistical mechanics is discussed. 

Recently " we propo.sed to describe the nonequilibrium alloy kinetics using the fundamental master equation for probability P of finding the occupation number set... 

The general problem of finding nonequilibrium solutions to Boltzman s equation is, as already mentioned above, an exceedingly difficult problem. Two tools that have proven invaluable in providing insight into nonequilibrium phenomena,... 

One shortcoming of Schild analysis is an overemphasized use of the control dose-response curve (i.e., the accuracy of every DR value depends on the accuracy of the control EC o value). An alternative method utilizes nonlinear regression of the Gaddum equation (with visualization of the data with a Clark plot [10], named for A. J. Clark). This method, unlike Schild analysis, does not emphasize control pECS0, thereby giving a more balanced estimate of antagonist affinity. This method, first described by Lew and Angus [11], is robust and theoretically more sound than Schild analysis. On the other hand, it is not as visual. Schild analysis is rapid and intuitive, and can be used to detect nonequilibrium steady states in the system that can corrupt... 

The third approach is called the thermodynamic theory of passive systems. It is based on the following postulates (1) The introduction of the notion of entropy is avoided for nonequilibrium states and the principle of local state is not assumed, (2) The inequality is replaced by an inequality expressing the fundamental property of passivity. This inequality follows from the second law of thermodynamics and the condition of thermodynamic stability. Further the inequality is known to have sense only for states of equilibrium, (3) The temperature is assumed to exist for non-equilibrium states, (4) As a consequence of the fundamental inequality the class of processes under consideration is limited to processes in which deviations from the equilibrium conditions are small. This enables full linearization of the constitutive equations. An important feature of this approach is the clear physical interpretation of all the quantities introduced. 

Figure 38. Classification of nonequilibrium fluctuations. (Reprinted from M. Asanuma and R. Aogaki, Non-equilibrium fluctuation theory on pitting dissolution. I. Derivation of dissolution current equations." J. Chem. Phys. 106,9938,1997. Copyright 1997, American Institute of Physics.)...

Couette Flow Simulation. MD typically simulate systems at thermodynamic equilibrium. For the simulation of systems undergoing flow various methods of nonequilibrium MD have been developed (Ifl iZ.). In all of these methods the viscosity Is calculated directly from the constitutive equation. 

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. 

Evans and Baranyai [51, 52] have explored what they describe as a nonlinear generalization of Prigogine s principle of minimum entropy production. In their theory the rate of (first) entropy production is equated to the rate of phase space compression. Since phase space is incompressible under Hamilton s equations of motion, which all real systems obey, the compression of phase space that occurs in nonequilibrium molecular dynamics (NEMD) simulations is purely an artifact of the non-Hamiltonian equations of motion that arise in implementing the Evans-Hoover thermostat [53, 54]. (See Section VIIIC for a critical discussion of the NEMD method.) While the NEMD method is a valid simulation approach in the linear regime, the phase space compression induced by the thermostat awaits physical interpretation even if it does turn out to be related to the rate of first entropy production, then the hurdle posed by Question (3) remains to be surmounted. 

The first question posed in the introduction, Question (3), makes the point that one cannot have a theory for the nonequilibrium state based on the first entropy or its rate of production. It ought to be clear that the steady state, which corresponds to the most likely flux, x(x, i), gives neither the maximum nor the minimum of Eq. (61), the rate of first entropy production. From that equation, the extreme rates of first entropy production occur when x = oo. Theories that invoke the Principle of Minimum Dissipation, [10-12, 32] or the Principle of... 

Perhaps the most common computer simulation method for nonequilibrium systems is the nonequilibrium molecular dynamics (NEMD) method [53, 88]. This typically consists of Hamilton s equations of motion augmented with an artificial force designed to mimic particular nonequilibrium fluxes, and a constraint force or thermostat designed to keep the kinetic energy or temperature constant. Here is given a brief derivation and critique of the main elements of that method. 

The first and second integrals have their coordinate systems centered on the catalytic C and noncatalytic N spheres, respectively. The local nonequilibrium average microscopic density field for species a is pa(r) = [Y = 5(r - ( )) The solution of the diffusion equation can be used to estimate this nonequilibrium density, and thus the velocity of the nanodimer can be computed. The simple model yields results in qualitative accord with the MPC dynamics simulations and shows how the nonequilibrium density field produced by reaction, in combination with the different interactions of the B particles with the noncatalytic sphere, leads to directed motion [117],... 

Many solvents do not possess the simple structure that allows their effects to be modeled by the Langevin equation or generalized Langevin equation used earlier to calculate the TS trajectory [58, 111, 112]. Instead, they must be described in atomistic detail if their effects on the effective free energies (i.e., the time-independent properties) and the solvent response (i.e., the nonequilibrium or time-dependent properties) associated with the... 

Monte Carlo heat flow simulation, nonequilibrium molecular dynamics, 73-74, 77-81 multiparticle collision dynamics hydrodynamic equations, 105-107 macroscopic laws and transport coefficients, 102-104 single-particle friction and diffusion, 114-118... 

Two thermal nonequilibrium equations (the difference between the enthalpy of a phase and the saturation enthalpy)... 

The convective diffusion equations presented above have been used to model tablet dissolution in flowing fluids and the penetration of targeted macro-molecular drugs into solid tumors [5], In comparison with the nonequilibrium thermodynamics approach described below, the convective diffusion equations have the advantage of theoretical rigor. However, their mathematical complexity dictates a numerical solution in all but the simplest cases. 

Nonequilibrium thermodynamics provides a second approach to combined convection and diffusion problems. The Kedem-Katchalsky equations, originally developed to describe combined convection and diffusion in membranes, form the basis of this approach [6,7] ... 

Kedem-Katchalsky equations Phenomenological equations for combined convection and diffusion, derived from nonequilibrium thermodynamics. See Eqs. (19) and (20). 

---