scholarly journals Velocity autocorrelation function and self-diffusion coefficient in large molecular dynamics models of liquid argon and water

2021 ◽  
Vol 13 (2) ◽  
pp. 149-156
Author(s):  
Yuri I. Naberukhin ◽  
◽  
Alexey V. Anikeenko ◽  
Vladimir P. Voloshin ◽  
◽  
...  

Autocorrelation function of the particle velocity Z(t) is calculated using the molecular dynamics method in the models of liquid argon and water. The large size of the models (more than a hundred thousand particles) allowed us to trace these functions up to 50 picoseconds in argon and up to 10 picoseconds in water, and to achieve a calculation accuracy sufficient for analytical analysis of their shape. The difference in the determination of the self-diffusion coefficient using Einstein's law and the integral of Z(t) (Green-Kubo integral) is analyzed and it is shown to be 3% at best when t is of the order of several picoseconds. The asymptote of the function Z(t) in argon is close to the power law αt–3/2 predicted by hydrodynamics, but with an amplitude that depends on the time interval under consideration. In water, the asymptote of Z(t) has nothing in common with that in argon: it has α < 0 and the exponent is close to -5/2, and not to -3/2.

2018 ◽  
Vol 63 (4) ◽  
pp. 317 ◽  
Author(s):  
N. P. Malomuzh ◽  
K. S. Shakun ◽  
A. A. Kuznetsova

Long-time tails of the molecular velocity autocorrelation function (VACF) in liquid argon at temperatures higher and lower than the spinodal temperature have been analyzed. By considering the time dependence of the VACF, the self-diffusion and shear viscosity coefficients, and the Maxwell relaxation time are determined, as well as their changes when crossing the spinodal. It is shown that the characteristic changes in the temperature dependences of the indicated kinetic coefficients allow the spinodal position to be determined with a high accuracy. A possibility toapply the proposed method to other low-molecular liquids is considered. As an example, nitrogen and oxygen are used, for which the averaged potential of intermolecular interaction has the Lennard-Jones form.


2016 ◽  
Vol 30 (11) ◽  
pp. 1650064 ◽  
Author(s):  
Dipendra Bhandari ◽  
N. P. Adhikari

Molecular dynamics study of diffusion of two krypton atoms in 300 SPC/E water molecules at temperatures 293, 303, 313, 323 and 333 K has been carried out. Self-diffusion coefficient of krypton and water along with their mutual diffusion coefficients are estimated. Self-diffusion coefficient for krypton is calculated by using Mean Square Displacement (MSD) method and Velocity Autocorrelation (VACF) method, while that for water is calculated by using MSD method only. The mutual diffusion coefficient is estimated by using the Darken’s relation. The diffusion coefficients are found to follow the Arrhenius behavior. The structural properties of the system have been estimated by the study of solute–solute, solvent–solvent, and solute–solvent Radial Distribution Function (RDF).


2018 ◽  
Vol 16 (2) ◽  
Author(s):  
Niraj Kumar ◽  
Narayan Prasad Adhikari

Molecular Dynamics simulation was performed using 2 xenon atoms as solute and 300 water molecules as solvent. We have studied the structural properties as well as transport property. As structural properties, we have determined the radial distribution function (RDF) of xenon-xenon, xenon-water, and water-water interactions. Study of RDF of xenon-xenon and oxygen-oxygen interactions of water shows that there is hydrophobic behavior of xenon in the presence of water. We have studied the self diffusion coefficient of xenon, water, and mutual diffusion coefficients of xenon in water. The self diffusion coefficient of xenon was estimated using both mean-squared displacement (MSD) and velocity autocorrelation function (VACF), while only MSD was used for water. The temperature dependence of the diffusion coefficient of xenon and water were found to follow the Arrhenius behavior. The activation energies obtained are 12.156 KJ/mole with MSD and 14.617 KJ/mole with VACF in the temperature range taken in this study.


2014 ◽  
Vol 28 (14) ◽  
pp. 1450084 ◽  
Author(s):  
Keshav Sharma ◽  
Narayan P. Adhikari

We have carried out the molecular dynamics (MD) simulation to study the structural properties and to estimate the diffusivity of molecular nitrogen ( N 2) gas (solute) in extended simple point charge model (SPC/E) water (solvent) with N 2 mole fraction of 0.018 at different temperatures. For the structural properties of the system, we have determined radial distribution function (RDF). The solute–solute, solute–solvent and solvent–solvent RDF have been evaluated. Self-diffusion coefficient of N 2 was estimated by evaluating mean-squared displacement (MSD) and velocity autocorrelation function (VACF) separately. The diffusion coefficients obtained from the two methods agree within 3%. The results are in agreement with the experimentally determined values within 10%. The self-diffusion coefficient of water ( H 2 O ) was also estimated by evaluating MSD. Mutual diffusion coefficient of the system have also been estimated invoking Darken's relation. The temperature dependance of the diffusion coefficients were found to follow Arrhenius relation.


Sign in / Sign up

Export Citation Format

Share Document