Molecular simulation of hydrogen storage and transport in cellulose

Research output: Contribution to journalArticle

Abstract

In this work we describe a computational workflow to model the sorption and transport of molecular hydrogen in cellulose frameworks. The work demonstrates the value of the molecular dynamics code, DL_POLY and Monte Carlo code, DL_MONTE sharing common input formats to enhance the compatibility of the codes, being supported by DL_FIELD. Structures generated using cellulose-builder were processed by DL_FIELD to generate input files for DL_POLY using the OPLS_2005 force field. After relaxation in molecular dynamics, structures were used for GCMC simulations in DL_MONTE before passing back to DL_POLY to evaluate transport properties at different levels of sorption. While no hydrogen sorption was seen in pure crystalline cellulose, increasing separation between layers did allow sorption. When slit-pores were sufficiently wide, interactions with the cellulose led to the volumetric density of adsorbed hydrogen exceeding vacuum density at accessible partial pressures as well as allowing diffusion through the system. These model systems can give useful insight into the behaviour of amorphous cellulose in future simulation and experiment.
LanguageEnglish
Pages1-10
Number of pages10
JournalMolecular Simulation
Early online date27 Mar 2019
DOIs
StatusE-pub ahead of print - 27 Mar 2019

Keywords

  • Molecular modelling
  • Molecular Dynamics
  • Monte Carlo
  • Cellulose
  • Hydrogen Storage

Cite this

@article{97ce9b1464364e36bb9085e72c954fb0,
title = "Molecular simulation of hydrogen storage and transport in cellulose",
abstract = "In this work we describe a computational workflow to model the sorption and transport of molecular hydrogen in cellulose frameworks. The work demonstrates the value of the molecular dynamics code, DL_POLY and Monte Carlo code, DL_MONTE sharing common input formats to enhance the compatibility of the codes, being supported by DL_FIELD. Structures generated using cellulose-builder were processed by DL_FIELD to generate input files for DL_POLY using the OPLS_2005 force field. After relaxation in molecular dynamics, structures were used for GCMC simulations in DL_MONTE before passing back to DL_POLY to evaluate transport properties at different levels of sorption. While no hydrogen sorption was seen in pure crystalline cellulose, increasing separation between layers did allow sorption. When slit-pores were sufficiently wide, interactions with the cellulose led to the volumetric density of adsorbed hydrogen exceeding vacuum density at accessible partial pressures as well as allowing diffusion through the system. These model systems can give useful insight into the behaviour of amorphous cellulose in future simulation and experiment.",
keywords = "Molecular modelling, Molecular Dynamics, Monte Carlo, Cellulose, Hydrogen Storage",
author = "Megan Stalker and Robert Grant and Yong, {Chin W.} and Leyorla Ohene-Yeboah and Timothy Mays and Stephen Parker",
year = "2019",
month = "3",
day = "27",
doi = "10.1080/08927022.2019.1593975",
language = "English",
pages = "1--10",
journal = "Molecular Simulation",
issn = "0892-7022",
publisher = "Taylor and Francis",

}

TY - JOUR

T1 - Molecular simulation of hydrogen storage and transport in cellulose

AU - Stalker, Megan

AU - Grant, Robert

AU - Yong, Chin W.

AU - Ohene-Yeboah, Leyorla

AU - Mays, Timothy

AU - Parker, Stephen

PY - 2019/3/27

Y1 - 2019/3/27

N2 - In this work we describe a computational workflow to model the sorption and transport of molecular hydrogen in cellulose frameworks. The work demonstrates the value of the molecular dynamics code, DL_POLY and Monte Carlo code, DL_MONTE sharing common input formats to enhance the compatibility of the codes, being supported by DL_FIELD. Structures generated using cellulose-builder were processed by DL_FIELD to generate input files for DL_POLY using the OPLS_2005 force field. After relaxation in molecular dynamics, structures were used for GCMC simulations in DL_MONTE before passing back to DL_POLY to evaluate transport properties at different levels of sorption. While no hydrogen sorption was seen in pure crystalline cellulose, increasing separation between layers did allow sorption. When slit-pores were sufficiently wide, interactions with the cellulose led to the volumetric density of adsorbed hydrogen exceeding vacuum density at accessible partial pressures as well as allowing diffusion through the system. These model systems can give useful insight into the behaviour of amorphous cellulose in future simulation and experiment.

AB - In this work we describe a computational workflow to model the sorption and transport of molecular hydrogen in cellulose frameworks. The work demonstrates the value of the molecular dynamics code, DL_POLY and Monte Carlo code, DL_MONTE sharing common input formats to enhance the compatibility of the codes, being supported by DL_FIELD. Structures generated using cellulose-builder were processed by DL_FIELD to generate input files for DL_POLY using the OPLS_2005 force field. After relaxation in molecular dynamics, structures were used for GCMC simulations in DL_MONTE before passing back to DL_POLY to evaluate transport properties at different levels of sorption. While no hydrogen sorption was seen in pure crystalline cellulose, increasing separation between layers did allow sorption. When slit-pores were sufficiently wide, interactions with the cellulose led to the volumetric density of adsorbed hydrogen exceeding vacuum density at accessible partial pressures as well as allowing diffusion through the system. These model systems can give useful insight into the behaviour of amorphous cellulose in future simulation and experiment.

KW - Molecular modelling

KW - Molecular Dynamics

KW - Monte Carlo

KW - Cellulose

KW - Hydrogen Storage

U2 - 10.1080/08927022.2019.1593975

DO - 10.1080/08927022.2019.1593975

M3 - Article

SP - 1

EP - 10

JO - Molecular Simulation

T2 - Molecular Simulation

JF - Molecular Simulation

SN - 0892-7022

ER -