Correlation functions nonequilibrium

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]

The present theory can be placed in some sort of perspective by dividing the nonequilibrium field into thermodynamics and statistical mechanics. As will become clearer later, the division between the two is fuzzy, but for the present purposes nonequilibrium thermodynamics will be considered that phenomenological theory that takes the existence of the transport coefficients and laws as axiomatic. Nonequilibrium statistical mechanics will be taken to be that field that deals with molecular-level (i.e., phase space) quantities such as probabilities and time correlation functions. The probability, fluctuations, and evolution of macrostates belong to the overlap of the two fields. 

Arguably a more practical approach to higher-order nonequilibrium states lies in statistical mechanics rather than in thermodynamics. The time correlation function gives the linear response to a time-varying field, and this appears in computational terms the most useful methodology, even if it may lack the... 

Sharp and Lohr proposed recently a somewhat different point of view on the relation between the electron spin relaxation and the PRE (126). They pointed out that the electron spin relaxation phenomena taking a nonequilibrium ensemble of electron spins (or a perturbed electron spin density operator) back to equilibrium, described in Eqs. (53) and (59) in terms of relaxation superoperators of the Redfield theory, are not really relevant for the PRE. In an NMR experiment, the electron spin density operator remains at, or very close to, thermal equilibrium. The pertinent electron spin relaxation involves instead the thermal decay of time correlation functions such as those given in Eq. (56). The authors show that the decay of the Gr(T) (r denotes the electron spin vector components) is composed of a sum of contributions... 

Nonequilibrium Aging State (NEAS). The system is initially prepared in a nonequilibrium state and put in contact with the sources. The system is then allowed to evolve alone but fails to reach thermal equilibrium in observable or laboratory time scales. In this case the system is in a nonstationary slowly relaxing nonequilibrium state called aging state and is characterized by a very small entropy production of the sources. In the aging state two-times correlations decay slower as the system becomes older. Two-time correlation functions depend on both times and not just on their difference. 

Let us make a final comment, concerning the violation of the Green-Kubo relation. There is a close connection between the breakdown of this fundamental prescription of nonequilibrium statistical physics and the breakdown of the agreement between the density and trajectory approach. We have seen that the CTRW theory, which rests on trajectories undergoing abrupt and unpredictable jumps, establishes the pdf time evolution on the basis of v /(f), whereas the density approach to GME, resting on the Liouville equation, either classical or quantum, and on the convenient contraction over the irrelevant degrees of freedom, eventually establishes the pdf time evolution on the basis of a correlation function, the correlation function in the dynamical case... 

If H (a,b) as provided by Eq. (3.79) is positive definite, we shall not find any theoretical difficulty in proposing the CFP as a calculation technique for this nonequilibrium process. Throughout this volume we shall consider stochastic variables A with vanishing mean value at equilibrium, that is, (A) = 0. The correlation function t) = A A t) / A A then shares the standard formal properties of the normalized equilibrium correlation functions ... 

While there is no unique criterion for choosing 4 E, the selection must lead to an accurate theory of solvation dynamics without invoking two-time many-point correlation functions. We have found that this goal can be achieved with a new theory for the nonequilibrium distribution function in which the renormalized solute-solvent interactions enter linearly. In this theory and are chosen such that the renormalized linear response theory accurately describes the essential solute-solvent static correlations that rule the equilibrium solvation both at t = 0 (when solvent is in equilibrium with the initial charge distribution of the solute) and at 1 = oc (when the solvent has reached equilibrium with the new solute charge distribution). ... 

Equation [116] is the central result of this section and is one of the most important equations in nonequilibrium statistical mechanics. It is from this point that one makes contact with the more familiar relations relating transport coefficients to time correlation functions under the guise of the so-called Green-Kubo formulas.i2.39,40... 

Recent years have seen the extensive application of computer simulation techniques to the study of condensed phases of matter. The two techniques of major importance are the Monte Carlo method and the method of molecular dynamics. Monte Carlo methods are ways of evaluating the partition function of a many-particle system through sampling the multidimensional integral that defines it, and can be used only for the study of equilibrium quantities such as thermodynamic properties and average local structure. Molecular dynamics methods solve Newton s classical equations of motion for a system of particles placed in a box with periodic boundary conditions, and can be used to study both equilibrium and nonequilibrium properties such as time correlation functions. 

In Section II we will review thermodynamics and the fluctuation-dissipation theorem for excess heat production based on the Boltzmann equilibrium distribution. We will also mention the nonequilibrium work relation by Jarzynski. In Section III, we will extend the fluctuation-dissipation theorem for the superstatisitcal equilibrium distribution. The fluctuation-dissipation theorem can be written as a superposition of correlation functions with different temperatures. When the decay constant of a correlation function depends on temperature, we can expect various behaviors in the excess heat. In Section IV, we will consider the case of the microcanonical equilibrium distribution. We will numerically show the breaking of nonergodic adiabatic invariant in the mixed phase space. In the last section, we will conclude and comment. 

