Specific Features of the Simulation of the Particle Diffusion Processes in Spatially Periodic Fields

Abstract

This paper is devoted to the studies of the specific features of the simulation of the particle diffusion processes in space – periodic potentials using Langevin equations. Different methods used for the presetting of initial conditions and their effect on the obtained solutions have been analyzed. It is shown that the system is nonequilibrium for all the methods of the presetting of initial conditions during a certain time interval of ttrm. This interval is increased as 1/γ with a decrease in the friction coefficient. A reasonable description of the transient processes of particle transport and diffusion requires a preliminary system thermalization procedure. A new method of the presetting of initial conditions that provides the most accurate description of equilibrium system has been suggested. It consists in the generation of the initial particle coordinates and velocities that correspond to the equilibrium distribution of harmonic oscillators with a specified temperature. The use of such initial conditions enables the computations with a good accuracy using no thermalization procedure at T < 0.1. The classic method of the determination of diffusion coefficients D as a limit limt→ꝏ (σ2/t) has been analyzed. It was shown that the use of it for computer-aided calculations is limited by the restricted computational time. It results in that the computation of D under certain conditions becomes impossible. A new method was suggested for the determination of the diffusion coefficient using the linear approximation of the dependence of dispersion on time. This approximation can only be possible after the kinetic temperature attains its stationary value. The suggested method requires several orders of magnitude less time in comparison to the classic method. As a result, it enables the computation of the diffusion coefficient even in the cases of total previous failure. The obtained data are of great importance for correct simulation computations of diffusion processes and for the appropriate physical interpretations of obtained data.