Kubo relation

This can be obtained also from the velocity according to the Green-Kubo relation ... 

Sindzinger and Gillan have calculated the thermal conductivity for NaCl and KCl melts as well as for sohds on the basis of MD simulations in Ml thermal equilibrium using the Green-Kubo relations (Table 17). In a single molten salt system, the local fluxes jz and of charge and energy... 

This result is a very stringent test of the present expression for the steady-state probability distribution, Eq. (160). There is one, and only one, exponent that is odd, linear in Xr, and that satisfies the Green-Kubo relation. 

Computer simulations of transport properties using Green-Kubo relations... 

T. Ihle and D. M. Kroll, Stochastic rotation dynamics. I. Formalism, Galilean invariance, and Green—Kubo relations, Phys. Rev. E 67, 066705 (2003). 

As A x was supposed stationary the integral is independent of time. The effect of the fluctuations is therefore to renormalize A0 by adding a constant term of order a2 to it. The added term is the integrated autocorrelation function of At. In particular, if one has a non-dissipative system described by A0, this additional term due to the fluctuations is usually dissipative. This relation between dissipation and the autocorrelation function of fluctuations is analogous to the Green-Kubo relation in many-body systems 510 but not identical to it, because there the fluctuations are internal, rather than added as a separate term as in (2.1). 

W. H. Miller The expression for the reaction rate (in terms of a flux-flux autocorrelation function) obtained by myself, Schwartz, and Tromp in 1983 is very similar (though not identical) to the one given earlier by Yamamoto. It is also an example of Green-Kubo relations. 

A further test of the pair potentials used is the translational diffusion coefficients, D. The coefficients in Tables II, III, and IV were all evaluated using the Einstein-Kubo relation ... 

To proceed it is necessary to evaluate this coefficient. One possible procedure is to use a measured value of the transport coefficient which is related to Cu(0) through a Kubo relation. Another possibility is to relate Xn to the moments of C (to). 

Thus, we find that a dilemma is created. We have to make a choice between two possible ways out. The first is that the Green-Kubo relation of Eq. (180) does not apply, and the second is that the Green-Kubo relation of Eq. (180) does apply, but the common sense expectation that a weak perturbation creates a finite current is violated. 

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

We note that in Ref. 87 the projection method was used to study the effect of friction on a process that, in the free case, is known to produce a Levy walk. As already discussed in Section VIII, the key point has been that the violation of the Green-Kubo relation must imply a different form of linear response, used later to justify the WS form of noncanonical equilibrium. It is worth mentioning that an interesting result of that article has been the following equation of motion ... 

In the molecular dynamics calculations the trajectories of methane molecules in the pore are followed using the equation of motion with appropriate temperature control. A diffuse reflection condition is applied at the pore wall. For the EMD simulations a collective transport coefficient obtained from autocorrelation of the fluctuating axial streaming velocity via a Green-Kubo relation [S]... 

In the earliest paper on FRs, a footnote remarked that it is possible to obtain the Green-Kubo relations from them if the fluctuations have a Gaussian distribution. ... 

When we evaluate the Green-Kubo relations for the transport coefficients we solve the equations of motion for the molecules. They are often modelled as rigid bodies. Therefore we review some of basic definition of rigid body dynamics [10]. The centres of mass of the molecules evolve according to the ordinary Newtonian equations of motion. The motion in angular space is more complicated. Three independent coordinates a, =(a,a,-2, ,3), i = 1, 2,. . N where N is the number of molecules, are needed to describe the orientation of a rigid body. (Note that a, is not a vector because it does not transform like a vector when the coordinate system is rotated.) The rate of change of a, is... 

Inserting this expression into the linear response relation (2.17) and comparing it to the Green-Kubo relation (3.3) gives the following thermal conductivity,... 

Comparison of the thermal conductivities of prolate (p) and oblate (o) nematic liquid crystals. The entries for zero field have been obtained by using the Green-Kubo relation (3.3). The entries for finite field have been obtained by applying the heat flow algorithm (3.5). Note that the EMD GK estimates and the NEMD estimates agree within the statistical error. 

One also finds that fixing the director generates a new equilibrium ensemble where the Green-Kubo relations for the viscosities are considerably simpler compared to the conventional canonical ensemble. They become linear functions of time correlation function integrals instead of rational functions. The reason for this is that all the thermodynamic forces are constants of motion and all the thermodynamic fluxes are zero mean fluctuating phase functions in the constrained ensemble. 

S. Sarman, Green-Kubo Relations for the Viscosity of Biaxial Nematic Liquid Crystals, J. Chem. Phys. 105 (1996) 4211 ... 

Silicon thin film thermal conductivities are predicted using equilibrium molecular dynamics and the Grccn-Kubo relation. Periodic boundary conditions are applied in the direetions parallel to the thin film surfaees (Fig. 5). Atoms near the surfaces of the thin film are subjeeted to the above-described repulsive potential in addition to the Stillinger-Weber potential [75]. Simulations were also performed adding to each surface four layers of atoms kept frozen at their crystallographic positions, in order to eompare the dependence of the predieted thermal eonduetivities on the surface boundary eonditions. We found that the thermal eonduetivities obtained using frozen atoms or the repulsive potential are identical within the statistical deviations, exeept for the in-plane thermal eonduetivity of films with thickness less than 10 nm [79]. Therefore, in the present study, we present only the predietions obtained with the repulsive potential. 

L. Bocquet and J. L. Barrat, Phys. Rev. E, 49,3079 (1994). Hydrodynamic Boundary Conditions, Correlation Functions and Kubo Relations for Confined Fluids. 

Liquid simulation studies have been essential in assessing the applicability of various fluctuation relations to real physical systems. These are important relations in nonequilibrium statistical mechanics that are valid far from equilibrium and can be used to derive Green-Kubo relations for transport coefficients.223,224 They show how thermodynamic irreversibility emerges from... 

Equilibrium MD simulations of self-diffusion coefficients, shear viscosity, and electrical conductivity for C mim][Cl] at different temperatures were carried out [82] The Green-Kubo relations were employed to evaluate the transport coefficients. Compared to experiment, the model underestimated the conductivity and self-diffusion, whereas the viscosity was over-predicted. These discrepancies were explained on the basis of the rigidity and lack of polarizability of the model [82], Despite this, the experimental trends with temperature were remarkably well reproduced. The simulations reproduced remarkably well the slope of the Walden plots obtained from experimental data and confirmed that temperature does not alter appreciably the extent of ion pairing [82],... 

The structure of a water/silicon interface was studied (Ursenbach et al, 1997), in addition to a water/copper interface (Halley et al, 1998) and a water/palladium interface (Klesingeza/., 1998). Finally, two studies have used Car-Parrinello simulations in conjunction with the Green-Kubo relations to calculate viscosities in liquid metals (Alfe and Gillan, 1998 Stadler et al, 1999). 

---