Quantifying Pore Width Effects on Diffusivity via a Novel 3D Stochastic Approach with Input from Atomistic Molecular Dynamics Simulations

Maria Apostolopoulou, Mirella S. Santos, Muhammad Hamza, Tai Bui, Ioannis G. Economou, Michail Stamatakis, Alberto Striolo

Research output: Contribution to journalArticle

Abstract

The increased production of unconventional hydrocarbons emphasizes the need to understand the transport of fluids through narrow pores. Although it is well-known that confinement affects fluids structure and transport, it is not yet possible to quantitatively predict properties such as diffusivity as a function of pore width in the range of 1-50 nm. Such pores are commonly found not only in shale rocks but also in a wide range of engineering materials, including catalysts. We propose here a novel and computationally efficient methodology to obtain accurate diffusion coefficient predictions as a function of pore width for pores carved out of common materials, such as silica, alumina, magnesium oxide, calcite, and muscovite. We implement atomistic molecular dynamics (MD) simulations to quantify fluid structure and transport within 5 nm-wide pores, with particular focus on the diffusion coefficient within different pore regions. We then use these data as input to a bespoke stochastic kinetic Monte Carlo (KMC) model, developed to predict fluid transport in mesopores. The KMC model is used to extrapolate the fluid diffusivity for pores of increasing width. We validate the approach against atomistic MD simulation results obtained for wider pores. When applied to supercritical methane in slit-shaped pores, our methodology yields data within 10% of the atomistic simulation results, with significant savings in computational time. The proposed methodology, which combines the advantages of MD and KMC simulations, is used to generate a digital library for the diffusivity of gases as a function of pore chemistry and pore width and could be relevant for a number of applications, from the prediction of hydrocarbon transport in shale rocks to the optimization of catalysts, when surface-fluid interactions impact transport.

Original languageEnglish
JournalJournal of Chemical Theory and Computation
DOIs
Publication statusAccepted/In press - 1 Jan 2019

Fingerprint

diffusivity
Molecular dynamics
molecular dynamics
porosity
Fluids
Computer simulation
simulation
Shale
Hydrocarbons
Kinetics
fluids
Rocks
Magnesium Oxide
Catalysts
Digital libraries
Aluminum Oxide
Calcium Carbonate
Magnesia
Calcite
Methane

ASJC Scopus subject areas

  • Computer Science Applications
  • Physical and Theoretical Chemistry

Cite this

Quantifying Pore Width Effects on Diffusivity via a Novel 3D Stochastic Approach with Input from Atomistic Molecular Dynamics Simulations. / Apostolopoulou, Maria; Santos, Mirella S.; Hamza, Muhammad; Bui, Tai; Economou, Ioannis G.; Stamatakis, Michail; Striolo, Alberto.

In: Journal of Chemical Theory and Computation, 01.01.2019.

Research output: Contribution to journalArticle

@article{3380a31abe264f43be4c75002fc83ebb,
title = "Quantifying Pore Width Effects on Diffusivity via a Novel 3D Stochastic Approach with Input from Atomistic Molecular Dynamics Simulations",
abstract = "The increased production of unconventional hydrocarbons emphasizes the need to understand the transport of fluids through narrow pores. Although it is well-known that confinement affects fluids structure and transport, it is not yet possible to quantitatively predict properties such as diffusivity as a function of pore width in the range of 1-50 nm. Such pores are commonly found not only in shale rocks but also in a wide range of engineering materials, including catalysts. We propose here a novel and computationally efficient methodology to obtain accurate diffusion coefficient predictions as a function of pore width for pores carved out of common materials, such as silica, alumina, magnesium oxide, calcite, and muscovite. We implement atomistic molecular dynamics (MD) simulations to quantify fluid structure and transport within 5 nm-wide pores, with particular focus on the diffusion coefficient within different pore regions. We then use these data as input to a bespoke stochastic kinetic Monte Carlo (KMC) model, developed to predict fluid transport in mesopores. The KMC model is used to extrapolate the fluid diffusivity for pores of increasing width. We validate the approach against atomistic MD simulation results obtained for wider pores. When applied to supercritical methane in slit-shaped pores, our methodology yields data within 10{\%} of the atomistic simulation results, with significant savings in computational time. The proposed methodology, which combines the advantages of MD and KMC simulations, is used to generate a digital library for the diffusivity of gases as a function of pore chemistry and pore width and could be relevant for a number of applications, from the prediction of hydrocarbon transport in shale rocks to the optimization of catalysts, when surface-fluid interactions impact transport.",
author = "Maria Apostolopoulou and Santos, {Mirella S.} and Muhammad Hamza and Tai Bui and Economou, {Ioannis G.} and Michail Stamatakis and Alberto Striolo",
year = "2019",
month = "1",
day = "1",
doi = "10.1021/acs.jctc.9b00776",
language = "English",
journal = "Journal of Chemical Theory and Computation",
issn = "1549-9618",
publisher = "American Chemical Society",

}

