scholarly journals Self-Diffusion Coefficients of Methane/n-Hexane Mixtures at High Pressures: An Evaluation of the Finite-Size Effect and a Comparison of Force Fields

2019 ◽  
Author(s):  
Thiago José Pinheiro dos Santos ◽  
Charlles Abreu ◽  
Bruno Horta ◽  
Frederico W. Tavares
Keyword(s):  
Molecular Dynamics ◽  
Diffusion Coefficients ◽  
High Pressures ◽  
Force Fields ◽  
Transport Coefficients ◽  
Finite Size ◽  
Finite Size Effect ◽  
Self Diffusion ◽  
Analytical Correction ◽  
Dynamics Simulations

Mass transport coefficients play an important role in process design and in compositional grading of oil reservoirs. As experimental measurements of these properties can be costly and hazardous, Molecular Dynamics simulations emerge as an alternative approach. In this work, we used Molecular Dynamics to calculate the self-diffusion coefficients of methane/n-hexane mixtures at different conditions, in both liquid and supercritical phases. We evaluated how the finite box size and the choice of the force field affect the calculated properties at high pressures. Results show a strong dependency between self-diffusion and the simulation box size. The Yeh-Hummer analytical correction [J. Phys. Chem. B, 108, 15873 (2004)] can attenuate this effect, but sometimes makes the results depart from experimental data due to issues concerning the force fields. We have also found that different all-atom and united-atom models can produce biased results due to caging effects and to different dihedral configurations of the n-alkane.

scholarly journals Self-Diffusion Coefficients of Methane/n-Hexane Mixtures at High Pressures: An Evaluation of the Finite-Size Effect and a Comparison of Force Fields

2019 ◽  
Author(s):  
Thiago José Pinheiro dos Santos ◽  
Charlles Abreu ◽  
Bruno Horta ◽  
Frederico W. Tavares
Keyword(s):  
Molecular Dynamics ◽  
Diffusion Coefficients ◽  
High Pressures ◽  
Force Fields ◽  
Transport Coefficients ◽  
Finite Size ◽  
Finite Size Effect ◽  
Self Diffusion ◽  
Analytical Correction ◽  
Dynamics Simulations

Mass transport coefficients play an important role in process design and in compositional grading of oil reservoirs. As experimental measurements of these properties can be costly and hazardous, Molecular Dynamics simulations emerge as an alternative approach. In this work, we used Molecular Dynamics to calculate the self-diffusion coefficients of methane/n-hexane mixtures at different conditions, in both liquid and supercritical phases. We evaluated how the finite box size and the choice of the force field affect the calculated properties at high pressures. Results show a strong dependency between self-diffusion and the simulation box size. The Yeh-Hummer analytical correction [J. Phys. Chem. B, 108, 15873 (2004)] can attenuate this effect, but sometimes makes the results depart from experimental data due to issues concerning the force fields. We have also found that different all-atom and united-atom models can produce biased results due to caging effects and to different dihedral configurations of the n-alkane.

scholarly journals Fick Diffusion Coefficients via Molecular Dynamics: An Alternative Approach in the Fourier Domain

2020 ◽  
Author(s):  
Thiago José Pinheiro dos Santos ◽  
Frederico W. Tavares ◽  
Charlles Abreu
Keyword(s):  
Molecular Dynamics ◽  
Diffusion Coefficients ◽  
High Pressures ◽  
Finite Size ◽  
Correlation Method ◽  
Finite Size Effect ◽  
New Approach ◽  
Molecular Systems ◽  
Alternative Approach ◽  
Wide Range

Mutual diffusion coefficient data are required for several systems of scientific and engineering interest to properly describe mass transport phenomena over a wide range of pressures, temperatures, and compositions. In this work, we calculated Fick diffusion coefficients for some CO2+n-alkane mixtures at high pressures using a new method, which we derived by introducing modifications to the Fourier Correlation Method (FCM) originally proposed by Nichols and Wheeler [I&EC Research, 54, 12156–12164 (2015)]. The modified FCM (mFCM) results were validated through comparisons with experimental data and with Fick coefficients calculated by employing well-established Molecular Dynamics methodologies. The new approach has some interesting advantages, such as providing Fick coefficients for molecular systems directly through a single equilibrium calculation, in contrast to traditional methods in which an extra calculation is needed to obtain the so-called thermodynamic factor. It is shown that the new approach considerably reduces the finite-size effect of the simulation box on the calculated diffusion coefficients, which are thus obtained in the thermodynamic limit.<br>

scholarly journals Fick Diffusion Coefficients via Molecular Dynamics: An Alternative Approach in the Fourier Domain

2020 ◽  
Author(s):  
Thiago José Pinheiro dos Santos ◽  
Frederico W. Tavares ◽  
Charlles Abreu
Keyword(s):  
Molecular Dynamics ◽  
Diffusion Coefficients ◽  
High Pressures ◽  
Finite Size ◽  
Correlation Method ◽  
Finite Size Effect ◽  
New Approach ◽  
Molecular Systems ◽  
Alternative Approach ◽  
Wide Range

Mutual diffusion coefficient data are required for several systems of scientific and engineering interest to properly describe mass transport phenomena over a wide range of pressures, temperatures, and compositions. In this work, we calculated Fick diffusion coefficients for some CO2+n-alkane mixtures at high pressures using a new method, which we derived by introducing modifications to the Fourier Correlation Method (FCM) originally proposed by Nichols and Wheeler [I&EC Research, 54, 12156–12164 (2015)]. The modified FCM (mFCM) results were validated through comparisons with experimental data and with Fick coefficients calculated by employing well-established Molecular Dynamics methodologies. The new approach has some interesting advantages, such as providing Fick coefficients for molecular systems directly through a single equilibrium calculation, in contrast to traditional methods in which an extra calculation is needed to obtain the so-called thermodynamic factor. It is shown that the new approach considerably reduces the finite-size effect of the simulation box on the calculated diffusion coefficients, which are thus obtained in the thermodynamic limit.<br>

