Nonequilibrium physics

Unlike the analysis of the structure and equilibrium concentrations of point defects, questions of defect motion take us directly to the heart of some of the most important unsolved questions in nonequilibrium physics. In particular, while there are a host of useful empirical constructs built around the idea of Arrhenius activated processes with rates determined generically by expressions of the form... 

In conclusion, we cannot escape the feeling that we live in an age of transition, an age that demands a better understanding of our environment. We must find and explore new resources, and achieve a less destructive coexistence with Nature. We cannot anticipate the outcome of this period of transition, but it is clear that science and especially nonequilibrium physics, is bound to play an increasingly important role in our effort to meet the challenge of understanding and reshaping our global environment. 

This new field is in a state of explosive development, if I can say so. A few years ago, when we studied a simple model, known since in the literature as the Brusselator, we had to do the mathematics by ourselves when we needed bifurcations, we had to adapt the mathematical tools in our own amateurish way. Today, this field is in full blossom. I think that this is an interesting example of interaction between physics and mathematics. In fact, the new developments in nonequilibrium physics, coupled with advances in modern dynamics, had a revigorizing influence on nonlinear mathematics, which we could compare to the progress in mathematics induced earlier by other fields like relativity or quantum theory. 

Classical nucleation theory has been successfully applied to the study of ice nucleation in pure water using a new equation of state for supercooled water [42]. Application of the theory to atmospheric solution droplets via eq 3 would require that the parameters C,Tmp (iw) be known however, elucidation of the equilibrium and nonequilibrium physical chemistry of these droplets in the atmospherically relevant thermodynamic regimes represents a formidable challenge in current atmospheric research [43, 44]. 

Actual crystal planes tend to be incomplete and imperfect in many ways. Nonequilibrium surface stresses may be relieved by surface imperfections such as overgrowths, incomplete planes, steps, and dislocations (see below) as illustrated in Fig. VII-5 [98, 99]. The distribution of such features depends on the past history of the material, including the presence of adsorbing impurities [100]. Finally, for sufficiently small crystals (1-10 nm in dimension), quantum-mechanical effects may alter various physical (e.g., optical) properties [101]. 

Holian B L 1996 The character of the nonequilibrium steady state beautiful formalism meets ugly reality Monte Carlo and Molecular Dynamics of Condensed Matter Systems, vol 49, ed K Binder and G Ciccotti (Bologna Italian Physical Society) pp 791-822... 

The volume of substance in a composite material that exists in a nonequilibrium state due to its proximity to an interface has been termed an interphase [1]. The interphase is a zone of distinct composition and properties formed by chemical or physical processes such as interdiffusion of mutually soluble components or chemical interaction between reactive species. 

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 this review we put less emphasis on the physics and chemistry of surface processes, for which we refer the reader to recent reviews of adsorption-desorption kinetics which are contained in two books [2,3] with chapters by the present authors where further references to earher work can be found. These articles also discuss relevant experimental techniques employed in the study of surface kinetics and appropriate methods of data analysis. Here we give details of how to set up models under basically two different kinetic conditions, namely (/) when the adsorbate remains in quasi-equihbrium during the relevant processes, in which case nonequilibrium thermodynamics provides the needed framework, and (n) when surface nonequilibrium effects become important and nonequilibrium statistical mechanics becomes the appropriate vehicle. For both approaches we will restrict ourselves to systems for which appropriate lattice gas models can be set up. Further associated theoretical reviews are by Lombardo and Bell [4] with emphasis on Monte Carlo simulations, by Brivio and Grimley [5] on dynamics, and by Persson [6] on the lattice gas model. 

An electrochemical reaction is said to be polarized or retarded when it is limited by various physical and chemical factors. In other words, the reduction in potential difference in volts due to net current flow between the two electrodes of the corrosion cell is termed polarization. Thus, the corrosion cell is in a state of nonequilibrium due to this polarization. Figure 4-415 is a schematic illustration of a Daniel cell. The potential difference (emf) between zinc and copper electrodes is about one volt. Upon allowing current to flow through the external resistance, the potential difference falls below one volt. As the current is increased, the voltage continues to drop and upon completely short circuiting (R = 0, therefore maximum flow of current) the potential difference falls toward about zero. This phenomenon can be plotted as a polarization diagram shown in Figure 4-416. 

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. 

Introduction.—Statistical physics deals with the relation between the macroscopic laws that describe the internal state of a system and the dynamics of the interactions of its microscopic constituents. The derivation of the nonequilibrium macroscopic laws, such as those of hydrodynamics, from the microscopic laws has not been developed as generally as in the equilibrium case (the derivation of thermodynamic relations by equilibrium statistical mechanics). The microscopic analysis of nonequilibrium phenomena, however, has achieved a considerable degree of success for the particular case of dilute gases. In this case, the kinetic theory, or transport theory, allows one to relate the transport of matter or of energy, for example (as in diffusion, or heat flow, respectively), to the mechanics of the molecules that make up the system. 

Considering the similarity between Figs. 1 and 2, the electrode potential E and the anodic dissolution current J in Fig. 2 correspond to the control parameter ft and the physical variable x in Fig. 1, respectively. Then it can be said that the equilibrium solution of J changes the value from J - 0 to J > 0 at the critical pitting potential pit. Therefore the critical pitting potential corresponds to the bifurcation point. From these points of view, corrosion should be classified as one of the nonequilibrium and nonlinear phenomena in complex systems, similar to other phenomena such as chaos. 

In all these treatments, nonequilibrium fluctuation plays the most important role. This is defined as the fluctuation of a physical quantity that deviates from the standard state determined by the nonzero flux in a nonequilibrium state. Such fluctuation has a kind of symmetry in that the area average is equal to zero although the flux changes locally. Therefore, macroscopically, such fluctuation does not affect the flux itself. This means that the flux must be determined a priori and is indifferent to the fluctuations. 

The newly formed equilibrium, however, is broken easily and incessantly by the thermal motion of solution particles. Since the electrode system is not in Nemstian equilibrium at the potential, such a breakdown (nonequilibrium fluctuation) produces pitting dissolution. The physical quantities related to the dissolution fluctuate on one side of the electrostatic equilibrium, that is, the fluctuations take place toward the direction in which the reaction proceeds. 

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.)...

Figure 40. Plot of the fluctuation-diffusion current J vs. iwr.91 id is the slope of the fluctuation-diffusion current given by Eq. (115). Solid and dotted lines correspond to the theoretical and experimental results, respectively. (NiCljJ = 0.1 mol nT3. [NsCl] = 7 mol m 3. V = 0.1 V, T= 300 K. (Reprinted from M. Asanuma and R. Aogaki, Nonequilibrium fluctuation theory on pitting dissolution, n. Determination of surface coverage of nickel passive film, J. Chem. Phys. 106, 9938, 1997, Fig. 8. Copyright 1997, American Institute of Physics.)...

---