TY - JOUR

T1 - Quantifying Pore Width Effects on Diffusivity via a Novel 3D Stochastic Approach with Input from Atomistic Molecular Dynamics Simulations

AU - Apostolopoulou, Maria

AU - Santos, Mirella S.

AU - Hamza, Muhammad

AU - Bui, Tai

AU - Economou, Ioannis G.

AU - Stamatakis, Michail

AU - Striolo, Alberto

PY - 2019/1/1

Y1 - 2019/1/1

N2 - The increased production of unconventional hydrocarbons emphasizes the need to understand the transport of fluids through narrow pores. Although it is well-known that confinement affects fluids structure and transport, it is not yet possible to quantitatively predict properties such as diffusivity as a function of pore width in the range of 1-50 nm. Such pores are commonly found not only in shale rocks but also in a wide range of engineering materials, including catalysts. We propose here a novel and computationally efficient methodology to obtain accurate diffusion coefficient predictions as a function of pore width for pores carved out of common materials, such as silica, alumina, magnesium oxide, calcite, and muscovite. We implement atomistic molecular dynamics (MD) simulations to quantify fluid structure and transport within 5 nm-wide pores, with particular focus on the diffusion coefficient within different pore regions. We then use these data as input to a bespoke stochastic kinetic Monte Carlo (KMC) model, developed to predict fluid transport in mesopores. The KMC model is used to extrapolate the fluid diffusivity for pores of increasing width. We validate the approach against atomistic MD simulation results obtained for wider pores. When applied to supercritical methane in slit-shaped pores, our methodology yields data within 10% of the atomistic simulation results, with significant savings in computational time. The proposed methodology, which combines the advantages of MD and KMC simulations, is used to generate a digital library for the diffusivity of gases as a function of pore chemistry and pore width and could be relevant for a number of applications, from the prediction of hydrocarbon transport in shale rocks to the optimization of catalysts, when surface-fluid interactions impact transport.

AB - The increased production of unconventional hydrocarbons emphasizes the need to understand the transport of fluids through narrow pores. Although it is well-known that confinement affects fluids structure and transport, it is not yet possible to quantitatively predict properties such as diffusivity as a function of pore width in the range of 1-50 nm. Such pores are commonly found not only in shale rocks but also in a wide range of engineering materials, including catalysts. We propose here a novel and computationally efficient methodology to obtain accurate diffusion coefficient predictions as a function of pore width for pores carved out of common materials, such as silica, alumina, magnesium oxide, calcite, and muscovite. We implement atomistic molecular dynamics (MD) simulations to quantify fluid structure and transport within 5 nm-wide pores, with particular focus on the diffusion coefficient within different pore regions. We then use these data as input to a bespoke stochastic kinetic Monte Carlo (KMC) model, developed to predict fluid transport in mesopores. The KMC model is used to extrapolate the fluid diffusivity for pores of increasing width. We validate the approach against atomistic MD simulation results obtained for wider pores. When applied to supercritical methane in slit-shaped pores, our methodology yields data within 10% of the atomistic simulation results, with significant savings in computational time. The proposed methodology, which combines the advantages of MD and KMC simulations, is used to generate a digital library for the diffusivity of gases as a function of pore chemistry and pore width and could be relevant for a number of applications, from the prediction of hydrocarbon transport in shale rocks to the optimization of catalysts, when surface-fluid interactions impact transport.

UR - http://www.scopus.com/inward/record.url?scp=85074712978&partnerID=8YFLogxK

UR - http://www.scopus.com/inward/citedby.url?scp=85074712978&partnerID=8YFLogxK

U2 - 10.1021/acs.jctc.9b00776

DO - 10.1021/acs.jctc.9b00776

M3 - Article

AN - SCOPUS:85074712978

JO - Journal of Chemical Theory and Computation

JF - Journal of Chemical Theory and Computation

SN - 1549-9618

ER -