Molecular-Dynamics Analysis of the Structural Properties of Silica during Cooling

2008 ◽  
Vol 139 ◽  
pp. 101-106 ◽  
Author(s):  
Byoung Min Lee ◽  
Shinji Munetoh ◽  
Teruaki Motooka ◽  
Yeo Wan Yun ◽  
Kyu Mann Lee
Keyword(s):  
Molecular Dynamics ◽  
Glass Transition ◽  
Transition Temperature ◽  
Structural Properties ◽  
Diffusion Coefficients ◽  
The Self ◽  
Dynamics Analysis ◽  
Self Diffusion ◽  
Dynamics Simulations ◽  
Potential Parameters

The structural properties of SiO2 liquid during cooling have been investigated by molecular dynamics simulations. The interatomic forces acting on the particles are calculated by the modified Tersoff potential parameters. The glass transition temperature and structural properties of the resulting SiO2 system at various temperatures have been investigated. The fivefold coordinations of Si and threefold coordinations of O atoms were observed, and the coordination defects of system decrease with decreasing temperature up to 17 % at 300 K. The self-diffusion coefficients for Si and O atoms drop to almost zero below 3000 K. The structures were distorted at high temperatures, but very stable atomic network persisted up to high temperature in the liquid state.

A Molecular Dynamics Simulation of Molten (Li-Rb)Cl Implying the Chemla Effect of Mobilities

1980 ◽  
Vol 35 (5) ◽  
pp. 493-499 ◽  
Author(s):  
Isao Okada ◽  
Ryuzo Takagi ◽  
Kazutaka Kawamura
Keyword(s):  
Molecular Dynamics ◽  
Molecular Dynamics Simulation ◽  
Diffusion Coefficients ◽  
Strong Correlation ◽  
Pair Correlation ◽  
The Self ◽  
Dynamics Simulation ◽  
Self Diffusion ◽  
Pure Salt ◽  
Dynamics Simulations

Abstract A new transport property, the self-exchange velocity (SEV) of neighbouring unlike ions, has been evaluated from molecular dynamics simulations of molten LiCl, RbCl and LiRbCl2 at 1100 K and the mixture at 750 K. From the increase of the SEV's in the order Rb+ (pure salt) <Li+ (mixture) < Rb+ (mixture) < Li+ (pure salt), it is conjectured that there is a strong correlation between the SEV’s and the internal mobilities. An interpretation of the Chemla effect in its dependence on temperature is given. The pair correlation functions and the self-diffusion coefficients are also calculated and discussed.

scholarly journals Molecular-Dynamics Simulation of Self-Diffusion of Molecular Hydrogen in X-Type Zeolite

2013 ◽  
Vol 2013 ◽  
pp. 1-7 ◽  
Author(s):  
Xiaoming Du
Keyword(s):  
Molecular Dynamics ◽  
Diffusion Coefficients ◽  
Hydrogen Diffusion ◽  
Mean Free Path ◽  
The Self ◽  
Dynamics Simulation ◽  
Pulse Field ◽  
Self Diffusion ◽  
Hydrogen Molecules ◽  
Dynamics Simulations

The self-diffusion of hydrogen in NaX zeolite has been studied by molecular-dynamics simulations for various temperatures and pressures. The results indicate that in the temperature range of 77–293 K and the pressure range of 10–2700 kPa, the self-diffusion coefficients are found to range from 1.61 × 10−9 m2·s−1to 3.66 × 10−8 m2·s−1which are in good agreement with the experimental values from the quasielastic neutron scattering (QENS) and pulse field gradients nuclear magnetic resonance (PFG NMR) measurements. The self-diffusion coefficients decrease with increasing pressure due to packing of sorbate-sorbate molecules which causes frequent collusion among hydrogen molecules in pores and increase with increasing temperature because increasing the kinetic energy of the gas molecules enlarges the mean free path of gas molecule. The activated energy for hydrogen diffusion determined from the simulation is pressure-dependent.

Molecular dynamics approach for crystal structures of methane A and B

2017 ◽  
Vol 28 (04) ◽  
pp. 1750048 ◽  
Author(s):  
César G. Galván ◽  
José M. Cabrera-Trujillo ◽  
Ivonne J. Hernández-Hernández ◽  
Luis A. Pérez
Keyword(s):  
Molecular Dynamics ◽  
High Pressures ◽  
Force Fields ◽  
Recent Experimental Data ◽  
Reactive Force ◽  
Carbon Structures ◽  
Reactive Force Fields ◽  
Temperature And Pressure ◽  
Dynamics Simulations ◽  
Good Agreement

The carbon structures of phases A and B of methane are investigated through classical molecular dynamics simulations using optimized potentials for liquid simulations all-atom force fields as well as ReaxFF reactive force fields. Both final thermodynamic states were obtained by the proper ramping of temperature and pressure through well-known regions of methane’s phase diagram using the isothermal–isobaric (NPT) ensemble. Our calculated structures are in good agreement with very recent experimental data. The knowledge of these phases is the basis for the study of methane at high pressures.

Sign in / Sign up

Export Citation Format

Share Document