Cutoff Radius Effect of Water Configuration Using the Wolf Method

Author(s):  
Kazuaki Takahashi ◽  
Tetsu Narumi ◽  
Kenji Yasuoka

Molecular dynamics simulation has been applied for water to compare the Wolf method to the IPS method and the Ewald sum by evaluating the diffusion coefficient and liquid structure. In our previous study, we applied the IPS method for bulk water and found notable deviation of the radial distribution function g(r). The Wolf method gives a good estimation for the potential energy and the self-diffusion coefficient at a cutoff radius, rc, greater than 2.2 nm while avoiding the notable deviation of g(r) which appeared in the case of IPS. The distance dependent Kirkwood factor Gk(r) was also calculated, and the truncation of a long-range interaction of the cutofflike method (such as cutoff with or without the switch function and the reaction field) show serious shortcomings for dipole-dipole correlations in bulk water systems. This was observed by comparing the shape to that of the Ewald sum. Gk(r) of the cutofflike method greatly deviates from that of the Ewald sum. However, the discrepancy of Gk(r) for the Wolf method was found to be much less than that of other typical cutoff-like methods. We conclude that the Wolf method is an adequately accurate technique for estimating transport coefficients and the liquid structure of water in a homogeneous system at long cutoff distances.

Author(s):  
Lawrence Novak

Rate-based models suitable for equipment or transport-reaction modeling require a capability for predicting transport coefficients over a sufficient range of temperature and pressure. This paper demonstrates a relatively simple novel approach to correlate and estimate transport coefficients for pure components over the entire fluid region.The use of Chapman-Enskog transport coefficients for reducing self-diffusion coefficient and viscosity to dimensionless form results in relatively simple mathematical relationships between component dimensionless transport coefficients and residual entropy over the entire fluid region. Dimensionless self-diffusion coefficients and viscosities were calculated from extensive molecular dynamics simulation data and experimental data on argon, methane, ethylene, ethane, propane, and n-decane. These dimensionless transport coefficients were plotted against dimensionless residual entropy calculated from highly accurate reference equations of state.Based on experimental data, the new scaling model introduced here shows promise as: (1) an equation of state-based transport coefficient correlation over the entire fluid region (liquid, gas, and critical fluid), (2) a component transport coefficient correlation for testing transport data consistency, and (3) a component transport coefficient correlation for interpolation and extrapolation of self-diffusion coefficient and viscosity.


2019 ◽  
Vol 21 (35) ◽  
pp. 19163-19171 ◽  
Author(s):  
Mingbing Wu ◽  
Wei Wei ◽  
Xiaowei Liu ◽  
Kang Liu ◽  
Song Li

Water confined in nanochannels can be stretched with the variation of external pressure, leading to the more disordered microstructure and higher diffusion coefficient than bulk water.


Molecules ◽  
2021 ◽  
Vol 26 (13) ◽  
pp. 4030
Author(s):  
Gengbiao Chen ◽  
Zhiwen Liu

The diffusion behavior of fluid water in nanochannels with hydroxylation of silica gel and silanization of different modified chain lengths was simulated by the equilibrium molecular dynamics method. The diffusion coefficient of fluid water was calculated by the Einstein method and the Green–Kubo method, so as to analyze the change rule between the modification degree of nanochannels and the diffusion coefficient of fluid water. The results showed that the diffusion coefficient of fluid water increased with the length of the modified chain. The average diffusion coefficient of fluid water in the hydroxylated nanochannels was 8.01% of the bulk water diffusion coefficient, and the diffusion coefficients of fluid water in the –(CH2)3CH3, –(CH2)7CH3, and –(CH2)11CH3 nanochannels were 44.10%, 49.72%, and 53.80% of the diffusion coefficients of bulk water, respectively. In the above four wall characteristic models, the diffusion coefficients in the z direction were smaller than those in the other directions. However, with an increase in the silylation degree, the increased self-diffusion coefficient due to the surface effect could basically offset the decreased self-diffusion coefficient owing to the scale effect. In the four nanochannels, when the local diffusion coefficient of fluid water was in the range of 8 Å close to the wall, Dz was greater than Dxy, and beyond the range of 8 Å of the wall, the Dz was smaller than Dxy.


2021 ◽  
Author(s):  
Kazushi Fujimoto ◽  
Tetsuro Nagai ◽  
Tsuyoshi Yamaguchi

<div>The position-dependent diffusion coefficient along with free energy profile are important parameters needed to study mass transport in heterogeneous systems such as biological and polymer membranes, and molecular dynamics (MD) calculation is a popular tool to obtain them. Among many methodologies, the Marrink-Berendsen (MB) method is often employed to calculate the position-dependent diffusion coefficient, in which the autocorrelation function of the force on a fixed molecule is related to the friction on the molecule. However, the diffusion coefficient is shown to be affected by the period of the removal of the center-of-mass velocity, which is necessary when performing MD calculations using the Ewald method for Coulombic interaction. We have clarified theoretically in this study how this operation affects the diffusion coefficient calculated by the MB method, and the theoretical predictions are proven by MD calculations. Therefore, we succeeded in providing guidance on how to select an appropriate the period of the removal of the center-of-mass velocity in estimating the position-dependent diffusion coefficient by the MB method. This guideline is applicable also to the Woolf-Roux method.</div>


2021 ◽  
Vol 1026 ◽  
pp. 39-48
Author(s):  
Han Bing He ◽  
Yu Si Wang ◽  
Ze Xiang Luo ◽  
Jing Zeng

