Transport and Helfand moments in the Lennard-Jones fluid. I. Shear viscosity.

The authors propose a new method, the Helfand-moment method, to compute the shear viscosity by equilibrium molecular dynamics in periodic systems. In this method, the shear viscosity is written as an Einstein-type relation in terms of the variance of the so-called Helfand moment. This quantity is modified in order to satisfy systems with periodic boundary conditions usually considered in molecular dynamics. They calculate the shear viscosity in the Lennard-Jones fluid near the triple point thanks to this new technique. They show that the results of the Helfand-moment method are in excellent agreement with the results of the standard Green-Kubo method.

Transport and Helfand moments in the Lennard-Jones fluid. II. Thermal conductivity.

The thermal conductivity is calculated with the Helfand-moment method in the Lennard-Jones fluid near the triple point. The Helfand moment of thermal conductivity is here derived for molecular dynamics with periodic boundary conditions. Thermal conductivity is given by a generalized Einstein relation with this Helfand moment. The authors compute thermal conductivity by this new method and compare it with their own values obtained by the standard Green-Kubo method. The agreement is excellent.

Viscosity in molecular dynamics with periodic boundary conditions.

We report a study of viscosity by the method of Helfand moment in systems with periodic boundary conditions. We propose a new definition of Helfand moment which takes into account the minimum image convention used in molecular dynamics with periodic boundary conditions. Our Helfand-moment method is equivalent to the method based on the Green-Kubo formula and is not affected by ambiguities due to the periodic boundary conditions. Moreover, in hard-ball systems, our method is equivalent to that developed by Alder, Gass, and Wainwright [J. Chem. Phys. 53, 3813 (1970)]. We apply and verify our method in a fluid composed of N> or =2 hard disks in elastic collisions. We show that the viscosity coefficients already take values in good agreement with Enskog's theory for N=2 hard disks in a hexagonal geometry.

Viscosity in the escape-rate formalism.

We apply the escape-rate formalism to compute the shear viscosity in terms of the chaotic properties of the underlying microscopic dynamics. A first-passage problem is set up for the escape of the Helfand moment associated with viscosity out of an interval delimited by absorbing boundaries. At the microscopic level of description, the absorbing boundaries generate a fractal repeller. The fractal dimensions of this repeller are directly related to the shear viscosity and the Lyapunov exponent, which allows us to compute its values. We apply this method to the Bunimovich-Spohn minimal model of viscosity which is composed of two hard disks in elastic collision on a torus. These values are in excellent agreement with the values obtained by other methods such as the Green-Kubo and Einstein-Helfand formulas.