TY - JOUR
T1 - A process-level perspective of the impact of molecular force fields on the computational screening of MOFs for carbon capture
AU - Cleeton, C
AU - de Oliveira, FL
AU - Neumann, RF
AU - Farmahini, AH
AU - Luan, BQ
AU - Steiner, M
AU - Sarkisov, L
PY - 2023/9/1
Y1 - 2023/9/1
N2 - The question we pose in this study is to what extent the ranking of metal organic frameworks (MOFs) for adsorption-based carbon capture, and the selection of top performers identified in Pressure Swing Adsorption (PSA) process modelling, depends on the choice of the commonly available forcefields. To answer this question, we first generated distributions of CO
2 and N
2 adsorption isotherms via molecular simulation in 690 MOFs using six typical forcefields: the UFF or Dreiding sets of Lennard-Jones parameters, in combination with partial charges derived from ab initio calculations or by charge equilibration schemes. We then conducted a systematic uncertainty quantification study using PSA process-level modelling. We observe that: (i) the ranking of MOFs significantly depends on the choice of forcefield; (ii) partial charge assignment is the prevailing source of uncertainty, and that charge equilibration schemes produce results which are inconsistent with ab initio-derived charges; (iii) the choice of Lennard-Jones parameters is a considerable source of uncertainty. Our work highlights that is not really possible to obtain material rankings with high resolution using a single molecular modelling approach and that, as a minimum, some uncertainty should be estimated for the performance of MOFs shortlisted using high throughput computational screening workflows. Future prospects for computational screening studies are also discussed.
AB - The question we pose in this study is to what extent the ranking of metal organic frameworks (MOFs) for adsorption-based carbon capture, and the selection of top performers identified in Pressure Swing Adsorption (PSA) process modelling, depends on the choice of the commonly available forcefields. To answer this question, we first generated distributions of CO
2 and N
2 adsorption isotherms via molecular simulation in 690 MOFs using six typical forcefields: the UFF or Dreiding sets of Lennard-Jones parameters, in combination with partial charges derived from ab initio calculations or by charge equilibration schemes. We then conducted a systematic uncertainty quantification study using PSA process-level modelling. We observe that: (i) the ranking of MOFs significantly depends on the choice of forcefield; (ii) partial charge assignment is the prevailing source of uncertainty, and that charge equilibration schemes produce results which are inconsistent with ab initio-derived charges; (iii) the choice of Lennard-Jones parameters is a considerable source of uncertainty. Our work highlights that is not really possible to obtain material rankings with high resolution using a single molecular modelling approach and that, as a minimum, some uncertainty should be estimated for the performance of MOFs shortlisted using high throughput computational screening workflows. Future prospects for computational screening studies are also discussed.
UR - https://www.webofscience.com/api/gateway?GWVersion=2&SrcApp=pure_starter&SrcAuth=WosAPI&KeyUT=WOS:001039203700001&DestLinkType=FullRecord&DestApp=WOS
UR - http://www.scopus.com/inward/record.url?scp=85167521388&partnerID=8YFLogxK
UR - https://www.mendeley.com/catalogue/f024a55b-d559-365b-a055-141a8bb75633/
U2 - 10.1039/d3ee00858d
DO - 10.1039/d3ee00858d
M3 - Article
SN - 1754-5692
VL - 16
SP - 3899
EP - 3918
JO - Energy & Environmental Science
JF - Energy & Environmental Science
IS - 9
ER -