The effect of different additives Ni, Fe, Cu on the structure and properties of electrolyte system 78% Na3AlF6- -9.5%AlF3-5.0%CaF2-7.5%Al2O3 at 1200K and 1.01Mpa was studied by molecular dynamics method. The radial distribution function, coordination number, diffusion coefficient, conductivity, and viscosity of the system were discussed in detail. The results demonstrated that the order of the self-diffusion coefficient of ions in the electrolyte system is: Na+ > F- > O2- > Ca2+ >Al3+. The addition of Ni and Fe connected the free aluminum composite ion groups in the system through fluorine bridges, which enhanced the interaction between Al3+ and Al3+. The addition of Cu weakened the interaction between Al3 + and Al3+ and the F-. The interaction between Al3+ and Na+, [AlF7]4- ionic groups might appeared in the melt system. After adding NiO, Fe2O3, and Cu, the electrical conductivity of the system increased, and the viscosity decreased. The research work revealed the influence of Ni, Fe, Cu on the ion existence form, mobility, inter-ion interaction and diffusion mechanism of cryolite molten salt system, which has important guiding significance for aluminum electrolysis production.


2018 ◽  
Vol 8 (10) ◽  
pp. 1874 ◽  
Author(s):  
Jie Wu ◽  
Jia Wang ◽  
Haiou Ni ◽  
Guimin Lu ◽  
Jianguo Yu

Molten chloride salts are the main components in liquid metal batteries, high-temperature heat storage materials, heat transfer mediums, and metal electrolytes. In this paper, interest is centered on the influence of the LiCl component and temperature on the local structure and transport properties of the molten LiCl-NaCl-KCl system over the temperature range of 900 K to 1200 K. The liquid structure and properties have been studied across the full composition range by molecular dynamics (MD) simulation of a sufficient length to collect reliable values, such as the partial radial distribution function, angular distribution functions, coordination numbers distribution, density, self-diffusion coefficient, ionic conductivity, and shear viscosity. Densities obtained from simulations were underestimated by an average 5.7% of the experimental values. Shear viscosities and ionic conductivity were in good agreement with the experimental data. The association of all ion pairs (except for Li-Li and Cl-Cl) was weakened by an increasing LiCl concentration. Ion clusters were formed in liquids with increasing temperatures. The self-diffusion coefficients and ionic conductivity showed positive dependences on both LiCl concentration and temperature, however, the shear viscosity was the opposite. By analyzing the hydrodynamic radii of each ion and the coordination stability of cation-anion pairs, it was speculated that ion clusters could be the cation-anion coordinated structure and affected the macro properties.


2020 ◽  
Author(s):  
Tianyu Zhu ◽  
Troy Van Voorhis

<p>The dipole moment of a single water molecule in liquid water has been a critical concept for understanding water’s dielectric properties. In this work, we investigate the dipole moment of liquid water through a self-attractive Hartree (SAH) decomposition of total electron density computed by density functional theory, on water clusters sampled from ab initio molecular dynamics simulation of bulk water. By adjusting one parameter that controls the degree of density localization, we reveal two distinct pictures of water dipoles that are consistent with bulk dielectric properties: a localized picture with smaller and less polarizable monomer dipoles, and a delocalized picture with larger and more polarizable monomer dipoles. We further uncover that the collective dipole-dipole correlation is stronger in the localized picture and is key to connecting individual dipoles with bulk dielectric properties. Based on these findings, we suggest considering both individual and collective dipole behaviors when studying the dipole moment of liquid water, and propose new design strategies for developing water models.</p>


2020 ◽  
Author(s):  
Tianyu Zhu ◽  
Troy Van Voorhis

<p>The dipole moment of a single water molecule in liquid water has been a critical concept for understanding water’s dielectric properties. In this work, we investigate the dipole moment of liquid water through a self-attractive Hartree (SAH) decomposition of total electron density computed by density functional theory, on water clusters sampled from ab initio molecular dynamics simulation of bulk water. By adjusting one parameter that controls the degree of density localization, we reveal two distinct pictures of water dipoles that are consistent with bulk dielectric properties: a localized picture with smaller and less polarizable monomer dipoles, and a delocalized picture with larger and more polarizable monomer dipoles. We further uncover that the collective dipole-dipole correlation is stronger in the localized picture and is key to connecting individual dipoles with bulk dielectric properties. Based on these findings, we suggest considering both individual and collective dipole behaviors when studying the dipole moment of liquid water, and propose new design strategies for developing water models.</p>


Author(s):  
Mohammad Moulod ◽  
Gisuk Hwang

Fundamental understanding of the water in graphene is crucial to optimally design and operate the sustainable energy, water desalination, and bio-medical systems. A numerous atomic-scale studies have been reported, primarily articulating the surface interactions (interatomic potentials) between the water and graphene. However, a systematic comparative study among the various interatomic potentials is rare, especially for the water transport confined in the graphene nanostructure. In this study, the effects of different interatomic potentials and gap sizes on water self-diffusivity are investigated using the molecular dynamics simulation at T = 300 K. The water is confined in the rigid graphene nanogap with the various gap sizes Lz = 0.7 to 4.17 nm, using SPC/E and TIP3P water models. The water self-diffusivity is calculated using the mean squared displacement approach. It is found that the water self-diffusivity in the confined region is lower than that of the bulk water, and it decreases as the gap size decreases and the surface energy increases. Also, the water self-diffusivity nearly linearly decreases with the increasing surface energy to reach the bulk water self-diffusivity at zero surface energy. The obtained results provide a roadmap to fundamentally understand the water transport properties in the graphene geometries and surface interactions.


Sign in / Sign up

Export Citation Format

Share Document