Universal dynamic exponent at the liquid-gas transition from molecular dynamics

The liquid-gas system is expected to exhibit distinct dynamic behavior in the fluid's critical region (model H). We present molecular dynamics simulations of a Lennard-Jones fluid model starting from specially designed, near-equilibrium, initial conditions. By following the fluid's relaxation towards equilibrium, we calculate the requisite transport coefficients in the critical region. The results yield the scaling behavior of the thermal diffusion coefficient D(T) approximately xi(-1.023+/-0.018) (xi is the correlation length) and a nonconventional divergent heat conductivity, all of which are in accord with mode-coupling and renormalization group predictions, as well as some experimental data.