The nonequilibrium solvation function iS (Z), which is directly observable (e.g. by monitoring dynamic line shifts as in Fig. 15.2), is seen to be equal in the linear response approximation to the time correlation function, C(Z), of equilibrium fluctuations in the solvent response potential at the position of the solute ion. This provides a route for generalizing the continuum dielectric response theory of Section 15.2 and also a convenient numerical tool that we discuss further in the next section. 

Fig. 15.3 The nonequilibrium solvation function S(t) (full lines) and the solvation correlation functions C(i) for a model solute ion of diameter 3.1 A in acetonitrile computed with the positive solute (dotted line) and neutral solute (dashed line). (From M. Maroncelli, J. Chem. Phys. 94, 2084 (1991).)...

In the spirit of the discussion of Section II, we may consider the time evolution of the correlation functions of 6/i (X ) rather than its nonequilibrium average. These correlation functions will satisfy the same kinetic equations as the 5F (1, t). This is the approach we take here. We now outline a method for explicitly constructing K p(l] z). 

We discussed some aspects of the responses of chemical systems, linear or nonlinear, to perturbations on several earlier occasions. The first was the responses of the chemical species in a reaction mechanism (a network) in a nonequilibrium stable stationary state to a pulse in concentration of one species. We referred to this approach as the pulse method (see chapter 5 for theory and chapter 6 for experiments). Second, we studied the time series of the responses of concentrations to repeated random perturbations, the formulation of correlation functions from such measurements, and the construction of the correlation metric (see chapter 7 for theory and chapter 8 for experiments). Third, in the investigation of oscillatory chemical reactions we showed that the responses of a chemical system in a stable stationary state close to a Hopf bifurcation are related to the category of the oscillatory reaction and to the role of the essential species in the system (see chapter 11 for theory and experiments). In each of these cases the responses yield important information about the reaction pathway and the reaction mechanism. 

The two previous secfions were devoted to modeling quantum resonances by means of effective Hamiltonians. From the mathematical point of view we have used two principal tools projection operators that permit to focus on a few states of interest and analytic continuation that allows to uncover the complex energies. Because the time-dependent Schrodinger equation is formally equivalent to the Liouville equation, it is attractive to try to solve the Liouville equation using the same tools and thus establishing a link between the dynamics and the nonequilibrium thermodynamics. For that purpose we will briefly recall the definition of the correlation functions which are similar to the survival and transition amplitudes of quantum mechanics. Then two models of regression of a fluctuation and of a chemical kinetic equation including a transition state will be presented. 

The nonequilibrium concentration fluctuations can be measured experimentally by dynamic light scattering. Fluctuating hydrodynamics predicts that the time-dependent correlations function C(k,t) ofthe scattered light is given by... 

Dynamic Density Functional Theory (DDFT), Fig. 1 Illustration of the local equilibrium approximation involved in the development of the DDFT. The left-hand side illustrates the nonequilibrium evolution of the density p(r t) thin lines) up to time t thick line). For the time evolution, the equal-time correlation function g(r, r t) is... 

The local equilibrium approximation for the two-point correlation function involved in the development of the DDET has two issues. Eirst, it is not a priori clear when it is justifiable to approximate the nonequilibrium correlations by equilibrium correlations. It has been shown that there are cases in which this approximation breaks down, in particular in driven steady-state systems. Second, the equilibrium sum rule in... 

Here T is the local-equilibrium temperature. In extended irreversible thermodynamics fluxes are independent variables. The kinetic temperature associated to the three spatial directions of along the flow, along the velocity gradient, and perpendicular to the previous to directions may be different from each other. To define temperature from the entropy is the most fundamental definition, and the nonequilibrium temperature may come from the derivative of a nonequilibrium entropy du/dS) -p. Effective nonequilibrium temperature may be defined from the fluctuation-dissipation theorem relating response function and correlation function. 

In the next four sections, we discuss the four principal types of application of molecular dynamics. Section 3 very briefly describes the problem of the approach to equilibrium. Section 4 deals with the evaluation of equilibrium thermodynamic functions through a discussion of the dynamical equation of state. In Section 5, we consider the evaluation of equilibrium time correlation functions, detailing the application of the combined Monte Carlo-molecular-dynamics method to the time correlation functions for self-diffusion. Section 6 deals with nonequilibrium molecular dynamics and in particular with a calculation for self-diffusion. 

---