Critical anomaly and finite size scaling of the self-diffusion coefficient for Lennard Jones fluids by non-equilibrium molecular dynamic simulation  

作　　者：Ahmed Asad 吴江涛 

机构地区：[1]State Key Laboratory of Multiphase Flow in Power Engineering, Xi'an Jiaotong University

出　　处：《Chinese Physics B》2011年第10期362-367,共6页中国物理B（英文版）

基　　金：supported by the National Natural Science Foundation of China (Grant No. 51076128);the National High Technology Research and Development Program of China (Grant No. 2009AA05Z107)

摘　　要：We use non-equilibrium molecular dynamics simulations to calculate the self-diffusion coefficient, D, of a Lennard Jones fluid over a wide density and temperature range. The change in self-diffusion coefficient with temperature decreases by increasing density. For density ρ＊ = ρσ3 = 0.84 we observe a peak at the value of the self-diffusion coefficient and the critical temperature T＊ = kT/ε = 1.25. The value of the self-diffusion coefficient strongly depends on system size. The data of the self-diffusion coefficient are fitted to a simple analytic relation based on hydrodynamic arguments. This correction scales as N-α, where α is an adjustable parameter and N is the number of particles. It is observed that the values of a 〈 1 provide quite a good correction to the simulation data. The system size dependence is very strong for lower densities, but it is not as strong for higher densities. The self-diffusion coefficient calculated with non-equilibrium molecular dynamic simulations at different temperatures and densities is in good agreement with other calculations fronl the literature.We use non-equilibrium molecular dynamics simulations to calculate the self-diffusion coefficient, D, of a Lennard Jones fluid over a wide density and temperature range. The change in self-diffusion coefficient with temperature decreases by increasing density. For density ρ＊ = ρσ3 = 0.84 we observe a peak at the value of the self-diffusion coefficient and the critical temperature T＊ = kT/ε = 1.25. The value of the self-diffusion coefficient strongly depends on system size. The data of the self-diffusion coefficient are fitted to a simple analytic relation based on hydrodynamic arguments. This correction scales as N-α, where α is an adjustable parameter and N is the number of particles. It is observed that the values of a 〈 1 provide quite a good correction to the simulation data. The system size dependence is very strong for lower densities, but it is not as strong for higher densities. The self-diffusion coefficient calculated with non-equilibrium molecular dynamic simulations at different temperatures and densities is in good agreement with other calculations fronl the literature.

关 键 词：self-diffusion coefficient non-equilibrium molecular dynamic simulation Lennard Jonesfluid critical dynamics 

分 类 号：O35[理学—流体力学]